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

相关文章

Java的栈与队列实现代码解析

《Java的栈与队列实现代码解析》栈是常见的线性数据结构,栈的特点是以先进后出的形式,后进先出,先进后出,分为栈底和栈顶,栈应用于内存的分配,表达式求值,存储临时的数据和方法的调用等,本文给大家介绍J... 目录栈的概念(Stack)栈的实现代码队列(Queue)模拟实现队列(双链表实现)循环队列(循环数组

C++如何通过Qt反射机制实现数据类序列化

《C++如何通过Qt反射机制实现数据类序列化》在C++工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作,所以本文就来聊聊C++如何通过Qt反射机制实现数据类序列化吧... 目录设计预期设计思路代码实现使用方法在 C++ 工程中经常需要使用数据类,并对数据类进行存储、打印、调试等操作。由于数据类

Python实现图片分割的多种方法总结

《Python实现图片分割的多种方法总结》图片分割是图像处理中的一个重要任务,它的目标是将图像划分为多个区域或者对象,本文为大家整理了一些常用的分割方法,大家可以根据需求自行选择... 目录1. 基于传统图像处理的分割方法(1) 使用固定阈值分割图片(2) 自适应阈值分割(3) 使用图像边缘检测分割(4)

Android实现在线预览office文档的示例详解

《Android实现在线预览office文档的示例详解》在移动端展示在线Office文档(如Word、Excel、PPT)是一项常见需求,这篇文章为大家重点介绍了两种方案的实现方法,希望对大家有一定的... 目录一、项目概述二、相关技术知识三、实现思路3.1 方案一:WebView + Office Onl

C# foreach 循环中获取索引的实现方式

《C#foreach循环中获取索引的实现方式》:本文主要介绍C#foreach循环中获取索引的实现方式,本文给大家介绍的非常详细,对大家的学习或工作具有一定的参考借鉴价值,需要的朋友参考下吧... 目录一、手动维护索引变量二、LINQ Select + 元组解构三、扩展方法封装索引四、使用 for 循环替代

Spring Security+JWT如何实现前后端分离权限控制

《SpringSecurity+JWT如何实现前后端分离权限控制》本篇将手把手教你用SpringSecurity+JWT搭建一套完整的登录认证与权限控制体系,具有很好的参考价值,希望对大家... 目录Spring Security+JWT实现前后端分离权限控制实战一、为什么要用 JWT?二、JWT 基本结构

Java实现优雅日期处理的方案详解

《Java实现优雅日期处理的方案详解》在我们的日常工作中,需要经常处理各种格式,各种类似的的日期或者时间,下面我们就来看看如何使用java处理这样的日期问题吧,感兴趣的小伙伴可以跟随小编一起学习一下... 目录前言一、日期的坑1.1 日期格式化陷阱1.2 时区转换二、优雅方案的进阶之路2.1 线程安全重构2

Android实现两台手机屏幕共享和远程控制功能

《Android实现两台手机屏幕共享和远程控制功能》在远程协助、在线教学、技术支持等多种场景下,实时获得另一部移动设备的屏幕画面,并对其进行操作,具有极高的应用价值,本项目旨在实现两台Android手... 目录一、项目概述二、相关知识2.1 MediaProjection API2.2 Socket 网络

使用Python实现图像LBP特征提取的操作方法

《使用Python实现图像LBP特征提取的操作方法》LBP特征叫做局部二值模式,常用于纹理特征提取,并在纹理分类中具有较强的区分能力,本文给大家介绍了如何使用Python实现图像LBP特征提取的操作方... 目录一、LBP特征介绍二、LBP特征描述三、一些改进版本的LBP1.圆形LBP算子2.旋转不变的LB

Redis消息队列实现异步秒杀功能

《Redis消息队列实现异步秒杀功能》在高并发场景下,为了提高秒杀业务的性能,可将部分工作交给Redis处理,并通过异步方式执行,Redis提供了多种数据结构来实现消息队列,总结三种,本文详细介绍Re... 目录1 Redis消息队列1.1 List 结构1.2 Pub/Sub 模式1.3 Stream 结