python连续小波变换CWT

赵孝正

已于 2022-10-06 07:41:55 修改

阅读量1.2k

点赞数 1

分类专栏：小波分析文章标签： python 开发语言 numpy

于 2022-09-28 19:16:59 首次发布

原文链接：https://blog.csdn.net/weixin_44566230/article/details/120958282?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522166435364116782425185501%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fall%255Fv2.%2522%257D&request_id=1664353641167824251

版权

小波分析专栏收录该内容

23 篇文章 0 订阅 ¥299.90 ¥99.00

订阅专栏

超级会员免费看

// An highlighted block
t = np.arange(0, len(x_f)/fs, 1/fs)
wavename = 'morl'
totalscal = 2000
fc = pywt.central_frequency(wavename)
cparam = 2 * fc * totalscal
scales = cparam / np.arange(totalscal, 1, -1)

[cwtmatr, frequencies] = pywt.cwt(x_f, scales, wavename, 1.0 / fs)
# fig
plt.figure(figsize=(20, 10))
plt.contourf(t, frequencies, abs(cwtmatr))
plt.ylabel('hz')
plt.xlabel('s')
plt.subplots_adjust(hspace=0.4)
plt.show()

连续小波变换中可用的小波：
**A wide range of continous wavelets are now available. These include the following:**
	Gaussian wavelets (gaus1…``gaus8``)
	Mexican hat wavelet (mexh)
	Morlet wavelet (morl)
	Complex Gaussian wavelets (cgau1…``cgau8``)
	Shannon wavelet (shan)
	Frequency B-Spline wavelet (fbsp)
	Complex Morlet wavelet (cmor)