program1_8.m

此为matlab的代码。包括数值o45算法及o23算法

%MATLAB Program1_8
A1=[0 1;-2 -3];B1=[0;1];C1=[1 0];D1=[0];
A2=[0 1;-12 -7];B2=[0;1];C2=[2 1];D2=[0];
[num1,den1]=ss2tf(A1,B1,C1,D1)
G1=tf(num1,den1)
[num2,den2]=ss2tf(A2,B2,C2,D2)
G2=tf(num2,den2)
%并联连接
G11=G1+G2
num11=[1,6,15,16];
den11=[1 10 35 50 24];
[A,B,C,D]=tf2ss(num11,den11)
%串联连接
G22=G2*G1
num22=[8.882e-016,1,2];
den22=[1,10,35,50,24];
[A,B,C,D]=tf2ss(num22,den22)
%负反馈连接
G33=G1*(inv(1+G2*G1))
num33=[8.882e-016,1,10,35,50,24];
den33=[1,13,67,176,249,180,52];
[A,B,C,D]=tf2ss(num33,den33)