m0722.m

来自「MATLAB仿真技术与实例应用教程」· M 代码 · 共 10 行

[z,p,k]=ellipap(4,5,45);
[A,B,C,D]=zp2ss(z,p,k);
[At1,Bt1,Ct1,Dt1]=lp2lp(A,B,C,D,10*pi);
[num1,den1]=ss2tf(At1,Bt1,Ct1,Dt1);
figure;
freqs(num1,den1);
pause(10);
[At2,Bt2,Ct2,Dt2]=lp2hp(A,B,C,D,60*pi);
[num2,den2]=ss2tf(At2,Bt2,Ct2,Dt2);
freqs(num2,den2);