Wave Walker DSP

DSP Algorithms for RF Systems

Trending

Buy the Book!

DSP for Beginners: Simple Explanations for Complex Numbers! The second edition includes a new chapter on complex sinusoids.

How to Create Additive White Gaussian Noise (AWGN)
June 1, 2022

Table of Contents

Introduction

Creating AWGN is useful in realistic simulations of DSP systems. For example, creating AWGN at proper power levels is useful in simulating bit error rates. If you’re writing a thesis or dissertation on DSP, or you plan to, you’re going to have to generate noise! Read on to find out how to do it correctly.

Check out these other blogs on DSP:

What is Noise?

Noise is what corrupts your signal. Noise is different than interference, which is another signal conflicting with your signal of interest.

However, adding up many types of interference produces Gaussian noise through the central limit theorem. Noise is the addition of all types of interference from cellular radio, AM and FM radio, broadcast TV. Noise also comes from places like the sun and lightning strikes.

Have you ever seen an image of an infrared camera? Physical bodies have a non-zero Kelvin temperature and therefore radiate which is what the infrared camera receives. This also contributes to the noise at your receiver.

The coolest type of noise is cosmic noise. Your radio is even affected by the noise generated from a quasar. I don’t know what that means but it sounds awesome!

Additive White Gaussian Noise (AWGN)

The central limit theorem allows the Gaussian distribution to be used as the model for AWGN. The Gaussian distribution is often referred to as the normal distribution. The noise in this blog is generated, or simulated, using NumPy’s random.normal() function.

Noise Power (or Variance)

The variance of AWGN in the real world will depend on the noise power, however the simplest way to start is to generate complex noise with a variance of 1 and then scale the noise to the desired power.

The variance \sigma^2 of the noise w[n] is given by

(1)   \begin{equation*}\sigma^2 = \mathbb{E} \left\{ (w[n] - \mu)(w[n] - \mu)^* \right\}\end{equation*}

where \mu is the mean of w[n],

(2)   \begin{equation*}\mu = \sum_{n=0}^{N-1} w[n].\end{equation*}

The mean \mu of AWGN is assumed to be zero, therefore the noise variance (or power) in (1) is

(3)   \begin{equation*}\sigma^2 = \mathbb{E} \left\{ w[n] w^*[n] \right\}\end{equation*}

Creating Real and Imaginary Noise

The time domain of AWGN [proakis2011, p.231] is Gaussian noise for both the real and imaginary components. The noise w[n] is constituted by its real and imaginary components,

(4)   \begin{equation*}w[n] = w_{R}[n] + jw_{I}[n].\end{equation*}

Substituting (4) into (3),

(5)   \begin{equation*}\begin{split}\sigma^2 & = \mathbb{E} \left\{ \left(w_{R}[n] + jw_{I}[n] \right) \left(w_{R}[n] + jw_{I}[n] \right)^* \right\} \\& = \mathbb{E} \left\{ w_{R}[n]^2 + w_{I}[n]^2 \right\} \\& = \mathbb{E} \left\{ w_{R}[n]^2 \right\} + \mathbb{E} \left\{ w_{I}[n]^2 \right\}.\end{split}\end{equation*}

The noise power for the real and imaginary terms are therefore

(6)   \begin{equation*}\sigma_{R}^2 = \mathbb{E} \left\{ w_{R}[n]^2 \right\},\end{equation*}

(7)   \begin{equation*}\sigma_{I}^2 = \mathbb{E} \left\{ w_{I}[n]^2 \right\}.\end{equation*}

The real and imaginary noise each contribute half of the power to the complex noise power. To generate complex noise with variance 1 then

(8)   \begin{equation*}\sigma_{R}^2 = \frac{1}{2},\end{equation*}

(9)   \begin{equation*}\sigma_{I}^2 = \frac{1}{2}.\end{equation*}

NumPy’s normal.random function requires the standard deviation \sigma of the Gaussian noise, which is

