Skip to main content

Gaussian random variable and its PDF in MATLAB


Home / Wireless Communication / Gaussian random variable and its PDF



What exactly are Gaussian Random Variable and its probability distribution function (PDF) are


 The practical communication system is modeled as 

y = x + n;

Where y=received  signal 

x= transmitted signal 

n= noise


What is the significance of the Gaussian Random Variable?  

We know, especially for wireless communication, whenever we transmit a signal from transmitter to receiver, there will be some additive white Gaussian noise to the signal when we receive it from the receiver. The additive white Gaussian noise has some properties, like zero mean and a specific standard deviation. We learn later what exactly they mean, what Deviations are, and the relation of the Gaussian random variable with it. Here, the word "random" is used because noise is always unexpected in the communication system. We can't predict it before the transmission of the signal. But we can draw its probability distribution function (PDF) from several experiments or values. 



What exactly is Gaussian Random Variable PDF is

PDF of Gaussian random variable is defined as


Here, Ïƒ = Standard Deviation of random variable samples

μ = mean of random variable samples

In the above figure probability distribution function of the Gaussian random variable is shown. Students often need clarification with the title of the x label and y label. x tag defines the variation of the standard deviation value of Gaussian noise collected from large samples or populations or many experiments. After getting the standard Deviation of noise,e we plot the probability of standard deviations derived from large samples. 


MATLAB Code for gaussian random variable and its PDF

clear;

close all;


% Number of samples to generate

n = 100000;


% Generate Gaussian distribution (Standard Normal Distribution)

gaussian_values = randn(1, n);  % Standard normal distribution (mean = 0, std = 1)


% Calculate mean and standard deviation of the Gaussian values

mu = mean(gaussian_values);

sigma = std(gaussian_values);


% Calculate the range for 1, 2, and 3 standard deviations

range_1sigma = sum(gaussian_values >= (mu - sigma) & gaussian_values <= (mu + sigma)) / n * 100;  % Percentage within 1 standard deviation

range_2sigma = sum(gaussian_values >= (mu - 2*sigma) & gaussian_values <= (mu + 2*sigma)) / n * 100;  % Percentage within 2 standard deviations

range_3sigma = sum(gaussian_values >= (mu - 3*sigma) & gaussian_values <= (mu + 3*sigma)) / n * 100;  % Percentage within 3 standard deviations


% Display the results

fprintf('Percentage of values within 1 standard deviation: %.2f%%\n', range_1sigma);

fprintf('Percentage of values within 2 standard deviations: %.2f%%\n', range_2sigma);

fprintf('Percentage of values within 3 standard deviations: %.2f%%\n', range_3sigma);


% Plotting the Gaussian distribution

figure;

histogram(gaussian_values, 30, 'Normalization', 'pdf');  % Normalized to show probability density

title('Gaussian Distribution (Standard Normal)');

xlabel('Value');

ylabel('Probability Density');


Output

Percentage of values within 1 standard deviation: 68.36%

Percentage of values within 2 standard deviations: 95.40%

Percentage of values within 3 standard deviations: 99.73%









Copy the aforementioned MATLAB Code from here003


Real-world mathematical examples to understand mean and standard Deviation

Mean of a Random Variable

As we have mentioned above, noise is random in a communication system. So, we take hundreds of values of that parameter and draw a PDF. For example, we have received ten random variables, i.e., X1, X2, X3, X4,..., X9, and X10. Then we calculate its mean or average. That is also meaningful.


Xmean = (X1 + X2 + X3+... +X8 +X9 +X10)/10


Standard Deviation of a Random Variable

In electronic communication, the standard signal deviation tells us how the signal varies over time. For example, we measure a signal in different time instants, from a different position, or at another aspect. Then we can calculate the standard Deviation to see how the signal varies. That value also matters for electronic devices. Similarly, we calculate the standard deviation value from many samples in the case of a Gaussian random variable. For example, in a class, marks obtained in math by students are as follows:

Student 1: 92 out of 100

Student 2: 85 out of 100

Student 3: 74 out of 100

Student 4: 70 out of 100

Student 5: 60 out of 100

