多项式最小二乘法
基于一元二次方程进行推理
方程式是:f(x)=ax^2+bx+c 其中点 N(x_i,y_i) 表示多个点信息,拟合一个二次方程;其中原理为:
∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c -y_i )^2=min〗
对各个系数求导
∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c -y_i )^2=min〗
对a求导:
∑_(i=1)N▒〖(a2 〖x_i〗4+2a〖x_i〗2 (bx_i+c -y_i )+ (bx_i+c -y_i )^2 )=K〗
dK/da= ∑_(i=1)N▒〖(2a〖x_i〗4+2〖x_i〗^2 (bx_i+c –y_i ))=0〗
dK/da= 2〖x_i〗^2 ∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c –y_i )=0〗
对b求导:
∑_(i=1)N▒〖(b2 〖x_i〗^2+2bx_i (a〖x_i〗^2+c -y_i )+ (a〖x_i〗^2+c -y_i )^2 )=K〗
dK/db= ∑_(i=1)N▒〖(2b〖x_i〗2+2x_i (a〖x_i〗^2+c -y_i ))=0〗
dK/db= 2x_i ∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c –y_i )=0〗
对c求导:
∑_(i=1)N▒〖(c2+2c(a〖x_i〗^2+bx_i -y_i )+ (a〖x_i〗^2+bx_i -y_i )^2 )=K〗
dK/dc= ∑_(i=1)N▒〖(2c+2(a〖x_i〗2+bx_i -y_i ))=0〗
dK/dc= 2∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c –y_i )=0〗
对求导公式进行整理及变换
2〖x_i〗^2 ∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c –y_i )=0〗
2x_i ∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c –y_i )=0〗
2∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c –y_i )=0〗
变换如下:
2〖x_i〗^2 ∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c)=〗 2〖x_i〗^2 ∑_(i=1)^N▒y_i
2x_i ∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c )=2x_i ∑_(i=1)^N▒y_i 〗
2∑_(i=1)N▒〖(a〖x_i〗2+bx_i+c )=2∑_(i=1)^N▒y_i 〗