第一次写博客，最近老师布置的地震工程学的大作业，要求编写几种不同的数值积分方法求解地震反应谱，到目前只做了前两个小题，以下是用matlab写的几个不同积分法的函数，可以比较几种积分法的结果，与大家分享，不足之处请读者慷慨指出。
http://s9/middle/87d82d08gc0ff6b0071d8&690


（1）精细积分法

function [Umax,Vmax,Amax]= JXJF(t,g,ksi,T0)
dt=0.02;
du0=0;
u0=0;
Y=[du0;u0];
%YY(1,1)=Y(1,1)
%YY(2,1)=Y(2,1)
du(1)=du0;
u(1)=u0;
LT=length(T0);
lt=length(t)-1;
for j=1:LT
w(j)=6.28/T0(j);
%计算迭代参数
%首先定义矩阵A,注意力向量Fe（t）是变化的，不在循环外定义。
A=[-2*ksi*w(j),-w(j)^2;1,