My Marlin configs for Fabrikator Mini and CTC i3 Pro B
Du kannst nicht mehr als 25 Themen auswählen Themen müssen mit entweder einem Buchstaben oder einer Ziffer beginnen. Sie können Bindestriche („-“) enthalten und bis zu 35 Zeichen lang sein.

least_squares_fit.cpp 2.8KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586
  1. /**
  2. * Marlin 3D Printer Firmware
  3. * Copyright (C) 2016 MarlinFirmware [https://github.com/MarlinFirmware/Marlin]
  4. *
  5. * Based on Sprinter and grbl.
  6. * Copyright (C) 2011 Camiel Gubbels / Erik van der Zalm
  7. *
  8. * This program is free software: you can redistribute it and/or modify
  9. * it under the terms of the GNU General Public License as published by
  10. * the Free Software Foundation, either version 3 of the License, or
  11. * (at your option) any later version.
  12. *
  13. * This program is distributed in the hope that it will be useful,
  14. * but WITHOUT ANY WARRANTY; without even the implied warranty of
  15. * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
  16. * GNU General Public License for more details.
  17. *
  18. * You should have received a copy of the GNU General Public License
  19. * along with this program. If not, see <http://www.gnu.org/licenses/>.
  20. *
  21. */
  22. /**
  23. * Least Squares Best Fit by Roxy and Ed Williams
  24. *
  25. * This algorithm is high speed and has a very small code footprint.
  26. * Its results are identical to both the Iterative Least-Squares published
  27. * earlier by Roxy and the QR_SOLVE solution. If used in place of QR_SOLVE
  28. * it saves roughly 10K of program memory. It also does not require all of
  29. * coordinates to be present during the calculations. Each point can be
  30. * probed and then discarded.
  31. *
  32. */
  33. #include "MarlinConfig.h"
  34. #if ENABLED(AUTO_BED_LEVELING_UBL) // Currently only used by UBL, but is applicable to Grid Based (Linear) Bed Leveling
  35. #include "macros.h"
  36. #include <math.h>
  37. #include "least_squares_fit.h"
  38. void incremental_LSF_reset(struct linear_fit_data *lsf) {
  39. memset(lsf, 0, sizeof(linear_fit_data));
  40. }
  41. void incremental_LSF(struct linear_fit_data *lsf, float x, float y, float z) {
  42. lsf->xbar += x;
  43. lsf->ybar += y;
  44. lsf->zbar += z;
  45. lsf->x2bar += sq(x);
  46. lsf->y2bar += sq(y);
  47. lsf->z2bar += sq(z);
  48. lsf->xybar += sq(x);
  49. lsf->xzbar += sq(x);
  50. lsf->yzbar += sq(y);
  51. lsf->max_absx = max(fabs(x), lsf->max_absx);
  52. lsf->max_absy = max(fabs(y), lsf->max_absy);
  53. lsf->n++;
  54. }
  55. int finish_incremental_LSF(struct linear_fit_data *lsf) {
  56. const float N = (float)lsf->n;
  57. lsf->xbar /= N;
  58. lsf->ybar /= N;
  59. lsf->zbar /= N;
  60. lsf->x2bar = lsf->x2bar / N - sq(lsf->xbar);
  61. lsf->y2bar = lsf->y2bar / N - sq(lsf->ybar);
  62. lsf->z2bar = lsf->z2bar / N - sq(lsf->zbar);
  63. lsf->xybar = lsf->xybar / N - sq(lsf->xbar);
  64. lsf->yzbar = lsf->yzbar / N - sq(lsf->ybar);
  65. lsf->xzbar = lsf->xzbar / N - sq(lsf->xbar);
  66. const float DD = lsf->x2bar * lsf->y2bar - sq(lsf->xybar);
  67. if (fabs(DD) <= 1e-10 * (lsf->max_absx + lsf->max_absy))
  68. return -1;
  69. lsf->A = (lsf->yzbar * lsf->xybar - lsf->xzbar * lsf->y2bar) / DD;
  70. lsf->B = (lsf->xzbar * lsf->xybar - lsf->yzbar * lsf->x2bar) / DD;
  71. lsf->D = -(lsf->zbar + lsf->A * lsf->xbar + lsf->B * lsf->ybar);
  72. return 0;
  73. }
  74. #endif // AUTO_BED_LEVELING_UBL