《插值方法课件.ppt》由会员分享，可在线阅读，更多相关《插值方法课件.ppt（81页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第1页，此课件共81页哦第2页，此课件共81页哦第3页，此课件共81页哦第4页，此课件共81页哦第5页，此课件共81页哦第6页，此课件共81页哦第7页，此课件共81页哦 拉格朗日插值公式的基本思想是:把Pn(x)的构造问题转化为n+1个插值基函数Li(x)（i=0,1,n)的构造。插值函数Pn(x)是这n1插值基函数的线性组合，其组合系数就是对应点上的函数值。knknijjjijnkKknYxxxxxLYxP)()()(000 第8页，此课件共81页哦拉格朗日插值的算法：Input X,(Xi,Yi),I=1,1,n0S,1kFOR K=1 TO n t=1 FOR j=1 TO n jk .
2、F.T.t*(X-Xj)/(Xk-Xj)t S+t*Yk SOutput S第9页，此课件共81页哦绘制拉格朗日插值基函数图形绘制拉格朗日插值基函数图形输入插值节点输入插值节点xi(0,n)输入插值点的个数输入插值点的个数 mh=(xm-x0)/(m-1)计算计算m个插值点上的插值基函数个插值点上的插值基函数l(i,k)和和 ll(k)=l(i,k)for k=0 to m xx(k)=x0+k*h ll(k)=0 for i=0 to n l(i,k)=1 for j=0 to n if ji l(i,k)=l(i,k)*(xx(k)-x(j)/(x(i)-x(j)end if next j
3、 ll(k)=ll(k)+l(i,k)next Inext k 确定坐标确定坐标 绘制图形绘制图形第10页，此课件共81页哦牛顿插值公式的基本思想是牛顿插值公式的基本思想是:把把P Pn n(x)(x)设设计为递推形式。计为递推形式。P Pn n(x)(x)P Pn-1n-1(x)+a(x)+an n(x-x(x-x0 0)(x-x)(x-xn-1n-1)=a =a0 0+a+a1 1(x-x(x-x0 0)+a)+a2 2(x-x(x-x0 0)(x-x)(x-x1 1)+)+a +an n(x-x(x-x0 0)(x-x)(x-xn-1n-1)满足满足P Pn n(x(xi i)=f(x)
4、=f(xi i)=Y)=Yi i i=0,1,2,ni=0,1,2,n 称为牛顿插值公式。称为牛顿插值公式。第11页，此课件共81页哦 为了简化上式，引进记号：为了简化上式，引进记号：0 0(x)=1 i i(x)=(x-xi-1i-1)i-1i-1(x)=(x-xx-x0 0)(x-x)(x-x1 1)(x-x)(x-xi-1i-1)i=1,2 i=1,2，n n P Pn n(x)(x)a a0 0 0 0(x)+a+a1 1 1 1(x)+a+an n n n(x)0 0(x),1 1(x),n n(x)称为牛顿插值以称为牛顿插值以x0 0,x1 1,xn n为插值节点的为插值节点的基函
5、数基函数。如果增加一个新节点，只需增加一个新项如果增加一个新节点，只需增加一个新项 a an n n n(x)第12页，此课件共81页哦其系数为各阶差商：),.,(.),(101001011nnxxxfaxxfxxyya第13页，此课件共81页哦不保留各阶差商的中间结果：For j=1 to n(计算各阶差商）For i=0 to n-j(计算J阶差商 中的每个差商）(Y(i+1)-Y(i)/(x(j+i)-x(i)Y(i)保留各阶差商的中间结果：For j=1 to n(计算各阶差商）For i=j to n(计算J阶差商 中的每个差商）y(j,i)=(y(j-1,i)-y(j-1,i-1)
6、/(x(i)-x(i-j)差商表差商表X f(x)一阶差商 二阶差商 三阶差商X0 f(x0)X1 f(x1)f(x0,x1)X2 f(x2)f(x1,x2)f(x0,x1,x2)X3 f(x3)f(x2,x3)f(x1,x2,x3)f(x0,x1,x2,x3)第14页，此课件共81页哦牛顿插值的算法(不保留各阶差商的中间结果）：Input X,(Xi,Yi),I=0,1,2,ny0y,1tFOR j=1 TO n t=t*(x-xj-1)FOR i=0 TO n-j yi=(yi+1-yi)/(xj+i-xi)y+t*Y0 yOutput y第15页，此课件共81页哦input 请输入的值：
7、to xxuse sj n=reccount()dime x(n),y(n)k=1do while.not.eof()x(k)=ax y(k)=ay k=k+1 skipenddot=1s=y(1)j=1do while j=n t=t*(xx-x(j)i=1 do while i=n-j y(i)=(y(i+1)-y(i)/(x(j+i)-x(i)i=i+1 enddo s=s+t*y(1)j=j+1enddo?F(+str(xx,6,4)+)=+str(s,10,6)return*牛顿插值的程序：插值点数据存放在数据库文件SJ中第16页，此课件共81页哦牛顿插值的算法(保留并输出各阶差商）
8、：Input X,(Xi,Y0,i),I=0,1,2,nY0,iy,1tFOR j=1 TO n t=t*(x-xj-1)FOR i=j TO n y(j,i)=(y(j-1,i)-y(j-1,i-1)/(x(i)-x(i-j)y+t*Y(j,j)yOutput y(j,i),x,y第17页，此课件共81页哦 埃特金插值公式的基本思想是：平面上的两个点可以连成一埃特金插值公式的基本思想是：平面上的两个点可以连成一条线，对应一个线性函数；把线性函数看作形式点，经线性条线，对应一个线性函数；把线性函数看作形式点，经线性组合，可以构成二次函数；把二次函数看作形式点，经线性组合，可以构成二次函数；把二
9、次函数看作形式点，经线性组合，可以构成三次函数。因次埃特金插值公式也叫组合，可以构成三次函数。因次埃特金插值公式也叫逐步线逐步线性插值性插值。过两点（过两点（X0,Y0),(X1,Y1）的一次线性插值函数为）的一次线性插值函数为:)()()(101001011,0 xfxxxxxfxxxxxP第18页，此课件共81页哦过两点（X0,Y0),(X2,Y2）的一次线性插值函数为:)()()(202002022,0 xfxxxxxfxxxxxP过（X1,P0,1(x),(X2,P0,2(x)两形式上的点的线性插值函数为:)()()(2,01211,02122,1,0 xPxxxxxPxxxxxP第1
10、9页，此课件共81页哦 对于一般情况：过两个形式点（对于一般情况：过两个形式点（Xn-1,P0,1,2,n-1（x）和）和（Xn,P0,1,n-2,n(x)）的线性插值函数为）的线性插值函数为:)()()(,2,1,0111,1,01,1,0 xPxxxxxPxxxxxnnnnnnnnnnP第20页，此课件共81页哦 Y(i)=(x-x(i)/(x(k-1)-x(i)*y(k-1)+(x-x(k-1)/(x(i)-x(k-1)*y(i)或 y(i)=y(k-1)+(y(i)-y(k-1)/(x(i)-x(k-1)*(x-x(k-1)埃特金插值表埃特金插值表X f(x)X0 f(x0)X1 f(
11、x1)P0,1(x)X2 f(x2)P0,2(x)P0,1,2(x)X3 f(x3)P0,3(x)P0,1,3(x)P0,1,2,3(x)第21页，此课件共81页哦埃特金插值的算法：Input X,(Xi,Yi),I=0,1,2,nFOR k=1 TO n FOR i=k TO n y(i)=y(k-1)+(y(i)-y(k-1)/(x(i)-x(k-1)*(x-x(k-1)Output y(n)第22页，此课件共81页哦 定理1 有惟一的n次多项式Pn(x)，满足条件：Pn(xi)=yi（i=0,1,n)第23页，此课件共81页哦 设设Pn(x)是在点是在点X0,X1,，Xn处关于处关于f(
12、x)f(x)的插值多项式，当的插值多项式，当X X XiXi时，时，f(x)f(x)与与 Pn(x)的偏差的偏差Rn(x)f(x)f(x)Pn(x)，称为插值余项。，称为插值余项。定理定理 2 若若f(x)f(x)在在包含插值节点包含插值节点X X0,X,X1,，X Xn的区间的区间a,ba,b上上n+1n+1次可微分，则对任意次可微分，则对任意x,xx,x a,b,a,b,有与有与x x有关的有关的(a(a b)b)存在，存在，使得：使得：R Rn(x)=(x)=f(x)f(x)Pn(x)其中其中(x)=(x-x0)(x-x1)(x-xn)()!1()()1(nfnx第24页，此课件共81页
13、哦随着插值次数的增高,插值余项的变化因被插值函数的不同而异。例如对于被插值函数 在两个端点x=5和x=-5的附近插值函数P10(x)与原函数 偏离很远。而在-2，2范围内，P10(x)能较好地逼近f(x)。这说明插值多项式这说明插值多项式P P1010(x)(x)对对f(x)f(x)的近似效果并不是随着插值次的近似效果并不是随着插值次数的增高而增高的数的增高而增高的。随着插值次数的增高，在段点处插值函数与原函数偏离原大的现象称为高次插值的龙格现象高次插值的龙格现象。高次插值的高次插值的SinSin现象现象说明的是插值次数越高插值多项式Pn(x)对函数Sin(x)的逼近效果越好。当插值函数高到一
14、定的次数后，插值函数的曲线和原函数的曲线就完全重叠。211)(xxf211)(xxf第25页，此课件共81页哦输入插值节点输入插值节点xi(0,n)计算插值节点上的原函数计算插值节点上的原函数yi(0,n)输入插值点的个数输入插值点的个数mh=(xn-x0)/(m-1)计算插值点上的插值基函数计算插值点上的插值基函数l、插值函数、插值函数ll(k)和原函数和原函数yy(k)for k=0 to m xx(k)=x0+i*h ll(k)=0 for i=0 to n l=1 for j=0 to n if ji l=l*(xx(k)-x(j)/(x(i)-x(j)end if next j ll
15、(k)=ll(k)+l*y(i)next I yy(k)=1/(1+xx(k)*xx(k)next k确定坐标确定坐标绘制图形绘制图形第26页，此课件共81页哦 为了避免龙格现象的发生，我们很自然地会想到把为了避免龙格现象的发生，我们很自然地会想到把区间区间5，5等分为等分为10个小区间，在每一个小个小区间，在每一个小区间内应用低次插之。但因为每个小区间只有两区间内应用低次插之。但因为每个小区间只有两个端点，按照我们已知的方法，得到的将是一个个端点，按照我们已知的方法，得到的将是一个分段线性插值函数，这显然不能满足对任意一个分段线性插值函数，这显然不能满足对任意一个曲线函数的逼近要求。要想在每
16、一个区间得到一曲线函数的逼近要求。要想在每一个区间得到一个曲线函数，不仅要利用插值节点处的函数值，个曲线函数，不仅要利用插值节点处的函数值，还要用到插值节点处的导数值，这就是还要用到插值节点处的导数值，这就是分段三次分段三次埃尔米特埃尔米特（Hermite）插值问题。）插值问题。第27页，此课件共81页哦 已知已知xi,f(xi),f(xi)(i=0,1,2,n)，求三次插值函数求三次插值函数H(x)H(x)满足：满足：H(xi)=f(xi)，H(xi)=f (xi)(i=0,1,2,n)在任意的子区间在任意的子区间Xi,Xi+1,i(0,1,2,n-1)上上构造插值函数构造插值函数 H(x)
17、=f(xi)h1(x)+f(xi)h2(x)+f (xi)h3(x)+f (xi+1)h4(x)hk(x)(k=1,2,3,4)称为称为埃尔米特插值的插值基函数埃尔米特插值的插值基函数。根据插值函数根据插值函数H(x)所满足的条件，可以确定插值基函数，从所满足的条件，可以确定插值基函数，从而确定插值函数。而确定插值函数。第28页，此课件共81页哦 x,x1,2,3,4)k(x)h(x),hH(x),.,1,0(),()(),()(:H(x)1iikk上所满足的插值条件点在插值节（构可以导出的函数结及所满足的插值条件根据nixfxHxfxHiiii第29页，此课件共81页哦 0)(xh)(xh)
18、(xh 1)(xh)(xf)(x)h(xf )(x)h(xf)(x)hf(x)(x)hf(x)(xH0)(xh)(xh)(xh 1)(xh)f(x)(x)h(xf )(x)h(xf)(x)hf(x)(x)hf(x)H(xxxiiiiii1iiii1iiiii4i3i2i1ii41ii3ii21ii1iii42134321时当第30页，此课件共81页哦 0)(xh)(xh)(xh 1)(xh)(xf)(x)h(xf )(x)h(xf)(x)hf(x)(x)hf(x)(xH0)(xh)(xh)(xh 1)(xh)f(x)(x)h(xf )(x)h(xf)(x)hf(x)(x)hf(x)H(xxx1
19、i1i1i1i1i1i1i1ii1i1i1ii1i1i41i31i11i21i1i41i1i3i1i21i1i1i1i1i32144321时当第31页，此课件共81页哦(x)h(x)h0)(xh)(xh)(xh)(xh 0)(xh)(xh)(xh)(xh 0)(xh)(xh)(xh)(xh 1)(xh)(xh)(xh)(xh :kk1i4i41i3i31i2i21i1i1i1i3i21i11ii31i2i144定可以确在插值点上的取值，就确定了得到第32页，此课件共81页哦31i1i31i1iii21iii121iiii1ii121i1121i111i1i11i1)(x3xa )(x2b 0)
20、x-x)bx2(a)x-b(x)(xh 1)x-x)bx(a)(xh ba,0)(xh1)(xh)x-xbx)(a(x)h(x)h)x-x(x)h(x)hx0)(xh)(xh 1iiixxx（可以确定，由（次多项式，我们令为三因子，又因为含有（或者说的二重根，为，可知由第33页，此课件共81页哦 21111)(21()(iiiiiixxxxxxxxxh21112)(21()(iiiiiixxxxxxxxxh第34页，此课件共81页哦 2114)()(iiiixxxxxxxh2113)()(iiiixxxxxxxh第35页，此课件共81页哦分段三次埃尔米特插值的算法：输入插值节点输入插值节点xi
21、(i=0,n)计算插值节点上的原函数计算插值节点上的原函数f(xi)和其和其输入插值点的个数输入插值点的个数mfor i=0 to n-1 h=(xi+1-xi)/(m-1)for k=0 to m x=xi+h*k ll(k)=h1f(xi)+h2f(xi+1)+h3f(xi)+h4f(xi+1)yy(k)=1/(1+x2)next k 绘制曲线绘制曲线 next i 21111)(21()(iiiiiixxxxxxxxxh21112)(21()(iiiiiixxxxxxxxxh2113)()(iiiixxxxxxxh2114)()(iiiixxxxxxxh第36页，此课件共81页哦 “样条
22、样条”是指在工业设计过程中为了描绘出光滑是指在工业设计过程中为了描绘出光滑的外形曲线所用的一种工具，即一个具有弹性的细的外形曲线所用的一种工具，即一个具有弹性的细长木条。长木条。绘图员首先把一串有序点标在平面上，让样条依绘图员首先把一串有序点标在平面上，让样条依次通过这串点，并在每个点上压一块重的压铁，然次通过这串点，并在每个点上压一块重的压铁，然后用铅笔沿着这样一根样条绘出一条光滑的曲线。后用铅笔沿着这样一根样条绘出一条光滑的曲线。第37页，此课件共81页哦 人们由此得到启发：如果建立起样条的数学模型，并人们由此得到启发：如果建立起样条的数学模型，并编成程序存放在计算机中，那么只要输入一串有
23、序点编成程序存放在计算机中，那么只要输入一串有序点列的坐标，再经过一些简单的计算，绘图机便能划出列的坐标，再经过一些简单的计算，绘图机便能划出一条光滑的曲线。一条光滑的曲线。这样的数学模型经简化后为分段三次多项式，在节这样的数学模型经简化后为分段三次多项式，在节点处左右两端的曲线的切线和曲率是连续的，这种点处左右两端的曲线的切线和曲率是连续的，这种模型称为三次样条函数。模型称为三次样条函数。第38页，此课件共81页哦 定义：给定定义：给定a,ba,b的划分：的划分：a=xa=x0 0 xx1 1xxn n=b,=b,如果函数如果函数S(x)S(x)在在区间区间a,ba,b上满足以下条件：上满足
24、以下条件：1.1.在每一个子区间（在每一个子区间（X Xi i,X,Xi+1i+1）(i=1,2,n-1)(i=1,2,n-1)上上S(x)S(x)是三是三次多项式；次多项式；2.S(x)2.S(x)在区间在区间a,ba,b上具有连续的二阶导数；上具有连续的二阶导数；3.S(x)=Y3.S(x)=Yi i(i=0,1,2,n),S(i=0,1,2,n),S(x(x0 0)=y)=y 0 0,S,S(x(xn n)=y)=y n n我们称我们称S(x)S(x)为三次样条函数。为三次样条函数。)()()()()(413211xsmxsmxsyxsyxsiiii第39页，此课件共81页哦其中其中21
25、111)(21()(iiiiiixxxxxxxxxs21112)(21()(iiiiiixxxxxxxxxs2113)()(iiiixxxxxxxs2114)()(iiiixxxxxxxs第40页，此课件共81页哦第41页，此课件共81页哦其中其中第42页，此课件共81页哦)(2)(42 )(2)(2)(22)()(2(2)(2(1 )(1)(21()(21()(13113 12113121321221111iiiiiiiiiiiiiiiiiiiiiiiiiiiixxhxxhxxxxhxxhsxxxxxxhhsxxxxhhxxhxxhxxxxxxxxxs第43页，此课件共81页哦)(2)(42
26、 )(2)(2)(22)()(2(2)(2(1 )(1)(21()(21()(1313 2213221322121112iiiiiiiiiiiiiiiiiiiiiiiiiiiixxhxxhxxxxhxxhsxxxxxxhhsxxxxhhxxhxxhxxxxxxxxxs第44页，此课件共81页哦)()(22 )()()(2)()(2)(1)()(1 )()(12112 3121232122113iiiiiiiiiiiiiiiiiixxxxhxxxxxxhxsxxxxxxhxsxxxxhxxxxxxxs第45页，此课件共81页哦)()(22 )()()(2)()(2)(1)()(1 )()(121
27、2 412242122114iiiiiiiiiiiiiiiiiixxxxhxxxxxxhxsxxxxxxhxsxxxxhxxxxxxxs第46页，此课件共81页哦002004210200302202030100210300104103002101000011i0i10ym00h1)(xs10)(h1)(xs0)(xs0)h-(hh2 )x-(xh)x-(xh2)(xs)(xsm)(xsm)(xsy)(xsy y)(xsxxx x,xxx,x得时当上，在区间xx第47页，此课件共81页哦nnn4n3n2n1n4nn31-nn2nn11-nnnnn1i1-nin1-nym1)(xs0)(xs0)(
28、xs0)(xs)(xsm)(xsm)(xsy)(xsy y)(xsxxx x,xxx,x得时当上，同样在区间第48页，此课件共81页哦1000120020-1 01101201 400111201 32011001301 22001011301 1 411 01 211 0-1 11i0i10-142y6y6)(xs4)()(2(h2)(xs2)()(2(h2)(xs6)(2)(4(h2)(xs6)(2)(4(h2)(xs)(xsm)(xsm)(xsy)(xsy)(xsxx,xxx,xxx131mhmhhhhxxxxhxxxxhxxhxxhxxhxx中，时，即在区间当第49页，此课件共81页哦
29、21112211211 02111211 411121211 32121111311 22111121311 1 421 11 221 11 21i1i21124y6y6)(xs2)()(2(h2)(xs4)()(2(h2)(xs6)(2)(4(h2)(xs6)(2)(4(h2)(xs)(xsm)(xsm)(xsy)(xsy)(xsxx,xxx,xxx131mhmhhhhxxxxhxxxxhxxhxxhxxhxx中，时，即在区间当第50页，此课件共81页哦)y-yy-y(312)y-(y6)y-(y62)411(22-4-y6y6-42y6y6)(xs)(xs20012112211010000
30、12012212111000211122112110001200201-1 hhhmmhhhhhmhhmhmhhmhmhmhhhmhmhhh由第51页，此课件共81页哦)y-yy-y(32)y-yy-y(31200110111210021001010110102001211221101000hhhhhhhhmhhhmmhhhhhhhhhhmmhhhhhm式两边同乘第52页，此课件共81页哦110010101i1-i1-ii1001001101112100210010101-1-,)y-yy-y(32hhhhhhhhhhhhhhhhhhhhhhmhhhmmhhh令第53页，此课件共81页哦得到方
31、程组）（）（令）（）（则上式简化为1-n1,2,.,i )y-yy-y-1(3 )y-yy-y-1(3)y-yy-y-1(3 2-1ii1ii1-i1-iiii1121001111121001121101hhhhhhmmm第54页，此课件共81页哦1-nn1-n121-n2-n1-n2-n232-n23221212110100 2-1 2-1.2-12-1nnnnnnymmmmmmmmmmmmmym）（）（）（）（得到方程组第55页，此课件共81页哦第56页，此课件共81页哦第57页，此课件共81页哦第58页，此课件共81页哦第59页，此课件共81页哦 例：要在程控铣床上加工直升飞机的旋转机翼
32、，给出了例：要在程控铣床上加工直升飞机的旋转机翼，给出了1818个点的外形轮廓，以及头部外形的轮廓为半径个点的外形轮廓，以及头部外形的轮廓为半径R=6.92mmR=6.92mm，TanTan 0.305,B1(0.52,5.288),B2(2.6,-3.615).0.305,B1(0.52,5.288),B2(2.6,-3.615).第60页，此课件共81页哦 第61页，此课件共81页哦 解解:要在程控铣床上加工工件要在程控铣床上加工工件,必须计算出整个外形必须计算出整个外形曲线足够密的点的坐标值。根据所给条件，头部圆弧曲线足够密的点的坐标值。根据所给条件，头部圆弧B B1 1B B2 2可由
33、圆的方程直接计算出点的坐标，其余部分必可由圆的方程直接计算出点的坐标，其余部分必须根据给出的点，利用插值方法计算加密点的坐标。须根据给出的点，利用插值方法计算加密点的坐标。第62页，此课件共81页哦现在先讨论头部圆弧现在先讨论头部圆弧B1B2的计算方法。根据所给条的计算方法。根据所给条件，可求得圆心坐标（件，可求得圆心坐标（X0，Y0)如下如下:第63页，此课件共81页哦13.5396.6196.92x-0.3016.619-6.92x6.926.619-x47.886)619.6(0)619.6(886.472.0186yy0186.2)619.6(886.470186.2)619.6(88
34、6.47)619.6(886.470186.2)619.6(886.47)0186.2(2122212221222xxxyxyxyxy时当第64页，此课件共81页哦 下面主要讨论如何根据表函数，应用插值多项式，下面主要讨论如何根据表函数，应用插值多项式，计算曲线上点的坐标。由于上下轮廓处理方法是计算曲线上点的坐标。由于上下轮廓处理方法是一样的，因此我们只讨论上轮廓的计算。一样的，因此我们只讨论上轮廓的计算。我们己知上轮廓线我们己知上轮廓线1818个点的坐标，加上与头部连个点的坐标，加上与头部连接的接的B B1 1，点共有，点共有1919个点。为了在整个区间加密函数个点。为了在整个区间加密函数值
35、，是不是必须按上面所讲的值，是不是必须按上面所讲的1919个点做一个个点做一个1818次次的插值多项式的插值多项式?这取决于对插值函数的性质分析。这取决于对插值函数的性质分析。在解决实际问题中，关键的一点就是要善于对具在解决实际问题中，关键的一点就是要善于对具体问题进行具体的析。为此，先作出差商表以便体问题进行具体的析。为此，先作出差商表以便分析。分析。第65页，此课件共81页哦 从差商表可以看到，除插值区间两端外，四阶差商从差商表可以看到，除插值区间两端外，四阶差商均近似于均近似于0 0，而在区间，而在区间28.65,50728.65,507中，三阶差商也中，三阶差商也近似于近似于0 0，这
36、说明被插值函数并不是一条高次曲线。，这说明被插值函数并不是一条高次曲线。实际上，三阶差商也近似于实际上，三阶差商也近似于0 0的部分，曲线与二次的部分，曲线与二次抛抛线近似。如果用高次多项式近似，不仅计算量抛抛线近似。如果用高次多项式近似，不仅计算量增大，效果反而不好。增大，效果反而不好。因此，对这种情况我们一般不采用高阶插值，也就是说因此，对这种情况我们一般不采用高阶插值，也就是说不去做不去做1818次插值多项式，而是采用分段插值的办法。次插值多项式，而是采用分段插值的办法。第66页，此课件共81页哦 即根据被插值函数的不同变化趋势，把整个插值区间即根据被插值函数的不同变化趋势，把整个插值区
37、间分成若干段，分别采用低次插值多项式做各段的近似函分成若干段，分别采用低次插值多项式做各段的近似函数，然后再由各段插值多项计算各段加密点的坐标。数，然后再由各段插值多项计算各段加密点的坐标。根据本例给出的曲线图形和差商表的情况，这个问根据本例给出的曲线图形和差商表的情况，这个问题可将区间分成三段考虑，即题可将区间分成三段考虑，即0.52,28.65,28.650.52,28.65,28.65，507,507507,507，520520。第67页，此课件共81页哦 现在先考虑区间现在先考虑区间28.65,507 28.65,507。由于三阶差商近似于由于三阶差商近似于0 0，故可考虑用二次插值，
38、即每相，故可考虑用二次插值，即每相邻三个点邻三个点X Xk k,X,Xk+1k+1,X,Xk+2k+2为插值点做一个二次插值多项式，为插值点做一个二次插值多项式，作为作为XXk k,X,Xk+2k+2 区间上函数区间上函数Y=f(x)Y=f(x)的近似值。的近似值。这样处理能否满足精度要求，还要估计误差。我们以这样处理能否满足精度要求，还要估计误差。我们以区间区间x4x4，x6x6为例进行分析，以插值点为例进行分析，以插值点 x x4 4=28.65=28.65，X X5 5=39.62=39.62，x x6 6=50.65=50.65，确定二次插值多项式，由，确定二次插值多项式，由Newdo
39、nNewdon插值多项式，有插值多项式，有:第68页，此课件共81页哦第69页，此课件共81页哦基本上能满足加工精度要求。从差商表看到，在区间基本上能满足加工精度要求。从差商表看到，在区间xx4 4，x x1717 中三阶差商都比中三阶差商都比0.0000230.000023小。因此，不必每个区间都小。因此，不必每个区间都进行分析，即可认为在这一段采用二次插值是合适的进行分析，即可认为在这一段采用二次插值是合适的。由。由NewdonNewdon插值公式得：插值公式得：这里这里K=4K=4，5 5，1616。用这个公式计算区间。用这个公式计算区间XX4 4,X,X1717 上的函数值上的函数值，
40、即用，即用Y=PY=P2 2(x)(x)的值作为的值作为Y=f(x)Y=f(x)的近似值。这就解决了区间的近似值。这就解决了区间X4,X17这一段曲线点的加密计算问题。这一段曲线点的加密计算问题。第70页，此课件共81页哦 下面讨论与头部连接部分的计算，也就是区间下面讨论与头部连接部分的计算，也就是区间xx0 0，x x4 4 这一段。这一段。当区间当区间x=xx=x0 0时曲线与圆弧时曲线与圆弧B Bl lB B2 2相连接，通常除了要求函相连接，通常除了要求函数值相等外，还必须要求在数值相等外，还必须要求在B B1 1处相切，即导数值必处相切，即导数值必须相等。因此，我们将区间须相等。因此
41、，我们将区间xx0 0，x x4 4 分成两段分成两段xx0 0，x x2 2 和和xx2 2，x x4 4 讨论。讨论。首先讨论首先讨论xx0 0，x x2 2 区间。在区间区间。在区间xx0 0，x x2 2 上构造插值多上构造插值多项式项式P(x)P(x)使其满足条件使其满足条件:第71页，此课件共81页哦第72页，此课件共81页哦第73页，此课件共81页哦第74页，此课件共81页哦第75页，此课件共81页哦第76页，此课件共81页哦通过上面分段插值的方法，可以把整个机翼截面曲通过上面分段插值的方法，可以把整个机翼截面曲线加密点的坐标计算出来。具体计算只要根据各段线加密点的坐标计算出来。具体计算只要根据各段所得到的让算公式拍成程序，在计算机上进行计算所得到的让算公式拍成程序，在计算机上进行计算即可。即可。第77页，此课件共81页哦 至于计算结果是否满足要求，还应通过实践检验。至于计算结果是否满足要求，还应通过实践检验。分段插值只能保证插值曲线的连续性，不能保证曲线分段插值只能保证插值曲线的连续性，不能保证曲线的光滑性，即导数不连续。的光滑性，即导数不连续。如果对精度要求高，就可以采用样条函数。如果对精度要求高，就可以采用样条函数。第78页，此课件共81页哦第79页，此课件共81页哦第80页，此课件共81页哦感谢大家观看感谢大家观看第81页，此课件共81页哦