| 查看: 535 | 回复: 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时有初值。 若能切上程序代码的话重谢！ |
回复此楼» 猜你喜欢
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
-
投稿Elsevier的杂志（返修），总是在选择OA和subscription界面被踢皮球
已经有6人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有4人回复
-
求个博导看看
已经有18人回复
-
自荐读博
已经有6人回复
-
青基代表作，AAAI之类的A会的special track在国内认可度高吗？还是归为workshop之流？
已经有3人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有6人回复
-
上海工程技术大学张培磊教授团队招收博士生
已经有4人回复
-
临港实验室与上科大联培博士招生1名
已经有9人回复
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 145890
- 红花: 1374
- 帖子: 93089
- 在线: 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
peterflyer
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 145890
- 红花: 1374
- 帖子: 93089
- 在线: 7693.9小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
3楼2015-11-08 22:49:23