scipy
신호 평활화

Savitzky-Golay 필터 사용

시끄러운 신호가 주어진다 :

import numpy as np
import matplotlib.pyplot as plt
np.random.seed(1)

x = np.linspace(0,2*np.pi,100)
y = np.sin(x) + np.random.random(100) * 0.2

plt.plot(x,y)
plt.show()

scipy.signal.savgol_filter() 메서드를 사용하여 Savitzky-Golay 필터를 사용하여 부드럽게 처리 할 수 ​​있습니다.

import scipy.signal
import numpy as np
import matplotlib.pyplot as plt
np.random.seed(1)

x = np.linspace(0,2*np.pi,100)
y = np.sin(x) + np.random.random(100) * 0.2
yhat = scipy.signal.savgol_filter(y, 51, 3) # window size 51, polynomial order 3

plt.plot(x,y)
plt.plot(x,yhat, color='red')
plt.show()