小波函数绘制（Python）

import pywt
import matplotlib.pylab as plt
plt.rcParams.update({"font.size": 12, "font.weight": 'bold', "lines.linewidth": 1.5})
# Plot discrete Mayer wavelet function
phi, psi, x = pywt.Wavelet('dmey').wavefun(level = 5)
plt.figure(figsize=(4, 3), dpi=100, layout='constrained')
plt.plot(x, psi, '-b', label = 'Wavelet')
# plt.plot(x, phi, '-b', label = 'Scaling')
plt.xlim([20, 40])
#plt.legend(frameon = False, loc = 'upper right')
plt.title('Meyer Wavelet', fontsize=14, fontweight='bold')
plt.xticks([])
plt.yticks([])
# plt.savefig('scaldmey.png', dpi=400, bbox_inches='tight', pad_inches=0.02)
plt.show()

# Coiflet wavelet
phi, psi, x = pywt.Wavelet('coif3').wavefun(level = 5)


plt.figure(figsize=(4, 3), dpi=100, layout='constrained')
# plt.plot(x, psi, '-b', label = 'Wavelet')
plt.plot(x, phi, '-b', label = 'Scaling')
#plt.text(0.005, 1.14, "(c)", fontdict = font)
#plt.legend(frameon = False, loc = 'upper right')
plt.title('Coif3: Scaling Function', fontsize=14, fontweight='bold')
plt.xlim([0, 16])
# plt.savefig('scalcoif3.png', dpi=400, bbox_inches='tight', pad_inches=0.02)
plt.show()

# Daubechies wavelet
phi, psi, x = pywt.Wavelet('db6').wavefun(level = 5)


plt.figure(figsize=(4, 3), dpi=100, layout='constrained')
# plt.plot(x, psi, '-b', label = 'Wavelet')
plt.plot(x, phi, '-b',label = 'Scaling')
#plt.text(0.005, 1.14, "(c)", fontdict = font)
#plt.legend(frameon = False, loc = 'upper right')
plt.title('Db6: Scaling Function', fontsize=14, fontweight='bold')
plt.xlim([0, 11])
plt.savefig('scaldb6.png', dpi=400, bbox_inches='tight', pad_inches=0.02)


plt.show()

# Symlet wavelet
phi, psi, x = pywt.Wavelet('sym10').wavefun(level = 5)


plt.figure(figsize=(4, 3), dpi=100, layout='constrained')
plt.plot(x, psi, '-b', label = 'Wavelet')
# plt.plot(x, phi, '-b', label = 'Scaling')
#plt.text(0.005, 1.14, "(c)", fontdict = font)
#plt.legend(frameon = False, loc = 'upper right')
plt.title('Symlet Wavelet', fontsize=14, fontweight='bold')
plt.xlim([4, 16])
# plt.savefig('scalsym10.png', dpi=400, bbox_inches='tight', pad_inches=0.02)
plt.show()

phid, psid, xd = pywt.Wavelet('db6').wavefun(level = 5) # Daubechies wavelet
phic, psic, xc = pywt.Wavelet('coif3').wavefun(level = 5) # Coiflet wavelet
phis, psis, xs = pywt.Wavelet('sym10').wavefun(level = 5) # Symlet wavelet
phim, psim, xm = pywt.Wavelet('dmey').wavefun(level = 5) # Meyer wavelet


fig, ax = plt.subplots(2, 2)
ax[0, 0].plot(xd, psid, '-b')
ax[0, 1].plot(xc, psic, '-b')
ax[1, 0].plot(xs, psis, '-b')
ax[1, 1].plot(xm, psim, '-b')


plt.show()

工学博士，担任《Mechanical System and Signal Processing》审稿专家，担任《中国电机工程学报》优秀审稿专家，《控制与决策》，《系统工程与电子技术》，《电力系统保护与控制》，《宇航学报》等EI期刊审稿专家，擅长领域：现代信号处理，机器学习，深度学习，数字孪生，时间序列分析，设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。
知乎学术咨询：https://www.zhihu.com/consult/people/792359672131756032?isMe=1
工学博士，担任《Mechanical System and Signal Processing》审稿专家，担任《中国电机工程学报》优秀审稿专家，《控制与决策》，《系统工程与电子技术》，《电力系统保护与控制》，《宇航学报》等EI期刊审稿专家，擅长领域：现代信号处理，机器学习，深度学习，数字孪生，时间序列分析，设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。