版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

>论坛更新日志 (5376)
>基金申请 (1143)
>虫友互识 (749)
>论文投稿 (425)
>文献求助 (381)
>导师招生 (310)
>考博 (271)
>休闲灌水 (183)
>硕博家园 (157)
>博后之家 (132)
>健康生活 (83)
>考研 (70)
>留学DIY (50)
>公派出国 (48)
>绿色求助(高悬赏) (41)
>学术会议 (40)

返回列表

glgjw

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 252.4
帖子: 232
在线: 93.2小时
虫号: 2115194
注册: 2012-11-08
性别: GG
专业: 信息处理方法与技术

[求助] MATLAB求解二阶常微分方程组出错已有1人参与

MATLAB求解一个二自由度振动微分方程组，主程序如下：
>> 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 10],[-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(5))-(k1+4.48e6/(6-y(5)))*y(5)-(c1+c2)*y(6)+c2*y(4)+4.48e6*y(6)/(6-y(5)))/m1;y(4);(4.48e6*v0*t/(6-y(5))-m2*10+c2*y(3)+4.48e6*y(5)/(6-y(5))-c2*y(4)-4.48e6*y(6)/(6-y(5)))/m2]

运行后出现如下错误：
??? Attempted to access y(5); index out of bounds because numel(y)=4.

