Updated comments to remove reference to a particular version of TS38.212 and to better reflect the operation of main_FAR.m.

Author: robmaunder

PBCH_decoder.m

@@ -1,6 +1,6 @@
 function a_hat = PBCH_decoder(f_tilde, A, L, min_sum, a_tilde)
 % PCBH_DECODER Polar decoder for the Public Broadcast Channel (PBCH) of 3GPP New Radio, as
-% defined in Section 7.1 of TS38.212 V1.2.1. Implements the Cyclic Redudancy
+% defined in Section 7.1 of TS38.212. Implements the Cyclic Redudancy
 % Check (CRC) attachment of Section 7.1.3, the channel coding of Section 7.1.4
 % and the rate matching of Section 7.1.5. Note that this code does not
 % implement the payload generation of Section 7.1.1 or the scrambling of
@@ -12,7 +12,7 @@
 %   f_tilde should be a real row vector comprising 864 Logarithmic
 %   Likelihood Ratios (LLRS), each having a value obtained as LLR =
 %   ln(P(bit=0)/P(bit=1)). The first LLR corresponds to f_0 from Section 
-%   7.1.5 of TS38.212 V1.2.1, while the last LLR corresponds to 
+%   7.1.5 of TS38.212, while the last LLR corresponds to 
 %   f_E-1.
 %
 %   A should be 32. It specifies the number of bits in the
@@ -39,12 +39,12 @@
 %   several higher-layer parameters. If scrambling is implemented
 %   externally to this code, then a_tilde should pertain to the scrambled
 %   bit values. The first input bit corresponds to a'_0 from 
-%   Section 7.1.3 of TS38.212 V1.2.1, while the last input bit corresponds 
+%   Section 7.1.3 of TS38.212, while the last input bit corresponds 
 %   to a'_A-1.
 %
 %   a_hat will be a binary row vector comprising 32 bits, each 
 %   having the value 0 or 1. The first output bit corresponds to a'_0 from 
-%   Section 7.1.3 of TS38.212 V1.2.1, while the last output bit corresponds 
+%   Section 7.1.3 of TS38.212, while the last output bit corresponds 
 %   to a'_A-1.
 %
 %

PBCH_encoder.m

@@ -1,6 +1,6 @@
 function f = PBCH_encoder(a, E)
 % PBCH_ENCODER Polar encoder for the Public Broadcast Channel (PBCH) of 3GPP New Radio, as
-% defined in Section 7.1 of TS38.212 V1.2.1. Implements the Cyclic Redudancy
+% defined in Section 7.1 of TS38.212. Implements the Cyclic Redudancy
 % Check (CRC) attachment of Section 7.1.3, the channel coding of Section 7.1.4
 % and the rate matching of Section 7.1.5. Note that this code does not
 % implement the payload generation of Section 7.1.1 or the scrambling of
@@ -10,7 +10,7 @@
 %
 %   a should be a binary row vector comprising 32 bits, each
 %   having the value 0 or 1. The first input bit corresponds to a'_0 from 
-%   Section 7.1.3 of TS38.212 V1.2.1, while the last input bit corresponds 
+%   Section 7.1.3 of TS38.212, while the last input bit corresponds 
 %   to a'_A-1.
 %
 %   E should be 864. It specifies the number of bits in the
@@ -19,7 +19,7 @@
 %
 %   f will be a binary row vector comprising 864 bits, each having
 %   the value 0 or 1. The first output bit corresponds to f_0 from Section 
-%   7.1.5 of TS38.212 V1.2.1, while the last output bit corresponds to 
+%   7.1.5 of TS38.212, while the last output bit corresponds to 
 %   f_E-1.
 %
 %   See also PBCH_DECODER

PDCCH_decoder.m

@@ -1,6 +1,6 @@
 function a_hat = PDCCH_decoder(f_tilde, A, L, min_sum, RNTI)
 % PDCCH_DECODER Polar decoder for the Physical Downlink Control Channel (PDCCH) of 3GPP New
-% Radio, as defined in Section 7.3 of TS38.212 V1.2.1. Implements the zero-
+% Radio, as defined in Section 7.3 of TS38.212. Implements the zero-
 % padding to increase the length of short payloads to 12 bits of Section 7.3.1,
 % the Cyclic Redudancy Check (CRC) attachment of Section 7.3.2, the channel
 % coding of Section 7.3.3 and the rate matching of Section 7.3.4. Note that
@@ -13,7 +13,7 @@
 %   f_tilde should be a real row vector comprising E number of Logarithmic
 %   Likelihood Ratios (LLRS), each having a value obtained as LLR =
 %   ln(P(bit=0)/P(bit=1)). The first LLR corresponds to f_0 from Section 
-%   7.3.4 of TS38.212 V1.2.1, while the last LLR corresponds to 
+%   7.3.4 of TS38.212, while the last LLR corresponds to 
 %   f_E-1.
 %
 %   A should be an integer scalar. It specifies the number of bits in the
