q803.m

<matlab在电子信息中的应用>一书中个章节的配套源代码。

%《MATLAB在电子信息课程中的应用》第八章例8.3程序q803
% 用LTI模型作信号流图的计算
% 电子工业出版社出版  陈怀琛 吴大正 高西全合著 2001年10月
% 
% 本题与第六章例6.20相同，解法不同，第六章用符号运算，本章用LTI对象模型来完成

clear
k=input('用什么模型？传递函数-键入1，零极增益-键入2,状态空间-键入3,k= ');
switch k
case 1
   G1=tf([1,0],[1,1]), 		% 定义lti对象，全用tf形式
   G2=tf(3,[1,2]),   G3=tf([1,4],[1,5,6])
case 2
   G1=zpk(0,-1,1), 		% 定义lti对象,有两个用zpk形式
   G2=zpk([ ],-2,3),   G3=tf([1,4],[1,5,6])
otherwise
	G1=ss(zpk(0,-1,1)), 		% 定义lti对象,有两个用zpk形式
   G2=zpk([ ],-2,3),   G3=tf([1,4],[1,5,6])

end
R(3,2)=G1;					% 采用字符矩阵，第一条赋值语句右端必须是字符变量					
R(2,1)=1;R(2,3)=-1;R(2,5)=-1;	% 列出连接矩阵
R(4,3)=1;R(4,1)=1;R(4,5)=-1;
R(5,4)=G2;
R(6,3)=1;R(6,5)=1;R(6,7)=-1;
R(7,6)=G3;
R(8,7)=5;
R(:,8)=zeros(8,1);	% 加一个全零列，补成方阵
P=[1;0;0;0;0;0;0;0];	
I=eye(size(R));
W=(I-R)\P;			% 求出完整的传递矩阵
if k==3 W, 	% 若要状态空间，输出W
else W8=W(8),	% 若非状态空间，输出x8的传递函数为W中的第八行
end