pid.py 2.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475
  1. import numpy as np
  2. from numbers import Number
  3. from openpilot.common.numpy_fast import clip, interp
  4. class PIDController:
  5. def __init__(self, k_p, k_i, k_f=0., k_d=0., pos_limit=1e308, neg_limit=-1e308, rate=100):
  6. self._k_p = k_p
  7. self._k_i = k_i
  8. self._k_d = k_d
  9. self.k_f = k_f # feedforward gain
  10. if isinstance(self._k_p, Number):
  11. self._k_p = [[0], [self._k_p]]
  12. if isinstance(self._k_i, Number):
  13. self._k_i = [[0], [self._k_i]]
  14. if isinstance(self._k_d, Number):
  15. self._k_d = [[0], [self._k_d]]
  16. self.pos_limit = pos_limit
  17. self.neg_limit = neg_limit
  18. self.i_unwind_rate = 0.3 / rate
  19. self.i_rate = 1.0 / rate
  20. self.speed = 0.0
  21. self.reset()
  22. @property
  23. def k_p(self):
  24. return interp(self.speed, self._k_p[0], self._k_p[1])
  25. @property
  26. def k_i(self):
  27. return interp(self.speed, self._k_i[0], self._k_i[1])
  28. @property
  29. def k_d(self):
  30. return interp(self.speed, self._k_d[0], self._k_d[1])
  31. @property
  32. def error_integral(self):
  33. return self.i/self.k_i
  34. def reset(self):
  35. self.p = 0.0
  36. self.i = 0.0
  37. self.d = 0.0
  38. self.f = 0.0
  39. self.control = 0
  40. def update(self, error, error_rate=0.0, speed=0.0, override=False, feedforward=0., freeze_integrator=False):
  41. self.speed = speed
  42. self.p = float(error) * self.k_p
  43. self.f = feedforward * self.k_f
  44. self.d = error_rate * self.k_d
  45. if override:
  46. self.i -= self.i_unwind_rate * float(np.sign(self.i))
  47. else:
  48. i = self.i + error * self.k_i * self.i_rate
  49. control = self.p + i + self.d + self.f
  50. # Update when changing i will move the control away from the limits
  51. # or when i will move towards the sign of the error
  52. if ((error >= 0 and (control <= self.pos_limit or i < 0.0)) or
  53. (error <= 0 and (control >= self.neg_limit or i > 0.0))) and \
  54. not freeze_integrator:
  55. self.i = i
  56. control = self.p + self.i + self.d + self.f
  57. self.control = clip(control, self.neg_limit, self.pos_limit)
  58. return self.control