C语言平滑性算法,优化了的过关键点的光滑曲线拟合算法

这个是我一个数学老师(教授，数学高手，经常自己做算法)给我的例子，用于多个离散点拟合光滑曲线的，他优化了追赶法，这个例子适用于闭合和不闭合两种情况。当时由于工程情况，写的急，代码不好看，但是很好用。为了方便传递参数，我做了一个链表，用时候根据自己情况可以修改，核心算法不动即可。

class CFoldPoint

{public:

double X;    double Y;

};

typedef CTypedPtrList CFoldPointList;

typedef CArray CDoubleArray;

三个函数，SPLine 调用另外两个。用时候直接调用SPLine函数，入口pList是已知离散点链表，pDestList是生成的点的链表。SM是在两个点中间插入点的数目，continue=0是采样点无规律，要求生成闭合曲线。1是采样点x坐标连续 2是y连续

void ZG(CDoubleArray *A,CDoubleArray *B,CDoubleArray *C,CDoubleArray *G,int &LOGI)

{

//追赶法

register long I;

int N;

N=A->GetSize();

if(LOGI==0)

{

(*C)[0]=(*C)[0]/(*B)[0];

for(I=1;I   {

(*B)[I]=(*B)[I]-(*A)[I]*(*C)[I-1];

(*C)[I]=(*C)[I]/(*B)[I];

}

(*A)[0]=0.;

(*C)[N-1]=0.;

LOGI=1;

}

(*G)[0]=(*G)[0]/(*B)[0];

for(I=1;I  {

(*G)[I]=((*G)[I]-(*A)[I]*(*G)[I-1])/(*B)[I];

}

for(I=N-2;I>-1;I–)//DO 30 I=N-1,1,-1

{