xianboqi.m

运用MATLAB设计一个LMS自适应滤波器

%陷波器实验 
%输入信号：x  正弦信号：s   噪声信号： v  输出信号：e
%输入信号功率：Px  正弦信号功率：Ps  噪声信号功率：Pv

clear;
close all;

N=20;                           %FIR滤波器长度
M=400;                          %信号长度
a=1;                            %噪声幅值
W0=0.1*pi;                      %正弦信号频率

v=a*[rand(1,M)-rand(1,M)];     
Pv=a^2/3;
Ps=10*Pv;
s=sin(W0*(1:M));   
A=sqrt(M*Ps/(s*s'));
s=A*s;
x=s+v;
Px=Ps+Pv;
u=1/(10*N*Px*1.1);
x1=[0,x(1:M-1)];               %延时信号
[w,y]=lms(x1,x,u,N);           %调用LMS算法
e=x-y;
figure(1);
subplot(211);plot(s);title('正弦信号s[n]');
subplot(212);plot(x);title('正弦加噪信号x[n]');

% 求系统的频率响应
fftn=1024;
H=fft(w,fftn);
n=[1:fftn/2]./(fftn/2);
figure(2);
plot(n,1-abs(H([1:fftn/2]))),title('幅频响应');grid on;

figure(3);
plot(v,'g-');
hold on;
plot(e,'m-');
legend('白噪声v[n]','系统输出信号e[n]');