24小时热门版块排行榜

>论坛更新日志 (3462)
>文献求助 (332)
>虫友互识 (257)
>导师招生 (251)
>考博 (130)
>休闲灌水 (92)
>招聘信息布告栏 (75)
>硕博家园 (73)
>博后之家 (72)
>基金申请 (57)
>论文投稿 (55)
>论文道贺祈福 (46)
>考研 (45)
>公派出国 (42)
>找工作 (38)
>教师之家 (32)

返回列表

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)

回复此楼