拉格朗日插值法 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

相关文章

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi

让树莓派智能语音助手实现定时提醒功能

最初的时候是想直接在rasa 的chatbot上实现,因为rasa本身是带有remindschedule模块的。不过经过一番折腾后,忽然发现,chatbot上实现的定时,语音助手不一定会有响应。因为,我目前语音助手的代码设置了长时间无应答会结束对话,这样一来,chatbot定时提醒的触发就不会被语音助手获悉。那怎么让语音助手也具有定时提醒功能呢? 我最后选择的方法是用threading.Time

Android实现任意版本设置默认的锁屏壁纸和桌面壁纸(两张壁纸可不一致)

客户有些需求需要设置默认壁纸和锁屏壁纸  在默认情况下 这两个壁纸是相同的  如果需要默认的锁屏壁纸和桌面壁纸不一样 需要额外修改 Android13实现 替换默认桌面壁纸: 将图片文件替换frameworks/base/core/res/res/drawable-nodpi/default_wallpaper.*  (注意不能是bmp格式) 替换默认锁屏壁纸: 将图片资源放入vendo

C#实战|大乐透选号器[6]:实现实时显示已选择的红蓝球数量

哈喽,你好啊,我是雷工。 关于大乐透选号器在前面已经记录了5篇笔记,这是第6篇; 接下来实现实时显示当前选中红球数量,蓝球数量; 以下为练习笔记。 01 效果演示 当选择和取消选择红球或蓝球时,在对应的位置显示实时已选择的红球、蓝球的数量; 02 标签名称 分别设置Label标签名称为:lblRedCount、lblBlueCount

科研绘图系列:R语言扩展物种堆积图(Extended Stacked Barplot)

介绍 R语言的扩展物种堆积图是一种数据可视化工具,它不仅展示了物种的堆积结果,还整合了不同样本分组之间的差异性分析结果。这种图形表示方法能够直观地比较不同物种在各个分组中的显著性差异,为研究者提供了一种有效的数据解读方式。 加载R包 knitr::opts_chunk$set(warning = F, message = F)library(tidyverse)library(phyl

透彻!驯服大型语言模型(LLMs)的五种方法,及具体方法选择思路

引言 随着时间的发展,大型语言模型不再停留在演示阶段而是逐步面向生产系统的应用,随着人们期望的不断增加,目标也发生了巨大的变化。在短短的几个月的时间里,人们对大模型的认识已经从对其zero-shot能力感到惊讶,转变为考虑改进模型质量、提高模型可用性。 「大语言模型(LLMs)其实就是利用高容量的模型架构(例如Transformer)对海量的、多种多样的数据分布进行建模得到,它包含了大量的先验

Kubernetes PodSecurityPolicy:PSP能实现的5种主要安全策略

Kubernetes PodSecurityPolicy:PSP能实现的5种主要安全策略 1. 特权模式限制2. 宿主机资源隔离3. 用户和组管理4. 权限提升控制5. SELinux配置 💖The Begin💖点点关注,收藏不迷路💖 Kubernetes的PodSecurityPolicy(PSP)是一个关键的安全特性,它在Pod创建之前实施安全策略,确保P

工厂ERP管理系统实现源码(JAVA)

工厂进销存管理系统是一个集采购管理、仓库管理、生产管理和销售管理于一体的综合解决方案。该系统旨在帮助企业优化流程、提高效率、降低成本,并实时掌握各环节的运营状况。 在采购管理方面,系统能够处理采购订单、供应商管理和采购入库等流程,确保采购过程的透明和高效。仓库管理方面,实现库存的精准管理,包括入库、出库、盘点等操作,确保库存数据的准确性和实时性。 生产管理模块则涵盖了生产计划制定、物料需求计划、