拉格朗日插值法 C语言实现

C语言实现拉格朗日插值法

最新推荐文章于 2025-10-27 09:22:13 发布

原创

最新推荐文章于 2025-10-27 09:22:13 发布 · 2.7w 阅读

60 ·

CC 4.0 BY-SA版权

文章标签：

#c 计算方法数值计算

这篇博客介绍了如何用C语言实现拉格朗日插值法，用于求解数值近似值。作者KDF5000在2013年4月15日编写了这段代码，通过结构体存储插值节点，并提供了输入插值点和计算指定x值的近似值的功能。

/*
*作者：KDF5000
*功能：利用拉格朗日插值法求解近似值
*时间：2013.4.15
*/
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
//存放插值节点
struct Data{
	double x;
	double y;
	struct Data *next;
};
/****************************************************
*LagrangeInsert()
*功能：拉格朗日插值法
*****************************************************/
double LagrangeInsert(struct Data *header,double x)
{
	Data *pi,*pj,*p;
	pi=pj=header->next;
	double temp1,temp2;
	temp1=0;                //记录内循环的积
	temp2=1;				//记录外循环的和
	while(pi!=NULL)
	{
		while(pj!=NULL)
		{
			if(pi!=pj)
				temp2 *=(x-pj->x)/(pi->x-pj->x);
			pj = pj->next;
		}
		temp1 +=temp2*pi->y;
		temp2=1;
		pj = header->next;
		pi = pi->next;
	}
	return temp1;   //返回计算结果
}

void main()
{
	Data *header = (Data *)malloc(sizeof(Data));
	char str[20];
	Data *p,*newData;
	char strx[20],stry[20];
	double x;

	p=header;
	p->x=0;
	p->y=0;
	p->next=NULL;

	//输出提示信息
	printf("***************