matlab开发-谐波激励。应用一些数值方法研究SDOF在时间域中的强迫振动。