说明：利用newmark法求动力响应的程序，可以参考。
说明：用matlab编写的newmark－belta法求解结构动力响应
说明：NEWMARK方法求解结构动力学源程序，可求解较大自由度时的响应。
说明：该计算程序可以对无固定约束碰撞振动系统在激励力作用 下的系统响应进行计算，模拟系统动态响应，并自动判断颤振发生情 况，记录不同参数组合（阻尼比，质量比，恢复系数，间隙，频率比） 情况下系统的响应数据，输出发生颤振的参数组合数据。
说明：本程序主要是针对一九杆结构进行的动力学分析，用newmark法，威尔逊法，状态空间法，子空间迭代法等求解了结构动态响应，并给出了图像显示
说明：Wilson法求解动力方程，可以画出节点的位移响应。