python插值变换地震波间隔：

采用插值法变换地震波间隔，把原本不等间距的地震动记录插值方法变换成等间距的地震动记录：

def interpolation(acc,time,timestep):
#插值方法是两点之间插值。acc是原时程，time是原非等间距时间，timestep是处理过后的时间间距
    t=0
    timeidx=0
    acc2idx=0
    acc2=np.zeros([int(time[-1]/timestep)+1,1])
    while t< time[-1]:

        if t == time[timeidx]:
            print(t)
            acc2[acc2idx]=acc[timeidx]
        else:
            while not ((t>time[timeidx] and t< time[timeidx+1]) or (t == time[timeidx])):
                if t<time[timeidx]:
                    timeidx=timeidx-1
                if t>time[timeidx+1]:
                    timeidx=timeidx+1
                a1=acc[timeidx]
                a2=acc[timeidx+1]
                t1=time[timeidx]
                t2=time[timeidx+1]
                acc2[acc2idx]=(a1*t2-a1*t+a2*t-a2*t1)/(t2-t1)#两点间插值
                if t== time[timeidx+1]:
                    timeidx=timeidx+1
        t=t+timestep
        acc2idx=acc2idx+1   
        timeseries=np.array(range(len(acc2)))*timestep
    return acc2,timeseries