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

本文主要是介绍拉格朗日插值法 C语言实现，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！

  /* 
 *作者：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("*******************************************\n");  
     printf("使用说明：\n1.用户输入插值点，每一行输入一组：x y;\n2.输入换行表示输入结束。\n");  
     printf("*******************************************\n");  
     printf("x    y\n");  
   
     //接收用户输入知道第一次输入非换行为止  
     memset(str,0,sizeof(str));  
     while(strlen(str)==0)  
         gets(str);  
     //数据输入完毕，输入换行结束输入  
     while(strlen(str)!=0)  
     {  
         newData = (Data *)malloc(sizeof(Data));  
         sscanf(str,"%s%s",strx,stry);     //获取输入的前两个字符串 第一个为x,第二个为y  
         newData->x = strtod(strx,NULL);   //将输入转换成浮点数  
         newData->y = strtod(stry,NULL);     
         newData->next=NULL;  
         p->next=newData;  
         p = p->next;  
         gets(str);  
     }  
     printf("请输入要计算的x值：");  
     scanf("%lf",&x);  
     printf("L(%f) = %f\n",x,LagrangeInsert(header,0.20));  
     return ;  
 

[cpp]  view plain copy  
 print ? 
 }  

  顶 2    踩  

这篇关于拉格朗日插值法 C语言实现的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！

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

相关文章

MySQL双主搭建+keepalived高可用的实现

Java实现文件图片的预览和下载功能

使用Sentinel自定义返回和实现区分来源方式

Java实现时间与字符串互相转换详解

opencv图像处理之指纹验证的实现

Springboot处理跨域的实现方式(附Demo)

Spring Boot 3.4.3 基于 Spring WebFlux 实现 SSE 功能(代码示例)

基于SpringBoot实现文件秒传功能

SpringBoot日志配置SLF4J和Logback的方法实现

Python如何使用slots实现节省内存和性能优化