text1.c

拉格拉日算法

float Lagrange(float px[], float py[], int n, float x)
{//px,py:插值点(Xi,Yi) n:插值点个数  x:待计算的函数点
    float y = 0;
    for(int k=0; k<n; k++)//k控制Lagrange基函数序列
    {
        float tmp = 1;//tmp表示Lagrange基函数
        for(int i=0;i<n;i++)
        {
            if(i!=k)
            {
                tmp = tmp * (x-px[i])/(px[k]-px[i]);  //key step
            }
        }
        y = y + py[k]*tmp;
    }
    return y;
}

int main()
{
    float x;//插值
    float px[10];//已知(x0,y0),(x1,y1)
    float py[10];
    int n;//输入已知插值组数

    cout<<"输入插值组数:"<<endl;
    cin>>n;
    cout<<"输入"<<n<<"组已知插值数（X,Y）"<<endl;
    for(int i=0; i<n; i++)
    {
        cin>>px[i]>>py[i];
    }
    cout<<"输入插值:"<<endl;
    cin>>x;

    cout<<"Lagrange插值结果:"<<Lagrange(px, py, n, x)<<endl;

    system("pause");
    return 0;
}