如何为cos(theta)设置Matlab ODE函数？二阶微分方程实现疑问

首先明确二阶微分方程转一阶方程组的正确逻辑：
对于方程 d²θ/dt² = cos(θ)，我们需要定义状态向量：

• Y(1) = θ(t)（当前角度）
• Y(2) = dθ/dt（当前角速度）

对应的一阶方程组为：

1. dY(1)/dt = Y(2)（角度的导数是角速度）
2. dY(2)/dt = cos(Y(1))（角速度的导数是角加速度，对应原方程右侧）

你的代码存在几处问题，修正后的版本如下：

function dYdt = osciODE(t,Y)
    dYdt = zeros(2,1);                      
    dYdt(1) = Y(2);                % dθ/dt = 角速度（即Y(2)）
    dYdt(2) = cos(Y(1));           % d²θ/dt² = cos(θ)，θ对应状态变量Y(1)
end

关键修正说明： #

• 移除多余的theta参数：θ是随时间变化的状态变量，不是固定输入参数，直接用Y(1)表示当前时刻的θ即可
• 替换未定义的theta1为Y(1)：Y(1)就是我们定义的θ(t)，对应方程中的角度变量
• 纠正对dYdt(1)的误解：dYdt(1)是Y(1)的一阶导数，也就是dθ/dt；二阶导对应dYdt(2)，即Y(2)的导数

调用该函数时，需要传入初始条件（初始角度+初始角速度），比如：

Y0 = [pi/4; 0];  % 初始角度π/4，初始角速度0
tspan = [0 10];  % 求解时间范围
[t,Y] = ode45(@osciODE, tspan, Y0);

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