rkthree.m

来自「运用MATLAB对控制系统进行仿真操作」· M 代码 · 共 11 行

function y=RKthree(h,A,B,C,U,N,n)
K=zeros(N,3);
X=zeros(N,1);
H=[0,1/3,2/3];
 for i=1:n
    K(:,1)=h*(A*X+B*U);
    K(:,2)=h*(A*(X+H(2)*K(:,1))+B*U);
    K(:,3)=h*(A*(X+H(3)*K(:,2))+B*U);
    X=X+1/4*(K(:,1)+3*K(:,2));
     y(i)=C*X 
 end