%《數(shù)字信號(hào)處理教程——MATLAB釋義與實(shí)現(xiàn)》第六章例6.6.1程序hc661
% 正余弦信號(hào)產(chǎn)生器的狀態(tài)空間分析
% 電子工業(yè)出版社出版  陳懷琛編著 2004年9月
%
syms  z  w0                         % 定義符號(hào)變量
A=[2*cos(w0), -1;1,0]; B=[1;0];		% 給A,B,C,D賦值
C=[sin(w0),0;cos(w0),-1];D=[0;1];
H=C*inv(z*eye(2)-A)*B+D;			% 調(diào)用由狀態(tài)方程求系統(tǒng)函數(shù)的公式 H=C*(zI-A)\B+D
pretty(simple(H))					% 用美觀的形式簡(jiǎn)化結(jié)果