%《數字信號處理教程——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))