pinlvxiangying.m

自己编写的自适应notch fileter算法

b=[1,-2*cos(pi*0.5),1];
a=[1,-2*(1-0.021)*cos(pi*0.5),1-2*0.021];
[H,w]=freqz(b,a,100);
magh=abs(H);
subplot(2,1,1);
plot(w/pi,magh);
grid;
xlabel('频率单位pi');
ylabel('幅度');
title('幅度响应');
phah=angle(H);
subplot(212);
plot(w/pi,phah/pi);
xlabel('频率单位pi');
ylabel('相位单位pi');
title('相位相应');
grid;

%b=[1,-2*cos(pi/4),1];
%a=[1,-(1+0.9825)*cos(pi/4),0.9825];