hc611.m

书籍中matlab第二部分代码

%《数字信号处理教程——MATLAB释义与实现》第四章例6.1.1程序hc611
% 用信号流图法计算滤波器的传递函数
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear,
syms q                    % 规定q为符号变量
Q(2,1)=q;Q(1,4)=-0.1;    % 給矩阵Q赋值，第一个分量必须含符号
Q(3,1)=2;Q(3,2)=1;Q(3,4)=0.3;
Q(4,3)=q;
Q(5,3)=0.5;Q(5,4)=1;
Q(5,5)=0;        % 給Q右下角分量赋值，使全矩阵未赋值的分量全置零
P(1,1)=1;        % 給矩阵P赋值
P(5,1)=0;        % 給P右下角分量赋值，使全矩阵未赋值的分量全置零
W=inv(eye(5)-Q)*P   % 按公式（6.2.7）求系统函数
pretty(W(5))    % 求以y=x5为输出的系统函数