rktwo.m

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

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