MATLAB的答案

在MATLAB中求解三阶ODE

158次查看(最近30天)
显示旧的评论
塞吉奥Manzetti
塞吉奥Manzetti 2018年2月9日
  • 0
编辑： 塞吉奥Manzetti2018年2月15日
嗨，Matlab是非常有关这个命令的：
信谊a h Y(x) g x B E T
D3Y = diff(Y, 3)
eqn = a。* d3y -0.5 * x ^ 2 * y ==（abs（y））
D2Y = diff(Y, 2)
dy = diff（y）
cond1 = Y(0) == 1;
cond2 = DY(0) == 0;
cond3 = D2Y(0) == 0
Y(x) = dsolve(eqn, cond1, cond2, cond3)
乳胶(Y (x))
解它有极限吗?谢谢

登录回答这个问题。

接受的答案

卡兰吉尔
卡兰吉尔 2018年2月12日
你没有收到警告吗?如果你拿到了，警告清楚了吗?
警告:不能找到明确的解决方案。
>溶解(第201行)
Y（x）=
[空符号]
试着用数字解题 ODE45. 或类似的。
2的评论
显示＆nbsp1旧评论
卡兰吉尔
卡兰吉尔 2018年2月13日
你说的“其他方法”是什么意思?

登录评论。

更多答案（3）

塞吉奥Manzetti
塞吉奥Manzetti 2018年2月14日
Wolfram Alpha，它解决了它没有任何问题。
塞吉奥Manzetti
塞吉奥Manzetti 2018年2月15日
编辑：塞吉奥Manzetti 2018年2月15日
我在wolfram上尝试了这个，相当于这个:
信谊a h Y(x) g x B E T
D3Y = diff(Y, 3)
Y = x^2*Y = Y
D2Y = diff(Y, 2)
dy = diff（y）
cond1 = Y(0) == 1;
cond2 = DY(0) == 0;
cond3 = D2Y(0) == 1;
Y(x) = dsolve(eqn, cond1, cond2, cond3)
乳胶(Y (x))
我得到了一个结果，Z = 1/3*(exp(x) + 2*exp(-x/2)*cos((√(3)*x)/2))，然而，结果现在是不可见的，因为std计算时间超过了。
3评论
显示＆nbsp2旧的评论
Torsten.
Torsten. 2018年2月15日
．．．我很惊讶这个解决方案不依赖于“a”。

登录评论。

塞吉奥Manzetti
塞吉奥Manzetti 2018年2月15日
无论如何，ABS（Y）没有通过任何一种方法屈服的结果，而前者Y只在Wolfram中产生结果。
2的评论
显示＆nbsp1旧评论
塞吉奥Manzetti
塞吉奥Manzetti 2018年2月15日
是的，我知道这一点。
托斯顿，有没有其他解决办法
D3y - x^2y = ay，其中a是某个常数?

登录评论。

类别

标签

社区寻宝

在MATLAB中心找到宝藏，并发现社区如何可以帮助你!

开始狩猎！