主題
Search

最小二乘擬合--垂直偏移


LeastSquaresOffsets

在實踐中,通常最小化的是來自直線(多項式、曲面、超平面等)的垂直偏移,而不是垂直偏移。 這為自變數 X 提供了一個擬合函式,用於估計給定 xy (通常是實驗者想要的),允許簡單地納入資料點沿 x 軸和 y 軸的不確定性,並且還為擬合引數提供了比基於垂直偏移的擬合獲得的更簡單的解析形式。

對於一組 n 個點,使用點 (x_i,y_i) 的未平方垂直距離 d_i 的最佳擬合線的殘差由下式給出

 R__|_=sum_(i=1)^nd_i.
(1)

由於從直線 y=a+bx 到點 i 的垂直距離由下式給出

 d_i=(|y_i-(a+bx_i)|)/(sqrt(1+b^2)),
(2)

要最小化的函式是

 R__|_=sum_(i=1)^n(|y_i-(a+bx_i)|)/(sqrt(1+b^2)).
(3)

不幸的是,由於絕對值函式不具有連續導數,因此最小化 R__|_ 不利於解析解。 然而,如果垂直距離的平方

 R__|_^2=sum_(i=1)^n([y_i-(a+bx_i)]^2)/(1+b^2)
(4)

被最小化,則問題可以以閉合形式解決。 當滿足以下條件時,R__|_^2 達到最小值

 (partialR__|_^2)/(partiala)=2/(1+b^2)sum_(i=1)^n[y_i-(a+bx_i)](-1)=0
(5)

 (partialR__|_^2)/(partialb)=2/(1+b^2)sum_(i=1)^n[y_i-(a+bx_i)](-x_i)+sum_(i=1)^n([y_i-(a+bx_i)]^2(-1)(2b))/((1+b^2)^2)=0.
(6)

前者給出

a=(sum_(i=1)^(n)y_i-bsum_(i=1)^(n)x_i)/n
(7)
=y^_-bx^_,
(8)

後者給出

 (1+b^2)sum_(i=1)^n[y_i-(a+bx_i)]x_i+bsum_(i=1)^n[y_i-(a+bx_i)]^2=0.
(9)

但是

[y-(a+bx)]^2=y^2-2(a+bx)y+(a+bx)^2
(10)
=y^2-2ay-2bxy+a^2+2abx+b^2x^2,
(11)

所以 (10) 變為

(1+b^2)(sum_(i=1)^(n)x_iy_i-asum_(i=1)^(n)x_i-bsum_(i=1)^(n)x_i^2)+b(sum_(i=1)^(n)y_i^2-2asum_(i=1)^(n)y_i-2bsum_(i=1)^(n)x_iy_i+a^2sum_(i=1)^(n)1+2absum_(i=1)^(n)x_i+b^2sum_(i=1)^(n)x_i^2)=0
(12)
[(1+b^2)(-b)+b(b^2)]sum_(i=1)^(n)x_i^2+[(1+b^2)-2b^2]sum_(i=1)^(n)x_iy_i+bsum_(i=1)^(n)y_i^2+[-a(1+b^2)+2ab^2]sum_(i=1)^(n)x_i-2absum_(i=1)^(n)y_i+ba^2sum_(i=1)^(n)1=0
(13)
-bsum_(i=1)^(n)x_i^2+(1-b^2)sum_(i=1)^(n)x_iy_i+bsum_(i=1)^(n)y_i^2+a(b^2-1)sum_(i=1)^(n)x_i-2absum_(i=1)^(n)y_i+ba^2n=0.
(14)

將 (◇) 代入 (14) 然後給出

 -bsum_(i=1)^nx_i^2+(1-b^2)sum_(i=1)^nx_iy_i+bsum_(i=1)^ny_i^2+1/n(b^2-1)(sum_(i=1)^ny_i-bsum_(i=1)^nx_i)sum_(i=1)^nx_i-2/n(sum_(i=1)^ny_i-bsum_(i=1)^nx_i)bsum_(i=1)^ny_i+b/n(sum_(i=1)^ny_i-bsum_(i=1)^nx_i)^2 
=0
(15)

經過相當多的代數運算後,結果是

 b^2+(sum_(i=1)^(n)y_i^2-sum_(i=1)^(n)x_i^2+1/n[(sum_(i=1)^(n)x_i)^2-(sum_(i=1)^(n)y_i)^2])/(1/nsum_(i=1)^(n)x_isum_(i=1)^(n)y_i-sum_(i=1)^(n)x_iy_i)b-1=0.
(16)

所以定義

B=1/2([sum_(i=1)^ny_i^2-1/n(sum_(i=1)^ny_i)^2]-[sum_(i=1)^nx_i^2-1/n(sum_(i=1)^nx_i)^2])/(1/nsum_(i=1)^nx_isum_(i=1)^ny_i-sum_(i=1)^nx_iy_i)
(17)
=1/2((sum_(i=1)^ny_i^2-ny^_^2)-(sum_(i=1)^nx_i^2-nx^_^2))/(nx^_y^_-sum_(i=1)^nx_iy_i),
(18)

二次公式給出

 b=-B+/-sqrt(B^2+1),
(19)

其中 a 使用 (◇) 找到。 請注意公式中最佳擬合引數的相當笨拙的形式。 此外,對於二階或更高階多項式最小化 R__|_^2 會導致具有更高階的多項式方程,因此這種公式無法擴充套件。


另請參閱

最小二乘擬合

使用 探索

參考文獻

Sardelis, D. 和 Valahas, T. "最小二乘擬合-垂直偏移。" http://library.wolfram.com/infocenter/MathSource/5292/

請引用為

Weisstein, Eric W. "最小二乘擬合--垂直偏移。" 來自 Web 資源。 https://mathworld.tw/LeastSquaresFittingPerpendicularOffsets.html

主題分類