This model, doc_gray_code, shows a communications link using Gray-coded 8-PSK modulation. Gray coding is a technique often used in multilevel modulation schemes to minimize the bit error rate by ordering modulation symbols so that the binary representations of adjacent symbols differ by only one bit.
The example model includes these blocks:
The Random Integer Generator block serves as the source, producing a sequence of integers.
The Integer to Bit Converter block converts each integer into a corresponding binary representation.
The AWGN Channel block adds white Gaussian noise to the modulated data.
The M-PSK Demodulator Baseband block demodulates the corrupted data.
The Bit to Integer Converter block converts each binary representation to a corresponding integer.
One copy of the Error Rate
Calculation block (labeled
Error Rate Calculation1 in
this model) compares the demodulated integer data with the original
source data, yielding symbol error statistics. The output of the Error
Rate Calculation block is a three-element vector containing
the calculated error rate, the number of errors observed, and the
amount of data processed.
Another copy of the Error Rate Calculation library
Error Rate Calculation2 in this
model) compares the demodulated binary data with the binary representations
of the source data, yielding bit error statistics.
In this model, the M-PSK Modulator Baseband block:
Accepts binary-valued inputs that represent integers between 0 and M − 1, where M is the alphabet size
Maps binary representations to constellation points using a Gray-coded ordering
Produces unit-magnitude complex phasor outputs, with evenly spaced phases between 0 and 2π(M − 1)/M
The table indicates which binary representations in the input correspond to which phasors in the output. The second column of the table is an intermediate representation that the block uses in its computations.
|Modulator Input||Gray-Coded Ordering||Modulator Output|
The table below sorts the first two columns of the table above, according to the output values. This sorting makes it clearer that the overall effect of this subsystem is a Gray code mapping, as shown in the figure below. Notice that the numbers in the second column of the table below appear in counterclockwise order in the figure.
|Modulator Output||Modulator Input|
You can analyze the data that the example produces to compare theoretical performance with simulation performance.
The theoretical symbol error probability of MPSK is
erfc is the complementary error function, Es/N0 is
the ratio of energy in a symbol to noise power spectral density, and M is
the number of symbols.
To determine the bit error probability, the symbol error probability, PE, needs to be converted to its bit error equivalent. There is no general formula for the symbol to bit error conversion. Upper and lower limits are nevertheless easy to establish. The actual bit error probability, Pb, can be shown to be bounded by
The lower limit corresponds to the case where the symbols have undergone Gray coding. The upper limit corresponds to the case of pure binary coding.
To test the Gray code modulation scheme in this model, simulate
the graycode model for a range of Eb/N0 values.
If you want to study bit error rates but not symbol error rates, then
you can use the
user interface as described in BERTool.
The rest of this section studies both the bit and symbol error
rates and hence does not use
Because increasing the value of Eb/N0 lowers the number of errors produced, the length of each simulation must be increased to ensure that the statistics of the errors remain stable.
sim command to
run a Simulink® simulation
from the MATLAB® command
window, the following code generates data for symbol error rate and
bit error rate curves. It considers
values in the range 0 dB to 12 dB, in steps of 2 dB.
M = 8; Tsym = 0.2; BERVec = ; SERVec = ; EbNoVec = [0:2:12]; for n = 1:length(EbNoVec); EbNo = EbNoVec(n); sim('doc_gray_code') ; SERVec(n,:) = graySER; BERVec(n,:) = grayBER; end;
After simulating for the full set of Eb/N0 values, you can plot the results using these commands:
semilogy( EbNoVec,SERVec(:,1), 'o', EbNoVec, BERVec(:,1), '*' ); legend ( 'Symbol error rate', 'Bit error rate' ); xlabel ( 'Eb/No (dB)' ); ylabel( 'Error Probability' ); title ( 'Symbol and Bit Error Probability' );
As a further exercise, using data obtained from
berawgn, you can plot the theoretical curves
on the same axes with the simulation results. You can also compare
Gray coding with pure binary coding, by modifying the M-PSK
Modulator Baseband and M-PSK Demodulator Baseband blocks
so that their Constellation ordering parameters
Binary instead of