| 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]); |