%《MATLAB及其在理工課程中的應(yīng)用指南》第九章例9－3－4程序
% 由狀態(tài)方程求傳遞函數(shù)
% 西安電子科技大學(xué)出版社出版  陳懷琛編著 1999年10月初版，2004年11月第二版
%
clear, 
disp('輸入狀態(tài)方程系數(shù)矩陣A,B,C,D')
disp('輸入A=（書(shū)上取 [0.9501,0.8913,0.8214,0.9218;0.2311,0.7621,0.4447,0.7382;...')
disp('0.6068,0.4565,0.6154,0.1763;0.4860,0.0185,0.7919,0.4057] ')
A=input('A= ')
B=input('B= （書(shū)上取 [0.9355; 0.9169; 0.4103; 0.8936] )'), 
C=input('C= （書(shū)上取 [0.0579, 0.3529, 0.8132, 0.0099] ) '), 
D=input('D= （書(shū)上取 0 )')
g=poly(eig(A)); n=length(A);	        % 分母多項(xiàng)式系數(shù)g(i)由特征根求得
P(:,:,1)=eye(n);
f(1)=D*g(1);						    % 分子多項(xiàng)式系數(shù)f(i)由遞推求得
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