(10)   \begin{equation*}\begin{split}\sigma & = \sqrt{ \sigma^2 }\\& = \sqrt{ (1/2)^2 } \\& = 1/\sqrt{2} \\& = \sqrt{2}/2.\end{split}\end{equation*}

The noise for both real and imaginary channels is generated by:

noiseReal = np.random.normal(0,np.sqrt(2)/2,numSamples)
noiseImag = np.random.normal(0,np.sqrt(2)/2,numSamples)

which is the combined into the complex signal

noise = noiseReal + 1j*noiseImag

The variance of the complex noise can be checked by:

print(‘noise power = ‘ + str(np.mean(np.abs(noise)**2)))

An example of AWGN with variance and power 1.
An example of AWGN with variance and power 1.

Scaling the Noise Power

The above Python code creates complex AWGN with a variance (or power) of 1. However, the desired noise power will likely need to be a different value. The noise can be scaled to a different noise power P_n through the equation

(11)   \begin{equation*}w_{scaled}[n] = \sqrt{P_n} w[n]\end{equation*}

which is implemented in Python by:

noiseScaled = np.sqrt(Pn) * noise

Double-check the scaled noise power in Python:

print(‘scaled noise power = ‘ + str(np.mean(np.abs(noiseScaled)**2)))

An example of AWGN with power 0.1.
An example of AWGN with power 0.1.

Conclusion

Create complex noise by starting with real and imaginary noise with power 0.5 and then adding them together. The result is complex noise with power 1. Then scale the complex noise to the desired power.

Check out these other blogs on DSP:

Leave a Reply

God, the Lord, is my strength; He makes my feet like the deer's; He makes me tread on my high places. Habakkuk 3:19
For everything there is a season, and a time for every matter under heaven. A time to cast away stones, and a time to gather stones together. A time to embrace, and a time to refrain from embracing. Ecclesiastes 3:1,5
The earth was without form and void, and darkness was over the face of the deep. And the Spirit of God was hovering over the face of the waters. Genesis 1:2
Behold, I am toward God as you are; I too was pinched off from a piece of clay. Job 33:6
Enter His gates with thanksgiving, and His courts with praise! Give thanks to Him; bless His name! Psalm 100:4
Lift up your hands to the holy place and bless the Lord! Psalm 134:2
Blessed is the man who trusts in the Lord, whose trust is the Lord. He is like a tree planted by water, that sends out its roots by the stream, and does not fear when heat comes, for its leaves remain green, and is not anxious in the year of drought, for it does not cease to bear fruit. Jeremiah 17:7-8
He said to him, “You shall love the Lord your God with all your heart and with all your soul and with all your mind. This is the great and first commandment. And a second is like it: You shall love your neighbor as yourself. On these two commandments depend all the Law and the Prophets.” Matthew 22:37-39
Then He said to me, “Prophesy over these bones, and say to them, O dry bones, hear the word of the Lord. Thus says the Lord God to these bones: Behold, I will cause breath to enter you, and you shall live." Ezekiel 37:4-5
Riches do not profit in the day of wrath, but righteousness delivers from death. Proverbs 11:4
The angel of the Lord appeared to him in a flame of fire out of the midst of a bush. He looked, and behold, the bush was burning, yet it was not consumed. And Moses said, “I will turn aside to see this great sight, why the bush is not burned.” When the Lord saw that he turned aside to see, God called to him out of the bush, “Moses, Moses!” And he said, “Here I am.” Exodus 3:2-3
Daniel answered and said: “Blessed be the name of God forever and ever, to whom belong wisdom and might. He changes times and seasons; He removes kings and sets up kings; He gives wisdom to the wise and knowledge to those who have understanding." Daniel 2:20-21
Previous slide
Next slide

This website participates in the Amazon Associates program. As an Amazon Associate I earn from qualifying purchases.

© 2021-2024 Wave Walker DSP