MATLAB程序精确法求解反应谱