Skip to main content

MATLAB Code for Channel Impulse Response


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

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

People are good at skipping over material they already know!

View Related Topics to







Contact Us

Name

Email *

Message *

Popular Posts

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

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

Fading : Slow & Fast and Large & Small Scale Fading (with MATLAB Code + Simulator)

📘 Overview 📘 LARGE SCALE FADING 📘 SMALL SCALE FADING 📘 SLOW FADING 📘 FAST FADING 🧮 MATLAB Codes 📚 Further Reading LARGE SCALE FADING The term 'Large scale fading' is used to describe variations in received signal power over a long distance, usually just considering shadowing.  Assume that a transmitter (say, a cell tower) and a receiver  (say, your smartphone) are in constant communication. Take into account the fact that you are in a moving vehicle. An obstacle, such as a tall building, comes between your cell tower and your vehicle's line of sight (LOS) path. Then you'll notice a decline in the power of your received signal on the spectrogram. Large-scale fading is the term for this type of phenomenon. SMALL SCALE FADING  Small scale fading is a term that describes rapid fluctuations in the received signal power on a small time scale. This includes multipath propagation effects as well as movement-induced Doppler fr...

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

Theoretical BER vs SNR for binary ASK, FSK, and PSK with MATLAB Code + Simulator

📘 Overview & Theory 🧮 MATLAB Codes 📚 Further Reading Theoretical BER vs SNR for Amplitude Shift Keying (ASK) The theoretical Bit Error Rate (BER) for binary ASK depends on how binary bits are mapped to signal amplitudes. For typical cases: If bits are mapped to 1 and -1, the BER is: BER = Q(√(2 × SNR)) If bits are mapped to 0 and 1, the BER becomes: BER = Q(√(SNR / 2)) Where: Q(x) is the Q-function: Q(x) = 0.5 × erfc(x / √2) SNR : Signal-to-Noise Ratio N₀ : Noise Power Spectral Density Understanding the Q-Function and BER for ASK Bit '0' transmits noise only Bit '1' transmits signal (1 + noise) Receiver decision threshold is 0.5 BER is given by: P b = Q(0.5 / σ) , where σ = √(N₀ / 2) Using SNR = (0.5)² / N₀, we get: BER = Q(√(SNR / 2)) Theoretical BER vs ...

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

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 carriers adding up constructively) Low (less fluctuation in amplitude) Why PAPR is High Subcarriers can add in phase, causing spikes DFT "pre-spreads" data, smoothing it Used in Wi-Fi, LTE downlink LTE uplink (as SC-FDMA) In OFDM, all subcarriers can...

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

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

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

📘 Overview 🧮 Simulator for calculating BER 🧮 MATLAB Codes for calculating theoretical BER 🧮 MATLAB Codes for calculating simulated BER 📚 Further Reading 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. Simulator for calculating BER vs SNR for binary ASK, FSK, and PSK Calculate BER for Binary ASK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary FSK Modulation Enter SNR (dB): Calculate BER Calculate BER for Binary PSK Modulation Enter SNR (dB): Calculate BER BER vs. SNR Curves MATLAB Code for Theoretical BER % The code is written by SalimWireless.Com clc; clear; close all; % SNR va...

MATLAB Code for ASK, FSK, and PSK (with Online Simulator)

📘 Overview & Theory 🧮 MATLAB Code for ASK 🧮 MATLAB Code for FSK 🧮 MATLAB Code for PSK 🧮 Simulator for binary ASK, FSK, and PSK Modulations 📚 Further Reading ASK, FSK & PSK HomePage MATLAB Code MATLAB Code for ASK Modulation and Demodulation % The code is written by SalimWireless.Com % Clear previous data and plots clc; clear all; close all; % Parameters Tb = 1; % Bit duration (s) fc = 10; % Carrier frequency (Hz) N_bits = 10; % Number of bits Fs = 100 * fc; % Sampling frequency (ensure at least 2*fc, more for better representation) Ts = 1/Fs; % Sampling interval samples_per_bit = Fs * Tb; % Number of samples per bit duration % Generate random binary data rng(10); % Set random seed for reproducibility binary_data = randi([0, 1], 1, N_bits); % Generate random binary data (0 or 1) % Initialize arrays for continuous signals t_overall = 0:Ts:(N_bits...