hc621.m

数字信号处理的matlab程序

%《数字信号处理教程——MATLAB释义与实现》第六章例6.2.1程序hc621
% 用信号流图法计算滤波器的传递函数
% 电子工业出版社出版  陈怀琛编著 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右下角分量赋值，使矩阵Q中未赋值的分量全置零
P(1,1)=1;        % 給矩阵P赋值
P(5,1)=0;        % 給P右下角分量赋值，使矩阵P未赋值的分量全置零
W=inv(eye(5)-Q)*P   % 按公式（6.2.7）求X=[x1;x2;x3;x4;y] 对U的系统函数
pretty(W(5))   	 % 求以y=x5为输出的系统函数