| Alexander Hansen | 46a755f | 2025-10-27 16:31:08 +0100 | [diff] [blame] | 1 | // SPDX-License-Identifier: Apache-2.0 |
| 2 | // SPDX-FileCopyrightText: Copyright 2017 Google Inc |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 3 | |
| James Feist | 22c257a | 2018-08-31 14:07:12 -0700 | [diff] [blame] | 4 | #include "pidcontroller.hpp" |
| Patrick Venture | da4a5dd | 2018-08-31 09:42:48 -0700 | [diff] [blame] | 5 | |
| 6 | #include "ec/pid.hpp" |
| 7 | |
| James Feist | 572c43d | 2019-01-31 15:52:22 -0800 | [diff] [blame] | 8 | #include <cmath> |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 9 | |
| Patrick Venture | a076487 | 2020-08-08 07:48:43 -0700 | [diff] [blame] | 10 | namespace pid_control |
| 11 | { |
| 12 | |
| Delphine CC Chiu | 9788963 | 2023-11-06 11:32:46 +0800 | [diff] [blame] | 13 | double PIDController::calPIDOutput(double setpt, double input, |
| 14 | ec::pid_info_t* info) |
| 15 | { |
| 16 | double output; |
| 17 | auto name = getID(); |
| 18 | |
| 19 | if (info->checkHysterWithSetpt) |
| 20 | { |
| 21 | // Over the hysteresis bounds, keep counting pid |
| 22 | if (input > (setpt + info->positiveHysteresis)) |
| 23 | { |
| 24 | // Calculate new output |
| 25 | output = ec::pid(info, input, setpt, &name); |
| 26 | |
| 27 | // this variable isn't actually used in this context, but we're |
| Manojkiran Eda | 7ca8887 | 2024-06-17 11:55:48 +0530 | [diff] [blame] | 28 | // setting it here in case somebody uses it later it's the correct |
| Delphine CC Chiu | 9788963 | 2023-11-06 11:32:46 +0800 | [diff] [blame] | 29 | // value |
| 30 | lastInput = input; |
| 31 | } |
| 32 | // Under the hysteresis bounds, initialize pid |
| 33 | else if (input < (setpt - info->negativeHysteresis)) |
| 34 | { |
| 35 | lastInput = setpt; |
| 36 | info->integral = 0; |
| 37 | output = 0; |
| 38 | } |
| 39 | // inside the hysteresis bounds, keep last output |
| 40 | else |
| 41 | { |
| 42 | lastInput = input; |
| 43 | output = info->lastOutput; |
| 44 | } |
| 45 | |
| 46 | info->lastOutput = output; |
| 47 | } |
| 48 | else |
| 49 | { |
| 50 | // if no hysteresis, maintain previous behavior |
| 51 | if (info->positiveHysteresis == 0 && info->negativeHysteresis == 0) |
| 52 | { |
| 53 | // Calculate new output |
| 54 | output = ec::pid(info, input, setpt, &name); |
| 55 | |
| 56 | // this variable isn't actually used in this context, but we're |
| Manojkiran Eda | 7ca8887 | 2024-06-17 11:55:48 +0530 | [diff] [blame] | 57 | // setting it here in case somebody uses it later it's the correct |
| Delphine CC Chiu | 9788963 | 2023-11-06 11:32:46 +0800 | [diff] [blame] | 58 | // value |
| 59 | lastInput = input; |
| 60 | } |
| 61 | else |
| 62 | { |
| 63 | // initialize if the value is not set (NAN) or abnormal (+INF or |
| 64 | // -INF) |
| 65 | if (!(std::isfinite(lastInput))) |
| 66 | { |
| 67 | lastInput = input; |
| 68 | } |
| 69 | |
| 70 | // if reading is outside of hysteresis bounds, use it for reading, |
| 71 | // otherwise use last reading without updating it first |
| 72 | else if ((input - lastInput) > info->positiveHysteresis) |
| 73 | { |
| 74 | lastInput = input; |
| 75 | } |
| 76 | else if ((lastInput - input) > info->negativeHysteresis) |
| 77 | { |
| 78 | lastInput = input; |
| 79 | } |
| 80 | |
| 81 | output = ec::pid(info, lastInput, setpt, &name); |
| 82 | } |
| 83 | } |
| 84 | |
| 85 | return output; |
| 86 | } |
| 87 | |
| James Feist | 22c257a | 2018-08-31 14:07:12 -0700 | [diff] [blame] | 88 | void PIDController::process(void) |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 89 | { |
| Patrick Venture | 5f59c0f | 2018-11-11 12:55:14 -0800 | [diff] [blame] | 90 | double input; |
| 91 | double setpt; |
| 92 | double output; |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 93 | |
| 94 | // Get setpt value |
| Patrick Venture | 563a356 | 2018-10-30 09:31:26 -0700 | [diff] [blame] | 95 | setpt = setptProc(); |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 96 | |
| 97 | // Get input value |
| Patrick Venture | 563a356 | 2018-10-30 09:31:26 -0700 | [diff] [blame] | 98 | input = inputProc(); |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 99 | |
| James Feist | 572c43d | 2019-01-31 15:52:22 -0800 | [diff] [blame] | 100 | auto info = getPIDInfo(); |
| 101 | |
| Delphine CC Chiu | 9788963 | 2023-11-06 11:32:46 +0800 | [diff] [blame] | 102 | // Calculate output value |
| 103 | output = calPIDOutput(setpt, input, info); |
| James Feist | 572c43d | 2019-01-31 15:52:22 -0800 | [diff] [blame] | 104 | |
| Delphine CC Chiu | 9788963 | 2023-11-06 11:32:46 +0800 | [diff] [blame] | 105 | info->lastOutput = output; |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 106 | |
| 107 | // Output new value |
| Patrick Venture | 563a356 | 2018-10-30 09:31:26 -0700 | [diff] [blame] | 108 | outputProc(output); |
| Patrick Venture | d801218 | 2018-03-08 08:21:38 -0800 | [diff] [blame] | 109 | |
| 110 | return; |
| 111 | } |
| Patrick Venture | a076487 | 2020-08-08 07:48:43 -0700 | [diff] [blame] | 112 | |
| 113 | } // namespace pid_control |