基于反应谱计算的基本原理,对于单自由度弹性体系在水平地震运动xg(t)作用下的动力平衡微分方程,分别使用gradient和diff进行计算,并与使用cumtrapz命令积分计算的精确计算反应谱积分公式和简化的反应谱计算积分公式值进行了比较,xg(t)分别选取了实测的爆炸震动数据和天然地震数据,经过计算发现,对于高频比较显著的爆炸震动波,用上述方法计算出的反应谱在低频部分有小的差别,对反应谱的整体特征影响不大.但对于低频比较显著的天然地震波,用cumtrapz积分计算的积分公式与用diff命令计算