拉格朗日插值法實(shí)驗(yàn)報(bào)告

1、拉格朗日插值法實(shí)驗(yàn)報(bào)告一、實(shí)驗(yàn)?zāi)康?、學(xué)習(xí)和掌握拉格朗日插值多項(xiàng)式。2、運(yùn)用拉格朗日插值多項(xiàng)式進(jìn)行計(jì)算。二、算法步驟函數(shù)定義double lgrr（int n,double X）輸入所求節(jié)點(diǎn)x和已知的節(jié)點(diǎn)數(shù)n輸出插值結(jié)果s步驟：1、輸入 n 對數(shù)（xi，yi）（ i=0-n）2、令 s=03、對于 i=0-n T=yi對 j=0- n 但 j! = i T=T*(X-xj)/(xi-xj)令 s=s+T4、輸出結(jié)果三流程圖四、原程序#include main()(double X;int n;double lgrr(int n,double X);printf(input the aim nu

2、mberX);printf(input the number of pairs of numbers n:);scanf(%lf%d,&X,&n);lgrr(n,X);return 0;double lgrr(int n,double X)(double s=0,t;double x81,y81;int i,j;printf(input n pairs of numbers:n);for(i=0;in;i+)scanf(%lf%lf,&xi,&yi);for(i=0;in;i+)(t=yi;for(j=0;jn;j+)(if(j! = i)t=t*(X-xj)/(xi-xj);s=s+t;pr

3、intf(s=%fn,s);return 0;作業(yè)用程序#includemain()(double X;int n;double lgrr(int n,double X);printf(input the aim numberX);printf(input the number of pairs of numbers n:);scanf(%lf%d,&X,&n);lgrr(n,X);return 0;double lgrr(int n,double X)(FILE *fp;double s=0,t,z;double x81,y81;int i,j,k;if(fp=fopen(lgrr,txt,

4、r+)=NULL) printf(can not open the file.); printf(input n pairs of numbers:n);for(i=0;in;i+)scanf(%lf%lf,&xi,&yi);for(k=0;k=20;k+)(for(i=0;in;i+)(z=X+.01*k;t=yi;for(j=0;jn;j+)(if(j! = i)t=t*(z-xj)/(xi-xj);s=s+t;printf(s=%fn,s);fprintf(fp,s=%fn,s);return 0;五、作業(yè)：3-27s=0.198670s=0.407105s=0.625307s=0.853266s=1.090966s=1.338384s=1.595493s=1.862263s=2.138663s=2.424661s=2.720224s=3.025321s=3.339921s=3.663993s=