数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法

正交多项式

上一节中我们学会了最佳一致逼近和平方逼近的方式求解函数逼近多项式,但是发现在求解法方程时,多元方程求解难度较大,如果我们能够将对应的法方法的矩阵直接变为只有对角线元素不为零的矩阵,将极大降低计算的运算量,因此我们引入了正交函数族的概念,从而得出正交多项式。下面简单介绍两种正交多项式勒让德多项式切比雪夫多项式,其他的正交多项式如:
L a g u e r r e 多 项 式 Laguerre多项式 Laguerre H e r m i t e 多 项 式 Hermite多项式 Hermite有兴趣的可以自行了解。

勒让德多项式(Legendre)

数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第1张图片
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第2张图片

勒让德多项式的性质

数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第3张图片

数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第4张图片

其中最为重要的就是性质一到三和五,第四条基本上不咋用。然后最重要的表格来了,其实也可以用性质三来推出来:
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第5张图片
整个题目来理解一下其中的妙处
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第6张图片

性质

数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第7张图片

切比雪夫多项式(Shebyshev)

和上面学勒让德多项式一样来学习切比雪夫多项式;
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第8张图片
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第9张图片
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第10张图片
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第11张图片
同样的,咱搞个实战试一试
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第12张图片

最小二乘法

定义;
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第13张图片
求函数 S ∗ ( x ) = ∑ j = 0 n α j ∗ φ j ( x ) S^*(x)=\sum\limits_{j=0}^n\alpha_j^*\varphi_j(x) S(x)=j=0nαjφj(x)的问题可归结为求它的系数 α 0 ∗ , . . . , α n ∗ \alpha_0^*,...,\alpha_n^* α0,...,αn使得多元函数
F ( α 0 , α 1 , . . . , α n ) ≜ ∑ i = 1 m [ f ( x i ) − ∑ j = 0 n α j φ j ( x j ) ] 2 \qquad\qquad\qquad\qquad\qquad\Large F(\alpha_0,\alpha_1,...,\alpha_n)\triangleq\sum\limits_{i=1}^m[f(x_i)-\sum\limits_{j=0}^n\alpha_j\varphi_j(x_j)]^2 F(α0,α1,...,αn)i=1m[f(xi)j=0nαjφj(xj)]2取得最小值。由于 F ( α 0 , α 1 , . . . , α n ) F(\alpha_0,\alpha_1,...,\alpha_n) F(α0,α1,...,αn)是关于 α 0 , α 1 , . . . , α n \large \alpha_0,\alpha_1,...,\alpha_n α0,α1,...,αn的二次函数,所以可以使用多元函数取得极值的必要条件得:
∂ F ∂ α k = 0 , ( k = 0 , 1 , 2 , . . . , n ) \qquad\qquad\qquad\qquad\qquad\frac{\large\partial F}{\large\partial \alpha_k}=0 ,(k=0,1,2,...,n) αkF=0,(k=0,1,2,...,n)则通过
F ( α 0 , α 1 , . . . , α n ) ≜ ∑ i = 1 m [ f ( x i ) − ∑ j = 0 n α j φ j ( x j ) ] 2 \qquad\qquad\qquad\qquad\large F(\alpha_0,\alpha_1,...,\alpha_n)\triangleq\sum\limits_{i=1}^m[f(x_i)-\sum\limits_{j=0}^n\alpha_j\varphi_j(x_j)]^2 F(α0,α1,...,αn)i=1m[f(xi)j=0nαjφj(xj)]2有:
∂ F ∂ α k = − 2 ∑ i = 1 n { [ f ( x i ) − ∑ j = 0 n α j φ j ( x i ) ] φ k ( x i ) ] } = 0 , ( k = 0 ∼ n ) \qquad\qquad\qquad\qquad\frac{\large\partial F}{\large\partial \alpha_k}=-2\sum\limits_{i=1}^n \lbrace [f(x_i)-\sum\limits_{j=0}^n\alpha_j\varphi_j(x_i)]\varphi_k(x_i)]\rbrace=0,(k=0 \sim n) αkF=2i=1n{[f(xi)j=0nαjφj(xi)]φk(xi)]}=0,(k=0n)
得到方程组:
∑ i = 1 m ( ∑ j = 0 n φ k ( x i ) φ j ( x i ) ) α j = ∑ i = 1 m f ( x i ) φ k ( x k ) \qquad\qquad\qquad\sum\limits_{i=1}^m(\sum\limits_{j=0}^n\varphi_k(x_i)\varphi_j(x_i))\alpha_j=\sum\limits_{i=1}^mf(x_i)\varphi_k(x_k) i=1m(j=0nφk(xi)φj(xi))αj=i=1mf(xi)φk(xk)在离散情况下, f ( x ) , g ( x ) f(x),g(x) f(x),g(x)关于点 x 1 , x 2 , x 3 , . . . , x m x_1,x_2,x_3,...,x_m x1,x2,x3,...,xm的内积为 ( f , g ) = ∑ i = 1 m f ( x i ) g ( x i ) (f,g)=\sum\limits_{i=1}^mf(x_i)g(x_i) (f,g)=i=1mf(xi)g(xi)也就是说可以将方程组改写为如下形式: ∑ j = 0 n ( φ k , φ j ) φ j = ( f , φ k ) , ( k = 0 ∼ n ) \sum\limits_{j=0}^n(\varphi_k,\varphi_j)\varphi_j=(f,\varphi_k),(k=0\sim n) j=0n(φk,φj)φj=(f,φk),(k=0n),再改写为矩阵形式:
数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第14张图片
也就得到了对应的线性拟合模型。

多项式拟合

在离散数据 ( x i , y i ) i = 1 − m (x_i,y_i)_{i=1-m} (xi,yi)i=1m的最小二乘拟合中,最简单,最常用的数学模型就是多项式,取 φ 0 ( x ) = 1 , φ 1 ( x ) = x , φ 2 ( x ) = x 2 , . . . , φ n ( x ) = x n \varphi_0(x)=1,\varphi_1(x)=x,\varphi_2(x)=x^2,...,\varphi_n(x)=x^n φ0(x)=1,φ1(x)=x,φ2(x)=x2,...,φn(x)=xn得到的拟合函数 f ( x ) = c 0 + c 1 x + c 2 x 2 + . . . + c n x n f(x)=c_0+c_1x+c_2x^2+...+c_nx^n f(x)=c0+c1x+c2x2+...+cnxn
称为多项式拟合。
加权的最小二乘法一般是在通常拟合函数前加一个加权函数乘以拟合部分,具体的可以自己百度一下,嘻嘻。

最小二乘法实战

数值计算(六)——函数逼近 (2)正交多项式多项式和最小二乘法_第15张图片
函数逼近我们将告一段落,后面我们将开启数值积分的版块,船新版本哦,等你来学!!!

欢乐的时光总是短暂的,让我们下一次再见!!!

good good study,day day up! (study hard, improve every day)

预知后事,请听下回分解!!!!

你可能感兴趣的:(数值计算,最小二乘法,数值计算,正交多项式,勒让德多项式,切比雪夫正交多项式)