Main Content

awgn

Add white Gaussian noise to signal

Description

Y = awgn(X,snr) adds white Gaussian noise to the vector signal X. This syntax assumes that the power of X is 0 dBW. For more information about additive white Gaussian noise, see What is AWGN?

example

Y = awgn(X,snr,signalpower) accepts an input signal power value in dBW. To measure the power of X before adding noise, specify signalpower as 'measured'. The 'measured' option does not generate the requested average SNR for repeated awgn function calls in a loop if the input signal power varies over time due to fading and the coherence time of the channel is larger than the input duration.

example

Y = awgn(X,snr,signalpower,randobject) additionally accepts a random number stream object to generate normal random noise samples. For information about producing repeatable noise samples, see Tips.

example

Y = awgn(X,snr,signalpower,seed) specifies a seed value for initializing the normal random number generator that is used to add white Gaussian noise to the input signal.

Y = awgn(___,powertype) specifies the signal and noise power type as 'dB' or 'linear' in addition to the input arguments in any of the previous syntaxes. For information on the relationships between SNR and other measures of the relative power of the noise, such as Es/N0, and Eb/N0, see AWGN Channel Noise Level.

[Y,var] = awgn(___) also returns the total noise variance used to produce random noise samples.

Examples

collapse all

Create a sawtooth wave.

t = (0:0.1:60)';
x = sawtooth(t);

Add white Gaussian noise and plot the results.

y = awgn(x,10,'measured');
plot(t,[x y])
legend('Original Signal','Signal with AWGN')

Figure contains an axes object. The axes object contains 2 objects of type line. These objects represent Original Signal, Signal with AWGN.

Transmit and receive data using a nonrectangular 16-ary constellation in the presence of Gaussian noise. Show the scatter plot of the noisy constellation and estimate the symbol error rate (SER) for two different SNRs.

Create a 16-QAM constellation based on the V.29 standard for telephone-line modems.

c = [-5 -5i 5 5i -3 -3-3i -3i 3-3i 3 3+3i 3i -3+3i -1 -1i 1 1i];
sigpower = pow2db(mean(abs(c).^2));
M = length(c);

Generate random symbols.

data = randi([0 M-1],2000,1);

Modulate the data by using the genqammod function. General QAM modulation is necessary because the custom constellation is not rectangular.

modData = genqammod(data,c);

Pass the signal through an AWGN channel with a 20 dB SNR.

rxSig = awgn(modData,20,sigpower);

Display a scatter plot of the received signal and the reference constellation c.

h = scatterplot(rxSig);
hold on
scatterplot(c,[],[],'r*',h)
grid
hold off

Figure Scatter Plot contains an axes object. The axes object with title Scatter plot, xlabel In-Phase, ylabel Quadrature contains 2 objects of type line. One or more of the lines displays its values using only markers This object represents Channel 1.

Demodulate the received signal by using the genqamdemod function. Determine the number of symbol errors and the SER.

demodData = genqamdemod(rxSig,c);
[numErrors,ser] = symerr(data,demodData)
numErrors = 4
ser = 0.0020

Repeat the transmission and demodulation process with an AWGN channel with a 10 dB SNR. Determine the SER for the reduced SNR. As expected, the performance degrades when the SNR is decreased.

rxSig = awgn(modData,10,sigpower);
demodData = genqamdemod(rxSig,c);
[numErrors,ser] = symerr(data,demodData)
numErrors = 457
ser = 0.2285

Generate random data symbols and the 4-PSK modulated signal.

M = 4;
k = log2(M);
snr = 3;
data = randi([0 M-1],2000,1);
x = pskmod(data,M);

Set the random number generator seed.

seed = 12345;

Generate repeatable random noise using the rng function before calling the awgn function.

rng(seed);
y = awgn(x,snr);

Compute the bit errors.

dataHat = pskdemod(y,M);
numErr1 = biterr(data,dataHat,k)
numErr1 = 309

Reset the random number generator seed.

rng(seed);

Demodulate the PSK signal and compute the bit errors.

y = awgn(x,snr);
dataHat = pskdemod(y,M);
numErr2 = biterr(data,dataHat,k)
numErr2 = 309

Compare numErr1 to numErr2. The errors are equal even after you reset the random number generator seed.

isequal(numErr1, numErr2)
ans = logical
   1

Generate white Gaussian noise addition results by using a RandStream object and the reset object function.

Specify the input signal power of as 0 dBW, add noise to produce an SNR of 10 dB, and use a local random stream. Add white Gaussian noise to sigin two times to produce sigout1 and sigout2. Use isequal to compare sigout1 to sigout2. The outputs are not equal when you do not reset the random stream.

