如何实时显示Matlab的ode45结果？

Question

假设我们有以下一阶 ODE 系统的求解器：

   % func.m
   function dydt = func(t,y)
   dydt = [y(2); (1-y(1)^2)*y(2)-y(1)];

和主要代码：

% solver.m
tspan=0:1:10;
[t,y] = ode45(@func,tspan,[2; 0]);

如何实时显示结果 y1(t) 和 y2(t)，对于 ode45 进行的每个时间步长 t (t=0,1,2,...,10)，无需等待整个代码完成？

Answer 1

应使用 OutputFcn ode 求解器选项。例如，要绘制解与时间的关系图，可以使用内置输出函数 odeplot：

options= odeset('OutputFcn',@odeplot);
[t,y] = ode45(@func,[0 200],[2; 0],options);

您可以使用自己的输出函数。这是一个例子：

myOutputFcn= @(t,y,flag)fprintf('t= %s y=  %s\n',mat2str(t),mat2str(y))*0;
options= odeset('OutputFcn',myOutputFcn);
[t,y] = ode45(@f,0:1:10,[2; 0],options);

如何实时显示Matlab的ode45结果？

How to display Matlab's ode45 results in real time?

matlab

real-time

ode