Skip to main content

MATLAB Code for Channel Impulse Response (with Simulator)


MATLAB Code for Channel Impulse Response (CIR)

MATLAB Script for Simulating CIR

This MATLAB script allows you to generate and visualize the channel impulse response (CIR). You can choose to create a 'random' multi-path channel or a near-'ideal' single-path channel to understand their distinct characteristics.


% User input for choosing the type of impulse response
response_type = input('Enter "random" for random channel impulse response or "ideal" for near-ideal impulse response: ', 's');

if strcmpi(response_type, 'random')
    % Parameters for random impulse response
    num_taps = input('Enter the number of taps: '); % Number of taps in the channel
    delay_spread = input('Enter the maximum delay spread in samples: '); % Maximum delay spread in samples
    channel_gain = input('Enter the overall channel gain: '); % Overall channel gain

    % Generate random tap delays
    tap_delays = randi(delay_spread, 1, num_taps);

    % Generate random complex gains for each tap
    tap_gains = (rand(1, num_taps) + 1i * rand(1, num_taps)) * channel_gain;

    % Generate impulse response
    channel_impulse_response = zeros(1, max(tap_delays) + 1);
    for i = 1:num_taps
        channel_impulse_response(tap_delays(i) + 1) = tap_gains(i);
    end
elseif strcmpi(response_type, 'ideal')
    % Parameters for near-ideal impulse response
    num_taps = 1; % Number of taps in the channel
    channel_gain = input('Enter the overall channel gain: '); % Overall channel gain

    % Generate impulse response
    channel_impulse_response = zeros(1, num_taps);
    channel_impulse_response(1) = channel_gain;
else
    error('Invalid input. Please enter either "random" or "ideal"');
end

% Plot impulse response
stem(0:length(channel_impulse_response)-1, abs(channel_impulse_response), 'filled');
xlabel('Time (samples)');
ylabel('Magnitude');
if strcmpi(response_type, 'random')
    title('Random Channel Impulse Response');
else
    title('Near-Ideal Channel Impulse Response');
end

Output Examples

Random Channel Impulse Response

Example Input:

Enter "random" for random channel impulse response or "ideal" for near-ideal impulse response: random
Enter the number of taps: 3
Enter the maximum delay spread in samples: 3
Enter the overall channel gain: 0.5
Plot of a randomly generated channel impulse response in MATLAB
Fig: Channel Impulse Response (Random Generation)

Ideal Channel Impulse Response

Example Input:

Enter "random" for random channel impulse response or "ideal" for near-ideal impulse response: ideal
Enter the overall channel gain: 0.8
Plot of an ideal channel impulse response in MATLAB showing a single path
Fig: Channel Impulse Response (Ideal Generation)

How to mitigate Channel Distortion caused by Multi-paths?

To mitigate channel distortion caused by multipath in wireless communication is crucial for ensuring reliable and high-quality signal transmission. Multipath distortion occurs when a transmitted signal takes multiple paths to reach the receiver, causing interference and signal degradation. Here are several techniques to mitigate this issue, including Equalization, OFDM, and Channel Coding.

Using an Adaptive Equalizer

An adaptive equalizer is a digital filter that can adjust its coefficients automatically to compensate for channel distortion. It is a powerful tool for mitigating the effects of multipath fading.

Block diagram of an adaptive equalizer used to mitigate channel distortion

Interactive Wireless Channel Simulator

Visualize how multipath interference shapes your signal.

1 Define Input Signal \(x[n]\)

The Unit Impulse is used to "probe" the channel. The output will show exactly how the channel behaves.

2 Design the Channel \(h[n]\)

Clean (100 dB)
\( h(t) = \sum a_i e^{j\theta_i} \delta(t-\tau_i) \)

3 Receiver Output \(y[n] = x[n] * h[n]\)

How it works: Each path in the channel creates a delayed and scaled version of the input signal. The receiver sees the sum of all these versions.

Advanced Channel Impulse Response Simulator

Effect of CIR / multi-paths on BPSK










Further Reading

  1. Channel Impulse Response (CIR)
  2. FFT Based Channel Estimation
  3. Impulse Response of an ARMA System in MATLAB
  4. Channel Matrix Gain

Contact Us

Name

Email *

Message *

Popular Posts

Rayleigh vs Rician Fading (with MATLAB + Simulator)

  In Rayleigh fading , the channel coefficients tend to have a Rayleigh distribution, which is characterized by a random phase and magnitude with an exponential distribution. This means the magnitude of the channel coefficient follows an exponential distribution with a mean of 1. In Rician fading , there is a dominant line-of-sight component in addition to the scattered components. The channel coefficients in Rician fading can indeed tend towards 1, especially when the line-of-sight component is strong. When the line-of-sight component dominates, the Rician fading channel behaves more deterministically, and the channel coefficients may tend towards the value of the line-of-sight component, which could be close to 1.   MATLAB Script clc; clear all; close all; % Define parameters numSamples = 1000; % Number of samples K_factor = 5; % K-factor for Rician fading SNR_dB = 20; % Signal-to-noise ratio (in dB) % Generate complex Gaussian random variable for Rayleigh fading channel h_r...

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

