求助！无法得到共振现象！极度疑惑！！！

pureriver

用S-function模块 + 相应的m文件搭了一个简单的单自由度欠阻尼二阶振动系统（wn=1000，Q=50），然后用signal generator 模块制造振幅为1的正弦激励（外界强迫力），观察输出（振子的位移），结果发现只有当激励频率较小时（<100Hz），输出才符合理论计算（振幅与wn的平方成反比，相位滞后基本为0），而再提高频率输出就与预想完全不符（比如，在谐振点上振幅比低频100Hz时小了2个数量级，相位滞后不是理论上的90度，而是达到了180，完全看不出高Q二阶系统幅频特性应有的峰值特性，相频特性也是如此）。试过各种solver 方法，情况毫无变化。因此怀疑是不是在仿真此类高Q系统的谐振特性时，solver所用的离散化积分方法会产生失真？（记得连续域向离散域映射时，对变换频率好像有限制）特来请教大牛。谢谢！！！