2013-12-11 41 views
1

測試matplotlib的detrend_linear功能與matplotlib的detrend_linear顯示消解趨勢罪(t)的偏差

import numpy as np 
import matplotlib.pyplot as plt 
from matplotlib.mlab import detrend_linear 

n = 1000 
t = np.arange(n) 
y = np.sin(8.*np.pi/n*t) # exactly 4 periods 

plt.plot(y, label='raw') 
plt.plot(detrend_linear(y), label='detrended') 
plt.legend(loc='best') 
plt.show() 

顯示明顯的偏差: enter image description here

我除外趨勢線是水平的,看看沒有什麼區別。 我的思維在哪裏出錯?

感謝您的任何建議,迪特里希

回答

2

有一個在正弦函數一個趨勢,因爲它不是在該域中對稱。試用4.5個週期,趨勢應該消失。

+0

非常感謝你 - 這是我需要的暗示,die!! – Dietrich