%《MATLAB及在電子信息課程中的應(yīng)用》第八章例8.1程序q801
% 用新舊工具箱作MIMO模型變換的對比
% 電子工業(yè)出版社出版  陳懷琛 吳大正 高西全合著 2001年10月初版，2003年7月第二版
%
clear
A=[1,-1,0;0,2,0;1,0,4];B=[1;0;-1];C=[2,0,0;1,2,3];D=[0;0];
% 方法一  用舊的控制系統(tǒng)工具箱命令
disp('對輸入1的系數(shù)向量及傳遞函數(shù)')
[f,g]=ss2tf(A,B,C,D),pause
printsys(f,g,'s'),pause
disp('對輸入1的零極增益參數(shù)')
[z,p,k]=ss2zp(A,B,C,D),pause
% 方法二  用LTI對象和新的控制系統(tǒng)工具箱命令提取參數(shù)
disp('方法二  用LTI對象和新的控制系統(tǒng)工具箱命令提取參數(shù)')
A=input('輸入狀態(tài)方程系數(shù)矩陣A=  ');		% 允許輸入其他狀態(tài)矩陣參數(shù)
if isempty(A) 
    A=[1,-1,0;0,2,0;1,0,4],B=[1;0;-1],C=[2,0,0;1,2,3],D=[0;0],
else B=input('輸入狀態(tài)方程系數(shù)矩陣B=  ')
    C=input('輸入狀態(tài)方程系數(shù)矩陣C=  ')
    D=input('輸入狀態(tài)方程系數(shù)矩陣D=  ')
end
sys=ss(A,B,C,D);,pause		
[f1,g1]=tfdata(sys),pause			% 轉(zhuǎn)換提取tf系數(shù)向量
disp('要知道f,g的內(nèi)容，鍵入f1{:.:},g1{:,:}')
f1{:,:},g1{:,:},pause			% 提取tf系數(shù)向量具體值
[z1,p1,k1]=zpkdata(sys),pause		% 轉(zhuǎn)換提取tf系數(shù)向量
disp('要知道z,p的內(nèi)容，鍵入z1{:.:},p1{:,:}')
z1{:,:},p1{:,:},pause			% 提取zpk系數(shù)向量具體值
%
disp('方法三  用LTI對象和新的控制系統(tǒng)工具箱命令再建新模型')
disp('系統(tǒng)的傳遞函數(shù)模型')
systf=tf(sys),pause			% 生成等價的tf對象的LTI模型
disp('系統(tǒng)的零極增益模型')
syszp=zpk(sys),pause			% 生成等價的zpk對象的LTI模型