小波阈值去噪python代码

以下是一个使用小波阈值去噪的Python代码示例：

import pywt
import numpy as np

# 加载含有噪声的信号
noisy_signal = np.array([...])

# 将信号分解为多个小波系数
coeffs = pywt.wavedec(noisy_signal, 'db1', level=6)

# 将噪声系数转换为标准差，并乘以指定的阈值
threshold = 0.5
sigma = np.median(np.abs(coeffs[-1])) / 0.6745
threshold_value = sigma * threshold

# 将所有小波系数转换为零，并保留一个子带用于重构
thresholded_coeffs = []
for i, coeff in enumerate(coeffs):
    if i == 0:
        thresholded_coeffs.append(coeff)
    else:
        thresholded_coeffs.append(pywt.threshold(coeff, threshold_value, 'soft'))
        
# 使用修改后的小波系数重构信号
denoised_signal = pywt.waverec(thresholded_coeffs, 'db1')

注意：这只是一个示例代码，实际应用中需要根据数据和需求进行调整和优化。