如何使用ODE45解决耦合差分方程在Matlab(How to solve coupled diff

我有两个微分方程：DA / DT = A（0.3 / A ^ 3 + 0.7）^ 1/2和Dτ/ DT = 1 /一个。 初始条件是T = 0; 一个分别= 1和τ= 0，。 我该如何解决在Matlab中的方程？ 我需要计算的，t和τ不同的值还积τVs的。 谢谢。

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]);

这是结果：

文章来源: How to solve coupled differential equation in matlab using ode45