本文所要你擬合的曲線公式：

未知引數：Pi

所求引數數目m=4

已知的點數目n=7

初始的引數設定 p[m] = { 1,1,1,1 }

X[i]設定為y： x[n] = {0,1,2,2,3,4,2};

Data是已知點的x：data[7] = {0,1,2,3,3,4,10};

其他保留為預設。

目標函式：

void ros(double* p, double* x, int m, int
 
 n, void* data){
    LVTYPE* data1 = (LVTYPE*)data;
    register int i;
    for (i = 0; i < n; i++){
        x[i] =(p[0]+p[1] * data1[i]+p[2]*data1[i]*data1[i]+p[3]*data1[i]*data1[i]*data1[i]);
    }
}

雅可比函式：

void jacros(double* p, double *jac, int m, int n, void
 
* data){
    LVTYPE* data1 = (LVTYPE*)data;
    register int i, j;
    for (i = 0, j = 0; i < n; i++){
        jac[j++] = 1;
        jac[j++] = data1[i];
        jac[j++] = data1[i]*2;
        jac[j++] = data1[i] * data1[i] * data1[i] * 3;
    }
}

呼叫：

int ret = dlevmar_der(ros, jacros, p, x, m, n, 10000
 
, NULL, info, NULL, NULL, (void*)data);

結果：

轉載自：http://www.cnblogs.com/jiafenggang/p/5008044.html