说明：变步长龙格库塔法，用于解决四阶微分方程组求解析解的问题，尤其适用于两自由度间隙碰撞问题的解决。
间隙碰撞 间隙-MATLAB 间隙-碰撞 龙格库塔-间隙 碰撞微分方程
说明：该计算程序可以对无固定约束碰撞振动系统在激励力作用 下的系统响应进行计算，模拟系统动态响应，并自动判断颤振发生情 况，记录不同参数组合（阻尼比，质量比，恢复系数，间隙，频率比） 情况下系统的响应数据，输出发生颤振的参数组合数据。
颤振 Chatter 间隙碰撞 振动响应 颤振-MATLAB-pudn
