拉格朗日插值 C++ 和 Matlab实现

拉格朗日插值 C++ 和 Matlab实现

拉格朗日插值法

概述

利用 n 个变量 $x_i$ 与函数值 $y_i$ 构造插值函数Ln的方法
构造关键
构造多项式 $l_k$
$x_k-1$
另 $l_k(x)=\frac{(x-x_{k-1})(x-x_{k+1})}{(x_k-x_{k-1})(x_k-x_{k+1})}$
$l_k$ 满足 $ l_k(x_k) = 1 , l_k(x_i) = 0 , i = 1,…,n ,i !=k$
解出 $$l_n(x)=y_i\times \sum _{k=0}^n{l}_x(x_i),i=1,...,n$$
注 ： 上述公式是我用mathjax写的，如有错误请联系我修正
敬请指正

详解

点此查看百科 拉格朗日插值法

代码

C++实现

结果和输入格式在注释中

/*
/
                              结果 
/
输入需要插值的数目 : 6
Xi : 0.4 0.55 0.65 0.8 0.9 1.05
Yi : 0.41075 0.57815 0.69675 0.88811 1.02652 1.25382
输入需要插值的x :
0.596
y : 0.631917
/
tips : 运行环境 DevC++ 5.9.2
/
*/ 
#include<iostream> 
#include<iomanip>
#include<cstring>

#define maxn 100

int n;
double x[maxn];
double y[maxn];
double ln[maxn];

using namespace std;

int main()
{
	double xi,res;
	res = 0;
	memset(x,0,sizeof(x));
	memset(y,0,sizeof(y));
	cout << "输入需要插值的数目 : " ; 
	cin >> n;
	if(n > maxn)
	{
		cout << "超出上限" <<endl; 
		cout << "请修改 maxn 再运行" <<endl;
		return 0; 
	}
	cout << "Xi : ";
	for(int i = 0;i < n;i++)
		cin >> x[i];
	cout << "Yi : ";
	for(int i = 0;i < n;i++)
		cin >> y[i];
	cout << "输入需要插值的x : " << endl;
	cin >> xi;
	for(int i = 0;i < n;i++)
	{
		ln[i] = 1;
		for(int j = 0;j < n;j++)
		{
			if(i == j)
				continue;
			ln[i] *= ((xi - x[j]) / (x[i] - x[j]));
		}
		res += (y[i] * ln[i]);
	}
	cout << "y : " << res << endl;
}

Matlab 实现

将该函数存为m文件

function res= Lagrange(x,y,xi)
    xl = length(x);
    yl = length(y);
    if(xl ~= yl)
        disp('向量长度不等')
        return
    end
    res = 0;
    for i = 1:xl
        ln = 1;
        for j = 1:yl
            if i == j
                continue;
            end
            ln = ln * ((xi - x(j))/(x(i) - x(j)));
        end
        res = res + (y(i) * ln);
    end
end

调用下面语句测试函数

x = [0.4,0.55,0.65,0.8,0.9,1.05]
y = [0.41075,0.57815,0.69675,0.88811,1.02652,1.25382]
xi = 0.596
Lagrange(x,y,xi)