| 查看: 713 | 回复: 1 | ||
zhuang1017铁虫 (小有名气)
|
[求助]
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) |
回复此楼» 猜你喜欢
深圳大学2026年秋博士招生-物理学-活性胶体方向-高永祥课题组
已经有18人回复
-
论物质与能量的统一模型及物理现象解释
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有172人回复
-
基于基元I统一理论的数学相关应用推导
已经有0人回复
-
基元I统一理论：宇宙本质、层级演化与修炼文明的本源规律
已经有1人回复
-
基元I理论下三大核心空间现象精准推导与细节解析
已经有0人回复
-
基于基元 I 统一理论的反重力理论推导
已经有0人回复
-
基于基元I统一理论的量子力学本源推导
已经有1人回复
-
推荐一款可以AI辅助写作的Latex编辑器SmartLatexEditor，超级好用，AI润色，全免费
已经有20人回复
-
【EI|Scopus 双检索】第六届智能机器人系统国际会议（ISoIRS 2026）
已经有0人回复
-
2026年第四届电动车与车辆工程国际会议（CEVVE 2026）
已经有0人回复
chendequan
铁虫 (小有名气)
- 应助: 152 (高中生)
- 金币: 918.6
- 红花: 17
- 帖子: 262
- 在线: 208.6小时
- 虫号: 994474
- 注册: 2010-04-11
- 专业: 结构工程
2楼2017-06-17 15:09:36