10个重要的算法C语言实现源代码(其2--牛顿插值)

最新推荐文章于 2025-01-05 00:15:00 发布

转载最新推荐文章于 2025-01-05 00:15:00 发布 · 1.2k 阅读

文章标签：

#算法 #语言 #float #c #input

C/C++ 专栏收录该内容

5 篇文章

订阅专栏

本文介绍了一种使用C语言实现的牛顿插值多项式算法，该算法适用于离散数据点的拟合问题。通过输入一系列离散的数据点，程序能够计算出任意指定点处的插值结果。

2.牛顿插值多项式，用于离散数据的拟合

#include<stdio.h> #include<conio.h>
#include<alloc.h>
void difference(float *x,float *y,int n)
{float *f;
int k,i;
f=(float *)malloc(n*sizeof(float));
for(k=1;k<=n;k++)
{ f[0]=y[k];
   for(i=0;i<k;i++)
      f[i+1]=(f[i]-y[i])/(x[k]-x[i]);
    y[k]=f[k];
}
return;
}
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");
difference(x,(float *)y,n);
printf("Input xx:");
scanf("%f",&xx);
yy=y[20];
for(i=n-1;i>=0;i--) yy=yy*(xx-x[i])+y[i];
printf("NewtonInter(%f)=%f",xx,yy);
getch();
}