发动机悬置系统强迫振动响应计算问题！求助！

hover

某发动机悬置系统，已知[K],[M],以及结构阻尼因子c，以及简谐激振力[F]sin(wt+β），
          即可得到系统的强迫振动微分方程。现想求系统的振动速度响应。
          我的思路是这样的，第一，先求解强迫振动微分方程，得到系统的震动位移响应；
                                          第二通过数值微分得到系统的震动速度响应。
          但算出来的结果很大，觉得不对，不知道数值微分是否可行，因为方程比较复杂，存在耦合，无法计算系统的解析解。
          请教各位高手，有什么方便可行的思路，向小弟传授一下！不胜感激呀！！！！！！！！！