Error in ==> doublemass at 3
dyy=[y(3);(-4.48e6*v0*t/(6-y(5))-(k1+4.48e6/(6-y(5)))*y(5)-(c1+c2)*y(6)+c2*y(4)+4.48e6*y(6)/(6-y(5)))/m1;y(4);(4.48e6*v0*t/(6-y(5))-m2*10+c2*y(3)+4.48e6*y(5)/(6-y(5))-c2*y(4)-4.48e6*y(6)/(6-y(5)))/m2
Error in ==> funfun\private\odearguments at 110
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.

Error in ==> ode45 at 173
[neq, tspan, ntspan, next, t0, tfinal, tdir, y0, f0, odeArgs, ...

请教各位大神，谢谢！

回复此楼

» 猜你喜欢

今天务委会开完了，明天出结果吗已经有21人回复
2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？已经有14人回复
在坚冰还盖着北海的时候，我看到了怒放的梅花。已经有9人回复
明天应该可查了！？已经有6人回复
有没有大神帮我看看基金代码已经有10人回复
放榜前的不淡定已经有19人回复
某些机构，以效率低为荣，以效率低作为存在感已经有10人回复
2026国自然函评费到账已经有23人回复
今日不放榜？网传国自然预计 8 月 27 日可查结果已经有22人回复
什么时候开奖？已经有13人回复

1楼 2016-02-23 17:13:29

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

wurongjun

专家顾问 (职业作家)

专家经验: +831
数学EPI: 9
应助: 791 (博后)
贵宾: 0.308
金币: 24609
散金: 310
红花: 75
帖子: 3004
在线: 881.4小时
虫号: 1368482
注册: 2011-08-14
性别: GG
专业: 计算数学与科学工程计算
管辖: 数学

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
glgjw: 金币+20, ★★★★★最佳答案 2016-02-24 16:47:31

初值的个数不对!
你的y有6个分量!

赞一下

回复此楼

» 本帖已获得的红花（最新10朵）

glgjw

善恶到头终有报,人间正道是沧桑.

2楼2016-02-23 18:20:57

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

glgjw

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 252.4
帖子: 232
在线: 93.2小时
虫号: 2115194
注册: 2012-11-08
性别: GG
专业: 信息处理方法与技术

引用回帖:

2楼: Originally posted by wurongjun at 2016-02-23 18:20:57
初值的个数不对!
你的y有6个分量!

但是我的y的编号是从3到6，y的数目是4个，不是6个啊。难道一定要写成y(1)到y(4)?

发自小木虫Android客户端

赞一下

回复此楼

3楼2016-02-23 18:56:06

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

glgjw

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 252.4
帖子: 232
在线: 93.2小时
虫号: 2115194
注册: 2012-11-08
性别: GG
专业: 信息处理方法与技术

送红花一朵

引用回帖:

2楼: Originally posted by wurongjun at 2016-02-23 18:20:57
初值的个数不对!
你的y有6个分量!

感谢wurongjun的回复，已算出结果。但结果似乎不太对，y(1)和y(2)我指定的初值是-0.0154和-0.0049，但计算出的结果显示是从0开始迭代计算的。如下：
y =

1.0e+010 *

-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.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.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.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.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.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.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.0000
-0.0000 0.0002 -0.0000 -0.0000
-0.0000 0.0003 -0.0000 -0.0000
-0.0000 0.0004 -0.0000 -0.0000
-0.0000 0.0005 -0.0000 -0.0000
-0.0000 0.0007 -0.0000 -0.0000
-0.0000 0.0009 -0.0000 -0.0000
-0.0000 0.0012 -0.0000 -0.0000
-0.0000 0.0016 -0.0000 -0.0000
-0.0000 0.0020 -0.0000 -0.0000
-0.0000 0.0026 -0.0000 -0.0000
-0.0000 0.0034 -0.0000 -0.0000
-0.0000 0.0042 -0.0000 -0.0000
-0.0000 0.0052 -0.0000 -0.0000
-0.0000 0.0068 -0.0000 -0.0000
-0.0000 0.0086 -0.0000 -0.0000
-0.0000 0.0108 -0.0001 -0.0000
-0.0000 0.0131 -0.0001 -0.0000
-0.0001 0.0176 -0.0001 -0.0000
-0.0001 0.0228 -0.0001 -0.0001
-0.0001 0.0288 -0.0001 -0.0001
-0.0001 0.0357 -0.0001 -0.0001
-0.0001 0.0455 -0.0001 -0.0001
-0.0002 0.0567 -0.0001 -0.0001
-0.0002 0.0695 -0.0001 -0.0001
-0.0002 0.0839 -0.0002 -0.0001
-0.0003 0.1000 -0.0002 -0.0001
-0.0003 0.1179 -0.0002 -0.0001
-0.0004 0.1378 -0.0002 -0.0001
-0.0004 0.1596 -0.0002 -0.0001
-0.0005 0.1833 -0.0003 -0.0001
-0.0006 0.2094 -0.0003 -0.0001
-0.0006 0.2380 -0.0003 -0.0001
-0.0007 0.2689 -0.0003 -0.0001
-0.0008 0.2923 -0.0004 -0.0001
-0.0008 0.3172 -0.0004 -0.0001
-0.0009 0.3436 -0.0004 -0.0001
-0.0010 0.3714 -0.0004 -0.0001
-0.0011 0.4005 -0.0004 -0.0001
-0.0011 0.4312 -0.0005 -0.0001
-0.0012 0.4635 -0.0005 -0.0001
-0.0013 0.4974 -0.0005 -0.0001
-0.0014 0.5388 -0.0005 -0.0001
-0.0015 0.5825 -0.0006 -0.0001
-0.0016 0.6289 -0.0006 -0.0001
-0.0018 0.6777 -0.0006 -0.0001
-0.0019 0.7341 -0.0006 -0.0002
-0.0021 0.7941 -0.0007 -0.0002
-0.0022 0.8582 -0.0007 -0.0002
-0.0024 0.9255 -0.0008 -0.0002
-0.0026 0.9896 -0.0008 -0.0002
-0.0027 1.0576 -0.0008 -0.0002
-0.0029 1.1307 -0.0009 -0.0002
-0.0031 1.2066 -0.0009 -0.0002
-0.0032 1.2536 -0.0009 -0.0002
-0.0034 1.3022 -0.0010 -0.0002
-0.0035 1.3524 -0.0010 -0.0002
-0.0036 1.4040 -0.0010 -0.0002
而且通过绘图命令如下
>> plot(t,y(:,1),t,y(:,2),t,y(:,3),t,y(:,4))
>> legend('1位移','1速度','2位移','2速度')
执行后，为何只出现y(2)和y(4)的图线？
请专家赐教。我会追加金币的。谢谢了！！

赞一下

回复此楼

4楼2016-02-24 09:07:06

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

glgjw

银虫 (小有名气)

应助: 2 (幼儿园)
金币: 252.4
帖子: 232
在线: 93.2小时
虫号: 2115194
注册: 2012-11-08
性别: GG
专业: 信息处理方法与技术

引用回帖:

2楼: Originally posted by wurongjun at 2016-02-23 18:20:57
初值的个数不对!
你的y有6个分量!

doublemass文件改后的内容是：
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]

赞一下

回复此楼

5楼2016-02-24 09:10:47

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

wurongjun

专家顾问 (职业作家)

专家经验: +831
数学EPI: 9
应助: 791 (博后)
贵宾: 0.308
金币: 24609
散金: 310
红花: 75
帖子: 3004
在线: 881.4小时
虫号: 1368482
注册: 2011-08-14
性别: GG
专业: 计算数学与科学工程计算
管辖: 数学

引用回帖:

4楼: Originally posted by glgjw at 2016-02-24 09:07:06
感谢wurongjun的回复，已算出结果。但结果似乎不太对，y(1)和y(2)我指定的初值是-0.0154和-0.0049，但计算出的结果显示是从0开始迭代计算的。如下：
y =

1.0e+010 *

-0.0000 -0.0000 0 ...

后面的覆盖了前面的!
你分开画四个图形就明白了!

赞一下

回复此楼

善恶到头终有报,人间正道是沧桑.

6楼2016-02-24 17:40:40

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

山东赵子明

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2274.7
散金: 10
帖子: 119
在线: 45.9小时
虫号: 3722737
注册: 2015-03-09
专业: 流体力学

您好在么，我有问题请教您

发自小木虫IOS客户端

赞一下

回复此楼

7楼2016-11-29 17:53:41

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

波波和倩倩

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 2056.8
散金: 38
帖子: 111
在线: 96.7小时
虫号: 4232453
注册: 2015-11-19
性别: GG
专业: 力学中的基本问题和方法

你好，可以请教下您吗？关于二阶常微分方程组的问题

发自小木虫Android客户端

赞一下

回复此楼

8楼2016-12-17 00:25:27

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

相关版块跳转我要订阅楼主 glgjw 的主题更新

返回列表

最具人气热帖推荐 [查看全部]		作者	回/看	最后发表

[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +12	木水思豆 2026-08-25	14/700	2026-08-25 22:35 by yys300
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6	ziyangfang 2026-08-25	9/450	2026-08-25 20:26 by huagongfeihu
[基金申请] 明天应该可查了！？ +6	chengyan1220 2026-08-23	6/300	2026-08-25 19:45 by zfd97
[基金申请] 有没有大神帮我看看基金代码 28+4	1234567wang 2026-08-24	10/500	2026-08-25 19:15 by lfy8008
[基金申请] 放榜前的不淡定 20+4	snowwithsea 2026-08-19	19/950	2026-08-25 17:57 by chick875
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9	yuleib84 2026-08-25	10/500	2026-08-25 17:14 by alexon
[基金申请] 2026国自然函评费到账 +20	羊腰板 2026-08-21	23/1150	2026-08-25 16:34 by zsna
[基金申请] 什么时候开奖？ +12	CrisMessi 2026-08-18	13/650	2026-08-25 14:40 by ZJTJZ
[基金申请] 人气不行了 +11	fansofjerry 2026-08-21	11/550	2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3	苏东坡二世 2026-08-23	3/150	2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5	process2012 2026-08-23	7/350	2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13	且听虎啸 2026-08-20	14/700	2026-08-25 09:10 by mrkang
[基金申请] 能否退出参与的面上项目解除限项 +21	koalala 2026-08-24	24/1200	2026-08-24 19:25 by 家与远方
[基金申请] 建议基金发布提前给出明确的时间点 +13	kulium 2026-08-21	16/800	2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13	大萍1987 2026-08-20	16/800	2026-08-24 10:23 by 太傻了
[基金申请] 科研孤儿太难了 +18	我4大白菜 2026-08-20	19/950	2026-08-24 09:47 by 凯恩广盛催化分�
[基金申请] 时间戳又变了 +13	wuchongjun 2026-08-20	19/950	2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4	哈哈蛤？ 2026-08-21	4/200	2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9	yuleib84 2026-08-19	11/550	2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4	longfie172 2026-08-20	4/200	2026-08-21 08:58 by mark mao

24小时热门版块排行榜

[求助] MATLAB求解二阶常微分方程组出错 已有1人参与

» 猜你喜欢

【答案】应助回帖

» 本帖已获得的红花（最新10朵）

[求助] MATLAB求解二阶常微分方程组出错已有1人参与