拉格朗日插值多项式，用于离散数据的拟合

最新推荐文章于 2024-09-11 20:00:00 发布

zyming0815

最新推荐文章于 2024-09-11 20:00:00 发布

阅读量1.8k

点赞数

本文介绍了一个使用C语言实现的拉格朗日插值算法程序。该程序允许用户输入一系列节点坐标，通过拉格朗日插值法计算并输出在指定位置的插值结果。代码中包含了内存分配、计算和释放等步骤。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include <stdio.h>
#include <conio.h>
#include <alloc.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;
}
main()
{ int i,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();
}