Обозначать сигнал с Пивавелет?
Мне нужно подавить сигнал. Я попытался снять шум с помощью savgol_filter, но это привело к потере сингулярностей в сигнале. Чтобы устранить шум и сохранить сингулярности, я попытался использовать вейвлет-преобразование, пороговое вейвлет-преобразование и обратное вейвлет-преобразование, но мне это не удалось. Кто-нибудь знает, как использовать шум вейвлет?
вот текстовый файл с сигнальными данными
import numpy as np
from matplotlib import pyplot as plt
from scipy.signal import savgol_filter
import pywt
def readSignal(nomFichier, N):
x=np.zeros((N), dtype=float)
y=np.zeros((N), dtype=float)
fichier=open(nomFichier,"r")
for k in range(N):
l=fichier.readline()
l=l.split(' \t\t\t ',1)
l2=l[1].split(' \n',1)
l[1]=l2[0]
x[k]=float(l[0])
y[k]=float(l[1])
fichier.close()
return x, y
nomFichier='front1.txt'
N=1509
x, y=readSignal(nomFichier, N)
#y=savgol_filter(y, 51, 3)
#cA, cD=pywt.dwt(y, 'db1')
#cA=pywt.thresholding.hard(cA, 400)
#y=pywt.idwt(cA,cD,'db1', 'sp1')
plt.plot(x,y)
plt.show(block=False)