%《數字信號處理教程——MATLAB釋義與實現》第六章例6.6.1程序hc661
% 正余弦信號產生器的狀態空間分析
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
syms  z  w0                         % 定義符號變量
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;			% 調用由狀態方程求系統函數的公式 H=C*(zI-A)\B+D
pretty(simple(H))					% 用美觀的形式簡化結果