首页 > 代码库 > 样条之埃特金(Aitken)逐步插值函数

样条之埃特金(Aitken)逐步插值函数

 核心代码:

//////////////////////////////////////////////////////////////////////// 埃特金逐步插值//////////////////////////////////////////////////////////////////////static float GetValueAitken(const void* valuesPtr, int stride, int n, float t, float eps){     int i,j,k,m,l;    float z,xx[10],yy[10];    // 初值    z = 0.0f;    // 特例处理    if (n < 1)     {        return(z);    }    if (n == 1)     {         z = YfGetFloatValue(valuesPtr, stride, 0);          return(z);    }    float xStep = 1.0f/(n - 1);    // 开始插值    m=10;    if (m > n)     {        m = n;    }    if (t <= 0.0f)     {        k = 1;    }    else if (t >= (n-1)*xStep)     {        k = n;    }    else    {         k = 1;         j = n;        while ((k-j != 1) && (k-j != -1))        {             l = (k+j)/2;            if (t < (l-1)*xStep)                 j = l;            else                 k = l;        }        if (fabs(t-((l-1)*xStep)) > fabs(t-(j-1)*xStep))         {            k = j;        }    }    j = 1;     l = 0;    for (i = 1; i <= m; i++)    {         k = k+j*l;        if ((k<1) || (k>n))        {             l = l+1;             j = -j;             k = k+j*l;        }        xx[i-1] = (k-1)*xStep;         yy[i-1] = YfGetFloatValue(valuesPtr, stride, k - 1);        l = l+1;         j = -j;    }    i = 0;    do    {         i = i+1;         z = yy[i];        for (j = 0; j <= i-1; j++)        {            z = yy[j]+(t-xx[j])*(yy[j]-z)/(xx[j]-xx[i]);        }        yy[i] = z;    }     while ((i != m-1) && (fabs(yy[i]-yy[i-1]) > eps));    return(z);}

切图:

 

 

相关软件的下载地址为:http://files.cnblogs.com/WhyEngine/TestSpline.zip

样条之埃特金(Aitken)逐步插值函数