num=input('輸入分子系數＝');
den=input('輸入分母系數＝');
[z,p,k]=tf2zp(num,den);
m=abs(p);
disp('零點在');disp(z);
disp('極點在');disp(p);
disp('增益常數');disp(k);
disp('極點半徑');disp(m);
sos=zp2sos(z,p,k);
disp('二階部分');disp(real(sos));
zplane(num,den);