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

MATLAB Code for BER performance of QPSK with BPSK, 4-QAM, 16-QAM, 64-QAM, 256-QAM, etc

📘 Overview 🧮 MATLAB Codes 🧮 Online Simulator for Calculating BER of M-ary PSK and QAM 🧮 QPSK vs BPSK and QAM: A Comparison of Modulation Schemes in Wireless Communication 🧮 Are QPSK and 4-PSK same? 📚 Further Reading   QPSK offers double the data rate of BPSK while maintaining a similar bit error rate at low SNR when Gray coding is used. It shares spectral efficiency with 4-QAM and can outperform 4-QAM or 16-QAM in very noisy channels. QPSK is widely used in practical wireless systems, often alongside QAM in adaptive modulation schemes [Read more...] What is the Gray Code? Gray Code: Gray code is a binary numeral system where two successive values differ in only one bit. This property is called the single-bit difference or unit distance code. It is also known as reflected binary code. Let's convert binary 111 to Gray code: Binary bits: B = 1 1 1 Apply the rule: G[0] = B[0] = 1...

FFT Butterfly Method Explained (with Simulations)

4-Point FFT Using Butterfly Method Given: x[n] = {0, 1, 2, 3} Step 1: Split into Even & Odd Even indices: x e = {x[0], x[2]} = {0, 2} Odd indices: x o = {x[1], x[3]} = {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 4 k O[k] X[k + 2] = E[k] - W 4 k O[k] Twiddle Factors (N=4): W 4 0 = 1, W 4 1 = -j Final Calculations: X[0] = E[0] + W 4 0 O[0] = 2 + (1)(4) = 6 X[2] = E[0] - W 4 0 O[0] = 2 - (1)(4) = -2 X[1] = E[1] + W 4 1 O[1] = -2 + (-j)(-2) = -2 + 2j X[3] = E[1] - W 4 1 O[1] = -2 - (-j)(-2) = -2 - 2j Final Answer: X[k] = {6, -2 + 2j, -2, -2 - 2j} 8-Point FFT Using Butterfly Method Given: x[n] = {0,1,2,3,4,5,6,7} Step 1: Split into Bit-Reversed Order To perform DIT-FFT, split the 8 points into pairs of two: Group A: {x[0], x[4]} = {0, 4}...

MATLAB Code for QPSK Modulation and Demodulation

📘 Overview 🧮 MATLAB Codes 🧮 Theory 🧮 BER performance of QPSK with BPSK, 4-QAM, 16-QAM, 64-QAM, 256-QAM, etc 📚 Further Reading QPSK Passband Signal Generation Spectral Efficiency in QPSK   Quadrature Phase Shift Keying (QPSK) is a digital modulation scheme that conveys two bits per symbol by changing the phase of the carrier signal. Each pair of bits is mapped to one of four possible phase shifts: 0°, 90°, 180°, or 270° 00  ===> 0 degree phase shift of carrier signal 01  ===> 90 degree 11  ===> 180 degree 10  ===> 270 degree   MATLAB Script clc; clear all; close all; clc; M = 4; data = randi([0 (M-1)], 1000, 1); Phase = 0; modData=pskmod(data,M,Phase); figure(1); scatterplot(modData); channelAWGN = 15; rxData2 = awgn(modData, channelAWGN); figure(2); scatterplot(rxData2); demodData = pskdemod(rxData2,M,Phase);   Result data 1 0 2 2 0 2 1 . . . modData -1.0...

Online Simulator for ASK, FSK, and PSK Signal Generation

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

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

OFDM Symbols and Subcarriers Explained

This article explains how OFDM (Orthogonal Frequency Division Multiplexing) symbols and subcarriers work. It covers modulation, mapping symbols to subcarriers, subcarrier frequency spacing, IFFT synthesis, cyclic prefix, and transmission. Step 1: Modulation First, modulate the input bitstream. For example, with 16-QAM , each group of 4 bits maps to one QAM symbol. Suppose we generate a sequence of QAM symbols: s0, s1, s2, s3, s4, s5, …, s63 Step 2: Mapping Symbols to Subcarriers Assume N sub = 8 subcarriers. Each OFDM symbol in the frequency domain contains 8 QAM symbols (one per subcarrier): Mapping (example) OFDM symbol 1 → s0, s1, s2, s3, s4, s5, s6, s7 OFDM symbol 2 → s8, s9, s10, s11, s12, s13, s14, s15 … OFDM sym...

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

DFTs-OFDM vs OFDM: Why DFT-Spread OFDM Reduces PAPR Effectively (with MATLAB Code)

Understanding PAPR in DFT-spread OFDM vs. Standard OFDM In modern wireless communications like 4G LTE and 5G NR, managing the Peak-to-Average Power Ratio (PAPR) is critical for hardware efficiency. While OFDM is the gold standard for high-speed data, its high PAPR poses significant challenges for mobile devices. This is where DFTs-OFDM (also known as SC-FDMA) comes in. DFT-spread OFDM (DFTs-OFDM) has lower Peak-to-Average Power Ratio (PAPR) because it "spreads" the data in the frequency domain before applying IFFT, making the time-domain signal behave more like a single-carrier signal rather than a multi-carrier one like OFDM. Deeper Explanation: Aspect OFDM DFTs-OFDM Signal Type Multi-carrier Single-carrier-like Process IFFT of QAM directly QAM → DFT → IFFT PAPR Level High (due to many...