大家好,欢迎来到IT知识分享网。
if __name__ == '__main__': A = __LoadData__('dataSource2.txt', 200, 1500) __Draw__(A) A =__SingalToHankel__(A[:20,0]) print(A)
附上部分代码:
if __name__ == '__main__': A = __LoadData__('dataSource2.txt',200, 1500) __Draw__(A) #A =__SingalToHankel__(A[0,:20]) #print(A) B,S = __getFFT__(A[0:1,:20]) #行是道号,列是到达检波器时间 print(B) plt.plot(B[0]) plt.show() x = range(20) __DrawFrequencyDomain__(A[1,0:20], x) plt.show()
后续将继续实现其他相关属性…
# 快速傅里叶变换 fft_A = fft(A) print(fft_A)
def Chebyshevmap1(x0, max_g): """ Chebyshev 映射 :param x0: 初值 :param max_g:最大迭代次数 :return: x的列表 """ x = x0 x_list = [] for i in max_g: x = math.cos(i * math.cos(x) (-1)) x_list.append(x) return x_list
原始波形
(加速度谱、位移谱和速度谱,不知道对不对)
免责声明:本站所有文章内容,图片,视频等均是来源于用户投稿和互联网及文摘转载整编而成,不代表本站观点,不承担相关法律责任。其著作权各归其原作者或其出版社所有。如发现本站有涉嫌抄袭侵权/违法违规的内容,侵犯到您的权益,请在线联系站长,一经查实,本站将立刻删除。 本文来自网络,若有侵权,请联系删除,如若转载,请注明出处:https://haidsoft.com/129406.html