%《MATLAB及在電子信息課程中的應用》第六章例6.17程序q617
% 線性控制系統的模型描述和轉換
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
% 由傳遞函數模型轉為其他模型
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)               % 打印出系統函數