@@ -31,11 +31,11 @@
 %   RNTI should be a binary row vector comprising 16 bits, each having the
 %   value 0 or 1. If this parameter is omitted, then ones(1,16) will be
 %   used for the RNTI. The first bit corresponds to x_rnti,0 from Section
-%   7.3.2 of TS38.212 V1.2.1, while the last bit corresponds to x_rnti,15.
+%   7.3.2 of TS38.212, while the last bit corresponds to x_rnti,15.
 %
 %   a_hat will be a binary row vector comprising A number of bits, each 
 %   having the value 0 or 1. The first output bit corresponds to a_0 from 
-%   Section 7.3.1 of TS38.212 V1.2.1, while the last output bit corresponds 
+%   Section 7.3.1 of TS38.212, while the last output bit corresponds 
 %   to a_A-1.
 %
 %   See also PDCCH_ENCODER

PDCCH_encoder.m

@@ -1,6 +1,6 @@
 function f = PDCCH_encoder(a, E, RNTI)
 % PDCCH_ENCODER Polar encoder for the Physical Downlink Control Channel (PDCCH) of 3GPP New
-% Radio, as defined in Section 7.3 of TS38.212 V1.2.1. Implements the zero-
+% Radio, as defined in Section 7.3 of TS38.212. Implements the zero-
 % padding to increase the length of short payloads to 12 bits of Section 7.3.1,
 % the Cyclic Redudancy Check (CRC) attachment of Section 7.3.2, the channel
 % coding of Section 7.3.3 and the rate matching of Section 7.3.4. Note that
@@ -11,7 +11,7 @@
 %
 %   a should be a binary row vector comprising A number of bits, each
 %   having the value 0 or 1. A should be in the range 1 to 140. The first
-%   input bit corresponds to a_0 from Section 7.3.1 of TS38.212 V1.2.1, while the
+%   input bit corresponds to a_0 from Section 7.3.1 of TS38.212, while the
 %   last input bit corresponds to a_A-1.
 %
 %   E should be an integer scalar. It specifies the number of bits in the
@@ -20,11 +20,11 @@
 %   RNTI should be a binary row vector comprising 16 bits, each having the
 %   value 0 or 1. If this parameter is omitted, then ones(1,16) will be
 %   used for the RNTI. The first bit corresponds to x_rnti,0 from Section
-%   7.3.2 of TS38.212 V1.2.1, while the last bit corresponds to x_rnti,15.
+%   7.3.2 of TS38.212, while the last bit corresponds to x_rnti,15.
 %
 %   f will be a binary row vector comprising E number of bits, each having
 %   the value 0 or 1. The first output bit corresponds to f_0 from Section 
-%   7.3.4 of TS38.212 V1.2.1, while the last output bit corresponds to 
+%   7.3.4 of TS38.212, while the last output bit corresponds to 
 %   f_E-1.
 %
 %   See also PDCCH_DECODER

PUCCH_decoder.m

@@ -1,7 +1,7 @@
 function a_hat = PUCCH_decoder(f_tilde, A, L, min_sum)
 % PUCCH_DECODER Polar decoder for the Physical Uplink Control Channel (PUCCH) and the
 % Physical Uplink Shared Channel (PUSCH) of 3GPP New Radio, as defined in
-% Section 6.3 of TS38.212 V1.2.1. Implements the code block segmentation and
+% Section 6.3 of TS38.212. Implements the code block segmentation and
 % Cyclic Redudancy Check (CRC) attachment of Sections 6.3.1.2.1 and 6.3.2.2.1,
 % the channel coding of Sections 6.3.1.3.1 and 6.3.2.3.2, the rate matching of
 % Sections 6.3.1.4.1 and 6.3.2.4.1, as well as the code block concatenation of
@@ -19,7 +19,7 @@
 %   Likelihood Ratios (LLRS), each having a value obtained as LLR =
 %   ln(P(bit=0)/P(bit=1)), where G should be no greater than 8192 if A<360
 %   and no greater than 16384 if A>=360. The first LLR corresponds to g_0 
-%   from Sections 6.3.1.5 and 6.3.2.5 of TS38.212 V1.2.1, while the last 
+%   from Sections 6.3.1.5 and 6.3.2.5 of TS38.212, while the last 
 %   LLR corresponds to g_G-1.
 %
 %   A should be an integer scalar. It specifies the number of bits in the
@@ -36,7 +36,7 @@
 %
 %   a_hat will be a binary row vector comprising A number of bits, each
 %   having the value 0 or 1. The first output bit corresponds to a_0 from 
-%   Sections 6.3.1.2 and 6.3.2.2 of TS38.212 V1.2.1, while the last output 
+%   Sections 6.3.1.2 and 6.3.2.2 of TS38.212, while the last output 
 %   bit corresponds to a_A-1.
 %
 %   See also PUCCH_ENCODER

