void myLMS_poly4(const std::vector<double> src_x, const std::vector<double> src_y, int size, std::vector<double>& dst_y)
{
double a, b, c, d, e;
//Mat A = Mat_<double>(5, 5);
//Mat B = Mat_<double>(5, 1);
//Mat C = Mat_<double>(5, 1);
Mat A(5, 5, CV_64F);
Mat B(5, 1, CV_64F);
Mat X(5, 1, CV_64F);
double sum = 0;
for (size_t i = 0; i < size; ++i)
{
sum += pow(src_x[i], 8);
}
A.at<double>(0, 0) = sum;
sum = 0;
for (size_t i = 0; i < size; ++i)
{
sum += pow(src_x[i], 7);
}
A.at<double>(0, 1) = A.at<double>(1, 0) = sum;
sum