如何使用ode45在MATLAB中求解耦合微分方程

Question

我有兩個微分方程：da / dt = a（.3 / a ^ 3 + .7）^ 1/2和dτ/ dt = 1 / a。 初始條件為t = 0； a = 1和τ= 0。 如何在Matlab中求解方程式？ 我需要計算a，t和τ的不同值，並繪制τvs a。 謝謝。

Answer 1

那很容易。

首先編寫一個函數來實現您的微分方程，並將其保存為與函數名稱相對應的文件名：

function dy = my_ode(t,y)
dy(1) = y(1)*(0.3/y(1)^3 + 0.)^(1/2); % a
dy(2) = 1/dy(1); % tau

然后在MATLAB中，使用您的函數調用ode45求解器

[t,y] = ode45(@my_ode,[0 10],[1; 0]);

結果如下：