%EX0601 寫出二階系統的狀態方程、傳遞函數和部分分式
zeta=0.707;wn=1;
A=[0 1;-wn^2 -2*zeta*wn];
B=[0;wn^2];
C=[1 0];
D=0;
G=ss(A,B,C,D)  %建立狀態方程模型
 num=1;
den=[1 1.414 1];
G=tf(num,den)  %得出傳遞函數
 z=roots(num)
p=roots(den)
 zpk(z,p,1)