AWGN Channel

Add white Gaussian noise to input signal

• 

Libraries:
Communications Toolbox / Channels

Description

The AWGN Channel block adds white Gaussian noise to the input signal. It inherits the sample time from the input signal.

Examples

LLR vs. Hard Decision Demodulation in Simulink

This model shows the improvement in BER performance when using log-likelihood ratio (LLR) instead of hard decision demodulation in a convolutionally coded communication link.

For a MATLAB® version of this example, see Log-Likelihood Ratio (LLR) Demodulation.

System Setup

This example model simulates a convolutionally coded communication system having one transmitter, an AWGN channel and three receivers. The convolutional encoder has a code rate of 1/2. The system employs a 16-QAM modulation. The modulated signal passes through an additive white Gaussian noise channel. The top receiver performs hard decision demodulation in conjunction with a Viterbi decoder that is set up to perform hard decision decoding. The second receiver has the demodulator configured to compute log-likelihood ratios (LLRs) that are then quantized using a 3-bit quantizer. It is well known that the quantization levels are dependent on noise variance for optimum performance [2]. The exact boundaries of the quantizer are empirically determined here. A Viterbi decoder that is set up for soft decision decoding processes these quantized values. The LLR values computed by the demodulator are multiplied by -1 to map them to the right quantizer index for use with Viterbi Decoder. To compute the LLR, the demodulator must be given the variance of noise as seen at its input. The third receiver includes a demodulator that computes LLRs which are processed by a Viterbi decoder that is set up in unquantized mode. The BER performance of each receiver is computed and displayed.

modelName = 'commLLRvsHD';
open_system(modelName);

System Simulation and Visualization

Simulate this system over a range of information bit Eb/No values. Adjust these Eb/No values for coded bits and multi-bit symbols to get noise variance values required for the AWGN block and Rectangular QAM Baseband Demodulator block. Collect BER results for each Eb/No value and visualize the results.

EbNo     = 2:0.5:8; % information rate Eb/No in dB
codeRate = 1/2;     % code rate of convolutional encoder
nBits    = 4;       % number of bits in a 16-QAM symbol
Pavg     = 10;      % average signal power of a 16-QAM modulated signal
snr      = EbNo - 10*log10(1/codeRate) + 10*log10(nBits); % SNR in dB
noiseVarVector = Pavg ./ (10.^(snr./10)); % noise variance

% Initialize variables for storing the BER results
ber_HD  = zeros(1,length(EbNo));
ber_SD  = zeros(1,length(EbNo));
ber_LLR = zeros(1, length(EbNo));

% Loop over all noiseVarVector values
for idx=1:length(noiseVarVector)
    noiseVar = noiseVarVector(idx); %#ok<NASGU>
    sim(modelName);
    % Collect BER results
    ber_HD(idx)  = BER_HD(1);
    ber_SD(idx)  = BER_SD(1);
    ber_LLR(idx) = BER_LLR(1);
end

% Perform curve fitting and plot the results
fitBER_HD  = real(berfit(EbNo,ber_HD));
fitBER_SD  = real(berfit(EbNo,ber_SD));
fitBER_LLR = real(berfit(EbNo,ber_LLR));
semilogy(EbNo,ber_HD,'r*', ...
    EbNo,ber_SD,'g*', ...
    EbNo,ber_LLR,'b*', ...
    EbNo,fitBER_HD,'r', ...
    EbNo,fitBER_SD,'g', ...
    EbNo,fitBER_LLR,'b');
legend('Hard Decision Decoding', ...
    'Soft Decision Decoding','Unquantized Decoding');
xlabel('Eb/No (dB)');
ylabel('BER');
title('LLR vs. Hard Decision Demodulation with Viterbi Decoding');
grid on;

To experiment with this system further, try different modulation types. This system uses a binary mapped modulation scheme for faster error collection but it is well known that Gray mapped signal constellation provides better BER performance. Experiment with various constellation ordering options in the modulator and demodulator blocks. Configure the demodulator block to compute approximate LLR to see the difference in the BER performance compared to hard decision demodulation and LLR. Try out a different range of Eb/No values. Finally, investigate different quantizer boundaries for your modulation scheme and Eb/No values.

Using Dataflow in Simulink

You can configure this example to use data-driven execution by setting the Domain parameter to dataflow for Dataflow Subsystem. With dataflow, blocks inside the domain, execute based on the availability of data as rather than the sample timing in Simulink®. Simulink automatically partitions the system into concurrent threads. This autopartitioning accelerates simulation and increases data throughput. To learn more about dataflow and how to run this example using multiple threads, see Multicore Simulation of Comparing Demodulation Types.

