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