q617.m

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

%《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)               % 打印出系统函数