python小波阈值去噪

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import math
import pywt 

####################一些参数和函数############
def sgn（num）:
    if（num > 0.0）:
        return 1.0
    elif（num == 0.0）:
        return 0.0
    else:
        return -1.0

###软硬阈值折衷法 a 参数
#read data
data = pd.read_csv（‘C:\\Users\\1701\Desktop\\wt03-0.2-1.csv‘）
#y_value为原信号
a=0.5
x1 = np.array（data[‘WIND_SPEED‘]）
y_values = np.array（data[‘REAL_POWER‘]）

#plt.subplot（211）
#plt.scatter（x1 y_values s=10）
#小波基的选取
w = pywt.Wavelet（‘db1‘）#选用db5小波
#ca3 cd3 cd2 cd1 = pywt.wavedec（y_values w）
maxlev = pywt.dwt_max_level（len（data） w）#最大分解级别，返回max_level。db.dec_lenx为小波的长度
coeffs = pywt.wavedec（y_values w mode=‘constant‘level=maxlev）#分解波
#recoeffs = pywt.w