hc661.m

来自「数字信号处理的matlab程序」· M 代码 · 共 14 行

%《数字信号处理教程——MATLAB释义与实现》第六章例6.6.1程序hc661
% 正余弦信号产生器的信号流图分析
% 电子工业出版社出版  陈怀琛编著 2004年9月
%
clear
syms q  w0                              % 定义符号变量
Q(2,1) = q; Q(1,4)=1;Q(3,2)=cos(w0);    % 信号流图联接矩阵赋值，第一项右端需有符号变量
Q(4,3)=2; Q(4,5) = -1;
Q(5,2)=q; Q(6,2) = sin(w0);
Q(7,1) =1; Q(7,3)= -1; Q(7,7)=0;
P(1,1)=1; P(7,1)=0;                     % 信号流图联接矩阵赋值
W=inv(eye(size(Q))-Q)*P                 % 用W=(I-Q)\P求系统函数
pretty(W(6)), 
pretty(W(7))