%《數字信號處理教程——MATLAB釋義與實現》第六章例6.5.4程序hc654
% 用信號流圖求格型梯形濾波器的傳遞函數
% 電子工業出版社出版  陳懷琛編著 2004年9月
%
clear, syms q           		% 對符號變量進行定義
k1=1/4; k2=1/2; k3=1/3; k0=1;	% 給網絡中的常數
C0=-0.2; C1=0.8; C2=1.5; C3=1;
Q(4,7)=q;		        % Q的第一個賦值元素為符號變量
Q(5,2)=1; Q(1,4)=-k3; Q(2,1)=1;	% 為聯接矩陣的非零元素賦值
Q(3,2)=k3; Q(3,4)=1; Q(5,8)=-k2;
Q(6,5)=1; Q(7,6)=k2; Q(7,8)=1;
Q(8,11)=q; Q(9,6)=1;
Q(9,12)=-k1; Q(10,9)=1;
Q(11,10)=k1; Q(11,12)=1; Q(12,10)=q;
Q(13,3)=C3; Q(13,7)=C2;
Q(13,11)=C1; Q(13,12)=C0;
Q(13,13)=0; P(1,1)=k0; P(13,1)=0;	%給右下角元素賦值，使未賦值的元素全為零
W=inv(eye(size(Q))-Q)*P		% 信號流圖方程解
pretty(W(13))		% 以美觀的形式顯示W中第13行