关于三种插值方法的可视化处理

通过程序求出插值函数的表达式是比较麻烦的，常用的方法是描出插值曲线上尽量密集的有限个采样点，并用这有限个采样点的连线，即折线，近似插值曲线。取点越密集，所得折线就越逼近理论上的插值曲线。本实验中将所取的点的横坐标存放于动态数组Xn中，通过插值方法计算得到的对应纵坐标存放于动态数组Yn中。

以Visual C++.Net 2005为例。

本实验将Lagrange插值、Newton插值和三次样条插值实现为一个C++类CInterpolation，并在Button单击事件中调用该类相应函数，得出插值结果并画出图像。CInterpolation类为

class CInterpolation
{
public:
    CInterpolation();//构造函数
CInterpolation(float *x1, float *y1, int n1);//结点横坐标、纵坐标、下标上限
~ CInterpolation();//析构函数
…………
…………
int n, N;//结点下标上限，采样点下标上限
float *x, *y, *X;//分别存放结点横坐标、结点纵坐标、采样点横坐标
float *p_H,*p_Alpha,*p_Beta,*p_a,*p_b,*p_c,*p_d,*p_m;//样条插值用到的公有指针，分别存放，，，，，，和
};

 其中，有参数的构造函数为

CInterpolation(float *x1, float *y1, int n1) 
{
//动态数组x1，y1中存放结点的横、纵坐标，n1是结点下标上限（即n1+1个结点）
    n=n1;
    N=x1[n]-x1[0];
    X=new float[N+1];
    x=new float[n+1];
    y=new float[n+1];
    for(int i=0;i<=n;i++)
    {
        x[i]=x1[i];
        y[i]=y1[i];
    }
    for(int i=0;i<=N;i++)
        X[i]=x[0]+i;
}

2.1 Lagrange插值

 

float l(float xv,int i)    //求插值基函数的值
{
    float t=1;
    for(int j=0;j<=n;j++)
        if(j!=i)
            t=t*(xv-x[j])/(x[i]-x[j]);
    return t;
}
调用函数l(float x,int i)，可求出
float p_l(float x)         //求在一个点的插值结果
{
    float t=0;
    for(int i=0;i<=n;i++)
        t+=y[i]*l(x,i);
    return t;
}
调用p_l(float x)可实现整个区间的插值
float *Lagrange()         //求整个插值区间上所有采样点的插值结果
{
    float *Y=new float[N+1];
    for(int k=0;k<=N;k++)
        Y[k]=p_l(x[0]+k*h);
    return Y;
}

2.2 Newton插值

float *f()
{
//该函数的返回值是一个长度为n＋1的动态数组，存放各阶差商
}
float w(float x, int i)
{
//该函数计算
}
在求的函数中调用*f()得到各阶差商，然后在循环中调用w(float x)可得出插值结果
float p_n(float x)
{
//该函数计算在一点的值
}
调用p_n(float x)可实现整个区间的插值
float *Newton()
{
//该函数计算出插值区间内所有点的值
}

 2.3三次样条插值

三次样条插值程序可分为以下四步编写：

1. 1. 计算结点间的步长；
   2. 利用产生三对角方程组的系数矩阵和常数向量；
   3. 通过求解三对角方程组，得出中间结点的导数