【计算方法实验】实验5:拉格朗日插值

本文介绍了一个使用C++实现的拉格朗日插值法程序,该程序定义了`Point`类来表示二维坐标点,并通过`Lagrange`类实现了拉格朗日插值法的基本功能。`Lagrange`类包括设置插值点、计算基函数值以及计算给定x值下的插值结果等方法。

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

#include<iostream>

using namespace std;

class Point{
public:
    double x,y;
    Point(double x, double y):x(x),y(y){}
};

class Lagrange{
public:

    double getPrimaryFunctionValue(double x, int k);
    double getValue(double x);
    void setPoints(Point* ps, int n);

private:

    Point* ps;
    int n;

};

void Lagrange::setPoints(Point* ps, int n){
    this->ps = ps;
    this->n = n;
}

double Lagrange::getPrimaryFunctionValue(double x, int k){

    double son = 1, mother = 1;

    for(int i = 0; i < n; i++){
        if(k != i){
            son *= (x - ps[i].x);
            mother *= (ps[k].x - ps[i].x);
        }
    }

    return son/mother;
}

double Lagrange::getValue(double x){

    double result = 0;

    for(int i = 0; i < n; i++){
        result += getPrimaryFunctionValue(x, i) * ps[i].y;
    }

    return result;
}

int main(){

    Lagrange lg;

    Point p[] = {Point(0.4,0.41075),Point(0.55,0.57815),Point(0.8,0.88811),Point(0.9,1.02652),Point(1,1.17520)};

    lg.setPoints(p, 5);

    cout<<lg.getValue(0.5)<<endl;
    cout<<lg.getValue(0.7)<<endl;
    cout<<lg.getValue(0.85)<<endl;

    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值