| 查看: 1235 | 回复: 4 | ||
[求助]
求一个微分方程组的数值解，总出错，请大神们帮忙看看，很急！！！很急！！！
|
|
请大家帮忙看一下，拜托了，真的很急，想求出y1,y2,y3,y4,而且画出图来，可是老出错，我都弄了好久了。 NDSolve[{y[1]''[ t] == [Subscript[Log, 2] (y[1][t] + y[2][t] + y[3][t] + y[4][t])/ y[1][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t]]*400^-2，y[2]''[ t] == [Subscript[Log, 2] (y[1][t] + y[2][t] + y[3][t] + y[4][t])/ y[2][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t]]*400^-2，y[3]''[ t] == [Subscript[Log, 2] (y[1][t] + y[2][t] + y[3][t] + y[4][t])/ y[3][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t]]*400^-2，y[4]''[ t] == [Subscript[Log, 2] (y[1][t] + y[2][t] + y[3][t] + y[4][t])/ y[4][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t]]*400^-2，Derivative[1][y[1]][0] == 0, Derivative[1][y[2]][0] == 0, Derivative[1][y[3]][0] == 0, Derivative[1][y[4]][0] == 0, y[1][0] == 300000000, y[2][0] == 270000000, y[3][0] == 310000000, y[4][0] == 290000000}, {y[1][t], y[2][t], y[3][t], y[4][t]}, t]，出错为Syntax::sntxf: "{" cannot be followed by "RowBox[{RowBox[{\(<< 14 >>\), ",", RowBox[{\(<< 1 >>\), "==", RowBox[{"", \(<< 9 >>\), ""}]}], ",", \(\(\(y[4]\)[0]\) == 290000000\)}], "}"}]". Syntax::tsntxi: "<<15>>,y[4][0]==290000000" is incomplete; more input is needed. Syntax::sntxi: "Incomplete expression; more input is needed \!\(\"\"\)." |
回复此楼» 猜你喜欢
面上项目filecode邪修
已经有6人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有14人回复
-
2026国自然放榜时间
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有22人回复
-
这样的filecode谁见过
已经有21人回复
-
filecode与中标关系的预测
已经有5人回复
-
fileCode有新解读？
已经有18人回复
-
关于filecode，很负责任的告诉大家
已经有7人回复
-
【2027博士申请】纳米药物递送方向
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于含参数的常微分方程求解
已经有0人回复
- 如何求偏微分方程组的数值解
已经有1人回复
- 帮忙：求解一组微分代数方程，奖励100金币
已经有7人回复
- 数学史上最富传奇色彩的未解猜想——ABC猜想，来自日本京都大学教授望月新一的证明
已经有15人回复
- 浙江大学电路原理及网络理论（内部课件）
已经有144人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 求大神---MATLAB-ode23tb-微分方程组--运行失败
已经有2人回复
- 牛顿到庞卡莱，ODE = NP？——常微分方程发展史和评述 转载
已经有12人回复
- 求解一个微分方程组
已经有0人回复
- 青基，三个建议资助，为什么挂了？求各位给点意见！
已经有13人回复
- 【求助】求稳定性的问题~~急~~
已经有8人回复
- 【书籍】《常用算法程序集》
已经有44人回复
- 【书籍】《精通Matlab综合辅导与指南》
已经有16人回复
emailwzh
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 5319.3
- 帖子: 271
- 在线: 115.8小时
- 虫号: 196480
- 注册: 2006-02-24
- 性别: GG
- 专业: 岩土与基础工程
【答案】应助回帖
感谢参与，应助指数 +1
|
DSolve[{y[1]''[t] == Log[2, ((y[1][t] + y[2][t] + y[3][t] + y[4][t])/y[1][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t])]*400^-2, y[2]''[t] == Log[2, ((y[1][t] + y[2][t] + y[3][t] + y[4][t])/y[2][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t])]*400^-2, y[3]''[t] == Log[2, ((y[1][t] + y[2][t] + y[3][t] + y[4][t])/y[3][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t])]*400^-2, y[4]''[t] == Log[2, ((y[1][t] + y[2][t] + y[3][t] + y[4][t])/y[4][t] - 2*10^-9*(y[1][t] + y[2][t] + y[3][t] + y[4][t]) - 1*y[1]'[t])]*400^-2, Derivative[1][y[1]][0] == 0, Derivative[1][y[2]][0] == 0, Derivative[1][y[3]][0] == 0, Derivative[1][y[4]][0] == 0, y[1][0] == 300000000, y[2][0] == 270000000, y[3][0] == 310000000, y[4][0] == 290000000}, {y[1][t], y[2][t], y[3][t], y[4][t]}, t] 改成这样试试 |
2楼2016-03-23 16:16:09
emailwzh
铁虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 5319.3
- 帖子: 271
- 在线: 115.8小时
- 虫号: 196480
- 注册: 2006-02-24
- 性别: GG
- 专业: 岩土与基础工程
3楼2016-03-23 16:24:29
4楼2016-03-23 17:06:35
5楼2016-03-23 17:08:29
20