| 查看: 524 | 回复: 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时有初值。 若能切上程序代码的话重谢！ |
回复此楼» 猜你喜欢
博士读完未来一定会好吗
已经有6人回复
-
小论文投稿
已经有3人回复
-
Bioresource Technology期刊，第一次返修的时候被退回好几次了
已经有9人回复
-
心脉受损
已经有3人回复
-
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗
已经有8人回复
-
申请2026年博士
已经有6人回复
-
请问哪里可以有青B申请的本子可以借鉴一下。
已经有5人回复
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 145691
- 红花: 1374
- 帖子: 93067
- 在线: 7693.7小时
- 虫号: 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 (院士)
- 金币: 145691
- 红花: 1374
- 帖子: 93067
- 在线: 7693.7小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
3楼2015-11-08 22:49:23