b分析平方逼近.pptx

《b分析平方逼近.pptx》由会员分享，可在线阅读，更多相关《b分析平方逼近.pptx（27页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、6-14 函数的最佳平方逼近 前面对离散数据，我们利用最小二乘法求拟合函数（多项式），本节对一些连续函数，当其表达式较复杂不易于计算和研究时，我们利用最小二乘法，求这些连续函数的近似函数（较简单的函数），称为函数f(x)在a,b上的最佳平方逼函数(x)。第第2 2页页/共共2727页页第1页/共27页6-24.1 基本方法 设f(x)在a,b上连续，i(x)(i=0,1,2,m)在a,b上线性无关，H=Span0,1,m为k(x)的集合，求(x)使：定义6.2连续情况下的内积定义为：(x x)为权函数）为权函数）第第3 3页页/共共2727页页第2页/共27页6-3基本方法（续）要求出满足（6

2、-10）的(x)，与离散情况完全类似，即要求k(x)满足正规方程组（6-5），当k(x)线性无关可求出唯一解是H中关于权函数(x)的唯一的最佳平方逼近多项式。若k(x)=xk(k=0,1,2,m)，此时H为k(x)所有线性组合生成的多项式集合，则(x)称为关于(x)的m次最佳平方逼近多项式或最小二乘逼近多项式。关于权函数(x)一般应给定，若没有特别标明则(x)1。第第4 4页页/共共2727页页第3页/共27页6-4最佳平方逼近多项式举例例例77求求f f(x x)=cos)=cos x x在在0,10,1上的一次最佳平方逼近多项式上的一次最佳平方逼近多项式 问题问题：如何求二：如何求二次、三

3、次最佳平次、三次最佳平方逼近多项式，方逼近多项式，可：可：（1 1）如上，如上，H H=1,=1,x x,x x2 2 即取即取 2 2(x x)=)=x x2 2（2 2）或如后面例，按三项推式构造正交多项式或如后面例，按三项推式构造正交多项式第第5 5页页/共共2727页页第4页/共27页6-5 4.2 利用正交多项式求最佳平方逼近多项式 从上节知道从上节知道 利用正交函数系可以简化最小二乘法的求利用正交函数系可以简化最小二乘法的求解，并提高解的精度，而正交多项式系，由于其计算简便解，并提高解的精度，而正交多项式系，由于其计算简便，是函数逼近的重要工具，后面一致逼近，积分也要用到，是函数逼

4、近的重要工具，后面一致逼近，积分也要用到正交多项式。正交多项式。定义定义6.36.3如果函数系如果函数系 0 0(x x),),1 1(x x),),mm(x x),),满足满足:则称此函数为区间则称此函数为区间 a a,b b 上关于权函数上关于权函数(x x)的正的正交函数系。特别地，若交函数系。特别地，若A Ak k=1(k=0,1,2,)=1(k=0,1,2,)，则称其为，则称其为标准正交函数系，当标准正交函数系，当 k k(x x)为多项式时，称为为多项式时，称为正交多正交多项式。项式。第第6 6页页/共共2727页页第5页/共27页6-6正交多项式举例第第7 7页页/共共2727页

5、页第6页/共27页6-7正交函数系性质 正交函数系具有以下性质：正交函数系具有以下性质：定理定理6.36.3 定理定理6.46.4 设设 k k(x x)(k=0,1,2,)(k=0,1,2,)是最高次项系数不为零的是最高次项系数不为零的k k次多项式，则次多项式，则 k k(x x)是是 a a,b b 上关于权函数上关于权函数(x x)的正交多的正交多项式系的充要条件是对任意至多项式系的充要条件是对任意至多k k1 1次的多项式次的多项式QQk k1 1(x x)，均有：，均有：区间区间 a a,b b 上关于权函数上关于权函数(x x)的正交函数系的正交函数系 0 0,1 1,n n是线

6、性无关的。是线性无关的。第第8 8页页/共共2727页页第7页/共27页6-8定理6.4的证明第第9 9页页/共共2727页页第8页/共27页6-9定理6.5证明类似于定理证明类似于定理6.26.2，略。，略。构造正交多项式的一般方法由以下定理给出：构造正交多项式的一般方法由以下定理给出：定定理理6 6.5 5第第1010页页/共共2727页页第9页/共27页6-10几种常用的正交多项式 下面介绍几种常用的正交多项式下面介绍几种常用的正交多项式:（一）勒让德（一）勒让德（LegendreLegendre）多项式）多项式 LegendreLegendre多项式的一般表示式为多项式的一般表示式为:

7、具体表达具体表达式为式为:第第1111页页/共共2727页页第10页/共27页6-11Legendre多项式性质（1 1）P Pk k(x x)是区间是区间 1,11,1上关于权函数上关于权函数(x x)1 1的正交的正交函数系，且函数系，且 第第1212页页/共共2727页页第11页/共27页6-12Legendre多项式性质（续1）通过变量变换由通过变量变换由LegendreLegendre多项式可以得到在任意区间多项式可以得到在任意区间 a a,b b 上关于权函数上关于权函数(x x)1 1的正交多项式系。的正交多项式系。第第1313页页/共共2727页页第12页/共27页6-13Le

