Scribd
Upload
27 views

%simulation Length %channel Length %number of Independent Trials

This document presents code for simulating channel equalization using an NLMS algorithm over multiple trials. The code generates training data, applies it to a channel, and adds noise. It then performs equalization using NLMS and averages the results over trials, plotting the final cascade impulse response.

Uploaded by

Oral Robertz
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
27 views

%simulation Length %channel Length %number of Independent Trials

This document presents code for simulating channel equalization using an NLMS algorithm over multiple trials. The code generates training data, applies it to a channel, and adds noise. It then performs equalization using NLMS and averages the results over trials, plotting the final cascade impulse response.

Uploaded by

Oral Robertz
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 2

clear all

close all

%simulation length
N = 1000;

%channel length
M = 5;

%number of independent trials


T = 100;

cascade_impulse_response = zeros(1,2*M-1);
for j = 1 : T
%training signal
u = randn(1,N);

%channel to be equalized
c = randn(M,1);
c = c / norm(c);

%channel output
z = filter(c,1,u);

%additive noise to the channel output


SNR = 30;
var_v = var(z) * 10^(-SNR/10);
v = var_v^0.5 * randn(1,N);

%input to the equalizer


x = z + v;

%NLMS channel equalization


w = zeros(M,1);
x_regressor = zeros(1,M);
step = 0.1;
epsilon = 10^(-6);
for k = 4 : N
x_regressor = [x(k) x_regressor(1:M-1)];
e = u(k-3) - x_regressor * w;
w = w + step * x_regressor' * e / (x_regressor *
x_regressor' + epsilon);
end
cascade_impulse_response = cascade_impulse_response +
conv(w,c)';
disp(j);
end
figure;
stem(cascade_impulse_response/T);
title('cascade channel-equalizer impulse response');
xlabel('taps');

You might also like