图44四阶龙格-库塔数值积分方法的仿真曲线 Cb=[000.250.5 KI=Ab*x+Bb*U K2=Ab*(X+h*K1/2)+Bb*U K3=Ab*(X+ h*K2/2)+Bb*U K4=Ab*(X+h*K3)+Bb*U X=X+h*(K1+2*K2+2*K3+K4)/6 Y=[Y, Cb*X t=[t, t(i)+h]: plot(t,Y) 当由已知系统开环传递函数求取系统开环状态方程时,除了可以利用ss函数直 接求取外,还可以使用状态方程相似变换函数sss进行求取,下面给出上例中应用 相似变换函数ss求取系统开环状态方程的基本用法。这里主要利用矩阵水平翻转 函数 flipper构造出反对角阵作为变换阵,从而得到相应的状态方程 num=2*1, 2: den= conv(conv(1, 0],[1, 3 ).1, 2, 2D) 1= tf(num, den); G2=SS(G1) plr(eye(4)); G=sS2ss(G2, T) 12127 图 4-4 四阶龙格-库塔数值积分方法的仿真曲线 Bb = [1; 0; 0;0]; Cb = [0 0 0.25 0.5]; Db = 0; for i = 1:n K1 = Ab * X + Bb*U; K2 = Ab * (X + h*K1/2)+Bb*U; K3 = Ab * (X + h*K2/2)+Bb*U; K4 = Ab * (X + h*K3)+Bb*U; X = X+ h*(K1 + 2*K2+ 2*K3 +K4)/6; Y = [ Y, Cb*X]; t = [t, t(i)+h]; end plot (t,Y) Grid 当由已知系统开环传递函数求取系统开环状态方程时，除了可以利用 ss 函数直 接求取外，还可以使用状态方程相似变换函数 ss2ss 进行求取，下面给出上例中应用 相似变换函数 ss2ss 求取系统开环状态方程的基本用法。这里主要利用矩阵水平翻转 函数 fliplr 构造出反对角阵作为变换阵，从而得到相应的状态方程。 num = 2 *[1,2]; den = conv(conv([1,0],[1,3]),[1,2,2]); G1= tf (num, den); G2 = ss(G1); T = fliplr (eye(4)); G = ss2ss(G2,T)