hc662.m

来自「这是数字信号处理——matlab释义与实现一书的源代码」· M 代码 · 共 10 行

%《数字信号处理教程——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))					% 用美观的形式简化结果