PUCCH_encoder.m

@@ -1,7 +1,7 @@
 function f = PUCCH_encoder(a, G)
 % PUCCH_ENCODER Polar encoder for the Physical Uplink Control Channel (PUCCH) and the
 % Physical Uplink Shared Channel (PUSCH) of 3GPP New Radio, as defined in
-% Section 6.3 of TS38.212 V1.2.1. Implements the code block segmentation and
+% Section 6.3 of TS38.212. Implements the code block segmentation and
 % Cyclic Redudancy Check (CRC) attachment of Sections 6.3.1.2.1 and 6.3.2.2.1,
 % the channel coding of Sections 6.3.1.3.1 and 6.3.2.3.2, the rate matching of
 % Sections 6.3.1.4.1 and 6.3.2.4.1, as well as the code block concatenation of
@@ -16,16 +16,16 @@
 %
 %   a should be a binary row vector comprising A number of bits, each
 %   having the value 0 or 1, where A should be in the range 12 to 1706. The first
-%   input bit corresponds to a_0 from Sections 6.3.1.2 and 6.3.2.2 of TS38.212
-%   V1.2.1, while the last input bit corresponds to a_A-1.
+%   input bit corresponds to a_0 from Sections 6.3.1.2 and 6.3.2.2 of TS38.212, 
+%   while the last input bit corresponds to a_A-1.
 %
 %   G should be an integer scalar. It specifies the number of bits in the
 %   encoded bit sequence, where G should be no greater than 8192 if A<360
 %   and no greater than 16384 if A>=360.
 %
 %   f will be a binary row vector comprising G number of bits, each having
 %   the value 0 or 1. The first output bit corresponds to g_0 from Sections 
-%   6.3.1.5 and 6.3.2.5 of TS38.212 V1.2.1, while the last output bit 
+%   6.3.1.5 and 6.3.2.5 of TS38.212, while the last output bit 
 %   corresponds to g_G-1.
 %
 %   See also PUCCH_DECODER

main_FAR.m

@@ -1,8 +1,6 @@
 function main_FAR(code, A, E, L, min_sum, target_false_alarms, seed)
-% PLOT_BLER_VS_SNR Plots Block Error Rate (BLER) versus Signal to Noise
-% Ratio (SNR) for polar codes.
-%   plot_BLER_vs_SNR(code, A, E, L, min_sum, target_block_errors, target_BLER, EsN0_start, EsN0_delta, seed)
-%   generates the plots.
+% main_FAR calculates False Alarm Rate (FAR) for polar codes.
+%   main_FAR(code, A, E, L, min_sum, target_false_alarms, seed)
 %
 %   code should be a string. This identifies which encoder and decoder
 %   functions to call. For example, if code is 'custom1', then the
@@ -30,25 +28,12 @@ function main_FAR(code, A, E, L, min_sum, target_false_alarms, seed)
 %   better error correction capability than the min-sum, but it has higher
 %   complexity.
 %
-%   target_block_errors should be an integer scalar. The simulation of each
-%   SNR for each coding rate will continue until this number of block
-%   errors have been observed. A value of 100 is sufficient to obtain
+%   target_false_alarms should be an integer scalar. The simulation of each
+%   encoded block length will continue until this number of false
+%   alarms have been observed. A value of 100 is sufficient to obtain
 %   smooth BLER plots for most values of A. Higher values will give
 %   smoother plots, at the cost of requiring longer simulations.
 %
-%   target_BLER should be a real scalar, in the range (0, 1). The
-%   simulation of each coding rate will continue until the BLER plot
-%   reaches this value.
-%
-%   EsN0_start should be a real row vector, having the same length as the
-%   vector of coding rates. Each value specifies the Es/N0 SNR to begin at
-%   for the simulation of the corresponding coding rate.
-%
-%   EsN0_delta should be a real scalar, having a value greater than 0.
-%   The Es/N0 SNR is incremented by this amount whenever
-%   target_block_errors number of block errors has been observed for the
-%   previous SNR. This continues until the BLER reaches target_BLER.
-%
 %   seed should be an integer scalar. This value is used to seed the random
 %   number generator, allowing identical results to be reproduced by using
 %   the same seed. When running parallel instances of this simulation,
@@ -68,9 +53,9 @@ function main_FAR(code, A, E, L, min_sum, target_false_alarms, seed)
 
 % Default values
 if nargin == 0
-    code = 'custom1';
-    A = 16;
-    E = round(A./[0.8333 0.7500 0.6666 0.5000 0.4000 0.3333 0.2500 0.2000 0.1666 0.1250]);
+    code = 'PUCCH';
+    A = 32;
+    E = [54 108 216 432 864 1728 3456 6912 13824];
     L = 1;
     min_sum = true;
     target_false_alarms = 10;