拉格朗日插值上机实验报告

《拉格朗日插值上机实验报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、计算机科学与技术学院姓名：程伟业学号：U200714790sonystyl一、题目：完成L■插值公式的编程，并由下表求f(x*Lx*=431.X1234567f(x)83-12-145二、程序设计与分析L■插值公式：f(x)u^•Lof(xi)Li(x)输入n,Xi，yi，x*i=0s=0计算Li(x)s=s+Li(x)yji=i+l输出s三、原代码与运行结果代码如下：/*完成拉格朗日插值公式编程*/#includefloatM(floatx[],floatyjnti,intn){floatm=l;intj;for(j=0;j<二n;j++)if(j

2、!二i)m=m*(y-x[j]);returnm;}/*构造函数M=(x-x[O])*(x-x[l])(x-x[i-l])(x-x[i+l])(x-x[n])*/voidmain(){floats=0,x;/*s用来记录结果，x是待估计的点*/intn,k,i=0;floata[20],b[20];/*用数组存储插值节点及其函数值*/printf("Pleaseentern:");scanf「％d”,&n);/*输入n,n为插值节点的个数*/k=n-l;printf("Pleaseenterx,y:");for(i=0;i<=k;i++)scanf(“％f,%

3、f“,&a[i],&b[i]);/*读入插值节点值*/printf(''Pleaseenterx*:");scanf(吆化&x);/*输入x*的值*/for(i=0;i<=k;i++)s=s+M(azx,i,k)/M(a,a[i]/i/k)*b[i];/*Largrange插值公式*/printf("Theresultis%f",s);}运行结杲：Pleaseentern:7Pleaseenterx,y:1,82.33,-14,25,-16.47.5Pleaseenterx*:4.31Theresultis1.512093