8、gendre多项式性质（续2）（2 2）LegendreLegendre多项式满足递推公式多项式满足递推公式:例例如如：第第1414页页/共共2727页页第13页/共27页6-14（二）第一类切比雪夫（Chebyshev）多项式 第一类第一类ChebyshevChebyshev多项式的一般表示式为：多项式的一般表示式为：令令x x=cos=cos，当，当x x在在 1,11,1上变化时，对应的上变化时，对应的 在在0,0,上变化，上变化，（6-126-12）可改写成：）可改写成：具体表达式为：具体表达式为：由上式容易看出，由上式容易看出，T Tn n(x x)是首是首项系数为项系数为2 2n

9、n-1-1的的n n次次多项式。多项式。第第1515页页/共共2727页页第14页/共27页6-15第一类Chebyshev多项式性质 第一类第一类ChebyshevChebyshev多项式多项式有以下性质：有以下性质：第第1616页页/共共2727页页第15页/共27页6-16第一类Chebyshev多项式性质（续）性质（性质（3 3），（），（4 4）由余弦函数性质即得。由余弦函数性质即得。第第1717页页/共共2727页页第16页/共27页6-17（三）拉盖尔（Laguerre）多项式 LaguerreLaguerre多项式多项式定义为：定义为：其具体表达式为：其具体表达式为：第第181

10、8页页/共共2727页页第17页/共27页6-18Laguerre多项式性质LaguerreLaguerre多项式有以下性质：多项式有以下性质：（1 1）L Ln n(x x)是区间是区间0,+0,+)上关于权函数上关于权函数(x x)=)=e ex x的正交的正交多项多项 式系，且：式系，且：（2 2）LaguerreLaguerre多项式满足递推公式：多项式满足递推公式：由定理由定理6.56.5可以逐步构造在区间可以逐步构造在区间 a a,b b 上关于权函数上关于权函数(x x)的正交多项式系的正交多项式系 n n(x x)进而求出满进而求出满足式（足式（6-86-8）的）的最佳平方最佳

11、平方逼近多项式：逼近多项式：也可直接利用已知的正交多项式作出满足式（也可直接利用已知的正交多项式作出满足式（6-86-8）的）的最佳平方最佳平方逼近多项式：逼近多项式：第第1919页页/共共2727页页第18页/共27页6-19最佳平方逼近多项式举例例例8 8 利用正交多项式求利用正交多项式求y=tg1x在区间在区间0,10,1上的上的最佳平方逼近一最佳平方逼近一次式。次式。第第2020页页/共共2727页页第19页/共27页6-20例8（续）第第2121页页/共共2727页页第20页/共27页6-21Laguerre多项式举例（例9）例9求求f f(x x)=)=SinSin x x在在0,

12、10,1上的二次最佳平方逼近多项上的二次最佳平方逼近多项式式解解:用构造正交多项式用构造正交多项式的方法，的方法，(x x)=1.)=1.第第2222页页/共共2727页页第21页/共27页6-22Laguerre多项式举例（例10）例10求求f f(x x)=)=e ex x 在在 1,11,1上的三次最佳平方逼近多项式上的三次最佳平方逼近多项式 第第2323页页/共共2727页页第22页/共27页6-23LaguerreLaguerre多项式举例（例多项式举例（例1010续）续）若区间不一样，如求若区间不一样，如求f f(x x)=)=e ex x 在在0,10,1上的二次最佳平上的二次最

13、佳平方逼近则需要变换区间，将方逼近则需要变换区间，将 1,11,1上的上的P Pn n(x x)变换为变换为0,10,1区区间上的正交多项式间上的正交多项式:第第2424页页/共共2727页页第23页/共27页6-24第六章结 束第第2525页页/共共2727页页第24页/共27页6-25上机练习题：不同拟合模型的比较已知观测数据如下表所示，按下述方案求最小二乘拟合函数，并求出偏差平方和Q，比较拟合曲线的优劣。方案I拟合函数取为如下形式的三次多项式：方案II用离散正交多项式求三次拟合多项式方案III用离散正交多项式求四次拟合多项式方案IV拟合函数取为如下形式的函数：第第2626页页/共共2727页页第25页/共27页6-26x00.20.61.01.31.61.71.81.9y02.54.05.73.52.01.02.03.5x2.22.32.52.62.93.13.43.84.1y4.07.07.59.910.911.913.513.011.9x4.44.74.84.95.05.15.3y9.06.54.01.50.02.55.0观测数据表第第2727页页/共共2727页页第26页/共27页第六章 函数逼近6-27感谢您的观赏！第27页/共27页