我未能解决4个方程微分方程。我给代码。两方程相同的代码工作。程序卡是dsolve (ode1, ode2, ode3 ode4)或dsolve(常微分方程)和从来没有完成。

10视图(30天)
信谊u (t) v (t) x (t) y (t)
ode1 = diff (u, t) = = 3 * u + 4 * v + 5 * x + 6 * y;
ode2 = diff (v, t) = = 4 * u + 3 * v + 9 * x + 8 * y;
t ode3 = diff (x) = = 3 * u + 9 * x + v + 4 * 6 * y;
ode4 = diff (y, t) = = 6 * u + 8 * v + 9 * x + 5 * y;
常微分方程= [ode1; ode2; ode3 ode4);
S = dsolve(常微分方程);
uSol (t) = S.u;
vSol (t) = S.v;
xSol (t) = S.x;
ySol (t) = S.y;
[uSol (t) vSol (t) xSol (t) ySol (t)) = dsolve(常微分方程)
cond1 = u (0) = = 0;
cond2 = v (0) = = 0;
cond3 = x (0) = = 0;
cond4 = y (0) = = 0;
气孔导度= [cond1;cond2;cond3;cond4];
[uSol (t) vSol (t) xSol (t) ySol (t)) = dsolve(诗赋,气孔导度);
ezplot (uSol)
持有
ezplot (vSol)
持有
ezplot (xSol)
持有
ezplot (ySol)
网格;
传奇(“uSol”,“vSol”,“xSol”,“ySol”)

接受的答案

捕鸟者
捕鸟者 2017年11月4日
你的方程 ode1 ode4 实际上是包含线性相关的条款
3 * u + 4 * v (ode1), 6 * u + 8 * v (ode4)
这种情况带来了无限多的系统解决方案,因此实际上是一个独特的解决方案不能被发现。万博 尤文图斯例如,如果你改变你的术语
6 * u + 7 * y
在ode4,您将看到,很快发现一个独特的解决方案。

答案(1)

沃尔特·罗伯森
沃尔特·罗伯森 2017年11月4日
解决方案是u (t) = 0, v (t) = 0, x (t) = 0, y (t) = 0。MATLAB是能够想出的解决方案如果你跳过到最后dsolve后初始化变量。
你有三个dsolve,前两个都一样除了他们如何输出结果。第三人解决了常微分方程,但边界条件相同,MATLAB是能够处理。

社区寻宝

找到宝藏在MATLAB中央,发现社区如何帮助你!

开始狩猎!