单自由度受迫振动

单自由度系统的受迫振动在简谐激振力0()cos()FtFt的作用下，系统的运动方程为：0cos()mxcxkxFt或22cos()nnxxxht其中2/ncm2/nkm0/hFm对系统参数取：m=1,k=9,0F=1并令c分别等于6、4.5、3、1.5、0.6、0或：n=3，h=1，/2,=3并令分别等于1、0.75、0.5、0.25、0.1、0带入参数：sin(3)69xtxx假定有初始条件：4x，2x画得Simulink程序框图：初始值设为4初始值设为2c/mk/mhSineWaveScope1sIntegrator11sIntegrator1Gain29Gain16GainAdd分别设置c/m等于6、4.5、3、1.5、0.6、0运行后分别得到以下仿真图（其中对c/m等于0.6和0的仿真图扩大了时间轴长度）：10.750.50.250.10本系统设计的频率比恒为/1n，通过改变阻尼比的大小来观测系统响应的变化。通过对比，系统的位移响应在初始很短的一段时间内随阻尼比的改变而变化的范围很小，主要原因在于系统有初始的速度、位移的影响；但在后续时间内，随着阻尼比的减小，系统的位移响应的振荡幅值会逐渐增大；对于0.1和0时系统的响应振幅先减小后增大，但0.1时后续的响应振幅接近稳定，而无阻尼时，系统响应振幅会逐步放大。附：若m文件打不开，可复制一下链接在MATLAB的commandwindow运行即可；bdcloseall;set_param(0,'CharacterEncoding','ISO-8859-1')