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

登录注册

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

24小时热门版块排行榜

返回列表

zhuang1017

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 758
散金: 50
帖子: 63
在线: 38.6小时
虫号: 4078324
注册: 2015-09-17
性别: GG
专业: 机械动力学

[求助] matlab 编程求助已有1人参与

各位前辈，matlab编程求助，关于碰撞振动的。
g(1)=0.02;
g(2)=0.02;
u=0.02;
odefunhdb = @(t,y)[y(2); (-2*g(1)*y(2)-y(1)+2*u*g(2)*(y(4)-y(2))+u*(y(3)-y(1))+cos(t));...
y(4); (-2*g(2)*(y(4)-y(2))-(y(3)-y(1)))];%把动力学方程中的二阶微分换成一阶
h = 0.05;  t = (0:h:100)'; y = zeros(length(t),4); y(1,: )=[0,0,0,0];%参数设置与初值
for n = 1:length(t)-1
   k1 = odefunhdb(t(n),y(n,: )');
   k2 = odefunhdb(t(n)+h/2,y(n,: )'+h/2*k1);
   k3 = odefunhdb(t(n)+h/2,y(n,: )'+h/2*k2);
   k4 = odefunhdb(t(n)+h,y(n,: )'+h*k3);
   y(n+1,: ) = y(n,: )+h/6*(k1+2*k2+2*k3+k4)';

end

plot(t,y(:,1))
这是没有发生碰撞时的动力学方程编程，现需要加入一个发生碰撞的条件，当在循环过程中满足y(3)-y(1)=0.8+1e-4时，有
y(2)+=0.9765*y(2)-+0.235*y(4)-
y(4)+=1.1765*y(2)--0.1765*y(4)-
其中+表示变化后的，-表示变化之前，然后用此刻的y(1),y(2)+,y(3),y(4)+作为初值继续在龙格库塔计算中循环，最后需要得到所有的y(1)

回复此楼

» 猜你喜欢

1楼 2017-06-15 10:42:12

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

chendequan

铁虫 (小有名气)

应助: 152 (高中生)
金币: 918.6
红花: 17
帖子: 262
在线: 208.6小时
虫号: 994474
注册: 2010-04-11
专业: 结构工程

【答案】应助回帖

感谢参与，应助指数 +1

内容已删除

赞一下(1人)

回复此楼

QQ：516477448，真心帮助解决MATLAB相关问题，提供详细资料，Word文档明确具体问题及要求，尽力而为！

2楼2017-06-17 15:09:36