%《MATLAB在電子信息課程中的應用》第六章例6.17程序q617
% 線性控制系統的模型描述和轉換
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月

% 由傳遞函數模型轉為其他模型
format compact
f=input('傳遞函數分子系數數組f= [f(1),..., f(nf)] =   書上取[2,-5,3]');
g=input('傳遞函數分母系數數組g= [g(1), g(2),..., g(ng)] =  書上取[2,3,5,9]  ');
printsys(f,g,'s')
disp('轉為零極增益模型')
z=roots(f)
p=roots(g)
k=f(min(find(f(:)~=0)))/g(1)
[z1,p1,k1]=tf2zp(f,g)
disp('轉為零極留數模型')
[r,p,h]=residue(f,g)
disp('轉為狀態空間模型')
[A,B,C,D]=tf2ss(f,g);
printsys(A,B,C,D)