原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验 一 .拉格朗日插值法 C 语言的实现
实验目的：
进一步熟悉拉格朗日插值法。
掌握编程语言字符处理程序的设计和调试技术。
实验要求：
已知：某些点的坐标以及点数。
输入：条件点数以及这些点的坐标 。
输出：根据给定的点求出其对应的拉格朗日插值多项式的值 。
程序流程：
输入已知点的个数；
分别输入已知点的 X 坐标；
分别输入已知点的 Y 坐标；
通过调用函数 lagrange 函数，来求某点所对应的函数值。
拉格朗日插值多项式如下：
L (x
n j
) ? ?n
k ?0
y l (x
k k j
) ? y
j
j ? 0,1,……n
(x ? x
)……(x-x
)(x-x ) …(x-x )
其中l
(x) ???0
k-1 k+1 n
k ? 0,1,……, n
k (x
k
? x )……(x
0 k
-x
k-1
)(x
k
-x
k+1
) …(x -x )
k n
开始↓
开始
↓
输入已知点个数 n↓
输入已知点个数 n
↓
输入已知点的 X 坐标以及输入已知点的 Y 坐标
↓
调用函数 lagrange 函数
↓
输出结果
输出结果
程 序 如 下 ： #include iostream #include conio.h #include malloc.h
float lagrange(float *x,float *y,float xx,int n) /*拉格朗日插值算法*/
{
int i,j;
float *a,yy=0.0; /*a 作为临时变量，记录拉格朗日插值多项式*/ a=(float *)malloc(n*sizeof(float));
for(i=0;i=n-1;i++)
{
a[i]=y[i]; for(j=0;j=n-1;j++)
if(j!=i) a[i]*=(xx-x[j])/(x[i]-x[j]); yy+=a[i];
}
free(a); return yy;
}
int main()
{
int i; int n;
float x[20],y[20],xx,yy;
printf(Input n:);
scanf(%d,n); if(n=20)
{
printf(Error!The value of n must in (0,20).); getch();return 1;
}
if(n=0)
{
printf(Error! The value of n must in (0,20).); getch(); return 1;
}
for(i=0;i=n-1;i++)
{
printf(x[%d]:,i);
scanf(%f,x[i]);
}
printf(\n); for(i=0;i=n-1;i++)
{
printf(y[%d]:,i);scanf(%f,y[i]);
}
printf(\n);
printf(Input xx:);
scanf(%f,xx); yy=lagrange(x,y,xx,n); printf(x=%f,y=%f\n,xx,yy); getch();
}
举例如下：
已知当 x=1，-1,2 时 f(x)=0，-3,4，求 f(1.5)的值。运行结果如下：
文档评论（0）