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

正交多项式

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

勒让德多项式(Legendre)

在这里插入图片描述
在这里插入图片描述

勒让德多项式的性质

在这里插入图片描述

在这里插入图片描述

其中最为重要的就是性质一到三和五,第四条基本上不咋用。然后最重要的表格来了,其实也可以用性质三来推出来:
在这里插入图片描述
整个题目来理解一下其中的妙处
在这里插入图片描述

性质

在这里插入图片描述

切比雪夫多项式(Shebyshev)

和上面学勒让德多项式一样来学习切比雪夫多项式;
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
同样的,咱搞个实战试一试
在这里插入图片描述

最小二乘法

定义;
在这里插入图片描述
求函数 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),再改写为矩阵形式:
在这里插入图片描述
也就得到了对应的线性拟合模型。

多项式拟合

在离散数据 ( 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
称为多项式拟合。
加权的最小二乘法一般是在通常拟合函数前加一个加权函数乘以拟合部分,具体的可以自己百度一下,嘻嘻。

最小二乘法实战

在这里插入图片描述
函数逼近我们将告一段落,后面我们将开启数值积分的版块,船新版本哦,等你来学!!!

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

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

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

  • 6
    点赞
  • 34
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值