求解方程组:matlab 中的 2 个微分,1 个二次方程

Solve system of equations: 2 differential, 1 quadratic in matlab

编写一个 Matlab(或其他)代码来解决系统的数值问题: w'(t)=dw(t)/dt;

w'(t)=3*w(t)*y(t),
y'(t)=8*w(t)*y(t),
t^2=9+w(t)+y(t)

我不知道如何使用 ode45 解决这个问题,因为它有 2 个解决方案。

为什么你需要用数字来解决这个问题?对于数值解,您至少需要一个初始条件,即 w(0), y0).

注意,通过比较前两个等式:8w'(t) = 3y'(t)

然后,推导第三个方程得到:

2t = w'(t)+y'(t)

这意味着:

8*3*2t = 8*3*w'(t)+8*3*y'(t)

48t = 8*3*w'(t)+8*8*w'(t)

48t = 88*w'(t)

6t = 11*w'(t)

因此 w'(0)=0y'(0)=0

因此,从第一个等式:w(0)*y(0)=0.

因为方程是对称的,所以你说的解有两种。假设w(0)=0,则从第三个等式,'y(0)=-9'。从 6t = 11*w'(t) 我们有 w(t)=(6/11)ty(t)=-9+(48/33)t.

另一个解决方案是y(t)=(6/11)tw(t)=-9+(48/33)t