利用MATLAB对Newmark-β法求解振动方程

1、原理如下图公式所示：其中：Newmark-β正确选择两个参数β和γ很重要；当γ>=1/2, β>=γ/2时，Newmark-β法是无条件稳定的。

3、算例演示：问题：下图所示：一无阻尼两自由度受迫振动系统，已知：m1=1kg、m2=2kg,1=k2=k3=1000N/m,系统在m2上作用有一阶跃激振力f(t)=10N，求解该系统的动力响应。

5、结果显示：在matlab命令窗口输入以上禅旄褡瘦命令，运行上述命令，就可以得到系统的位移、速度、加速度纽马克法数值解，我们用MATLAB画图工具，直观的显示出系统两泥撼费瘠质量块的位移、速度、加速图时间图，如下图所示：