q618.m

MATLAB及在电子信息课程中的应用例题索引与程序

%《MATLAB在电子信息课程中的应用》章例618程序
% 线性控制系统的模型描述和转换
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月

%　由状态空间转换为传递函数
clear, 
disp('输入状态方程系数矩阵A,B,C,D')
disp('输入A=（书上取[0.2844,0.5828,0.4329,0.5298;0.4692,0.4235,0.2259,0.6405;...')
disp('0.0648,0.5155,0.5798,0.2091;0.9883,0.3340,0.7604,0.3798] }')
A=input('A= ')
B=input('B=（书上取 [0.7833;0.6808;0.4611;0.5678] )') 
C=input('C= （书上取 [0.7942,0.0592,0.6029,0.0503] )'), 
D=input('D= （书上取 0 )')
g=poly(eig(A)); n=length(A);	% 分母多项式系数由特征根求得
P(:,:,1)=eye(n);
f(1)=D*g(1);						% 分子多项式系数由递推求得
f(2)=C*P(:,:,1 )*B+D*g(2);
for i=2:n 
   P(:,:,i )=A*P(:,:,i-1 )+g(i)*eye(n); 
   f(i+1)=C*P(:,:,i )*B+D*g(i+1); 
end
f, g