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

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4238)
>基金申请 (262)
>虫友互识 (179)
>文献求助 (178)
>导师招生 (116)
>考博 (90)
>硕博家园 (82)
>论文投稿 (55)
>找工作 (50)
>考研 (49)
>论文道贺祈福 (42)
>教师之家 (39)
>公派出国 (39)
>博后之家 (37)
>招聘信息布告栏 (29)
>休闲灌水 (26)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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, ...

请教各位大神，谢谢！

回复此楼

» 猜你喜欢

售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急已经有6人回复
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急已经有7人回复
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急已经有12人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有7人回复
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急已经有14人回复
面上意见出来了已经有18人回复
有多少人是今天查系统知道结果的？已经有19人回复
能否申诉？已经有4人回复
找导师已经有7人回复
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急已经有4人回复

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

已阅回复此楼关注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的回帖

查看全部 8 个回答

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个分量!

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的回帖

查看全部 8 个回答

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

[基金申请] 面上意见出来了 +9	黄鸟于飞Chao 2026-08-29	18/900	2026-08-30 16:47 by 黄鸟于飞Chao
[找工作] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	jCd0dEvKHShX 2026-08-29	4/200	2026-08-30 14:05 by GkCVyODlx4mR
[公派出国] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3	jCd0dEvKHShX 2026-08-29	5/250	2026-08-30 14:03 by GkCVyODlx4mR
[考博] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4	jCd0dEvKHShX 2026-08-29	6/300	2026-08-30 13:54 by GkCVyODlx4mR
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +6	gy1nBQXYQJqL 2026-08-29	7/350	2026-08-30 12:13 by l0VvVHGBGRLv
[博后之家] 售SCI文章，我:8O.5.5.1O.54,科目全,可十急 +8	ASdOkHsho7FD 2026-08-28	9/450	2026-08-30 11:40 by l0VvVHGBGRLv
[考博] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3	G6APbkg8SA6w 2026-08-29	4/200	2026-08-30 07:28 by ZPa0EcMwuECS
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3	gy1nBQXYQJqL 2026-08-29	4/200	2026-08-30 06:57 by ZPa0EcMwuECS
[考博] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +6	ASdOkHsho7FD 2026-08-28	10/500	2026-08-30 06:11 by ZPa0EcMwuECS
[教师之家] 售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急 +5	ASdOkHsho7FD 2026-08-28	8/400	2026-08-30 05:48 by ZPa0EcMwuECS
[基金申请] 国社科申报系统有变化 +3	kynobel 2026-08-23	3/150	2026-08-29 21:57 by 余韵清
[基金申请] 为什么到现在没收到通知？ +4	tannykie 2026-08-29	4/200	2026-08-29 17:14 by lmz0216
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +9	苏东坡二世 2026-08-23	9/450	2026-08-29 14:41 by hustersqt
[基金申请] 国自然评审意见 +13	wangmingqi 2026-08-28	19/950	2026-08-29 10:22 by Poppy1104
[基金申请] 2026年叶企孙基金 +4	bud_bud 2026-08-27	7/350	2026-08-29 07:23 by foolishmani
[文学芳草园] 梦想 +3	myrtle 2026-08-26	3/150	2026-08-27 10:01 by angelyueyi
[基金申请] 为什么国际(地区)合作与交流项目没有放榜？ 10+3	majunge000 2026-08-26	11/550	2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6	wittyboy 2026-08-26	14/700	2026-08-26 10:55 by wittyboy
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8	scalable 2026-08-24	8/400	2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9	图啦图啦 2026-08-24	10/500	2026-08-25 11:59 by 南海小哥

24小时热门版块排行榜

glgjw

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

» 猜你喜欢

glgjw

wurongjun

【答案】应助回帖

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

glgjw

glgjw

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