how can i write a matlab code for lms algorithm for audio signal?

lipsa tripathy

2017 2 26

0 个回答

2 次查看（30 天）

0 个投票

clc 
clear all
close all
[x,Fs] = audioread('test11.wav');
sound(x,Fs);
pause(10);
delay1 = 0.6; % delay
delay2=0.9;
delay3=0.6;
alpha1 = 0.65;% echo strength
alpha2=.40;
alpha3=.5; 
t=0:1/Fs:(length(x)-1)/Fs;
subplot(3,2,1)
plot(t,x)
title('original signal')
D1 = delay1*Fs; 
D2=delay2*Fs;
D3=delay3*Fs; 
y = zeros(size(x));
y(1:D1) = x(1:D1);
snr=4;
for i=D1+40000:length(x) 
    y(i) = x(i)+alpha1*x(i-D1)+alpha2*x(i-D2)+alpha3*x(i-D3);
 end 
   sound(y,Fs);
   t1=0:1/Fs:(length(y)-1)/Fs;
   subplot(3,2,2)
   plot(t1,y)
   title('original+echo');
   % ----------------------------------------------------------------- 
   load nearspeech
   load farspeech
   N=100000;
   L=1000;
   mu=0.00003;
   h1=fir1(2*L-1,0.5); 
   h=10*h1(L+1:2*L);
   % ---------------------------------------------------------------------
   wlms=zeros(L,1); 
   xlms=zeros(L,1); 
   for m=1:N
    xlms=[ y(m);xlms(1:L-1) ];
    ylms= h * xlms;
    elms(m)= d(m) - ylms;
    wlms= wlms + (mu * xlms * elms(m));
   end
   subplot(3,2,4)
   plot(elms,'m')
   title('Recovered Signal')

here only the first element of y is passing not the whole y signal.how can i make the whole y signal pass through the filter without any dimension missmatch?