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

最后

以上就是朴实牛排为你收集整理的python插值变换地震波间隔：的全部内容，希望文章能够帮你解决python插值变换地震波间隔：所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。