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

2024-05-11 11:18

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

  1. /* 
  2. *作者:KDF5000 
  3. *功能:利用拉格朗日插值法求解近似值 
  4. *时间:2013.4.15 
  5. */  
  6. #include <stdio.h>  
  7. #include <stdlib.h>  
  8. #include <string.h>  
  9. //存放插值节点  
  10. struct Data{  
  11.     double x;  
  12.     double y;  
  13.     struct Data *next;  
  14. };  
  15. /**************************************************** 
  16. *LagrangeInsert() 
  17. *功能:拉格朗日插值法 
  18. *****************************************************/  
  19. double LagrangeInsert(struct Data *header,double x)  
  20. {  
  21.     Data *pi,*pj,*p;  
  22.     pi=pj=header->next;  
  23.     double temp1,temp2;  
  24.     temp1=0;                //记录内循环的积  
  25.     temp2=1;                //记录外循环的和  
  26.     while(pi!=NULL)  
  27.     {  
  28.         while(pj!=NULL)  
  29.         {  
  30.             if(pi!=pj)  
  31.                 temp2 *=(x-pj->x)/(pi->x-pj->x);  
  32.             pj = pj->next;  
  33.         }  
  34.         temp1 +=temp2*pi->y;  
  35.         temp2=1;  
  36.         pj = header->next;  
  37.         pi = pi->next;  
  38.     }  
  39.     return temp1;   //返回计算结果  
  40. }  
  41.   
  42. void main()  
  43. {  
  44.     Data *header = (Data *)malloc(sizeof(Data));  
  45.     char str[20];  
  46.     Data *p,*newData;  
  47.     char strx[20],stry[20];  
  48.     double x;  
  49.   
  50.     p=header;  
  51.     p->x=0;  
  52.     p->y=0;  
  53.     p->next=NULL;  
  54.   
  55.     //输出提示信息  
  56.     printf("*******************************************\n");  
  57.     printf("使用说明:\n1.用户输入插值点,每一行输入一组:x y;\n2.输入换行表示输入结束。\n");  
  58.     printf("*******************************************\n");  
  59.     printf("x    y\n");  
  60.   
  61.     //接收用户输入知道第一次输入非换行为止  
  62.     memset(str,0,sizeof(str));  
  63.     while(strlen(str)==0)  
  64.         gets(str);  
  65.     //数据输入完毕,输入换行结束输入  
  66.     while(strlen(str)!=0)  
  67.     {  
  68.         newData = (Data *)malloc(sizeof(Data));  
  69.         sscanf(str,"%s%s",strx,stry);     //获取输入的前两个字符串 第一个为x,第二个为y  
  70.         newData->x = strtod(strx,NULL);   //将输入转换成浮点数  
  71.         newData->y = strtod(stry,NULL);     
  72.         newData->next=NULL;  
  73.         p->next=newData;  
  74.         p = p->next;  
  75.         gets(str);  
  76.     }  
  77.     printf("请输入要计算的x值:");  
  78.     scanf("%lf",&x);  
  79.     printf("L(%f) = %f\n",x,LagrangeInsert(header,0.20));  
  80.     return ;  
[cpp]  view plain copy
print ?
  1. }  

2

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



http://www.chinasem.cn/article/979383

相关文章

Vue项目的甘特图组件之dhtmlx-gantt使用教程和实现效果展示(推荐)

《Vue项目的甘特图组件之dhtmlx-gantt使用教程和实现效果展示(推荐)》文章介绍了如何使用dhtmlx-gantt组件来实现公司的甘特图需求,并提供了一个简单的Vue组件示例,文章还分享了一... 目录一、首先 npm 安装插件二、创建一个vue组件三、业务页面内 引用自定义组件:四、dhtmlx

Vue ElementUI中Upload组件批量上传的实现代码

《VueElementUI中Upload组件批量上传的实现代码》ElementUI中Upload组件批量上传通过获取upload组件的DOM、文件、上传地址和数据,封装uploadFiles方法,使... ElementUI中Upload组件如何批量上传首先就是upload组件 <el-upl

Docker部署Jenkins持续集成(CI)工具的实现

《Docker部署Jenkins持续集成(CI)工具的实现》Jenkins是一个流行的开源自动化工具,广泛应用于持续集成(CI)和持续交付(CD)的环境中,本文介绍了使用Docker部署Jenkins... 目录前言一、准备工作二、设置变量和目录结构三、配置 docker 权限和网络四、启动 Jenkins

Python3脚本实现Excel与TXT的智能转换

《Python3脚本实现Excel与TXT的智能转换》在数据处理的日常工作中,我们经常需要将Excel中的结构化数据转换为其他格式,本文将使用Python3实现Excel与TXT的智能转换,需要的可以... 目录场景应用:为什么需要这种转换技术解析:代码实现详解核心代码展示改进点说明实战演练:从Excel到

如何使用CSS3实现波浪式图片墙

《如何使用CSS3实现波浪式图片墙》:本文主要介绍了如何使用CSS3的transform属性和动画技巧实现波浪式图片墙,通过设置图片的垂直偏移量,并使用动画使其周期性地改变位置,可以创建出动态且具有波浪效果的图片墙,同时,还强调了响应式设计的重要性,以确保图片墙在不同设备上都能良好显示,详细内容请阅读本文,希望能对你有所帮助...

C# string转unicode字符的实现

《C#string转unicode字符的实现》本文主要介绍了C#string转unicode字符的实现,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随... 目录1. 获取字符串中每个字符的 Unicode 值示例代码:输出:2. 将 Unicode 值格式化

python安装whl包并解决依赖关系的实现

《python安装whl包并解决依赖关系的实现》本文主要介绍了python安装whl包并解决依赖关系的实现,文中通过图文示例介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面... 目录一、什么是whl文件?二、我们为什么需要使用whl文件来安装python库?三、我们应该去哪儿下

Python脚本实现图片文件批量命名

《Python脚本实现图片文件批量命名》这篇文章主要为大家详细介绍了一个用python第三方库pillow写的批量处理图片命名的脚本,文中的示例代码讲解详细,感兴趣的小伙伴可以了解下... 目录前言源码批量处理图片尺寸脚本源码GUI界面源码打包成.exe可执行文件前言本文介绍一个用python第三方库pi

Java中将异步调用转为同步的五种实现方法

《Java中将异步调用转为同步的五种实现方法》本文介绍了将异步调用转为同步阻塞模式的五种方法:wait/notify、ReentrantLock+Condition、Future、CountDownL... 目录异步与同步的核心区别方法一:使用wait/notify + synchronized代码示例关键

Nginx实现动态封禁IP的步骤指南

《Nginx实现动态封禁IP的步骤指南》在日常的生产环境中,网站可能会遭遇恶意请求、DDoS攻击或其他有害的访问行为,为了应对这些情况,动态封禁IP是一项十分重要的安全策略,本篇博客将介绍如何通过NG... 目录1、简述2、实现方式3、使用 fail2ban 动态封禁3.1 安装 fail2ban3.2 配