主题
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 会导致具有更高阶的多项式方程,因此这种公式无法扩展。


另请参阅

最小二乘拟合

使用 Wolfram|Alpha 探索

参考文献

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

请引用为

Weisstein, Eric W. "最小二乘拟合--垂直偏移。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/LeastSquaresFittingPerpendicularOffsets.html

主题分类