| 查看: 532 | 回复: 2 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
求助如何用matlab解微分方程数值解 已有1人参与
|
||
|
请教如何用matlab解如下方程，获得数值解（解析解是没希望了）： a*y''*y-[1+(y')^2]^1.5=0; 其中a=100000，为常数；y(x)是x的函数，y'和y''是y的一次和二次微分。 x=0时，y、y'和y''值都未知。但是x=无穷大时，y、y'和y''都为0。 试过ode45，但这个貌似要求x=0时有初值。 若能切上程序代码的话重谢！ |
回复此楼» 猜你喜欢
求个博导看看
已经有17人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
带资进组求博导收留
已经有11人回复
-
自荐读博
已经有5人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
求助院士们，这个如何合成呀
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 145878
- 红花: 1374
- 帖子: 93087
- 在线: 7693.9小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
3楼2015-11-08 22:49:23
peterflyer
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 145878
- 红花: 1374
- 帖子: 93087
- 在线: 7693.9小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
canbill: 金币+30, ★★★★★最佳答案, 没想到还有解析解，太感谢啦 2015-11-08 22:22:34
感谢参与，应助指数 +1
canbill: 金币+30, ★★★★★最佳答案, 没想到还有解析解，太感谢啦 2015-11-08 22:22:34
|
也不能完全说不能得到解析解的。实际上最后可得到积分形式的解析解。 令V=dy/dx,则y"=dV/dx=dV/dy*dy/dx=V*dV/dy 带入方程中得：a*V*dV/dy*y=[1+V^2]^(3/2) V=+-sqrt{1/Ln[C/y^(3/2)-1} 即dy/dx=+-sqrt{1/Ln[C/y^(3/2)-1} x=D+-Integral{dy/{sqrt{1/Ln[C/y^(3/2)-1} 其中C、D为积分常数。 [ 发自手机版 http://muchong.com/3g ] |
2楼2015-11-08 18:42:56