OpenCV中的加权线性最小二乘法

Question

OpenCV中的加权线性最小二乘法

6

OpenCV的cvSolve可以解决像这样的线性最小二乘问题：

// model: y = a1*x1 + a2*x2 + a3
CvMat *y = cvCreateMat(N, 1, CV_64FC1);
CvMat *X = cvCreateMat(N, 3, CV_64FC1);
CvMat *coeff = cvCreateMat(3, 1, CV_64FC1);

// fill vector y and matrix X
for (int i=0; i<N; ++i)
{
    cvmSet(y, i, 0, my_y_value(i) );
    cvmSet(X, i, 0, my_x1_value(i) );
    cvmSet(X, i, 1, my_x2_value(i) );
    cvmSet(X, i, 2, 1 );
} 

cvSolve(X, y, coeff, CV_SVD);
// now coeff contains a1, a2, a3

但是，我希望对我的数据点应用不同的权重。我该如何应用这些权重？

- Hendrik

1个回答

网页内容由stack overflow 提供, 点击上面的

可以查看英文原文，
原文链接

- Hendrik · Accepted Answer

我发现其实并不难：

for (int i=0; i<N; ++i)
{
    double w = weight(i);

    cvmSet(y, i, 0, w * my_y_value(i) );
    cvmSet(X, i, 0, w * my_x1_value(i) );
    cvmSet(X, i, 1, w * my_x2_value(i) );
    cvmSet(X, i, 2, w );
}

cvSolve(X, y, coeff, CV_SVD);

这个片段只是用权重w将线性方程式的左手边和右手边相乘。对于样本i的错误项实际上被w²乘以。