scipy.signal.detrend中的不準確性/隨機性

Question

爲什麼scipy.signal.detrend對同一數據的結果略有不同？此外，它似乎取決於關鍵字「線性」是否包含以得到不同的結果（默認情況下，消除趨勢是線性反正）

編輯：我知道這個誤差是非常小的，有些不準確預計，由於浮動點算術。奇怪的是，結果是不同爲相同的數據+功能。

from scipy.signal import detrend as scipy_detrend 
from pylab import * 

x = arange(10) 
y = arange(10, dtype='int64') 

subplot(211) 
plot(x, scipy_detrend(y, type="linear"), label='scipy detrend linear') 
plot(x, scipy_detrend(y), label='scipy detrend') 
plot(x, detrend(y, "linear"), label='pylab detrend') 

subplot(212) 
plot(x, scipy_detrend(y, type="linear"), label='scipy detrend linear') 
plot(x, scipy_detrend(y), label='scipy detrend') 
plot(x, detrend(y, "linear"), label='pylab detrend') 

show() 

注：紅線爲pylab.detrend，藍線爲scipy.signal.detrend with linear keyword，綠色只是scipy.signal.detrend。

Answer 1

這是浮點舍入錯誤。在一般情況下，浮點錯誤並不一定是可重複的整個運行相同的CPU，同樣的數據，以及相同的代碼上，因爲它可以通過程序之外的事件受到影響（除非採取特殊措施）：

Consistency of Floating-Point Results using the Intel® Compiler or Why doesn’t my application always give the same answer? - Dr. Martyn J. Corden, David Kreitzer

這似乎是一個常見問題

Answer 2

您的數據是arange（10），您的去結果結果是1e-15的順序，這意味着差異是由於浮點精度造成的。 （您的去結果比您的輸入小15個數量級）