埃尔米特插值及其代码

简介

Lagrange 插值、Newton 插值虽然构造比较简单，但仅要求他们的函数值与被插函数的函数值在节点上相同，但没有考虑到导数值。Hermite 插值 (切触插值) 是要构造一个插值函数，不但在给定的节点上取已知函数值，而且取已知导数值，使插值函数与被插函数的密合程度更好。

关于两点三次 Hermite 插值

即，Hermite插值通式为：

实例

代码

#include<iostream>
using namespace std;
double x[10];
double y[10];
double y_[10];
double d[10],e[10]={1},f[10],g[10]; //倒数 平方基数  


double hermiter(double a,int n)
{
	for(int i=0;i<10;i++) e[i]=1.0,d[i]=0.0;
	double re = 0;
	for(int i=0; i<n; i++) {
        for(int j=0; j<n; j++) {
            if(j!=i) {
                d[i]+=1/(x[i]-x[j]);
                e[i]*=(a-x[j])/(x[i]-x[j]);
            }
        }
    }
    for(int i=0; i<n; i++) {
        g[i]=(1-2*(a-x[i])*d[i])*(e[i]*e[i]);	//y[i]的系数 
        f[i]=(a-x[i])*(e[i]*e[i]);	//y[i]的导数的系数 
        re+=(y[i]*g[i]+ y_[i]*f[i]);
    }
	return re;
}	

int main()
{ 
	int n;
	cout<<"输入共有几个点"<<endl;
	cin>>n;
	cout<<"请先输入x的坐标值，再输入y的值，最后输入y的导数的值"<<endl;
	for(int i=0;i<n;i++)
	{
		cin>>x[i]>>y[i]>>y_[i];
	 }   
	double test;
	while(1)
	{	
		cout<<"输入要测试的数据"<<endl; 
		cin>>test;
		double re = hermiter(test,n);
		cout<<"结果是："<<re<<endl;	
	} 
	return 0;	
}

测试结果：