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

登录注册

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

24小时热门版块排行榜

返回列表

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

请教各位大神，谢谢！

回复此楼

» 猜你喜欢

请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
真诚求助：手里的省社科项目结项要求主持人一篇中文核心，有什么渠道能发核心吗已经有6人回复
孩子确诊有中度注意力缺陷已经有14人回复
三甲基碘化亚砜的氧化反应已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有5人回复
2025冷门绝学什么时候出结果已经有3人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有4人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复
AI论文写作工具：是科研加速器还是学术作弊器？已经有3人回复

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

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

wurongjun

专家顾问 (职业作家)

专家经验: +831
数学EPI: 9
应助: 791 (博后)
贵宾: 0.308
金币: 24607.5
散金: 310
红花: 75
帖子: 3004
在线: 881.2小时
虫号: 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
金币: 24607.5
散金: 310
红花: 75
帖子: 3004
在线: 881.2小时
虫号: 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的回帖