main.m

来自「动态系统仿真例子」· M 代码 · 共 33 行

% 系统仿真例子（开环 无控制时）
% 主程序

% 系统函数句柄 
SysFun=@SysSs;

% 状态初值
% x0=[0.144,0.8862];
x0=[2,5];

% Simulation Settings
H=0.1;
t0=0;
tf=20;
tspan=t0:H:tf;
StepNum=size(tspan,2);
u=0;
xArr=x0;
% Simulate
for i=1:StepNum
   xArr(i,:)=x0;
  
%  开环，无控制
   x0=RK4Step(SysFun,H,x0);
   
end
x=x0
plot(tspan,xArr(:,1));
plot(tspan,xArr(:,2));

xlabel('time/s');
ylabel('state x1');