Student 6: 66 out of 100

Student 7: 82 out of 100

Student 8: 63 out of 100

Student 9: 76 out of 100

Student 10: 59 out of 100


The average marks obtained by students are calculated as

=(92+85+74+70+60+66+82+63+76+59)/10

=72.7

The mean value is 72.7


Now, we'll calculate Standard Deviation,

Std or Ïƒ= sqrt{(1/(N-1) * Σ(Ni -N0)^2}

Here, N= total number of sample

Ni denotes the instantaneous value of  N

N0 denotes the mean of N

'sqrt' denotes 'square root' here


The standard Deviation for obtained marks by students is,

Std or Ïƒ =sqrt{1/(10-1) * Î£ (Ni -72.7)^2} 

(as here several samples or population is 10 & mean/avg. =72.7)

Or, Ïƒ = sqrt [1/9 * {(92-72.7)^2 + (85-72.7)^2 + (74-72.7)^2 + (70-72.7)^2 + (60-72.7)^2 + ... +(76-72.7)^2 + (59-72.7)^2}]

Or, Ïƒ = 10.57

Standard Deviation, in many cases defined as the notation Ïƒ (sigma). The standard Deviation (σ ) indicates how far a 'typical' observation deviates from the data's average or mean value, Î¼.



People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

Online Simulator for ASK, FSK, and PSK

Try our new Digital Signal Processing Simulator!   Start Simulator for binary ASK Modulation Message Bits (e.g. 1,0,1,0) Carrier Frequency (Hz) Sampling Frequency (Hz) Run Simulation Simulator for binary FSK Modulation Input Bits (e.g. 1,0,1,0) Freq for '1' (Hz) Freq for '0' (Hz) Sampling Rate (Hz) Visualize FSK Signal Simulator for BPSK Modulation ...

Constellation Diagrams of ASK, PSK, and FSK

📘 Overview of Energy per Bit (Eb / N0) 🧮 Online Simulator for constellation diagrams of ASK, FSK, and PSK 🧮 Theory behind Constellation Diagrams of ASK, FSK, and PSK 🧮 MATLAB Codes for Constellation Diagrams of ASK, FSK, and PSK 📚 Further Reading 📂 Other Topics on Constellation Diagrams of ASK, PSK, and FSK ... 🧮 Simulator for constellation diagrams of m-ary PSK 🧮 Simulator for constellation diagrams of m-ary QAM BASK (Binary ASK) Modulation: Transmits one of two signals: 0 or -√Eb, where Eb​ is the energy per bit. These signals represent binary 0 and 1.    BFSK (Binary FSK) Modulation: Transmits one of two signals: +√Eb​ ( On the y-axis, the phase shift of 90 degrees with respect to the x-axis, which is also termed phase offset ) or √Eb (on x-axis), where Eb​ is the energy per bit. These signals represent binary 0 and 1.  BPSK (Binary PSK) Modulation: Transmits one of two signals...

MATLAB Code for Rms Delay Spread

RMS delay spread is crucial when you need to know how much the signal is dispersed in time due to multipath propagation, the spread (variance) around the average. In high-data-rate systems like LTE, 5G, or Wi-Fi, even small time dispersions can cause ISI. RMS delay spread is directly related to the amount of ISI in such systems. RMS Delay Spread [↗] Delay Spread Calculator Enter delays (ns) separated by commas: Enter powers (dB) separated by commas: Calculate   The above calculator Converts Power to Linear Scale: It correctly converts the power values from decibels (dB) to a linear scale. Calculates Mean Delay: It accurately computes the mean excess delay, which is the first moment of the power delay profile. Calculates RMS Delay Spread: It correctly calculates the RMS delay spread, defined as the square root of the second central moment of the power delay profile.   MATLAB Code  clc...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...

📘 Overview of BER and SNR 🧮 Online Simulator for BER calculation of m-ary QAM and m-ary PSK 🧮 MATLAB Code for BER calculation of M-ary QAM, M-ary PSK, QPSK, BPSK, ... 📚 Further Reading 📂 View Other Topics on M-ary QAM, M-ary PSK, QPSK ... 🧮 Online Simulator for Constellation Diagram of m-ary QAM 🧮 Online Simulator for Constellation Diagram of m-ary PSK 🧮 MATLAB Code for BER calculation of ASK, FSK, and PSK 🧮 MATLAB Code for BER calculation of Alamouti Scheme 🧮 Different approaches to calculate BER vs SNR What is Bit Error Rate (BER)? The abbreviation BER stands for Bit Error Rate, which indicates how many corrupted bits are received (after the demodulation process) compared to the total number of bits sent in a communication process. BER = (number of bits received in error) / (total number of tran...

LDPC Encoding and Decoding Techniques

📘 Overview & Theory 🧮 LDPC Encoding Techniques 🧮 LDPC Decoding Techniques 📚 Further Reading 'LDPC' is the abbreviation for 'low density parity check'. LDPC code H matrix contains very few amount of 1's and mostly zeroes. LDPC codes are error correcting code. Using LDPC codes, channel capacities that are close to the theoretical Shannon limit can be achieved.  Low density parity check (LDPC) codes are linear error-correcting block code suitable for error correction in a large block sizes transmitted via very noisy channel. Applications requiring highly reliable information transport over bandwidth restrictions in the presence of noise are increasingly using LDPC codes. 1. LDPC Encoding Technique The proper form of H matrix is derived from the given matrix by doing multiple row operations as shown above. In the above, H is parity check matrix and G is generator matrix. If you consider matrix H as [-P' | I] then matrix G will be...

Comparisons among ASK, PSK, and FSK | And the definitions of each

📘 Comparisons among ASK, FSK, and PSK 🧮 Online Simulator for calculating Bandwidth of ASK, FSK, and PSK 🧮 MATLAB Code for BER vs. SNR Analysis of ASK, FSK, and PSK 📚 Further Reading 📂 View Other Topics on Comparisons among ASK, PSK, and FSK ... 🧮 Comparisons of Noise Sensitivity, Bandwidth, Complexity, etc. 🧮 MATLAB Code for Constellation Diagrams of ASK, FSK, and PSK 🧮 Online Simulator for ASK, FSK, and PSK Generation 🧮 Online Simulator for ASK, FSK, and PSK Constellation 🧮 Some Questions and Answers Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK    Comparisons among ASK, PSK, and FSK Comparison among ASK, FSK, and PSK Parameters ASK FSK PSK Variable Characteristics Amplitude Frequency ...

What is a Chirp Signal?

📘 Overview & Theory 🧮 MATLAB Code 📚 Further Reading   Chirp signals are often used to find target objects. In a chirp signal, the frequency varies with time. For up-chirp signals, frequency increases with time. Oppositely, for down-chirp signals, the frequency decreases with time. Advantages of a chirp signal over a single-toned signal Better resolution Better Security The wide bandwidth of a chirp signal allows for capturing more detailed info about the target or object In a chirp signal, pulse compression enhances resolution by concentrating the signal energy into a shorter duration of time It is less susceptible to noise  It improves signal to noise ratio Up-Chirp Signal A sinusoidal up-chirp signal is denoted as Where A is the amplitude of this signal             f0 is the starting frequency of the chirp at t=0             Î± is the chirp rate or the rate at which the frequency incre...

What are the main lobe and side lobes in Beamforming

    What are the main lobe and side lobes in Beamforming? You've probably noticed that in the diagram of  beamforming , there are two types of lobes in beamforming patterns. One is the main lobe, while the others are side lobes. We intend to communicate with receivers with a stronger directional path from the transmitter when we produce beams for wireless communication. We can also see side lobes in this scenario. These side lobes, on the other hand, are not necessary for effective communication. As a result, we take various procedures to remove those side lobes or to reduce the number of side lobes as much as feasible; otherwise, inter-symbol interference  occurs, and signal quality suffers. Figure: Illustration of Main Lobe and Side lobes, where the x-axis denotes the angle of arrival (AOA) and angle of departure (AOD), respectively, while, the y-axis denotes the gain/power in dB (decibel).     In the case of MIMO antennas, our major goal is to reduce int...