% Cleanup
close_system(modelName,0);
clear modelName EbNo codeRate nBits Pavg snr noiseVarVector ...
    ber_HD ber_SD ber_LLR idx noiseVar fitBER_HD fitBER_SD fitBER_LLR;

Selected Bibliography

[1] J. L. Massey, "Coding and Modulation in Digital Communications", Proc. Int. Zurich Seminar on Digital Communications, 1974

[2] J. A. Heller, I. M. Jacobs, "Viterbi Decoding for Satellite and Space Communication", IEEE® Trans. Comm. Tech. vol COM-19, October 1971

Gray-Coded M-PSK Modulation Error Rate in AWGN Channel Using Simulink

This example uses the doc_gray_code to compute bit error rates (BER) and symbol error rates (SER) for M-PSK modulation. The theoretical error rate performance of M-PSK modulation in AGWN is compared to the error rate performance for Gray-coded symbol mapping and to the error rate performance of binary-coded symbol mapping.

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 M-PSK Modulator Baseband block in the doc_gray_code model:

• Accepts binary-valued inputs that represent integers in the range [0, (M - 1], where M is the modulation order.

• Maps binary representations to constellation points using a Gray-coded ordering.

• Produces unit-magnitude complex phasor outputs, with evenly spaced phases in the range [0, (2 (M - 1) / M)].

The AWGN Channel block adds white Gaussian noise to the modulated data. The M-PSK Demodulator Baseband block demodulates the noisy data. The Bit to Integer Converter block converts each binary representation to a corresponding integer. Then two separate Error Rate Calculation blocks calculate the error rates of the demodulated data. The block labeled SER Calculation compares the integer data to compute the symbol error rate statistics and the block labeled BER Calculation compares the bits data to compute the bit error rate 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.

To reduce simulation run time and ensure that the statistics of the errors remain stable as the Eb/N0 ratio increases, the model is configured to run until 100 errors occur or until 1e8 bits have been transmitted.

The model initializes variables used to configure block parameters by using the PreLoadFcn callback function. For more information, see Model Callbacks (Simulink).

Produce Error Rate Curves

Compute the theoretical BER for nondifferential 8-PSK in AWGN over a range of Eb/N0 values by using the berawgn function. Simulate the doc_gray_code model with Gray-coded symbol mapping over the same range of Eb/N0 values.

Compare Gray coding with binary coding, by modifying the M-PSK Modulator Baseband and M-PSK Demodulator Baseband blocks to set the Constellation ordering parameter to Binary instead of Gray. Simulate the doc_gray_code model with binary-coded symbol mapping over the same range of Eb/N0 values.

Plot the results by using the semilogy function. The Gray-coded system achieves better error rate performance than the binary-coded system. Further, the Gray-coded error rate aligns with the theoretical error rate statistics.

Ports

Input

In — Input data signal

vector | matrix

Input data signal, specified as an N_S-by-1 vector or an N_S-by-N_C matrix. N_S represents the number of samples in the input signal. N_C represents the number of channels, as determined by the number of columns in the input signal matrix. Both N_S and N_C can be equal to 1.

The block adds frames of length-N_S Gaussian noise to each of the N_C channels, using a distinct random distribution per channel.

Data Types: double | single
Complex Number Support: Yes

Var — Variance of additive white Gaussian noise

positive scalar | vector

Variance of additive white Gaussian noise, specified as a positive scalar or a 1-by-N_C vector. N_C represents the number of channels, as determined by the number of columns in the input signal matrix. For more information, see Specifying the Variance Directly or Indirectly.

Dependencies

To enable this port, set Mode to Variance and Noise variance source to Input port.

Data Types: double

Output

Out — Output data signal

vector | matrix

Output data signal for the AWGN channel, returned as a vector or matrix. The datatype and dimensions of Out match those of the input signal, In.

Parameters

To edit block parameters interactively, use the Property Inspector. From the Simulink^® Toolstrip, on the Simulation tab, in the Prepare gallery, select Property Inspector.

Mode — Variance mode

Signal to noise ratio (SNR) (default) | Signal to noise ratio (Eb/No) | Signal to noise ratio (Es/No) | Variance

Variance mode, specified as Signal to noise ratio (Eb/No), Signal to noise ratio (Es/No), Signal to noise ratio (SNR), or Variance. For more information, see Relationship Among Eb/No, Es/No, and SNR Modes and Specifying the Variance Directly or Indirectly.

SNR (dB) — Ratio of signal power to noise power

10 (default) | scalar | vector

Ratio of signal power to noise power in decibels, specified as a scalar or vector.

Tunable: Yes

Dependencies

To enable this parameter, set Mode to Signal to noise ratio (SNR).

Eb/No (dB) — Ratio of information bit energy per symbol to noise power spectral density

10 (default) | scalar | vector

Ratio of information bit energy per symbol to noise power spectral density in decibels, specified as a scalar or vector. The information bit energy is the magnitude without channel coding.

Tunable: Yes

Dependencies

To enable this parameter, set Mode to Signal to noise ratio (Eb/No).

Es/No (dB) — Ratio of information symbol energy per symbol to noise power spectral density

10 (default) | scalar | vector

Ratio of information symbol energy per symbol to noise power spectral density in decibels, specified as a scalar or vector. The information symbol energy is the magnitude without channel coding.

Tunable: Yes

Dependencies

To enable this parameter, set Mode to Signal to noise ratio (Es/No).

Number of bits per symbol — Number of bits in each input symbol

1 (default) | scalar | vector

Number of bits in each input symbol, specified as a scalar or vector.

Dependencies

To enable this parameter, set Mode to Signal to noise ratio (Eb/No).

Input signal power, referenced to 1 ohm (watts) — Mean square power of input

1 (default) | scalar | vector

Mean square power of the input in watts, specified as a scalar or vector.

• When Mode is Signal to noise ratio (Eb/No) or Signal to noise ratio (Es/No), the parameter is the mean square power of the input symbols.

• When Mode is Signal to noise ratio (SNR), this parameter is the mean square power of the input samples.

Tunable: Yes

Dependencies

To enable this parameter, set Mode to Signal to noise ratio (Eb/No), Signal to noise ratio (Es/No), or Signal to noise ratio (SNR).

Samples per symbol — Samples per symbol

1 (default) | positive scalar | vector

Samples per symbol in samples, specified as a positive scalar or vector.

Dependencies

To enable this parameter, set Mode to Signal to noise ratio (Eb/No) or Signal to noise ratio (Es/No).

Noise variance source — Noise variance source

Parameter (default) | Input port

Noise variance source, specified as a Parameter or Input port. For more information, see Specifying the Variance Directly or Indirectly.

Tunable: Yes

Dependencies

To enable this parameter, set Mode to Variance.

Noise variance — Variance of white Gaussian noise

1 (default) | scalar | vector

Variance of the white Gaussian noise, specified as a scalar or vector. For more information, see Specifying the Variance Directly or Indirectly.

Tunable: Yes

Dependencies

To enable this port, set Mode to Variance and Noise variance source to Parameter.

Randomization

Random number source — Random number source

Global stream (default) | mt19937ar with seed

Random number source, specified as a Global stream or mt19937ar with seed.

Tunable: Yes

Initial seed — Noise generator initial seed

67 (default) | positive scalar | vector

Noise generator initial seed for the mt19937ar algorithm, specified as a positive scalar or a 1-by-N_C vector.

When the input signal is complex, the block creates random data as:

randData = randn(2*NS,NC) 
noise = randData(1:2:end) + 1i(randData(2:2:end))

N_S is the number of samples and N_C is the number of channels.

You can specify different seed values for each DLL build.

Tunable: Yes

Dependencies

To enable this parameter, set Random number source to mt19937ar with seed.

Simulate using — Type of simulation to run

Code generation (default) | Interpreted execution

Type of simulation to run, specified as Code generation or Interpreted execution.

• Code generation — Simulate the model by using generated C code. The first time you run a simulation, Simulink generates C code for the block. The model reuses the C code for subsequent simulations unless the model changes. This option requires additional startup time, but the speed of the subsequent simulations is faster than with the Interpreted execution option.

• Interpreted execution — Simulate the model by using the MATLAB^® interpreter. This option shortens startup time, but the speed of subsequent simulations is slower than with the Code generation option. In this mode, you can debug the source code of the block.

Block Characteristics

Data Types	double | single
Multidimensional Signals	no
Variable-Size Signals	no

Tips

• You can tune parameters in normal mode, accelerator mode, or rapid accelerator mode.

• Unless otherwise indicated, parameters are nontunable.

  • For nontunable parameters, when you use the Simulink Coder™ rapid simulation (RSIM) target to build an RSIM executable, you cannot change their values without recompiling the model.

  • If a parameter is tunable, you can change its value at any time. This is useful for Monte Carlo simulations in which you run the simulation multiple times (such as on multiple computers) with different amounts of noise.

Algorithms

Relationship Among Eb/No, Es/No, and SNR Modes

For uncoded complex input signals, the AWGN Channel block relates E_b/N_0, E_s/N₀, and SNR according to these equations:

E_s/N₀ = (T_sym/T_samp) × SNR

E_s/N₀ = E_b/N₀ + 10log₁₀(k) in dB

• E_s represents the signal energy in joules.

• E_b represents the bit energy in joules.

• N₀ represents the noise power spectral density in watts/Hz.

• T_sym/T_samp equals the value of the Samples per symbol parameter.

• T_sym is the symbol period of the signal, in seconds.

• T_samp represents the inherited sample time of the block, in seconds.

• k represents the number of information bits per input symbol, Number of bits per symbol.

For real signal inputs, the AWGN Channel block relates E_s/N₀ and SNR according to this equation:

E_s/N₀ = 0.5 (T_sym/T_samp) × SNR

Note

• All values of power assume a nominal impedance of 1 ohm.

• The equation for the real case differs from the corresponding equation for the complex case by a factor of 2. Specifically, the object uses a noise power spectral density of N₀/2 watts/Hz for real input signals, versus N₀ watts/Hz for complex signals.

For more information, see AWGN Channel Noise Level.

Specifying the Variance Directly or Indirectly

To directly specify the variance of the noise generated by AWGN Channel, specify the Mode as Variance. When you set Noise variance source to:

• Parameter — Specify the variance as a positive scalar or vector by using the Noise variance parameter.

• Input port — Specify the variance as a positive scalar or vector by using the Var block input. The variance input sampling rate must equal that of the input signal.

When you specify variance as a:

• Scalar, all signal channels are uncorrelated but share the same variance.

• Vector whose length is the number of channels in the input signal, each element represents the variance of the corresponding signal channel.

  Note

  If you apply complex input signals to the AWGN Channel block, then it adds complex zero-mean Gaussian noise with the calculated or specified variance. The variance for each quadrature component of the complex noise is half of the calculated or specified value.

To specify the variance indirectly, that is, to have the block calculate the variance, specify the Mode as:

• Signal to noise ratio (Eb/No), where the block calculates the variance from these quantities that you specify in the dialog box:

  • Eb/No (dB), the ratio of bit energy to noise power spectral density

  • Number of bits per symbol

  • Input signal power, referenced to 1 ohm (watts), the actual power of the symbols at the input of the block

  • Samples per symbol

• Signal to noise ratio (Es/No), where the block calculates the variance from these quantities that you specify in the dialog box:

  • Es/No (dB), the ratio of signal energy to noise power spectral density

  • Input signal power, referenced to 1 ohm (watts), the actual power of the symbols at the input of the block

  • Samples per symbol

• Signal to noise ratio (SNR), where the block calculates the variance from these quantities that you specify in the dialog box:

  • SNR (dB), the ratio of signal power to noise power

  • Input signal power, referenced to 1 ohm (watts), the actual power of the samples at the input of the block

Changing the symbol period in the AWGN Channel block affects the variance of the noise added per sample, which also causes a change in the final error rate.

$$NoiseVariance=\frac{SignalPower\times SymbolPeriod}{SampleTime\times {10}^{\frac{Es/No}{10}}}$$

Tip

Select the symbol period equal to the symbol period of the model. The value depends on what constitutes a symbol and what the oversampling applied to it is. For example, a symbol could have 3 bits and be oversampled by 4. For more information, see AWGN Channel Noise Level.

References

[1] Proakis, John G. Digital Communications. 4th Ed. McGraw-Hill, 2001.

Extended Capabilities

C/C++ Code Generation
Generate C and C++ code using Simulink® Coder™.

Version History

Introduced before R2006a

R2023a: Variable-size support

This support enables you to vary the length of input signal each time the block executes in your simulation.

R2023a: Random number generator update

The AWGN Channel block updates the random number generator (RNG) source permitting selection of Global stream or mt19937ar with seed.

The new RNG provides better random number statistics. For simulations created before R2023a, the block defaults to legacy mode with prior RNG ( Ziggurat method - V5 RANDN algorithm). In this mode, the block uses the Random Source block to generate noise. Random numbers are generated using the Ziggurat method (V5 RANDN algorithm). The block reuses the same initial seeds every time you rerun the simulation, so that this block outputs the same signal each time you run a simulation. To improve statistics, it is recommended that you run Update Model on pre-R2023a simulations that use the AWGN Channel block.

R2023a: Code Generation

Code generation is enabled by default. For more information, see Simulate using.