sy1.m

matlab实现陷波器

clear;close all;
N=40;
NS=400;
Fs=8000;
u=0.002;
v=rand(1,NS)-rand(1,NS);
s=2*sin(2*pi*[0:(NS-1)]*100/Fs);
figure,hold on,grid on
x=s+v;
%plot(x,'k-.')
delayx=[0,x(1:NS-1)];
[w,y,e]=lms(x,delayx,u,N);
plot(v,'g-')
plot(s,'b-')
plot(y,'r-')
%legend('加噪信号','噪声信号','理想信号','滤波后信号');
legend('噪声信号','理想信号','滤波后信号');
%m=fft(w);
[m,a]=freqz(w);
figure
%subplot(211),plot(a,20*log10(abs(m))),title('幅频响应');
%subplot(212),plot(a,20*log10(angle(m))),title('相频响应');
subplot(211),plot(a,abs(m)),title('幅频响应');grid on
subplot(212),plot(a,angle(m)),title('相频响应');grid on