fir.m

该文件为Z的一个子程序

%用窗函数进行设计,wp=0.2pi,ws=0.3pi,Rs=40dB,Rp=3dB.
wp=0.2*pi;
ws=0.3*pi;
wdelta=ws-wp;
N=ceil(8*pi/wdelta);       
wc=(0.2+0.3)*pi/2;
r=(N-1)/2;
if rem(N,2)~=0
    hdn(r+1)=wc/pi;
    for n=(0:r-1)&(r+1:N-1)
        hdn(n+1)=sin(wc*(n-r))/pi/(n-r);
    end
else
    for n=0:N-1
        hdn(n+1)=sin(wc*(n-r))/pi/(n-r);
    end
end
wn=input('wn=');     %输入窗函数
wn=wn';
h=hdn.*wn;
H=fft(h,512);
w=2*[0:511]/512;
subplot(2,1,1);plot(w,20*log10(abs(H)));
grid on;
xlabel('频率 w/pi');
ylabel('幅度 dB');
Haphase=unwrap(angle(H));
Haphase=Haphase*180/pi;
subplot(2,1,2);plot(w,Haphase);
grid on;
xlabel('频率 w/pi');
ylabel('相位 度');