振型分解反应谱法matlab