program1_9.m

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

%MATLLAB Program1_9
num1=[2];
den1=[1 1 0];
G1=tf(num1,den1)
num2=[1 2];
den2=[1 3];
G2=tf(num2,den2)
num3=[10];
den3=[1 2];
G3=tf(num3,den3)
num4=[5];
den4=[1 4];
G4=tf(num4,den4)
G11=G1*(inv(1+G2*G1));
G22=G11*G3;
G=G4*(inv(1+G22*G4));
num=[5 55 235 515 650 480 160];
den=[1 15 91 291 642 1416 2316 1328];
[A,B,C,D]=tf2ss(num,den)