PI_Controller.h 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176
  1. #ifndef _PI_Contrller_H__
  2. #define _PI_Contrller_H__
  3. #include "math/fast_math.h"
  4. #include "bsp/bsp.h"
  5. #if 0
  6. typedef struct {
  7. s16q5_t kp;
  8. s16q5_t ki;
  9. s16q5_t max;
  10. s16q5_t min;
  11. s16q5_t Ui;
  12. s16q14_t DT;
  13. }PI_Controller;
  14. static __INLINE void PI_Controller_Reset(PI_Controller *pi, s16q5_t init) {
  15. pi->Ui = init;
  16. }
  17. static __INLINE s16q5_t PI_Controller_run(PI_Controller *pi, s16q5_t err) {
  18. s16q5_t kp_err = S16_mul(err,pi->kp, 5);
  19. s16q5_t ki_err = S16_mul(err,pi->ki, 5);
  20. s16q5_t integral = S16_mul(ki_err, pi->DT, 14);
  21. pi->Ui = MATH_sat(pi->Ui + integral, pi->min, pi->max);
  22. s32q14_t out = pi->Ui + kp_err;
  23. return MATH_sat(out, pi->min, pi->max);
  24. }
  25. #else
  26. typedef struct {
  27. float kp;
  28. float ki;
  29. float kd;
  30. float max;
  31. float min;
  32. float Ui;
  33. float sat;
  34. float DT;
  35. bool is_sat;
  36. float last_err;
  37. float b_dstart;
  38. }PI_Controller;
  39. static __INLINE void PI_Controller_max(PI_Controller *pi, float max, float min) {
  40. pi->max = max;
  41. pi->min = min;
  42. }
  43. static __INLINE void PI_Controller_Reset(PI_Controller *pi, float init) {
  44. pi->Ui = (init);
  45. pi->sat = 0.0f;
  46. pi->is_sat = false;
  47. pi->b_dstart = false;
  48. }
  49. static __INLINE float PI_Controller_Run(PI_Controller *pi, float err) {
  50. float kp_err = (err) * pi->kp;
  51. float ki_err = (err) * pi->ki;
  52. float integral = ki_err * pi->DT;
  53. pi->Ui = MATH_sat(pi->Ui + integral, pi->min, pi->max);
  54. float out = pi->Ui + kp_err;
  55. float sat_out = (MATH_sat(out, pi->min, pi->max));
  56. if (out != sat_out) {
  57. pi->is_sat = true;
  58. }else {
  59. pi->is_sat = false;
  60. }
  61. return sat_out;
  62. }
  63. static __INLINE float PI_Controller_RunLimit(PI_Controller *pi, float err) {
  64. float kp_err = (err) * pi->kp;
  65. float ki_err = (err) * pi->ki;
  66. float integral = ki_err * pi->DT;
  67. pi->Ui = MATH_sat(pi->Ui + integral, pi->min, pi->max);
  68. float out = pi->Ui + kp_err;
  69. float sat_out = MATH_sat(out, pi->min, pi->max);
  70. return sat_out;
  71. }
  72. static __INLINE float PI_Controller_RunWithDiff(PI_Controller *pi, float err) {
  73. float kp_err = (err) * pi->kp;
  74. float ki_err = (err) * pi->ki;
  75. float integral = ki_err * pi->DT;
  76. if (!pi->b_dstart) {
  77. pi->last_err = err;
  78. pi->b_dstart = true;
  79. }
  80. pi->Ui = MATH_sat(pi->Ui + integral, pi->min, pi->max);
  81. float out = pi->Ui + kp_err + pi->kd * (err - pi->last_err);
  82. float sat_out = (MATH_sat(out, pi->min, pi->max));
  83. if (out != sat_out) {
  84. pi->is_sat = true;
  85. }else {
  86. pi->is_sat = false;
  87. }
  88. pi->last_err = err;
  89. return sat_out;
  90. }
  91. #if 0
  92. static __INLINE float PI_Controller_RunSat(PI_Controller *pi, float err) {
  93. float kp_err = (err) * pi->kp;
  94. float ki_err = (err) * pi->ki;
  95. float integral = ki_err * pi->DT;
  96. pi->Ui = pi->Ui + integral + pi->sat * pi->kd;
  97. float out = pi->Ui + kp_err ;
  98. float out_sat = MATH_sat(out, pi->min, pi->max);
  99. if (out != out_sat) {
  100. pi->is_sat = true;
  101. }else {
  102. pi->is_sat = false;
  103. }
  104. pi->sat = out_sat - out;
  105. return out_sat;
  106. }
  107. #endif
  108. static __INLINE float PI_Controller_RunSerial(PI_Controller *pi, float err) {
  109. float kp_err = (err) * pi->kp;//S16_mul(err,pi->kp, 5);
  110. float ki_err = (kp_err) * pi->ki;
  111. float integral = ki_err * pi->DT;
  112. pi->Ui = MATH_sat(pi->Ui + integral, pi->min, pi->max);
  113. float out = pi->Ui + kp_err;
  114. float sat_out = (MATH_sat(out, pi->min, pi->max));
  115. if (out != sat_out) {
  116. pi->is_sat = true;
  117. }else {
  118. pi->is_sat = false;
  119. }
  120. return sat_out;
  121. }
  122. static __INLINE float _fmod(float v, s32 m) {
  123. int v_i = (int)v;
  124. int m_i = v_i % m;
  125. return (v-v_i) + m_i;
  126. }
  127. typedef struct {
  128. float observer;
  129. float kp;
  130. float ki;
  131. float Ui;
  132. float out;
  133. s32 max_wp;
  134. bool ob_wp;
  135. float DT;
  136. }PLL_t;
  137. static __INLINE void PLL_Reset(PLL_t *pll) {
  138. pll->observer = 0.0f;
  139. pll->out = 0.0f;
  140. pll->ob_wp = false;
  141. }
  142. /*
  143. Transfer func:
  144. ki*s/(s*s + kp*s + ki)
  145. ki = (kp/2) * (kp/2)
  146. Wn = kp/2 ==> kp = 2 * Wn, Wn为PLL的带宽
  147. */
  148. static __INLINE float PLL_run(PLL_t *pll, float sample, float comp) {
  149. float observer = pll->observer - comp;
  150. if (comp < 0) {
  151. observer = -comp - pll->observer;
  152. }
  153. float delta = sample - observer;
  154. pll->observer = observer + (pll->out + pll->kp * delta) * pll->DT;
  155. pll->out += pll->ki * delta * pll->DT;
  156. return pll->out;
  157. }
  158. #endif
  159. #endif /*_PI_Contrller_H__*/