上一节中我们学会了最佳一致逼近和平方逼近的方式求解函数逼近多项式,但是发现在求解法方程时,多元方程求解难度较大,如果我们能够将对应的法方法的矩阵直接变为只有对角线元素不为零的矩阵,将极大降低计算的运算量,因此我们引入了正交函数族的概念,从而得出正交多项式。下面简单介绍两种正交多项式勒让德多项式和切比雪夫多项式,其他的正交多项式如:
L a g u e r r e 多 项 式 Laguerre多项式 Laguerre多项式, H e r m i t e 多 项 式 Hermite多项式 Hermite多项式有兴趣的可以自行了解。
其中最为重要的就是性质一到三和五,第四条基本上不咋用。然后最重要的表格来了,其实也可以用性质三来推出来:
整个题目来理解一下其中的妙处
和上面学勒让德多项式一样来学习切比雪夫多项式;
同样的,咱搞个实战试一试
定义;
求函数 S ∗ ( x ) = ∑ j = 0 n α j ∗ φ j ( x ) S^*(x)=\sum\limits_{j=0}^n\alpha_j^*\varphi_j(x) S∗(x)=j=0∑nα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=1∑m[f(xi)−j=0∑nα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) ∂αk∂F=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=1∑m[f(xi)−j=0∑nα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) ∂αk∂F=−2i=1∑n{[f(xi)−j=0∑nαjφj(xi)]φk(xi)]}=0,(k=0∼n)
得到方程组:
∑ 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=1∑m(j=0∑nφk(xi)φj(xi))αj=i=1∑mf(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=1∑mf(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=0∑n(φk,φj)φj=(f,φk),(k=0∼n),再改写为矩阵形式:
也就得到了对应的线性拟合模型。
在离散数据 ( x i , y i ) i = 1 − m (x_i,y_i)_{i=1-m} (xi,yi)i=1−m的最小二乘拟合中,最简单,最常用的数学模型就是多项式,取 φ 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
称为多项式拟合。
加权的最小二乘法一般是在通常拟合函数前加一个加权函数乘以拟合部分,具体的可以自己百度一下,嘻嘻。
函数逼近我们将告一段落,后面我们将开启数值积分的版块,船新版本哦,等你来学!!!