python程序精确法求解反应谱,傅里叶谱，功率谱

最新推荐文章于 2024-06-05 18:31:46 发布

踏梦拾青

最新推荐文章于 2024-06-05 18:31:46 发布

阅读量4.3k

点赞数 3

分类专栏：地震动处理

本文链接：https://blog.csdn.net/qq_36145163/article/details/82939888

版权

这篇博客介绍了如何使用Python进行地震动处理，包括傅里叶谱和功率谱的计算，以及反应谱的求解。内容来源于MATLAB代码的算法转换，详细解释可在链接中查看。

摘要由CSDN通过智能技术生成

主要用于地震动处理种的代码：

傅里叶谱

    def myfft(inp,delta): #输入的x为竖向的array
        from scipy.fftpack import fft,ifft
        inp=np.array(inp).reshape(-1,1)
        inpf=np.fft.fft(inp,axis=0) #默认的是对行向量fourier变换，所以此处要定义下轴
        inpf=inpf*2/inp.shape[0]
        inpf2=inpf[0:int((inp.shape[0]+1)/2)]#取一半
        n_points=inpf2.shape[0]#地震波点数
        frequency=1/delta
        fseries=np.linspace(0,frequency/2,n_points)
        #     plt.plot(fseries,inpf2)
        return inpf2,fseries

功率谱

    def mypowverspectra(inp,delta):
        from scipy.fftpack import fft,ifft
        inp=np.array(inp).reshape(-1,1)
        inpf=np.fft.fft(inp,axis=0) #默认的是对行向量fourier变换，所以此处要定义下轴
        inpf=inpf*2/inp.shape[0]
        inpf2=inpf[0:int((inp.shape[0]+1)/2)]#取一半
        inpf2=np.abs(inpf2)**2
        n_points=inpf2.shape[0]#地震波点数