| 查看: 2637 | 回复: 3 | |||
glgjw银虫 (小有名气)
|
[求助]
MATLAB求解二阶常微分方程初值问题 已有1人参与
|
|
主程序如下： >> global m1 m2 k1 c1 c2 v0 >> m1=4080;m2=12.5e3;k1=7.35e6;c1=17.56e3;c2=0.225e3;v0=3.2/60 >> [t,y]=ode45(@doublemass,[0 5],[-0.0154;-0.0049;0;0]) 其中的doublemass的函数M文件定义如下： function dyy=doublemass(t,y) global m1 m2 k1 c1 c2 v0 dyy=[y(3);(-4.48e6*v0*t/(6-y(1))-(k1+4.48e6/(6-y(1)))*y(1)-(c1+c2)*y(2)+c2*y(4)+4.48e6*y(2)/(6-y(1)))/m1;y(4);(4.48e6*v0*t/(6-y(1))-m2*10+c2*y(3)+4.48e6*y(1)/(6-y(1))-c2*y(4)-4.48e6*y(2)/(6-y(1)))/m2] 但计算出的结果却不对，前面两个初值不是前面设置好的-0.0154；-0.0049，而是 -0.0000 -0.0000，如下图： y = 1.0e+009 * -0.0000 -0.0000 0 0 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0001 -0.0000 -0.0000 -0.0000 0.0002 -0.0000 -0.0001 -0.0000 0.0002 -0.0000 -0.0001 -0.0000 0.0002 -0.0000 -0.0001 -0.0000 0.0003 -0.0000 -0.0001 -0.0000 0.0003 -0.0000 -0.0001 -0.0000 0.0004 -0.0000 -0.0001 -0.0000 0.0005 -0.0000 -0.0002 -0.0000 0.0005 -0.0000 -0.0002 -0.0000 0.0006 -0.0000 -0.0002 -0.0000 0.0006 -0.0000 -0.0002 -0.0000 0.0007 -0.0000 -0.0002 -0.0000 0.0007 -0.0000 -0.0002 -0.0000 0.0007 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0008 -0.0000 -0.0003 -0.0000 0.0010 -0.0001 -0.0003 -0.0000 0.0013 -0.0001 -0.0003 -0.0000 0.0017 -0.0001 -0.0004 -0.0000 0.0022 -0.0001 -0.0004 -0.0000 0.0030 -0.0001 -0.0004 -0.0000 0.0041 -0.0001 -0.0004 -0.0000 0.0054 -0.0001 -0.0004 -0.0000 0.0069 -0.0002 -0.0004 -0.0000 0.0094 -0.0002 -0.0004 -0.0001 0.0123 -0.0002 -0.0004 -0.0001 0.0157 -0.0002 -0.0004 -0.0001 0.0197 -0.0003 -0.0004 -0.0001 0.0261 -0.0003 -0.0004 -0.0001 0.0335 -0.0003 -0.0004 -0.0002 0.0420 -0.0003 -0.0004 -0.0002 0.0516 -0.0004 -0.0004 -0.0002 0.0678 -0.0004 -0.0004 -0.0003 0.0864 -0.0005 -0.0004 -0.0004 0.1076 -0.0005 -0.0004 -0.0004 0.1314 -0.0006 -0.0005 -0.0005 0.1585 -0.0007 -0.0005 -0.0006 0.1885 -0.0007 -0.0005 -0.0007 0.2214 -0.0008 -0.0005 -0.0008 0.2575 -0.0008 -0.0005 -0.0009 0.2967 -0.0009 -0.0005 -0.0010 0.3392 -0.0010 -0.0006 -0.0011 0.3851 -0.0010 -0.0006 -0.0013 0.4345 -0.0011 -0.0006 -0.0014 0.4875 -0.0012 -0.0006 -0.0016 0.5441 -0.0013 -0.0006 -0.0017 0.6045 -0.0013 -0.0006 -0.0019 0.6688 -0.0014 -0.0007 -0.0021 0.7372 -0.0015 -0.0007 -0.0023 0.8096 -0.0016 -0.0007 -0.0025 0.8863 -0.0017 -0.0007 -0.0027 0.9673 -0.0018 -0.0007 -0.0029 1.0528 -0.0019 -0.0007 -0.0032 1.1428 -0.0020 -0.0008 -0.0034 1.2376 -0.0021 -0.0008 -0.0037 1.3371 -0.0021 -0.0008 -0.0040 1.4416 -0.0023 -0.0008 -0.0043 1.5512 -0.0024 -0.0008 -0.0046 1.6659 -0.0025 -0.0009 -0.0049 1.7859 -0.0026 -0.0009 -0.0051 1.8795 -0.0027 -0.0009 -0.0054 1.9762 -0.0027 -0.0009 -0.0056 2.0761 -0.0028 -0.0009 -0.0059 2.1792 -0.0029 -0.0009 不知是何原因？ |
回复此楼» 猜你喜欢
关于filecode，很负责任的告诉大家
已经有5人回复
-
关于豆爷回答的JTJC与%2F数量
已经有7人回复
-
2026国自然放榜时间
已经有3人回复
-
国基金的申报应该改成非等额制，评价高的钱多评价低的钱少，但是增加资助率
已经有7人回复
-
好奇怪的filecode
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
wurongjun
专家顾问 (职业作家)
-
专家经验: +831 - 数学EPI: 9
- 应助: 791 (博后)
- 贵宾: 0.308
- 金币: 24609
- 散金: 310
- 红花: 75
- 帖子: 3004
- 在线: 881.4小时
- 虫号: 1368482
- 注册: 2011-08-14
- 性别: GG
- 专业: 计算数学与科学工程计算
- 管辖: 数学
2楼2016-02-24 17:47:24
|
本帖内容被屏蔽 |
3楼2016-02-25 10:00:29
glgjw
银虫 (小有名气)
- 应助: 2 (幼儿园)
- 金币: 252.4
- 帖子: 232
- 在线: 93.2小时
- 虫号: 2115194
- 注册: 2012-11-08
- 性别: GG
- 专业: 信息处理方法与技术
4楼2016-02-25 11:11:56
20