插插方法.pptVIP

多项式拟合 对于给定的一组数据(xi,yi),(i=1,2,…,N) ，求作m次式 使 最小. 对于多项拟合，ak(k=0,1, …,m)应满足 设有某实验数据如下： x 1.36 1.73 1.95 2.28 y 14.094 16.844 18.475 20.963 试按最小二乘法求一次多项式拟合以上数据. 分析 这是直线拟合问题，直接按照直线拟的最小二乘法去做即可. 解 记x1=1.36，x2=1.73，x3=1.95，x4=2.28 y1=14.094，y2=16.844，y3=18.475 y4=20.963 设一次拟合多项为y=a+bx，则正则方程组为 其中 将以上数据代入得 解得 a=3.937 4，b=7.462 6 从而y=3.937 4+7.462 6x 求形如y=aebx（a，b为常数且a0）的经验公式，使它能和下表数据相拟合： xi 1.00 1.25 1.50 1.75 2.00 yi 5.10 5.79 6.53 7.45 8.46 分析 经验公式y=aebx不是多项式，应设法将其变为多项式，本题可通过取对数的办法将y=aebx变为lny=lna+bx，若令y‘=lny，A=lna则有y‘=A+bx，通过直线拟合的最小二乘法求出后，再变回y=aebx即可. 例 已知x=0, 2, 3, 5对应的函数值为y=1, 3, 2, 5，作三次Newton插值多项式.如再增加x=6时的函数数值为6，作四次Newton插值多项式. 解 首先构造差商表 xi f(xi) 一阶差商 二阶差商 三阶差商 0 1 2 3 1 3 2 -1 -2/3 5 5 3/2 5/6 3/10 三次Newton插值多项式为 例 设f(x)=x4 ，试用插值余项定理写出以-1,0,1,2为插值节点的三次插值多项式. 解 设f(x)以-1,0,1,2为插值节点的三次Lagrange插值多项式为L3(x)，由插值余项定理有 因而 例 设 是以x0 ,x1 ,? ,xn为节点的Lagrange插值基函数，试证： (1) (2) (3) (4) 证 (1)设f(x)= 1，则f(x)以x0 ,x1 ,? ,xn为插值节点的n次Lagrange插值多项式为 由插值余项定理知 从而 Ln(x)=f(x) 即 (2)设f(x)= xj ，则f(x)以x0 ,x1 ,? ,xn为插值节点的n次Lagrange插值多项式为 由插值余项定理知 从而 Ln(x)=f(x) 即 (3)将(xk-x)j按二项式展开，得 代入左端，得 利用(2)的结论，有 当j=0时，由(1)的结论知 当j=1, 2, …, n时，由(2)的结论知 当j=n+1时，令f(x)=xn+1 ，有f (n+1) (x)=(n+1)! f(x)以x0 ,x1 ,? ,xn为插值节点的n次Lagrange插值多项式为 由插值余项定理知 从而 即 令x=0，有 例 设 f(x)?C2[a,b] ，且f(a)= f(b)=0 ，求证: 证 以a，b为插值节点进行线性插值，其线性插值多项式为 线性插值余项为 从而 由于|(x-a)(x-b)| 在(a+b)/2 处取最大值，故 分段线性插值 Runge现象 给定函数 取等距插值节点 建立10次插值多项式 分段线性插值 设在区间[a, b]上给定n+1个插值节点a=x0 x1 … xn=b和相应的函数值y0, y1, …, yn ，求作一个插值函数φ(x) ，具有性质 1°φ(xi) ＝yi　 (i=1,2,…n)； 2°φ(x)在每个小区间[xi, xi+1]（i=0, 1, …, n-1）上是线性函数. 例 设给出了cosx的函数表(0?x