222,17 → 222,22 |
double A, B; |
double esty; |
|
// printf("PointEstimate points: "); |
for(int i = 0; i < nrp; i++) |
{ |
// printf("[%lf,%lf], ", points[2*i], points[2*i+1]); |
accx += points[2*i]; |
accy += points[2*i+1]; |
accxy += points[2*i]*points[2*i+1]; |
accx2 += points[2*i]*points[2*i]; |
} |
// printf("\n"); |
|
A = (accx2*accy - accx*accxy)/(nrp*accx2 - accx*accx); |
B = (nrp*accxy - accx*accy)/(nrp*accx2 - accx*accx); |
|
// printf("Fiting function = %lf + %lf*x\n", A, B); |
|
esty = A + B*points[2*nrp]; |
|
if(DBGSIG) printf("PointEstimate(): A = %lf, B = %lf, estimate = %lf, real = %lf, error = %lf\n", A, B, esty, points[2*nrp+1], abs(esty - points[2*nrp+1])/points[2*nrp+1]); |