| 查看: 554 | 回复: 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时有初值。 若能切上程序代码的话重谢！ |
回复此楼» 猜你喜欢
一志愿浙大生物325分求调剂
已经有4人回复
-
0831一轮调剂失败求助
已经有8人回复
-
本科南方医科大学 一志愿985 药学学硕284分 求调剂
已经有5人回复
-
295分求调剂
已经有12人回复
-
调剂求收留
已经有25人回复
-
272分材料子求调剂
已经有49人回复
-
211本科材料化工求调剂
已经有16人回复
-
279学硕食品专业求调剂院校
已经有12人回复
-
一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣
已经有11人回复
-
279求调剂
已经有8人回复
peterflyer 
木虫之王 (文学泰斗)
peterflyer
- 数学EPI: 10
- 应助: 20282 (院士)
- 金币: 146246
- 红花: 1374
- 帖子: 93091
- 在线: 7694.3小时
- 虫号: 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 (院士)
- 金币: 146246
- 红花: 1374
- 帖子: 93091
- 在线: 7694.3小时
- 虫号: 1482829
- 注册: 2011-11-08
- 性别: GG
- 专业: 功能陶瓷
3楼2015-11-08 22:49:23
10