第二章代数差值.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 5.3 三次样条插值函数的误差估计 1. 如果f(x)?C[a,b]，且划分的网格比 * 其中 一致有界，则当 时，S(x)一致收敛于f(x). 2. 如果f(x)?C4[a,b]，且S(x)满足边界条件I、II, 则 即h?0时， 如果划分比?一致有界， 5.4 三次样条插值函数的程序设计 1. 三对角方程求解的追赶法 * #include stdio.h #include math.h int cetrd(double b[],int n,int m,double d[]) { int k,j; double s; if(m!=(3*n-2)) { printf(error\n); return(-2); } for(k=0;k=n-2;k++) { j=3*k; s=b[j]; if(fabs(s)+1.0==1.0) { printf(fail\n); return(0);} b[j+1]=b[j+1]/s; d[k]=d[k]/s; b[j+3]=b[j+3]-b[j+2]*b[j+1]; d[k+1]=d[k+1]-b[j+2]*d[k]; } s=b[3*n-3]; if(fabs(s)+1.0==1.0) { printf(fail\n);return(0); } d[n-1]=d[n-1]/s; for(k=n-2;k=0;k--) d[k]=d[k]-b[3*k+1]*d[k+1]; return(2); } * 2. 调用追赶法解方程组 #include stdio.h #include cetrd.c main() { int i; static double b[13]={13.0,12.0,11.0,10.0,9.0,8.0,7.0,6.0,5.0,4.0,3.0,2.0,1.0}; static double d[5]={3.0,0.0,-2.0,6.0,8.0}; i=crtrd(b,5,13,d); if(i=0;i=4;i++) printf(“x(%d)=%e\n”,i,d[i]); } 本节(§5)问题 * 一、何为三次样条函数？三次样条插值多项式有何优点？ 二、如何通过 Hermite 插值多项式求解三次样条插值函数？ 三、在构造三次样条插值函数时，如何使用边界条件？ 四、如何在计算机上实现三次样条函数的近似计算？ 五、上机实习任务 1.分段线性插值； 2.分段Hermite插值 3.三次样条插值：边界条件I、边界条件II 练 习 二 * 2-1 当 x=1，-1，2时，f(x) 分别为 0，-3，4，求f(x) 的二次插值多项式 p2(x) 。 2-2 设 li(x) 是以xk=x0+kh, k=0,1,2,3 为插值节点的3次插值基函数，求 。 2-3 设 l0(x), l1(x), …, ln(x) 是以 x0, x1, … ,xn 为节点的n次Lagrange插值基函数，求证 （1） （2） * 2-4 设 f(x)∈C2[a,b] ，且 f(a)=f(b)=0 ， 其中 证明 2-5 利用 在 x=100,121,144点的函数值，用插值方法求 的近似值，并由误差公式给出误差界，同时与实际误差作比较。 2-6 在 -4≤x≤-4 上给出f(x)=ex 的等距节点函数表， 利用与距离最近的三点作二次插值作为的e2 近似，使其误差不超过10-6，问函数表的步长应多少？ * 2-7 证明n阶差商有下列性质 (1) 若F(x)=f(x)+g(x) ，则 F[x0 ,x1 ,…,xn] =f [x0 ,x1 ,…,x n]+ g [x0 ,x1 ,…,x n] （2）若 f(x) ∈Pm (m次多项式), m≥ n 则 f[x0 ,x1 ,…,x n-1 ,x] ∈ Pm-n. 2-8 f(x)= x5+4x4+3x+1 ,求差商 f [20,21,..,25 ] 和 f[20,21,..,26] . 2-9 设 f(x)=x5+x