24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1144  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

672418539

新虫 (初入文坛)

[求助] 求一个微分方程组的数值解,总出错,请大神们帮忙看看,很急!!!很急!!!

请大家帮忙看一下,拜托了,真的很急,想求出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 \!\(\"\"\)."
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

672418539

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by emailwzh at 2016-03-23 16:24:29
NDSolve用法不正确!

哪里有错?看不出来啊,求指教。
5楼2016-03-23 17:08:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

emailwzh

铁虫 (小有名气)

【答案】应助回帖

感谢参与,应助指数 +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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

emailwzh

铁虫 (小有名气)

【答案】应助回帖

NDSolve用法不正确!
3楼2016-03-23 16:24:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

672418539

新虫 (初入文坛)

引用回帖:
2楼: Originally posted by emailwzh at 2016-03-23 16:16:09
DSolve[{y'' ==
   Log*400^-2,
  y'' ==
   Log*400^-2,
  y'' ==
   Log*400^-2,
  y'' ==
   Log*400^-2, Derivative == 0,
  Derivative == 0, Derivative == 0,
  Derivative == 0, y == 30 ...

亲,不对哦,没只是改了个形式而已,没有结果,也没有提醒错误
4楼2016-03-23 17:06:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专业383求调剂 +6 郭阳阳阳成 2026-04-03 6/300 2026-04-03 22:32 by 无际的草原
[考研] 311求调剂 +20 zchqwer 2026-04-01 22/1100 2026-04-03 22:09 by lglzsd
[考研] 278求调剂 +5 Yy7400 2026-04-03 5/250 2026-04-03 22:01 by 心上人啊啊啊
[考研] 一志愿郑州大学材料与化工085600,求调剂 +17 吃的不少 2026-04-02 17/850 2026-04-03 21:48 by qlm5820
[考研] 求调剂 +3 wos666 2026-04-03 3/150 2026-04-03 21:36 by lbsjt
[考研] 化学308分调剂 +21 你好明天你好 2026-03-30 22/1100 2026-04-03 21:32 by 百灵童888
[考研] 兽医调剂 +3 wh119216 2026-04-02 3/150 2026-04-03 19:34 by zrongyan
[考研] 335求调剂 +7 沈清璃 2026-04-03 7/350 2026-04-03 18:55 by lijunpoly
[考研] 315求调剂 +12 小羊小羊_ 2026-04-02 12/600 2026-04-03 18:22 by ls刘帅
[考研] 303求调剂 +3 一色清羽 2026-04-02 4/200 2026-04-03 10:22 by 蓝云思雨
[考研] 326求调剂 +10 崽崽仔 2026-04-02 10/500 2026-04-03 09:08 by 帕尔马拉特
[考研] 求调剂推荐 +3 南山南@ 2026-04-01 3/150 2026-04-02 12:09 by xiaoranmu
[考研] 一志愿北交大材料工程总分358 +8 cs0106 2026-04-01 9/450 2026-04-02 10:36 by 不吃魚的貓
[考研] 279求调剂 +7 莫xiao 2026-04-01 7/350 2026-04-01 22:05 by 客尔美德
[考研] 322求调剂 +8 三水sss 2026-04-01 8/400 2026-04-01 10:19 by 唐沐儿
[考研] 288资源与环境专硕求调剂,不限专业,有学上就行 +25 lllllos 2026-03-30 26/1300 2026-04-01 09:52 by 一只好果子?
[考研] 求调剂,一志愿北林食品与营养095500,301分,已过六级,有科研经历 +4 快乐储蓄罐 2026-03-31 4/200 2026-04-01 09:26 by JourneyLucky
[考研] 335求调剂 +3 321* 2026-03-31 4/200 2026-04-01 00:00 by 321*
[考研] 318求调剂 +10 陈晨79 2026-03-30 10/500 2026-03-31 17:37 by 544594351
[考研] 279求调剂 +4 蝶舞轻绕 2026-03-29 4/200 2026-03-29 09:45 by laoshidan
信息提示
请填处理意见