untitled.m

lms算法的程序可以直接运行

clear 
clc
t=0:1:8000-1;%采样频率fs=1000,信号频率800
s=5*sin(2*pi*t/20)+5*sin(pi*t/20)+5*cos(2*pi*t/30);
%d=-2*log(s);
n=5*rand(size(t));%n1=-20*log(n);
x=n;
wq=fft(x);
%plot(abs(wq))
%w=[0 0.5 ];
k=40;j=size(t,2);
c=1:k;
w=zeros(1,k);
u=0.0001;

for i=1:j-k;
    y(i)=w*x(i:i+k-1)';
    e(i)=x(i+k)-y(i);
    w=w+2*u*e(i)*x(i:i+k-1);
end 
%x=-20*log(n);
%figure(2)
subplot(2,1,1)
plot(t,x);

xlabel('噪声信号样本数');
ylabel('输入噪声幅度')
z=1:size(e,2);
%e=-20*log(e);
subplot(2,1,2);
%figure(2)
plot(z,e);
xlabel('迭代次数');
ylabel('误差信号幅度')
% % figure(2);
% % plot(c,w)
% % subplot(3,1,3);
% % 
% plot(z,y);
% title('滤波结果');