Punctured Convolutional Encoding
This model shows how to use the Convolutional Encoder and Viterbi Decoder blocks to simulate a punctured coding system. The complexity of a Viterbi decoder increases rapidly with the code rate. Puncturing is a technique that allows the encoding and decoding of higher rate codes using standard rate 1/2 encoders and decoders.
The example is somewhat similar to the one that appears in Soft-Decision Decoding, which shows convolutional coding without puncturing.
Structure of the Example
This example contains these blocks.
Bernoulli Binary Generator: Create a sequence of random bits to use as a message.
Convolutional Encoder: Encode the message using the convolutional encoder.
BPSK Modulator Baseband: Modulate the encoded message.
AWGN Channel: Pass the modulated signal through a noisy channel.
Error Rate Calculation: Compute the number of discrepancies between the original and recovered messages.
Open the example,
doc_punct_conv_code, by entering the
following at the MATLAB® command
Generating Random Data
The Bernoulli Binary Generator block produces the information source for this simulation. The block generates a frame of three random bits at each sample time. The Samples per frame parameter determines the number of rows of the output frame.
Convolutional Encoding with Puncturing
The Convolutional Encoder block encodes the data from the Bernoulli Binary Generator. This example uses the same code as described in Soft-Decision Decoding.
The puncture pattern is specified by the Puncture vector parameter in the mask. The puncture vector is a binary column vector. A 1 indicates that the bit in the corresponding position of the input vector is sent to the output vector, while a 0 indicates that the bit is removed.
For example, to create a rate 3/4 code from the rate 1/2, constraint length 7 convolutional code, the optimal puncture vector is [1 1 0 1 1 0].' (where the .' after the vector indicates the transpose). Bits in positions 1, 2, 4, and 5 are transmitted, while bits in positions 3 and 6 are removed. Now, for every 3 bits of input, the punctured code generates 4 bits of output (as opposed to the 6 bits produced before puncturing). This makes the rate 3/4.
In this example, the output from the Bernoulli Binary Generator is a column vector of length 3. Because the rate 1/2 Convolutional Encoder doubles the length of each vector, the length of the puncture vector must divide 6.
The AWGN Channel block simulates transmission over a noisy channel. The parameters for the block are set in the mask as follows:
The Mode parameter for this block is set to
Signal to noise ratio (Es/No).
The Es/No parameter is set to 2 dB. This value typically is changed from one simulation run to the next.
The preceding modulation block generates unit power signals so the Input signal power is set to 1 Watt.
The Symbol period is set to 0.75 seconds because the code has rate 3/4.
In this simulation, the Viterbi Decoder block is set to accept unquantized inputs. As a result, the simulation passes the channel output through a Simulink® Complex to Real-Imag block that extracts the real part of the complex samples.
Viterbi Decoding of Punctured Codes
The Viterbi Decoder block is configured to decode the same rate 1/2 code specified in the Convolutional Encoder block.
In this example, the decision type is set to
codes without puncturing, you would normally set the Traceback
depth for this code to a value close to 40. However, for decoding
punctured codes, a higher value is required to give the decoder enough data to
resolve the ambiguities introduced by the punctures.
Since the punctured bits are not transmitted, there is no information to indicate their values. As a result they are ignored in the decoding process.
The Puncture vector parameter indicates the locations of the punctures or the bits to ignore in the decoding process. Each 1 in the puncture vector indicates a transmitted bit while each 0 indicates a puncture or the bit to ignore in the input to the decoder.
In general, the two Puncture vector parameters in the Convolutional Encoder and Viterbi Decoder must be the same.
Calculating the Error Rate
The Error Rate Calculation block compares the decoded bits to the original source bits. The output of the Error Rate Calculation block is a three-element vector containing the calculated bit error rate (BER), the number of errors observed, and the number of bits processed.
In the mask for this block, the Receive delay parameter is set to 96, because the Traceback depth value of 96 in the Viterbi Decoder block creates a delay of 96. If there were other blocks in the model that created delays, the Receive delay would equal the sum of all the delays.
BER simulations typically run until a minimum number of errors have occurred, or until the simulation processes a maximum number of bits. The Error Rate Calculation block uses its Stop simulation mode to set these limits and to control the duration of the simulation.
Generating a bit error rate curve requires multiple simulations. You can perform
multiple simulations using the
sim command. Follow these
In the model window, remove the Display block and the line connected to its port.
In the AWGN Channel block, set the Es/No parameter to the variable name
In the Error Rate Calculation block, set Output data to
Workspaceand then set Variable name to
Save the model in your working directory under a different name, such as
Execute the following code, which runs the simulation multiple times and gathers results.
CodeRate = 0.75; EbNoVec = [2:.5:5]; EsNoVec = EbNoVec + 10*log10(CodeRate); BERVec = zeros(length(EsNoVec),3); for n=1:length(EsNoVec), EsNodB = EsNoVec(n); sim('my_commpunctcnvcod'); BERVec(n,:) = BER_Data; end
To confirm the validity of the results, compare them to an established performance bound. The bit error rate performance of a rate r = (n-1)/n punctured code is bounded above by the expression:
In this expression,
erfc denotes the complementary error
function, r is the code rate, and both
ωd are dependent on the particular code. For the rate 3/4
code of this example, dfree = 5,
ω5 = 42, ω6 = 201,
ω7 = 1492, and so on. See reference [ 1 ] for more
The following commands compute an approximation of this bound using the first
seven terms of the summation (the values used for
nerr are found
in Table 2 of reference [ 2 ]:
dist = [5:11]; nerr = [42 201 1492 10469 62935 379644 2253373]; CodeRate = 3/4; EbNo_dB = [2:.02:5]; EbNo = 10.0.^(EbNo_dB/10); arg = sqrt(CodeRate*EbNo'*dist); bound = nerr*(1/6)*erfc(arg)';
To plot the simulation and theoretical results in the same figure, use the commands below.
berfit(EbNoVec',BERVec(:,1)); % Curve-fitted simulation results hold on; semilogy(EbNo_dB,bound,'g'); % Theoretical results legend('Simulated BER','Fit for simulated BER',... 'Theoretical bound on BER')
In some cases, at the lower bit error rates, you might notice simulation results that appear to indicate error rates slightly above the bound. This can result from simulation variance (if fewer than 500 bit errors are observed) or from the finite traceback depth in the decoder.
 Yasuda, Y., K. Kashiki, and Y. Hirata, "High Rate Punctured Convolutional Codes for Soft Decision Viterbi Decoding," IEEE Transactions on Communications, Vol. COM-32, March, 1984, pp. 315-319.
 Begin, G., Haccoun, D., and Paquin, C., "Further results on High-Rate Punctured Convolutional Codes for Viterbi and Sequential Decoding," IEEE Transactions on Communications, Vol. 38, No. 11, November, 1990, p. 1923.