使用C++编写拉格朗日3次插值法

本文原文发表于本人博客 哔哔哔哔-使用C++编写拉格朗日3次插值法

拉格朗日插值法是由Joseph-Louis Lagrange提出的插值的多项式，是估算函数数值的重要工具之一

前言

接上文使用C++计算3次牛顿插值法

GitHub地址

https://github.com/gst-be/lagrange_polynomial

代码

#include <iostream>
using namespace std;
int main()
{double x[4],y[4],t;for(int i=0;i<=3;i++){cout<<"input x"<<i<<"=";cin>>x[i];cout<<"input y"<<i<<"=";cin>>y[i];}cout<<"input x=";cin>>t;//输入x0-x3,y0-y3,xdouble l1[4],l2[4];             //定义基函数分母和分子for(int i=0;i<=3;i++){l1[i]=l2[i]=1;          //初始化分母,分子为1}double l[4];for(int i=0;i<=3;i++){for(int j=0;j<=3;j++){if(i==j)continue;l1[i]*=t-x[j];     //计算基函数分子l2[i]*=x[i]-x[j];  //计算基函数分母}l[i]=l1[i]/l2[i];}double p=0;for(int i=0;i<=3;i++){p+=y[i]*l[i];              //计算P(x)}cout<<"P("<<t<<")="<<p<<endl;return 0;
}

解析

基函数分子,分母的计算

为方便处理，笔者在此处将分子和分母分开计算，分别命名为数组l1和数组l2
这里使用for循环来处理分子分母的计算
遇到i=j的情况时使用if-continue跳过该轮循环，如下所示：

if(i==j)
continue;

为了使l1和l2计算时的循环完整，需将l1,l2初始化为1

for(int i=0;i<=3;i++){l1[i]=l2[i]=1;}

本文已发表于哔哔哔哔(bebebe.be)

使用C++编写拉格朗日3次插值法相关推荐

matlab编写拉格朗日插值代码函数
要求:根据拉格朗日多项式插值法原理,设计算法流程并且编写拉格朗日插值代码函数. 代码如下: function[y]=lagrange(x0,y0,x) %建立一个函数名为lagrange的函数,输 ...
理解插值法（拉格朗日、牛顿插值法）
引言我们首先理解下插值法主要用来做什么事:插值法就是利用已知的点建立合适的插值函数 f ( x ) f(x) f(x) ,未知点 x i x_i xi 由插值函数 f ( x ) f(x) f(x ...
拉格朗日插值的优缺点_拉格朗日与牛顿插值法的比较
第 1 页共 7 页拉格朗日插值法与牛顿插值法的比较一. 背景在工程和科学研究中出现的函数是多种多样的.常常会遇到这样的情况:在某个实际问题中,虽然可以断定所考虑的函数 ) ( x f 在区 ...
MATLAB编写拉格朗日插值与龙格现象
插值与龙格现象在区间[-1,1]上对函数 ,选取不同的插值节点构造插值多项式,比较他们的误差. (1) 取等距节点,n=5,10,15,20. (2) 取节点 ,k=0,1,2,-n:分别取n=5, ...
拉格朗日插值法MATLAB实现（附代码、实例、详解）
第一部分:问题分析 (1)实验题目:拉格朗日插值算法具体实验要求:要求学生运用拉格朗日插值算法通过给定的平面上的n个数据点,计算拉格朗日多项式Pn(x)的值,并将其作为实际函数f(x)的估计值.用m ...
数值计算方法-拉格朗日插值法
要求: 编写拉格朗日插值函数,输入为(x,f,order1) 其中x代表节点数据,为n*2矩阵 f代表是否分段 oreder1:若分段,oreder1代表函数阶次,只能为1或者2 若不分段,则函数阶次 ...
使用拉格朗日多项式(Lagrangian polynomials)的插值法(python,数值积分)
第三十五篇拉格朗日多项式插值插值多项式首先考虑一个函数的推导,该函数精准地通过一系列np离散数据点.虽然有无限多的函数具备这个条件,但我们将专注于最简单的一个,一个n阶多项式,其中n = np− ...
C#实现拉格朗日、牛顿、Hermite插值
拉格朗日.牛顿插值法初始时需要至少两个参考坐标点,在参考点的基础之上构造插值函数y=f(x),然后由插值函数确定需要求解的x坐标的函数值.相较于拉格朗日插值,牛顿插值公式更加复杂,但是在增加一个参考 ...
分段线性插值法实验报告_插值方法_实验报告.doc
插值方法_实验报告肖建计科三班开课学院.实验室: 数统学院实验时间 : 2011 年 5 月 8 日课程名称数学实验实验项目名称插值方法实验项目类型验证演示综合设计其他指导教师李东 ...

使用C++编写拉格朗日3次插值法

前言

GitHub地址

代码

解析

基函数分子,分母的计算

使用C++编写拉格朗日3次插值法相关推荐

最新文章

热门文章