mykey1.m

基于变步长LMS在干扰抵消中的应用

clear all;
syorder=8;a=0.97;
r=50/100000;
umax=0.005;
umin=0.0005;
[signal,fs]=wavread('speech');
signal=signal';
n=length(signal);
d=zeros(n+syorder-1,syorder);
wk=zeros(syorder,n);
en=zeros(1,n);
uk=en;
uk(1)=umax;
noise=zeros(1,n);
noise=awgn(noise,5);
noise=[zeros(1,syorder-1),noise];
signal=[zeros(1,syorder-1),signal];
dn=signal+noise;
for i=1:n
    d(i+syorder-1,:)=noise(i+syorder-1:-1:i);
end
for i=1:(n-1)
    y(i)=d(i+syorder-1,:)*wk(:,i);
   en(i)=dn(i+syorder-1)-y(i);
   wk(:,i+1)=wk(:,i)+2*uk(i)*en(i)*d(i+syorder-1,:)';
   uk(i+1)=a*uk(i)+r*[abs(en(i)).^2];
   if uk(i+1)>umax
       uk(i+1)=umax;
   elseif uk(i+1)<umin
       uk(i+1)=umin;
   else
       uk(i+1)=uk(i+1);
   end
end