Theoretical vs. simulated BER vs. SNR for ASK, FSK, and PSK (MATLAB Code + Simulator)

📘 Overview 🧮 Simulator 💻 Theoretical Code 📊 Simulated Code 📚 Resources Overview BER vs. SNR denotes how many bits in error are received for a given signal-to-noise ratio, typically measured in dB. Common noise types in wireless systems: 🚀 1. Additive White Gaussian Noise (AWGN) 🌊 2. Rayleigh Fading AWGN adds random noise; Rayleigh fading attenuates the signal variably. A good SNR helps reduce these effects. Bit Error Rate (BER) Equations BER formulas for ASK, FSK, and PSK modulation schemes. ASK BER = 0.5 × erfc(0.5 × √SNR) FSK BER = 0.5 × erfc(√(SNR / 2)) PSK BER = 0.5 × erfc(√SNR) erfc / Q-function (Click here) Live BER S...

UGC NET Electronic Science Previous Year Question Papers

Home / Engineering & Other Exams / UGC NET 2022 PYQ 📥 Download UGC NET Electronics PDFs Complete collection of previous year question papers, answer keys and explanations for Subject Code 88. Start Downloading UGC-NET (Electronics Science, Subject code: 88) Subject_Code : 88; Department : Electronic Science; 📂 View All Question Papers Q. UGC Net Electronic Science Question Paper [June 2025] A. UGC Net Electronic Science Question Paper With Answer Key Download Pdf [June 2025] with full explanation Q. UGC Net Electronic Science Question Paper [December 2024] A. UGC Net Electronic Science Question Paper With Answer Key Download Pdf [December 2024] Q. UGC Net Electronic Science Question Paper [Aug 2024] A. UGC Net Electronic Scien...

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

MATLAB code for BER vs SNR for M-QAM, M-PSK, QPSk, BPSK, ...(with Online Simulator)

🧮 MATLAB Code for BPSK, M-ary PSK, and M-ary QAM Together 🧮 MATLAB Code for M-ary QAM 🧮 MATLAB Code for M-ary PSK 📚 Further Reading MATLAB Script for BER vs. SNR for M-QAM, M-PSK, QPSK, BPSK % Written by Salim Wireless clc; clear; close all; snr_db = -5:2:25; psk_orders = [2, 4, 8, 16, 32]; qam_orders = [4, 16, 64, 256]; ber_psk_results = zeros(length(psk_orders), length(snr_db)); ber_qam_results = zeros(length(qam_orders), length(snr_db)); for i = 1:length(psk_orders) ber_psk_results(i, :) = berawgn(snr_db, 'psk', psk_orders(i), 'nondiff'); end for i = 1:length(qam_orders) ber_qam_results(i, :) = berawgn(snr_db, 'qam', qam_orders(i)); end figure; semilogy(snr_db, ber_psk_results(1, :), 'o-', 'LineWidth', 1.5, 'DisplayName', 'BPSK'); hold on; for i = 2:length(psk_orders) semilogy(snr_db, ber_psk_results(i, :), 'o-', 'DisplayName', sprintf('%d-PSK', psk_orde...

Constellation Diagrams of ASK, PSK, and FSK (with MATLAB Code + Simulator)

Constellation Diagrams: ASK, FSK, and PSK Comprehensive guide to signal space representation, including interactive simulators and MATLAB implementations. 📘 Overview 🧮 Simulator ⚖️ Theory 📚 Resources Definitions Constellation Tool Key Points MATLAB Code 📂 Other Topics: M-ary PSK & QAM Diagrams ▼ 🧮 Simulator for M-ary PSK Constellation 🧮 Simulator for M-ary QAM Constellation 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...

MATLAB Code for Zero-Forcing (ZF) Beamforming in 4×4 MIMO Systems

MATLAB Code for Zero-Forcing (ZF) Beamforming in 4×4 MIMO Systems clc; clear; close all; %% Parameters Nt = 4; % Transmit antennas Nr = 4; % Receive antennas (must be >= Nt for ZFBF) numBits = 1e4; % Number of bits per stream SNRdB = 0; % SNR in dB numRuns = 100; % Number of independent runs for averaging %% Precompute noise standard deviation noiseSigma = 10^(-SNRdB / 20); %% Accumulator for total errors totalErrors = 0; for run = 1:numRuns % Generate random bits: [4 x 10000] bits = randi([0 1], Nt, numBits); % BPSK modulation: 0 → +1, 1 → -1 txSymbols = 1 - 2 * bits; % Rayleigh channel matrix: [4 x 4] H = (randn(Nr, Nt) + 1j * randn(Nr, Nt)) / sqrt(2); %% === Zero Forcing Beamforming at Transmitter === W_zf = pinv(H); % Precoding matrix: [Nt x Nr] txPrecoded = W_zf * txSymbols; % Apply ZF precoding % Normalize transmit power (optional but useful) txPrecoded = txPrecoded / sqrt(mean(abs(txPrecoded(:)).^2)); %% Channel transmission with AWGN noise = noiseSigma * (randn(...