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

Q-function in BER vs SNR Calculation (with Simulation)

Q-function in BER vs. SNR Calculation In digital communications and signal processing, the Q-function plays a significant role in predicting system reliability. It allows engineers to quantify the probability that Gaussian noise will exceed a specific threshold, causing a bit error. What is the Q-function? The Q-function is a mathematical function representing the tail probability of the standard normal (Gaussian) distribution. It is the complementary cumulative distribution function (CCDF) of a standard Gaussian distribution. Q(x) = (1 / √(2Ï€)) ∫â‚“∞ e^(-t² / 2) dt Q-Function Interactive Simulator Move the slider to see how the "Tail Probability" (the area in red) changes. This area represents the Probability of Error (BER) . Threshold Distance ( x ) — (Simulates Increasing SNR) x = 1.0 Q(x) = 0.1587 ...

Design of CMOS Flip-Flops (SR, D, JK)

Design of CMOS Flip-Flops (SR, D, JK) A flip-flop or latch is a circuit with two stable states, used to store state information. It is the basic storage element in sequential logic and a fundamental building block in digital electronics systems, including computers and communication devices. Flip-flops and latches act as data storage elements for states, pulse counting, and synchronization of variably-timed input signals to a reference clock. Flip-flops can be transparent/opaque (latches) or clocked (synchronous, edge-triggered). Latches are level-sensitive, while flip-flops are edge-sensitive. In sequential logic, the output depends on current inputs and previous states. Fig.1 shows a sequential circuit combining a combinational block and a memory element. ...

Pulse Width Modulation (PWM)

Pulse-width modulation (PWM), or pulse-duration modulation (PDM), is a method of controlling the average power delivered by an electrical signal.   Fig: An example of PWM in an idealized inductor driven by a blue line voltage source modulated as a series of sawtooth pulses, resulting in a red line current in the inductor.    Generating a PWM Signal The simplest way to generate a PWM signal is the intersection method, which requires only a sawtooth or a triangle waveform (easily generated using a simple oscillator) and a comparator. When the value of the reference signal is more than the modulation waveform, the PWM signal (magenta) is in the high state; otherwise, it is in the low state.      Duty cycle A low duty cycle equates to low power because the power is off for most of the time; the word duty cycle reflects the ratio of "on" time to the regular interval or "period" of time. The duty cycle is measured in percent, with 100% representing full o...

BER vs SNR for M-ary QAM, M-ary PSK, QPSK, BPSK, ...(MATLAB Code + Simulator)

Bit Error Rate (BER) & SNR Guide Analyze communication system performance with our interactive simulators and MATLAB tools. 📘 Theory 🧮 Simulators 💻 MATLAB Code 📚 Resources BER Definition SNR Formula BER Calculator MATLAB Comparison 📂 Explore M-ary QAM, PSK, and QPSK Topics ▼ 🧮 Constellation Simulator: M-ary QAM 🧮 Constellation Simulator: M-ary PSK 🧮 BER calculation for ASK, FSK, and PSK 🧮 Approaches to BER vs SNR What is Bit Error Rate (BER)? The BER indicates how many corrupted bits are received compared to the total number of bits sent. It is the primary figure of merit f...

FFT Butterfly Method Explained (with Example of 4-point DFT)

  FFT Using Butterfly Method Given: x[n] = {0, 1, 2, 3} Step 1: Split into Even & Odd Even indices: x e = {0, 2} Odd indices: x o = {1, 3} Step 2: 2-point DFT For any {a, b}: DFT = {a + b, a - b} Even Part: E = {0+2, 0-2} = {2, -2} Odd Part: O = {1+3, 1-3} = {4, -2} Step 3: Combine Using Butterfly X[k] = E[k] + W k O[k] X[k + N/2] = E[k] - W k O[k] For N = 4: W 0 = 1 W 1 = -j Final Calculations X[0] = 2 + 4 = 6 X[2] = 2 - 4 = -2 X[1] = -2 + (-j)(-2) = -2 + 2j X[3] = -2 - (-j)(-2) = -2 - 2j Final Answer: X[k] = {6, -2 + 2j, -2, -2 - 2j} Try Interactive Online Simulations Interactive FFT Online Simulator (For understanding Fundamentals)  Interactive FFT Online Simulator (Analyze .CSV, .MP3, .MP4, etc. Further Reading Fourier Transform OFDM Return to Fourier Transform Main Page →

Frequency Shift Keying (FSK) Modulation & Demodulation (with Simulation)

Frequency Shift Keying (FSK) Theoretical Foundations: Frequency Shift Keying (FSK) is a discrete frequency modulation scheme wherein the digital information is encoded via instantaneous shifts in the carrier signal's frequency. The fundamental implementation is Binary FSK (BFSK), which maps binary data onto two distinct, discrete spectral states. A binary '1' (the "mark" state) is represented by a carrier frequency \( f_1 \), while a binary '0' (the "space" state) corresponds to frequency \( f_2 \). Each symbol is sustained for a bit interval denoted by \( T_b \). FSK Transmitter Characterization: The mathematical model for the modulated BFSK output \( s(t) \) is defined as: \[ s(t) = \begin{cases} A_c \cos(2\pi f_1 t), & \text{for } m = 1 \\ A_c \cos(2\pi f_2 t), & \text{for } m = 0 \end{cases} \] ...

AM Modulation Online Simulator

Amplitude Modulation Simulator s AM (t) = A c [1 + k a m(t)] cos(ω c t) where, ω = 2πf & k a = Amplitude Sensitivity Modulation index, μ = k a A m Message Frequency (fm): Carrier Frequency (fc): Carrier Amplitude (Ac): Modulation Index (m = Am / Ac):

Online Simulator for ASK, FSK, and PSK

Interactive Digital Signal Processing (DSP) Tutorial and Simulator for ASK, FSK, and BPSK modulation techniques. Try our new Digital Signal Processing Simulator!   •   Interactive ASK, FSK, and BPSK tools updated for 2025. Start Now Digital Modulation Visualizer: ASK, FSK, & BPSK Simulator Learn and visualize binary modulation techniques (ASK, FSK, BPSK) in real-time with adjustable carrier and sampling parameters. Perfect for DSP students and engineers. 📡 ASK Simulator 📶 FSK Simulator 🎚️ BPSK Simulator 📚 More Topics ASK Modulator FSK Modulator BPSK Modulator More Topics 1. ASK (Amplitude Shift Keying) Simulat...