q801.m

MATLAB及在电子信息课程中的应用,电子工业出版社出版 陈怀琛 吴大正 高西全合著

%《MATLAB及在电子信息课程中的应用》第八章例8.1程序q801
% 用新旧工具箱作MIMO模型变换的对比
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 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];
% 方法一  用旧的控制系统工具箱命令
disp('对输入1的系数向量及传递函数')
[f,g]=ss2tf(A,B,C,D),pause
printsys(f,g,'s'),pause
disp('对输入1的零极增益参数')
[z,p,k]=ss2zp(A,B,C,D),pause
% 方法二  用LTI对象和新的控制系统工具箱命令提取参数
disp('方法二  用LTI对象和新的控制系统工具箱命令提取参数')
A=input('输入状态方程系数矩阵A=  ');		% 允许输入其他状态矩阵参数
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('输入状态方程系数矩阵B=  ')
    C=input('输入状态方程系数矩阵C=  ')
    D=input('输入状态方程系数矩阵D=  ')
end
sys=ss(A,B,C,D);,pause		
[f1,g1]=tfdata(sys),pause			% 转换提取tf系数向量
disp('要知道f,g的内容，键入f1{:.:},g1{:,:}')
f1{:,:},g1{:,:},pause			% 提取tf系数向量具体值
[z1,p1,k1]=zpkdata(sys),pause		% 转换提取tf系数向量
disp('要知道z,p的内容，键入z1{:.:},p1{:,:}')
z1{:,:},p1{:,:},pause			% 提取zpk系数向量具体值
%
disp('方法三  用LTI对象和新的控制系统工具箱命令再建新模型')
disp('系统的传递函数模型')
systf=tf(sys),pause			% 生成等价的tf对象的LTI模型
disp('系统的零极增益模型')
syszp=zpk(sys),pause			% 生成等价的zpk对象的LTI模型