Welcome to EDAboard.com

Welcome to our site! EDAboard.com is an international Electronic Discussion Forum focused on EDA software, circuits, schematics, books, theory, papers, asic, pld, 8051, DSP, Network, RF, Analog Design, PCB, Service Manuals... and a whole lot more! To participate you need to register. Registration is free. Click here to register now.

matlab code help needed

Status
Not open for further replies.

jenabi

Newbie level 4
Joined
Oct 1, 2010
Messages
6
Helped
0
Reputation
0
Reaction score
0
Trophy points
1,281
Activity points
1,334
hii
i need help below is the code that is for dynamically changing the tap length of adaptive filter....but when i run this in matlab i get error saying that d(n:-1:n-L+1) cannot have negative indices...can someone tell me how to debug this...and where am i going wrong

clc;
clear;

delta = 2;
n_bits=15;
t=5000;

x = 0.5*((rand(1,t)>0.5)-.5);

r = 1;
imp = [-4.4 0 -7 -12 -9 -14 -20 -16.5 -20 -22.5 -20];
response = (10.^(imp/10));
quant = round(response*(2^n_bits));
quant = quant/2^n_bits;

channel = [0.4 1 0.2 0 0.06 0.1 0.04 0.01 0.02 0.01 0.005 0.01];
channel = channel.*2;

b = (randn(1,5)).*exp(-(0:4));
c = [1 0 -.3 0 .5 0 -.1 0 0.1];

% Generate noise
snr = 20;
n=1/sqrt(2).*(randn(1,length(x)));
noise = 10^(-snr/20)*n;
noise = 0.25.*noise;

var = 1;
noise_low_SN=0.55*(randn(1,t*r)*sqrt(var));
noise_low_q=round(noise_low_SN*2^n_bits)/2^n_bits;
noise_high_SN=0.01*(randn(1,t*r)*sqrt(var));
noise_high_q=round(noise_high_SN*2^n_bits)/2^n_bits;
noise = [noise_high_q];

% input signal
input = (filter(c,1,x));
%input = (x);

pad = zeros(1,delta);
pad1 = ones(1,4);
input_LMS = round(input*2^n_bits)/2^n_bits;
input_q = round(input*2^n_bits)/2^n_bits; %Quantisation for file output
input_des=round(x*2^n_bits)/2^n_bits;

%adaptive filtering
L = 9; % filter length
P = 3; % sub-filter size
mu = 0.005; % step size ( mu < 1/(L*std(x)^2) )
W = zeros(1,L); % adaptive coefficients are initialised to zero

ASE = 50;

e_array = zeros(1,ASE);
e_sub_array = zeros(1,50);

pad = zeros(1,L);
x_pad = zeros(1,L-3);
x_delay = [pad x];

d = [x_pad, input_LMS];
y = zeros(1,length(x));

e = zeros(1,length(x));
y_sub = zeros(1,length(x));
e_sub = zeros(1,length(x));
coeff = zeros(1,length(x));

for n = L:length(x)
n
L
X = d(n:-1: n-L+1);
y(n) = round((W*X')*2^n_bits)/2^n_bits; % adaptive filter out

X_sub = d(n:-1:n-L+(P+1));
W_sub = W(1:1:(L-P));
y_sub(n) = round((W_sub*X_sub')*2^n_bits)/2^n_bits; % adaptive sub-filter out

% error signals
e(n) = (round((x_delay(n) - y(n))*2^n_bits))/2^n_bits;
e_sub(n) = (round((x_delay(n) - y_sub(n))*2^n_bits))/2^n_bits;
e_array(2:ASE) = e_array(1:ASE-1); % calculation of ASE
e_array(1) = (e(n))^2;
e_sub_array(2:ASE) = e_sub_array(1:ASE-1);
e_sub_array(1) = (e_sub(n))^2;

ASE_res = sum(e_array);
ASE_res_sub = sum(e_sub_array);

if (ASE_res <= ASE_res_sub) % evaluate Filter length
L = L + P; % to add or remove P taps

elseif (ASE_res >= ASE_res_sub)
L = L - P;
end
coeff(n) = (round(2*mu*e(n)*2^n_bits)/2^n_bits);
W = W + coeff(n)*X; % LMS update
W=(round(W*2^n_bits)/2^n_bits);
e_sq(n) = e(n)^2;
end;


thanks
 

Status
Not open for further replies.

Part and Inventory Search

Welcome to EDABoard.com

Sponsor

Top