%《MATLAB在電子信息課程中的應用》（第三版）第七章例7.24程序q724
% 正余弦產生器系統函數的計算機推導
% 電子工業出版社出版  陳懷琛 吳大正 高西全合著 2005年12月
%
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;		% 給Q陣右下角元素賦值可使其余未賦值的元素自動置零
P(1,1)=1; P(7,1)=0;
W=inv(eye(size(Q))-Q)*P
pretty(W(6))
pretty(W(7))