zimpinvar.m

用MATLAB实现数字信号处理的程序

%Impinvar函数：模拟滤波器变换成数字滤波器的脉冲响应不变法
Wp=2*pi*400;
Ws=2*pi*600;
Rp=0.3;
Rs=600;
Fs=1000;
%选择滤波器的最小阶数
[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs,'s');
%创建低通巴特沃斯滤波器
[Z,P,K]=cheb1ap(N,Rp);
[A,B,C,D]=zp2ss(Z,P,K);
%低通到低通的转变
[At,Bt,Ct,Dt]=lp2lp(A,B,C,D,Wn);
[num1,den1]=ss2tf(At,Bt,Ct,Dt);
%脉冲响应不变法把模拟滤波器转换为数字滤波器
[num2,den2]=impinvar(num1,den1,Fs)
%绘制幅度频率响应图线
[H,W]=freqz(num2,den2);
plot(W*Fs/2/pi,abs(H));grid;
xlabel('频率/Hz'),ylabel('幅度');