拉格朗日（Lagrange）插值法（C语言）

最新推荐文章于 2022-05-18 22:01:57 发布

原创最新推荐文章于 2022-05-18 22:01:57 发布 · 1.8w 阅读

113 ·

CC 4.0 BY-SA版权

文章标签：

#计算方法 #拉格朗日插值

计算方法专栏收录该内容

16 篇文章

订阅专栏

本文提供了一个使用C++实现的拉格朗日插值法示例，包括二次插值和五次插值的计算过程。通过具体的数据点展示了如何进行插值计算，并给出了插值结果。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题目
这里写图片描述
拉格朗日插值法流程图

(1)二次插值：

#include<stdio.h>
float lagelangri(float x[],float y[],float xx,int n)
{
    int i,j;
    float *a,yy=0;
    a=new float[n];
    for(i=0;i<=n-1;i++)
    {
        a[i]=y[i];
        for(j=0;j<=n-1;j++)
            if(j!=i)a[i]*=(xx-x[j])/(x[i]-x[j]);

        yy+=a[i];
    }
    delete a;
    return yy;
}

void main()
{
    float x[5]={-3.0,-1.0,1.0,2.0,3.0};
    float y[5]={1.0,1.5,2.0,2.0,1.0};
    float xx1=-2,xx2=0,xx3=2.75,yy1,yy2,yy3;
    yy1=lagelangri(x,y,xx1,3);
    yy2=lagelangri(x,y,xx2,3);
    yy3=lagelangri(x,y,xx3,3);
    printf("x1=%-20f,y1=%f\n",xx1,yy1);
    printf("x2=%-20f,y2=%f\n",xx2,yy2);
    printf("x3=%-20f,y3=%f\n",xx3,yy3);
}

(2)五次插值：

#include<stdio.h>
float lagelangri(float x[],float y[],float xx,int n)
{
    int i,j;
    float *a,yy=0;
    a=new float[n];
    for(i=0;i<=n-1;i++)
    {
        a[i]=y[i];
        for(j=0;j<=n-1;j++)
            if(j!=i)a[i]*=(xx-x[j])/(x[i]-x[j]);

        yy+=a[i];
    }
    delete a;
    return yy;
}

void main()
{
    float x[6]={0.30,0.42,0.50,0.58,0.66,0.72};
    float y[6]={1.04403,1.08462,1.11803,1.15603,1.19817,1.23223};
    float xx1=0.46,xx2=0.55,xx3=0.60,yy1,yy2,yy3;
    yy1=lagelangri(x,y,xx1,6);
    yy2=lagelangri(x,y,xx2,6);
    yy3=lagelangri(x,y,xx3,6);
    printf("x1=%-20f,y1=%f\n",xx1,yy1);
    printf("x2=%-20f,y2=%f\n",xx2,yy2);
    printf("x3=%-20f,y3=%f\n",xx3,yy3);
}