如何在Matlab的ode45中固定采样时间？若不可行有何替代方案？

为什么你的ode45代码没做到固定采样时间 #

ode45是自适应变步长的ODE求解器，哪怕你指定了tspan = 0:T:10，它内部还是会根据误差自动调整计算步长，最终只是把计算结果插值到你指定的时间点上输出。你调试时看到的可变采样时间，就是它内部实际使用的计算步长——这是ode45的固有特性，没法强行改成固定步长运行。

不过如果你的需求只是输出固定时间间隔的结果，那你的代码其实已经实现了：最终plot里的t就是0、0.5、1…10这些点，只是求解过程是变步长的而已。

要真正固定采样时间求解，用这些方法 #

如果必须让求解过程全程使用固定步长，Matlab里有两种可行方案：

1. 用固定步长ODE求解器 #

Matlab自带固定步长的求解器，比如ode4（经典四阶Runge-Kutta），只需要通过odeset设置FixedStep参数即可。修改后的代码如下：

T = 0.5;
tspan = 0:T:10;
y0 = [0.2,0.3];
opts = odeset('FixedStep', T); % 设定固定步长
[t, y] = ode4(@(t,y) odefcn(t,y), tspan, y0, opts);
plot(t, y, '-o'), grid on
xlabel('Time')

function dx = odefcn(t,x)
dx = zeros(2,1);
u = [2 -2]*x(1:2) + 1;
dx(1) = x(1) - 2*x(2);
dx(2) = - x(2) + u;
end

注意：固定步长求解器的精度和稳定性由步长决定，步长太大可能导致结果偏离真实值，需要根据你的方程特性调整。

2. 手动实现固定步长数值算法 #

如果想完全掌控求解过程，可以自己写固定步长的数值方法，比如最常用的四阶Runge-Kutta（RK4）算法。示例代码：

T = 0.5;
t = 0:T:10;
y = zeros(length(t), 2);
y(1,:) = [0.2, 0.3]; % 初始值

for i = 1:length(t)-1
    current_t = t(i);
    current_y = y(i,:)';
    % RK4四步计算
    k1 = odefcn(current_t, current_y);
    k2 = odefcn(current_t + T/2, current_y + T/2*k1);
    k3 = odefcn(current_t + T/2, current_y + T/2*k2);
    k4 = odefcn(current_t + T, current_y + T*k3);
    % 更新下一个时间点的值
    y(i+1,:) = (current_y + T/6*(k1 + 2*k2 + 2*k3 + k4))';
end

plot(t, y, '-o'), grid on
xlabel('Time')

function dx = odefcn(t,x)
dx = zeros(2,1);
u = [2 -2]*x(1:2) + 1;
dx(1) = x(1) - 2*x(2);
dx(2) = - x(2) + u;
end

这种方式完全由你控制每一步的计算，适合对求解流程有严格要求的场景。

内容的提问来源于stack exchange，提问作者u1997