2 * This file is part of Cleanflight.
4 * Cleanflight is free software: you can redistribute it and/or modify
5 * it under the terms of the GNU General Public License as published by
6 * the Free Software Foundation, either version 3 of the License, or
7 * (at your option) any later version.
9 * Cleanflight is distributed in the hope that it will be useful,
10 * but WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 * GNU General Public License for more details.
14 * You should have received a copy of the GNU General Public License
15 * along with Cleanflight. If not, see <http://www.gnu.org/licenses/>.
25 #include "build/debug.h"
27 #include "common/axis.h"
28 #include "common/filter.h"
29 #include "common/maths.h"
30 #include "common/utils.h"
31 #include "common/calibration.h"
33 #include "config/config_reset.h"
34 #include "config/parameter_group.h"
35 #include "config/parameter_group_ids.h"
37 #include "drivers/accgyro/accgyro.h"
38 #include "drivers/accgyro/accgyro_mpu.h"
39 #include "drivers/accgyro/accgyro_mpu6000.h"
40 #include "drivers/accgyro/accgyro_mpu6500.h"
41 #include "drivers/accgyro/accgyro_mpu9250.h"
43 #include "drivers/accgyro/accgyro_bmi088.h"
44 #include "drivers/accgyro/accgyro_bmi160.h"
45 #include "drivers/accgyro/accgyro_bmi270.h"
46 #include "drivers/accgyro/accgyro_icm20689.h"
47 #include "drivers/accgyro/accgyro_icm42605.h"
48 #include "drivers/accgyro/accgyro_lsm6dxx.h"
49 #include "drivers/accgyro/accgyro_fake.h"
50 #include "drivers/sensor.h"
52 #include "fc/config.h"
53 #include "fc/runtime_config.h"
54 #include "fc/settings.h"
56 #include "io/beeper.h"
58 #include "sensors/acceleration.h"
59 #include "sensors/battery.h"
60 #include "sensors/boardalignment.h"
61 #include "sensors/gyro.h"
62 #include "sensors/sensors.h"
64 #ifdef USE_HARDWARE_REVISION_DETECTION
65 #include "hardware_revision.h"
69 FASTRAM acc_t acc
; // acc access functions
71 STATIC_FASTRAM zeroCalibrationVector_t zeroCalibration
;
73 STATIC_FASTRAM
float accADC
[XYZ_AXIS_COUNT
];
75 STATIC_FASTRAM filter_t accFilter
[XYZ_AXIS_COUNT
];
76 STATIC_FASTRAM filterApplyFnPtr accSoftLpfFilterApplyFn
;
77 STATIC_FASTRAM
void *accSoftLpfFilter
[XYZ_AXIS_COUNT
];
79 static EXTENDED_FASTRAM pt1Filter_t accVibeFloorFilter
[XYZ_AXIS_COUNT
];
80 static EXTENDED_FASTRAM pt1Filter_t accVibeFilter
[XYZ_AXIS_COUNT
];
82 static EXTENDED_FASTRAM filterApplyFnPtr accNotchFilterApplyFn
;
83 static EXTENDED_FASTRAM
void *accNotchFilter
[XYZ_AXIS_COUNT
];
85 static EXTENDED_FASTRAM
float fAccZero
[XYZ_AXIS_COUNT
];
86 static EXTENDED_FASTRAM
float fAccGain
[XYZ_AXIS_COUNT
];
88 PG_REGISTER_WITH_RESET_FN(accelerometerConfig_t
, accelerometerConfig
, PG_ACCELEROMETER_CONFIG
, 5);
90 void pgResetFn_accelerometerConfig(accelerometerConfig_t
*instance
)
92 RESET_CONFIG_2(accelerometerConfig_t
, instance
,
93 .acc_hardware
= SETTING_ACC_HARDWARE_DEFAULT
,
94 .acc_lpf_hz
= SETTING_ACC_LPF_HZ_DEFAULT
,
95 .acc_notch_hz
= SETTING_ACC_NOTCH_HZ_DEFAULT
,
96 .acc_notch_cutoff
= SETTING_ACC_NOTCH_CUTOFF_DEFAULT
,
97 .acc_soft_lpf_type
= SETTING_ACC_LPF_TYPE_DEFAULT
99 RESET_CONFIG_2(flightDynamicsTrims_t
, &instance
->accZero
,
100 .raw
[X
] = SETTING_ACCZERO_X_DEFAULT
,
101 .raw
[Y
] = SETTING_ACCZERO_Y_DEFAULT
,
102 .raw
[Z
] = SETTING_ACCZERO_Z_DEFAULT
104 RESET_CONFIG_2(flightDynamicsTrims_t
, &instance
->accGain
,
105 .raw
[X
] = SETTING_ACCGAIN_X_DEFAULT
,
106 .raw
[Y
] = SETTING_ACCGAIN_Y_DEFAULT
,
107 .raw
[Z
] = SETTING_ACCGAIN_Z_DEFAULT
111 static void updateAccCoefficients(void) {
113 for (uint8_t i
= 0; i
< XYZ_AXIS_COUNT
; i
++) {
115 fAccZero
[i
] = (float)accelerometerConfig()->accZero
.raw
[i
];
117 fAccGain
[i
] = (float)accelerometerConfig()->accGain
.raw
[i
] / 4096.0f
;
122 static bool accDetect(accDev_t
*dev
, accelerationSensor_e accHardwareToUse
)
124 accelerationSensor_e accHardware
= ACC_NONE
;
126 dev
->accAlign
= ALIGN_DEFAULT
;
128 requestedSensors
[SENSOR_INDEX_ACC
] = accHardwareToUse
;
130 switch (accHardwareToUse
) {
134 #ifdef USE_IMU_MPU6000
136 if (mpu6000AccDetect(dev
)) {
137 accHardware
= ACC_MPU6000
;
140 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
141 if (accHardwareToUse
!= ACC_AUTODETECT
) {
147 #if defined(USE_IMU_MPU6500)
149 if (mpu6500AccDetect(dev
)) {
150 accHardware
= ACC_MPU6500
;
153 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
154 if (accHardwareToUse
!= ACC_AUTODETECT
) {
160 #if defined(USE_IMU_MPU9250)
162 if (mpu9250AccDetect(dev
)) {
163 accHardware
= ACC_MPU9250
;
166 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
167 if (accHardwareToUse
!= ACC_AUTODETECT
) {
173 #if defined(USE_IMU_BMI160)
175 if (bmi160AccDetect(dev
)) {
176 accHardware
= ACC_BMI160
;
179 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
180 if (accHardwareToUse
!= ACC_AUTODETECT
) {
186 #if defined(USE_IMU_BMI088)
188 if (bmi088AccDetect(dev
)) {
189 accHardware
= ACC_BMI088
;
192 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
193 if (accHardwareToUse
!= ACC_AUTODETECT
) {
199 #ifdef USE_IMU_ICM20689
201 if (icm20689AccDetect(dev
)) {
202 accHardware
= ACC_ICM20689
;
205 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
206 if (accHardwareToUse
!= ACC_AUTODETECT
) {
212 #ifdef USE_IMU_ICM42605
214 if (icm42605AccDetect(dev
)) {
215 accHardware
= ACC_ICM42605
;
218 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
219 if (accHardwareToUse
!= ACC_AUTODETECT
) {
225 #ifdef USE_IMU_BMI270
227 if (bmi270AccDetect(dev
)) {
228 accHardware
= ACC_BMI270
;
231 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
232 if (accHardwareToUse
!= ACC_AUTODETECT
) {
237 #ifdef USE_IMU_LSM6DXX
239 if (lsm6dAccDetect(dev
)) {
240 accHardware
= ACC_LSM6DXX
;
243 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
244 if (accHardwareToUse
!= ACC_AUTODETECT
) {
251 if (fakeAccDetect(dev
)) {
252 accHardware
= ACC_FAKE
;
255 /* If we are asked for a specific sensor - break out, otherwise - fall through and continue */
256 if (accHardwareToUse
!= ACC_AUTODETECT
) {
263 case ACC_NONE
: // disable ACC
264 accHardware
= ACC_NONE
;
268 if (accHardware
== ACC_NONE
) {
272 detectedSensors
[SENSOR_INDEX_ACC
] = accHardware
;
273 sensorsSet(SENSOR_ACC
);
277 bool accInit(uint32_t targetLooptime
)
279 memset(&acc
, 0, sizeof(acc
));
281 // Set inertial sensor tag (for dual-gyro selection)
283 acc
.dev
.imuSensorToUse
= gyroConfig()->gyro_to_use
; // Use the same selection from gyroConfig()
285 acc
.dev
.imuSensorToUse
= 0;
288 if (!accDetect(&acc
.dev
, accelerometerConfig()->acc_hardware
)) {
292 acc
.dev
.acc_1G
= 256; // set default
293 acc
.dev
.initFn(&acc
.dev
);
294 acc
.accTargetLooptime
= targetLooptime
;
295 acc
.accClipCount
= 0;
297 updateAccCoefficients();
299 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
300 acc
.extremes
[axis
].min
= 100;
301 acc
.extremes
[axis
].max
= -100;
307 static bool calibratedPosition
[6];
308 static float accSamples
[6][3];
310 uint8_t accGetCalibrationAxisFlags(void)
312 if (STATE(ACCELEROMETER_CALIBRATED
)) {
313 return 0x3F; // All 6 bits are set
316 static const uint8_t bitMap
[6] = { 0, 1, 3, 5, 2, 4 }; // A mapping of bits to match position indexes in Configurator
318 for (int i
= 0; i
< 6; i
++) {
319 if (calibratedPosition
[i
]) {
320 flags
|= (1 << bitMap
[i
]);
327 static int getPrimaryAxisIndex(float accADCData
[3])
329 // Work on a copy so we don't mess with accADC data
332 applySensorAlignment(sample
, accADCData
, acc
.dev
.accAlign
);
334 // Tolerate up to atan(1 / 1.5) = 33 deg tilt (in worst case 66 deg separation between points)
335 if ((ABS(sample
[Z
]) / 1.5f
) > ABS(sample
[X
]) && (ABS(sample
[Z
]) / 1.5f
) > ABS(sample
[Y
])) {
337 return (sample
[Z
] > 0) ? 0 : 1;
339 else if ((ABS(sample
[X
]) / 1.5f
) > ABS(sample
[Y
]) && (ABS(sample
[X
]) / 1.5f
) > ABS(sample
[Z
])) {
341 return (sample
[X
] > 0) ? 2 : 3;
343 else if ((ABS(sample
[Y
]) / 1.5f
) > ABS(sample
[X
]) && (ABS(sample
[Y
]) / 1.5f
) > ABS(sample
[Z
])) {
345 return (sample
[Y
] > 0) ? 4 : 5;
351 void accStartCalibration(void)
353 int positionIndex
= getPrimaryAxisIndex(accADC
);
355 // Fail if we can't detect the side
356 if (positionIndex
< 0) {
360 // Fail if we have accelerometer fully calibrated and are NOT starting with TOP-UP position
361 if (STATE(ACCELEROMETER_CALIBRATED
) && positionIndex
!= 0) {
365 // Top+up and first calibration cycle, reset everything
366 if (positionIndex
== 0) {
367 for (int axis
= 0; axis
< 6; axis
++) {
368 calibratedPosition
[axis
] = false;
369 accSamples
[axis
][X
] = 0;
370 accSamples
[axis
][Y
] = 0;
371 accSamples
[axis
][Z
] = 0;
374 DISABLE_STATE(ACCELEROMETER_CALIBRATED
);
377 // Tolerate 5% variance in accelerometer readings
378 zeroCalibrationStartV(&zeroCalibration
, CALIBRATING_ACC_TIME_MS
, acc
.dev
.acc_1G
* 0.05f
, true);
381 static bool allOrientationsHaveCalibrationDataCollected(void)
383 // Return true only if we have calibration data for all 6 positions
384 return calibratedPosition
[0] && calibratedPosition
[1] && calibratedPosition
[2] &&
385 calibratedPosition
[3] && calibratedPosition
[4] && calibratedPosition
[5];
388 bool accIsCalibrationComplete(void)
390 return zeroCalibrationIsCompleteV(&zeroCalibration
);
393 static void performAcclerationCalibration(void)
395 // Shortcut - no need to do any math if acceleromter is marked as calibrated
396 if (STATE(ACCELEROMETER_CALIBRATED
)) {
400 // If zero calibration logic is finished - no need to do anything
401 if (accIsCalibrationComplete()) {
406 int positionIndex
= getPrimaryAxisIndex(accADC
);
408 // Check if sample is usable
409 if (positionIndex
< 0) {
413 if (!calibratedPosition
[positionIndex
]) {
418 zeroCalibrationAddValueV(&zeroCalibration
, &v
);
420 if (zeroCalibrationIsCompleteV(&zeroCalibration
)) {
421 if (zeroCalibrationIsSuccessfulV(&zeroCalibration
)) {
422 zeroCalibrationGetZeroV(&zeroCalibration
, &v
);
424 accSamples
[positionIndex
][X
] = v
.v
[X
];
425 accSamples
[positionIndex
][Y
] = v
.v
[Y
];
426 accSamples
[positionIndex
][Z
] = v
.v
[Z
];
428 calibratedPosition
[positionIndex
] = true;
431 calibratedPosition
[positionIndex
] = false;
434 beeperConfirmationBeeps(2);
438 if (allOrientationsHaveCalibrationDataCollected()) {
439 sensorCalibrationState_t calState
;
441 bool calFailed
= false;
443 /* Calculate offset */
444 sensorCalibrationResetState(&calState
);
446 for (int axis
= 0; axis
< 6; axis
++) {
447 sensorCalibrationPushSampleForOffsetCalculation(&calState
, accSamples
[axis
]);
450 if (!sensorCalibrationSolveForOffset(&calState
, accTmp
)) {
457 accelerometerConfigMutable()->accZero
.raw
[X
] = lrintf(accTmp
[X
]);
458 accelerometerConfigMutable()->accZero
.raw
[Y
] = lrintf(accTmp
[Y
]);
459 accelerometerConfigMutable()->accZero
.raw
[Z
] = lrintf(accTmp
[Z
]);
461 /* Not we can offset our accumulated averages samples and calculate scale factors and calculate gains */
462 sensorCalibrationResetState(&calState
);
464 for (int axis
= 0; axis
< 6; axis
++) {
467 accSample
[X
] = accSamples
[axis
][X
] - accelerometerConfig()->accZero
.raw
[X
];
468 accSample
[Y
] = accSamples
[axis
][Y
] - accelerometerConfig()->accZero
.raw
[Y
];
469 accSample
[Z
] = accSamples
[axis
][Z
] - accelerometerConfig()->accZero
.raw
[Z
];
471 sensorCalibrationPushSampleForScaleCalculation(&calState
, axis
/ 2, accSample
, acc
.dev
.acc_1G
);
474 if (!sensorCalibrationSolveForScale(&calState
, accTmp
)) {
481 for (int axis
= 0; axis
< 3; axis
++) {
482 accelerometerConfigMutable()->accGain
.raw
[axis
] = lrintf(accTmp
[axis
] * 4096);
486 // If failed - don't save and also invalidate the calibration data for all positions
487 for (int axis
= 0; axis
< 6; axis
++) {
488 calibratedPosition
[axis
] = false;
492 // saveConfigAndNotify will trigger eepromREAD and in turn call back the accelerometer gain validation
493 // that will set ENABLE_STATE(ACCELEROMETER_CALIBRATED) if all is good
494 saveConfigAndNotify();
495 //Recompute all coeffs
496 updateAccCoefficients();
501 static void applyAccelerationZero(void)
503 float tmp
[XYZ_AXIS_COUNT
];
506 arm_sub_f32(accADC
, fAccZero
, tmp
, XYZ_AXIS_COUNT
);
508 arm_mult_f32(tmp
, fAccGain
, accADC
, XYZ_AXIS_COUNT
);
512 * Calculate measured acceleration in body frame in m/s^2
514 void accGetMeasuredAcceleration(fpVector3_t
*measuredAcc
)
516 arm_scale_f32(acc
.accADCf
, GRAVITY_CMSS
, measuredAcc
->v
, XYZ_AXIS_COUNT
);
522 const acc_extremes_t
* accGetMeasuredExtremes(void)
524 return (const acc_extremes_t
*)&acc
.extremes
;
527 float accGetMeasuredMaxG(void)
532 void resetGForceStats(void) {
535 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
536 acc
.extremes
[axis
].min
= 100;
537 acc
.extremes
[axis
].max
= -100;
544 if (ARMING_FLAG(SIMULATOR_MODE_HITL
)) {
545 //output: acc.accADCf
546 //unused: acc.dev.ADCRaw[], acc.accClipCount, acc.accVibeSq[]
550 if (!acc
.dev
.readFn(&acc
.dev
)) {
553 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
554 accADC
[axis
] = acc
.dev
.ADCRaw
[axis
];
555 DEBUG_SET(DEBUG_ACC
, axis
, accADC
[axis
]);
558 if (!ARMING_FLAG(SIMULATOR_MODE_SITL
)) {
559 performAcclerationCalibration();
560 applyAccelerationZero();
563 applySensorAlignment(accADC
, accADC
, acc
.dev
.accAlign
);
564 applyBoardAlignment(accADC
);
566 // Calculate acceleration readings in G's
567 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
568 acc
.accADCf
[axis
] = (float)accADC
[axis
] / acc
.dev
.acc_1G
;
571 // Before filtering check for clipping and vibration levels
572 if (fabsf(acc
.accADCf
[X
]) > ACC_CLIPPING_THRESHOLD_G
|| fabsf(acc
.accADCf
[Y
]) > ACC_CLIPPING_THRESHOLD_G
|| fabsf(acc
.accADCf
[Z
]) > ACC_CLIPPING_THRESHOLD_G
) {
573 acc
.isClipped
= true;
577 acc
.isClipped
= false;
580 // Calculate vibration levels
581 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
582 // filter accel at 5hz
583 const float accFloorFilt
= pt1FilterApply(&accVibeFloorFilter
[axis
], acc
.accADCf
[axis
]);
585 // calc difference from this sample and 5hz filtered value, square and filter at 2hz
586 const float accDiff
= acc
.accADCf
[axis
] - accFloorFilt
;
587 acc
.accVibeSq
[axis
] = pt1FilterApply(&accVibeFilter
[axis
], accDiff
* accDiff
);
588 acc
.accVibe
= fast_fsqrtf(acc
.accVibeSq
[X
] + acc
.accVibeSq
[Y
] + acc
.accVibeSq
[Z
]);
591 // Filter acceleration
592 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
593 acc
.accADCf
[axis
] = accSoftLpfFilterApplyFn(accSoftLpfFilter
[axis
], acc
.accADCf
[axis
]);
596 if (accelerometerConfig()->acc_notch_hz
) {
597 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
598 acc
.accADCf
[axis
] = accNotchFilterApplyFn(accNotchFilter
[axis
], acc
.accADCf
[axis
]);
604 // Record extremes: min/max for each axis and acceleration vector modulus
605 void updateAccExtremes(void)
607 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
608 if (acc
.accADCf
[axis
] < acc
.extremes
[axis
].min
) acc
.extremes
[axis
].min
= acc
.accADCf
[axis
];
609 if (acc
.accADCf
[axis
] > acc
.extremes
[axis
].max
) acc
.extremes
[axis
].max
= acc
.accADCf
[axis
];
612 float gforce
= calc_length_pythagorean_3D(acc
.accADCf
[X
], acc
.accADCf
[Y
], acc
.accADCf
[Z
]);
613 if (gforce
> acc
.maxG
) acc
.maxG
= gforce
;
616 void accGetVibrationLevels(fpVector3_t
*accVibeLevels
)
618 accVibeLevels
->x
= fast_fsqrtf(acc
.accVibeSq
[X
]);
619 accVibeLevels
->y
= fast_fsqrtf(acc
.accVibeSq
[Y
]);
620 accVibeLevels
->z
= fast_fsqrtf(acc
.accVibeSq
[Z
]);
623 float accGetVibrationLevel(void)
628 uint32_t accGetClipCount(void)
630 return acc
.accClipCount
;
633 bool accIsClipped(void)
635 return acc
.isClipped
;
638 void accSetCalibrationValues(void)
640 if (!ARMING_FLAG(SIMULATOR_MODE_SITL
) &&
641 ((accelerometerConfig()->accZero
.raw
[X
] == 0) && (accelerometerConfig()->accZero
.raw
[Y
] == 0) && (accelerometerConfig()->accZero
.raw
[Z
] == 0) &&
642 (accelerometerConfig()->accGain
.raw
[X
] == 4096) && (accelerometerConfig()->accGain
.raw
[Y
] == 4096) &&(accelerometerConfig()->accGain
.raw
[Z
] == 4096))) {
643 DISABLE_STATE(ACCELEROMETER_CALIBRATED
);
646 ENABLE_STATE(ACCELEROMETER_CALIBRATED
);
650 void accInitFilters(void)
652 accSoftLpfFilterApplyFn
= nullFilterApply
;
654 if (acc
.accTargetLooptime
&& accelerometerConfig()->acc_lpf_hz
) {
656 switch (accelerometerConfig()->acc_soft_lpf_type
)
659 accSoftLpfFilterApplyFn
= (filterApplyFnPtr
)pt1FilterApply
;
660 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
661 accSoftLpfFilter
[axis
] = &accFilter
[axis
].pt1
;
662 pt1FilterInit(accSoftLpfFilter
[axis
], accelerometerConfig()->acc_lpf_hz
, US2S(acc
.accTargetLooptime
));
666 accSoftLpfFilterApplyFn
= (filterApplyFnPtr
)biquadFilterApply
;
667 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
668 accSoftLpfFilter
[axis
] = &accFilter
[axis
].biquad
;
669 biquadFilterInitLPF(accSoftLpfFilter
[axis
], accelerometerConfig()->acc_lpf_hz
, acc
.accTargetLooptime
);
676 const float accDt
= US2S(acc
.accTargetLooptime
);
677 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
678 pt1FilterInit(&accVibeFloorFilter
[axis
], ACC_VIBE_FLOOR_FILT_HZ
, accDt
);
679 pt1FilterInit(&accVibeFilter
[axis
], ACC_VIBE_FILT_HZ
, accDt
);
682 STATIC_FASTRAM biquadFilter_t accFilterNotch
[XYZ_AXIS_COUNT
];
683 accNotchFilterApplyFn
= nullFilterApply
;
685 if (acc
.accTargetLooptime
&& accelerometerConfig()->acc_notch_hz
) {
686 accNotchFilterApplyFn
= (filterApplyFnPtr
)biquadFilterApply
;
687 for (int axis
= 0; axis
< XYZ_AXIS_COUNT
; axis
++) {
688 accNotchFilter
[axis
] = &accFilterNotch
[axis
];
689 biquadFilterInitNotch(accNotchFilter
[axis
], acc
.accTargetLooptime
, accelerometerConfig()->acc_notch_hz
, accelerometerConfig()->acc_notch_cutoff
);
695 bool accIsHealthy(void)