实验插值法.docVIP

实验二 插值法 实验目的： 1、掌握直接利用拉格郎日插值多项式计算函数在已知点的函数值；观察拉格郎日插值的龙格现象。 2、了解Hermite插值法、三次样条插值法原理，结合计算公式，确定函数值。 2、实验要求: PC机一台，C语言、PASCAL语言、Matlab任选 3、实验内容： 1) 用拉格郎日插值公式确定函数值；对函数f(x)进行拉格郎日插值，并对f(x)与插值多项式的曲线作比较。 已知函数表：（0.56160,0.82741）、（0.56280,0.82659）、（0.56401,0.82577）、（0.56521,0.82495）用三次拉格朗日插值多项式求x=0.5635时函数近似值。 2) 用Hermite插值公式，对函数f(x)进行三次样条插值。 已知函数表：（1,2）、（2,3）及对应导数值（1,0）、（2,－1），求出满足上述条件的两点三次插值多项式，及在1.5,1.7处的函数值。 4、题目： 插值法 5、原理： 拉格郎日插值原理： 已知函数y=f(x)在n+1个不同的点x0 ,x1 ,…,x2 上的函数值分别为y0 ,y1 ,…,yn ,求一个次数不超过n的多项式Pn (x),使其满足Pn (xi )=yi , (i=0,1,…,n),即n+1个不同的点可以唯一决定一个n次多项式。1. 插值基函数过n+1个不同的点分别决定n+1个n次插值基函数l0 (x),l1 (x),…,ln (X)每个插值基本多项式li (x)满足：(1) li (x)是n次多项式;(2) li (xi )=1,而在其它n个li (xk )=0 ,(k≠i)。由于li (xk )=0 ,(k≠i), 故有因子: (x-x0 )…(x-xi-1 )(x-xi+1 )…(x-xn )因其已经是n次多项式，故而仅相差一个常数因子。令:li (x)=a(x-x0 )…(x-xi-1 )(x-xi+1 )…(x-xn )由li (xi )=1,可以定出a, 进而得到： Hermite插值 Hermite插值是利用未知函数f(x)在插值节点上的函数值及导数值来构造插值多项式的,起其提法为:给定n+1个互异的节点x0,x1,……,xn上的函数值和导数值 求一个2n+1次多项式H2n+1(x)满足插值条件 H2n+1(xk)=yk H2n+1(xk)=yk k=0,1,2,……,n (13) 如上求出的H2n+1(x)称为2n+1次Hermite插值函数,它与被插函数 一般有更好的密合度.拉格朗日插 n + 1个点x0, x1, ... xn(x0 x1 ... xn)的函数f(x0), f(x1) , ... , f(xn)推出n次多項式p(x)，然后n次多項式p(x)求出任意的点xf(x)的算法。Hermite插值将待求的n次插值多项式Pn（x）改写为具有承袭性的形式，然后利用插值条件确定Pn（x）的待定系数，以求出所要的插值函数。#includeiostream #includemath.h #define N 4 using namespace std; double fun(double *x,double *y, int n,double p); void main() { int i; double a[N],p; double b[N]; cout please input xiangliang a= endl; for(i=0;iN;i++) { cina[i];} coutplease input xiangliang b= endl; for(i=0;iN;i++) {cinb[i];} coutplease input LagelangrichazhiJieDian p= endl; cinp; coutThe Answer= fun(a,b,N,p)endl; } double fun(double x[],double y[], int n,double p) { double z=0,s=1.0; int k,i=0; double L[N]; k=0; while(kn) { if(k==0) { for(i=1;in;i++) {s=s*(p-x[i])/(x[0]-x[i]); } L[0]=s*y[0]; k=k+1; } else { s=1.0; for(i=0;i=k-1;i++) { s=s*((p-x[i])/(x[k]-x[i])); } for(i=k+1;