S = RandStream('mt19937ar',Seed=5489);
sigin = sqrt(2)*sin(0:pi/8:6*pi);
sigout1 = awgn(sigin,10,0,S);
sigout2 = awgn(sigin,10,0,S);
isequal(sigout1,sigout2)
ans = logical
   0

Reset the random stream object, returning the object to its state prior to adding AWGN to sigout1. Add AWGN to sigin to produce sigout3, and then compare sigout1 to sigout3. The outputs are equal when you reset the random stream.

reset(S);
sigout3 = awgn(sigin,10,0,S);
isequal(sigout1,sigout3)
ans = logical
   1

Input Arguments

collapse all

Input signal, specified as a scalar, vector, or a dlarray (Deep Learning Toolbox) object or a 3-D numeric array. The power of the input signal is assumed to be 0 dBW. If X is complex, awgn adds a complex noise. For more information, see Array Support.

Data Types: double
Complex Number Support: Yes

Signal-to-noise ratio in dB, specified as:

  • a scalar if the input signal is a scalar or a vector.

  • a scalar or a vector if the input signal is a 3D array. For more information see, Array Support.

The function applies the same snr value to each channel. The columns of the input signal represent the different channels of a multichannel signal.

Data Types: double

Signal power in dBW, specified as a scalar or 'measured'.

  • Scalar — The value is used as the signal level of X to determine the noise level required to achieve the specified snr.

  • 'measured' — The signal level of X is computed to determine the noise level required to achieve the specified snr. You cannot use this value if the input signal is a dlarray or if snr is a vector.

If the input signal is a multichannel signal, the function calculates the signalpower value across all channels as a single value. It then uses the value to calculate the noise level for all the channels.

Data Types: double

Random number stream object, specified as a RandStream object. The state of the random stream object determines the sequence of numbers produced by the randn function. Configure the random stream object using the reset (RandStream) function and its properties. You cannot specify a random stream object, if the input signal is a dlarray.

For information about producing repeatable noise samples, see Tips.

Random number generator seed value, specified as a scalar. You cannot specify a seed, if the input signal is a dlarray.

Data Types: double

Signal power unit, specified as 'dB' or 'linear'.

  • When powertype is 'dB', snr is measured in dB and signalpower is measured in dBW.

  • When powertype is 'linear', the snr is measured as a ratio and signalpower is measured in watts assuming a reference load of 1 ohms.

To set the powertype argument, you must also set snr and signalpower.

Output Arguments

collapse all

Output signal, returned as a scalar, vector, or a dlarray (Deep Learning Toolbox) object. The returned output signal is the input signal with added white Gaussian noise. If the input signal, X, is a dlarray, then this output, Y, is also a dlarray.

Total noise variance, returned as a positive scalar or a vector when snr is a vector. The function uses the noise variance to generate random noise samples.

More About

collapse all

What is AWGN?

Additive white Gaussian noise (AWGN) is a simple noise model that represents electron motion in the RF front end of a receiver. As the name implies, the noise gets added to the signal. The noise is called white because it is spectrally flat across the entire sampling bandwidth. Analogously, the color white contains equal spectral power levels at all frequencies of the visible light spectrum. The noise is Gaussian because its amplitude can be modeled with a normal probability distribution.

The AWGN channel is often used to model a satellite communications channel, since that channel typically does not suffer from common terrestrial impairments like fading, multipath, and interference. An AWGN channel serves as a good starting point for the analysis of terrestrial wireless links because it establishes a best-case bound on the bit error rate performance of a terrestrial link.

Array Support

In awgn function, the input signal X can be a numeric array, dlarray (Deep Learning Toolbox), or gpuArray (Parallel Computing Toolbox).

  • X — The input signal can be an array of up to three dimensions, specified as Ns-by-Nc-by-Nb array.

    Ns is the number of samples, Nc is the number of channels, and Nb is the number of batches.

  • snr can be a scalar or a vector of size Nb.

Tips

  • For information on the relationships between SNR and other measures of the relative power of the noise, such as Es/N0, and Eb/N0, see AWGN Channel Noise Level.

  • To generate repeatable white Gaussian noise samples, do one of the following:

    • Use rng(seed) before calling the awgn function to generate repeatable random noise.

    • Provide a static seed value as an input to awgn except when the input is a dlarray object.

    • Use the reset (RandStream) function on the randobject before passing it as an input to awgn except when the input is a dlarray object.

    • Provide randobject in a known state as an input to awgn except when the input is a dlarray object. For more information, see RandStream.

Extended Capabilities

Version History

Introduced before R2006a

expand all