《正交多项式ppt课件.ppt》由会员分享,可在线阅读,更多相关《正交多项式ppt课件.ppt(26页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、v主要内容主要内容:v正交多项式的构造;v常用的多项式;v一致逼近的基本概念;v最佳一致逼近多项式;v均方逼近的基本概念;v最佳均方逼近多项式;v最小二乘曲线拟合的基本概念;v用正交多项式作最小二乘曲线拟合。第六章 函数逼近与拟合(FunctionApproximationandInterpolation)v函数逼近问题:函数逼近问题:在实际应用中常需为解析式子比较复杂的函数寻找一个多项式来近似代替它,并要求其误差在某种度量意义下最小。v曲线拟合问题:曲线拟合问题:在实际应用中,往往并不需要多项式通过给定的数据点,而只要求用多项式近似代替列表函数时,其误差在某种度量意义下最小。第六章 函数逼近
2、与拟合(FunctionApproximationandInterpolation)第六章 函数逼近与拟合(FunctionApproximationandInterpolation)6. 1 6. 1 正交多项式正交多项式() , ,0,1,0,( )d( ) 0,)0()()(jmbanjmnxxaxbxmxnxx对于定义在区间上的一个函数系,如果其中任何两个函数在此区间上的积分为零,而他们之中每个函数自乘的积分不等于零,即,则称此函数系为在此区间上关于权正交函数的函数系。( )d)1( )bnnaxxxx规范的正交函数当时称之为;当此函数系中的每一个函数均为多项式时称之为系正交多项式(系
3、)。 ,0,coscosd,02 ,0cossind00,sinsind,0mnmnmnmnmxxxnxnnmnmnxxx xmxxnxx xx如在区间上的正交函数系1,cos ,sin ,cos2 , ,con,sin,中,6. 16. 1正交多项式正交多项式(Orthogonal(Orthogonal Multinomial)Multinomial)6. 16. 1正交多项式正交多项式(Orthogonal Multinomial)01011( )( )( )( ) , ( )0,11()( )( ),1,2jjjjjjQ xa bQ xQ xQxQ xQxjxxxj:为了构造在给定区间上
4、关于权函数1的正交多项式系,,可用递推方法正交多项式的构,造221( )( )d0,1,1,2d ,0,1bajjjjjbjajjdddQxxxjjxxjdQ其中,正交多项式的构造1110002111,( )()()( )( )Q xQjxxxdxQ xd如 01000001110102222212ddd() d( )( )( )(ddd) d)bbaabbaabbaabbaaxxx xxxx xxxQxQxQxddddddxQxxxx其中, ,, Example 6.11,14在区间上构造正交多项式。011102211()(),1,2d0,1,1,2d ,0(,1)( )( )( )( )(
5、 )( )jjjjjbjajjjjbjajjQ xQ xQxQ xQxxxjddxxjjxjddQxQx解:由递推方法,其中,Example 6.1111410042222000011 dd35d1 d43/485()( )8( )bbaaxxxxQxQxQxxxddxx解:,11141141111012222111102259d() d82565() dd59/25638,9/25683/464( )( )( )553511()()()886443( )( )( )2bbjaadQxQxQ xQ xQxxxxxxxxxxxxxxdxjdQd,为次多项式cos , 0cos()cos( arc
6、)( )cosnnxxnT xn切比:设则称为(第一类雪夫多项式切比雪夫) 阶多项式。12210,1, 1,1( )( )( )( )1101d010kmnT xTx Txmnkxxxxmnmn2:,在区间上关于权函数是正交多性质1切比雪夫多项式项式系,即,1. 正交多项式(OrthogonalMultinomial)0111( )( )(0,1,12,1,2)( )( )( )kjjjT xT xT xTxT xTxkxxj:,具有下列三项递推关系,性质2切比雪夫多项式11( )(2(2)nnnnnnT xT xnnxx: 阶是一个次代数多项式, 且其最性质3切比雪夫多高次幂项的系数为。项式
7、首一(为多项式)切比雪夫多项式 1,1 1,121,1,2,221cos,( )(1),2,2nkknTkknnkknnnxxxnT: 阶在区间上满足1, 且在区间上有性质4切比雪个零点即夫多项式切比雪夫多项式 1,1,0,1,cos,0,( )11,kknkknnkTxknxnnn: 阶在区性质5切比雪夫多项式间上有个峰值点即11111111 1,12max0max0 ,1max( )( )( )( )( )(2)nxxnxnnnnnnP xT xxxP xPnx :在区间上所有的次中,对零的偏差最小,即定理6-1首一且多项式切比雪夫多项式0110( )( 1,11)( )(2)( )nnf
8、 xf xT xTaxxaTa切比雪夫级:定义在区间上的函数数可以展开成01 1,1 1,1( )( )(2)1kkkf xf xa Tax当函数在上有连续导数时,在上是绝对、一切比雪夫级数致收敛的,且1. 正交多项式(OrthogonalMultinomial)cos121021d12d0,1,( )( )(cos )coskkxaxxf x T xfkk 其中Example 6.2切比雪夫级数(1)用计算函数的近似值。11210( )( )(cos 1,141221cosd0,1,(4)1) kkkkkaf xf xT xT xxxxxkxaa解:设函数切比雪夫级中,数,则其2210211
9、2122( )( )(01,2,1221cosd0,1,41)( )kkknkkknxT xf xTxxkxkaaTxaa当为奇数时,为奇数,,其中,Example 6.220121021cosd( )41kxxaT xxa解: 如果要求精确到小数点后第九位, 1211210-133( )d0322121cosd41232cos11.703283757342f xxfx xaffxx利用三点高斯 切比雪夫求积公式(第七章中)可求得Example 6.212222121(21)cosd41233221 cos10.1483581213242 xaxxx解:1422142421(881)cosd4
10、123332881 cos1322420.001921449xxx xxa Example 6.2164226461221(3248181)cosd4123333323248181 cos1222420.000009965xxxx xxa 解:80.000000028a 类似求出最后一个不为零的系数Example 6.20111( )( )(0,1,)( )( )(12,1,2)kjjjkxT xT xT xTxT xTxxj解:,具有下列三项递推关系,切比雪夫多项式234210321432543654253642221243288( )( )( )( )( )( )( )( )( )( )(
11、121)( )( )( )6205232481)8(1xxxxxxxxxxT xT xT xT xT xT xT xT xT xT xT xT xxxxT xT xxTxxxExample 6.202122426641cos420.8516418780.1483581210.0019214490.0000099650.9999994720.3084232530( )( )( )(.0158499130.0003188)(70)kkkf xTxT xTxxxaT xaxx解: , 22 , , ( )( )21)12(f xxtxtxttffFa bbabaa bbaba:对于定义在一般区间上的
12、函数,可以先作如下变换即在区间上变注化时在区间上变化,可利用切比雪夫级数 ( )kkT xx:与之间存在线性变换关系,所以对于一个函数的逼近多项式,可利用来找一个次数较低的线性变换关系切比新的近似多项式,雪夫多项式切且满足相同的比雪夫多项式精度要求。T0 = 11 = T0T1 = xx = T1T2 = 2x2-1x2 = (T0+T2)/2T3 = 4x3-3xx3 = (3T1+T3)/4T4 = 8x4-8x2+1x4 = (3T0+4T2+T4)/8 ( )kkT xx与之间存在的线性变换关系1. 正交多项式(OrthogonalMultinomial)Example 6.3(2)用
13、来降低逼近多切比雪夫多项式项式的次数。23456 1,101111126241200.0038( )720( )( )xxf xf xE xexxexxxxxex 解:设函数,如果对其在 点处作幂级数有限项展开作近似代替,即此五次多项式替代的误差为01021302413511()(3)42411(34)(105)1)9219 0(2TTTTTTTTTf xTTT现利用切比雪夫多项式和幂函数的线性变换关系Example 6.30123450123812171317116419248384192192081217131764192483( )(4)8f xf xTTTTTTTTTT进一步近似解:
14、1,11110.005819219200.0( )nT x在区间上此时由省去的最后两项增加的误差即用三次多项式替代的总误差为0.0096,仍不超过1。2233(21)81217131764192483(841(38239320868)3 )3 4( )48xxxxfxxxx线性变换关系2, 1d(1) 2!( )1d1nnnnnnP xxnx勒让德多项式的构造勒让德(L:在区间上定义egendre函数称为阶)多项式110,1, 10d22( )( )(11)1( )mknxmkP xPnmPxnnxx:,在区间上关于权函数是正交函数系,即,性质1勒让德多项式( )nnxnP性质2勒让:当为偶/
15、奇数时, 阶为偶德多项式/奇数。1. 正交多项式(OrthogonalMultinomial)1101( )( )( )( )( )(0,1,1(1)(21),1)( ),2, 1 1nnnnkP xnnnP xP xPxP xPxnkxxxnnnP:,具有下列递推关系,,显然阶是一个次代数多项式,且性质3勒让德多在区间上项式勒有个不同的让德多项式实零点。22(2 ),2 ( )d(1) (2 )( )d,(1 1)nnnnnnnnnnnnxP xxP x性质4勒让德多项式勒让德首一多项式!: 阶最高次项的系数为!阶,!且在区间上与零的平方误差最小。勒让德多项式(LegendreMultinomial)第六章习题第六章习题pp192193,习题6:1. ; 4. ; 5.;8.; 9.;11. ; 15. 。 Please wait for a while!I will be back soon!