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, Î¼.





Contact Us

Name

Email *

Message *

Popular Posts

Hybrid Beamforming | Page 1

Beamforming Techniques Hybrid Beamforming... Page 1 | Page 2 | Hybrid Beamforming: Hybrid beam formation was developed to address some of the limitations of digital pre-coding approaches. Every antenna element is connected to an RF chain in digital pre-coding (beam forming) method. We also know that each RF chain is in charge of providing a separate data stream between the transmitter and the receiver. We know that a larger number of independent data streams leads to higher data rates. It has a spatial multiplexing feature for MIMO. As a result, we may assume that switching from MIMO to massive MIMO will benefit us more in terms of spatial multiplexing in massive MIMO, where each antenna is coupled to a single RF chain. We'll proceed with a definition of hybrid beam forming. Overview of hybrid beam forming with example: Unlike digital beam forming, more than one antenna element is connected to a single RF chain in hybr...

UGC NET Electronic Science Previous Year Question Papers with Solutions

Home / Engineering & Other Exams / UGC NET 2026 PYQ ⬇️ Download Papers and Solutions 📋 Exam Pattern 💡 Preparation Tips ❓ FAQs 📊 Exam Highlights: Electronic Science (88) Feature Details Junior Research Fellowship (JRF) ₹37,000 + HRA per month Eligibility M.Sc/M.Tech in Electronics (55%) Validity of Certificate JRF (3 Years) | Lectureship (Lifetime) 📥 Download UGC NET Electronics PDFs Complete collection of previous year question papers, answer keys and explanations for Subject Code 88. Start Downloading 📂 View All Question Papers June 2025 - Question Paper Download PDF June 2025 - Solved Paper + Explanation ...

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...

Comparisons among ASK, PSK, and FSK (with MATLAB + Simulator)

Modulation ASK, FSK & PSK Constellation MATLAB Simulink MATLAB Code Comparisons among ASK, PSK, and FSK 📘 Comparisons among ASK, FSK, and PSK 🧮 Online Simulator Bandwidth 🧮 MATLAB Code BER Analysis 📚 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 Comparisons among ASK, PSK, and FSK Comparison among ASK, FSK, and PSK Parameters ASK FSK PSK Variable Characteristics Amplitude ...

How Windowing Affects Your Periodogram

The windowed periodogram is a widely used technique for estimating the Power Spectral Density (PSD) of a signal. It enhances the classical periodogram by mitigating spectral leakage through the application of a windowing function. This technique is essential in signal processing for accurate frequency-domain analysis.   Power Spectral Density (PSD) The PSD characterizes how the power of a signal is distributed across different frequency components. For a discrete-time signal, the PSD is defined as the Fourier Transform of the signal’s autocorrelation function: S x (f) = FT{R x (Ï„)} Here, R x (Ï„)}is the autocorrelation function. FT : Fourier Transform   Classical Periodogram The periodogram is a non-parametric PSD estimation method based on the Discrete Fourier Transform (DFT): P x (f) = \(\frac{1}{N}\) X(f) 2 Here: X(f): DFT of the signal x(n) N: Signal length However, the classical periodogram suffers from spectral leakage due to abrupt truncation of the ...

What is - 3dB Frequency Response? And applications

📘 Overview & Theory 📘 Application of -3dB Frequency Response 🧮 MATLAB Codes 🧮 Online Digital Filter Simulator 📚 Further Reading Filters What is -3dB Frequency Response?   Remember, for most passband filters, the magnitude response typically remains close to the peak value within the passband, varying by no more than 3 dB. This is a standard characteristic in filter design. The term '-3dB frequency response' indicates that power has decreased to 50% of its maximum or that signal voltage has reduced to 0.707 of its peak value. Specifically, The -3dB comes from either 10 Log (0.5) {in the case of power} or 20 Log (0.707) {in the case of amplitude} . Viewing the signal in the frequency domain is helpful. In electronic amplifiers, the -3 dB limit is commonly used to define the passband. It shows whether the signal remains approximately flat across the passband. For example, in pulse shapi...

Amplitude Shift Keying (ASK) Modulation & Demodulation (with Simulation)

Amplitude Shift Keying (ASK): Signal Analysis and Characterization Theoretical Overview: Amplitude Shift Keying (ASK) represents a primary digital modulation technique wherein information is encoded through discrete variations in the carrier signal's instantaneous amplitude. In a Binary ASK (BASK) framework, the modulation process maps binary data onto two distinct amplitude levels. Specifically, the binary '1' (mark) is conveyed by a sinusoidal carrier with amplitude A c and frequency f c over a bit interval T b , while the binary '0' (space) is represented by a null signal state. This particular signaling method is widely recognized as On-Off Keying (OOK) . It is technically realized by gating a carrier oscillator with a unipolar baseband sequence, effectively performing a product modulation that shifts the baseband spectrum to the carrier frequency. ASK Transmitter Architecture: ...