计算机数学基础(下)数值部分辅导.doc

《计算机数学基础(下)》数值部分辅导(3) 中央电大 冯泰 第11章 函数插值与最小二乘拟合 一、重点内容 1. 函数插值 已知函数f(x)的n个函数值yk=f(xk), k=0,1,2,…,n。构造一个多项式P(x)，使得 P(xk)=yk。P(x)就是插值多项式，f(x)就是被插函数，xk就是插值节点。误差R(x)=f (x)－P(x)。 2. 拉格朗日多项式 用n次多项式 Pn(x)=y0l0+y1l1+…+ynln= 其中基函数 ， 当n=1时，线性插值 P1(x)=yk lk (x)+yk+1 lk+1 (x) 其中基函数 , 。 当n=2时，得到二次多项式，就是二次插值。 拉格朗日插值多项式的余项为 其中 注意：过n+1个互异点，所得的多项式应该是次数不超过n的多项式。 3. 均差与牛顿插值多项式 函数值与自变量的差商就是均差， 一阶均差 (或记作f [x0,x1]); 二阶均差 (或记作f [x0,x1,x2]) 均差有两条常用性质：(1)均差用函数值的线性组合表示；(2)均差与插值节点顺序无关。 用均差为系数构造多项式，就是牛顿插值多项式 Nn(x)= f(x0)+f(x0,x1)(x-x0)+f(x0,x1,x2)(x-x0)(x-x1)+ …+f(x0,x1,x2,…,xn)(x-x0)(x-x1)(x-x2)…(x-xn-1) 牛顿插值多项式的余项为 Rn(x)=f(x)－Nn(x) =f(x,x0,x1,x2,…,xn)(x－x0)(x－x1)(x－x2)…(x－xn－1)(x－xn) 4. 分段线性插值 已知n+1个互异节点x,x1,…,xn,构造一个分段一次的多项式P(x)，且满足：(1)P(x)在[a ,b]上连续; (2) P(xk)=yk(k=0,1,2,…,n); (3)P(x)在[xk ,xk+1]上是线性函数。 分段线性插值函数 其中lk(x)(k=0,1,2,…,n)是分段线性插值基函数. li(x)= ln(x)= 5. 三次样条插值函数 其中S((xk)=mk(k=0,1,2,…,n), hk=xk+1－xk(k=0,1,2,…,n－1)…,mn满足的方程组是 (*) 其中： , (k=1,2,…,n－1) (1) 当已知S((x0)=y(0 ,S((xn)=(y(n时，(*)式中 (0=1, (n=1, . (2) 当已知S((x0)=y(0=m0, S((xn)=y(n=mn时，(*)式化为 6. 最小二乘法 用((x)拟合数据(xk,yk) (k=1,2,…,n)，使得误差的平方和 为最小，求((x)的方法，称为最小二乘法。 (1) 直线拟合 若，a0,a1满足法方程组 (2) 二次多项式拟合 若满足法方程组 二、实例 例1 已知函数y=f(x)的观察数据为 xk －2 0 4 5 yk 5 1 －3 1 试构造拉格朗日多项式Pn (x)，并计算P(－1)。 [只给4对数据，求得的多项式不超过3次] 解 先构造基函数 所求三次多项式为 P3(x)= ＝＋－＋ ＝ P3(－1)＝ 例2 已知函数y=f(x)的数据如表中第1，2列。计算它的各阶均差。 解 依据均差计算公式，结果列表中。 k Xk f(xk) 一阶均差 二阶均差 三阶均差 四阶均差 0 0.40 0.410 75 1 0.55 0.578 15 1.116 00 2 0.65 0.696 75 1.168 00 0.280 00 3 0.80 0.888 11 1.275 73 0.358 93 0.197 33 4 0.90 1.201 52 1.384 10 0.433 48 0.213 00 0.031 34 计算公式为 一阶均差 二阶均差 三阶均差 四阶均差 例3 设是n+1个互异的插值节点，是拉格朗日插值基函数，证明： (1) (2) 证明 (1) Pn(x