m0723.m

来自「matlab的经典讲义」· M 代码 · 共 10 行

Wp=300*2*pi;Ws=400*2*pi;
Rp=0.3;As=80;
Fs=1000;
[N,Wn]=cheb1ord(Wp,Ws,Rp,As,'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,10);
freqs(num2,den2);