pid.c 2.2 KB

12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364656667686970717273
  1. /*
  2. * Copyright (c) 2022 OpenLuat & AirM2M
  3. *
  4. * Permission is hereby granted, free of charge, to any person obtaining a copy of
  5. * this software and associated documentation files (the "Software"), to deal in
  6. * the Software without restriction, including without limitation the rights to
  7. * use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
  8. * the Software, and to permit persons to whom the Software is furnished to do so,
  9. * subject to the following conditions:
  10. *
  11. * The above copyright notice and this permission notice shall be included in all
  12. * copies or substantial portions of the Software.
  13. *
  14. * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
  15. * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS
  16. * FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR
  17. * COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER
  18. * IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN
  19. * CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
  20. */
  21. #include "pid.h"
  22. VAR_TYPE PID_IncrementalCal(PID_VarType *PIDVar)
  23. {
  24. VAR_TYPE out;
  25. VAR_TYPE ep, ei, ed;
  26. PIDVar->e0 = PIDVar->Target - PIDVar->Measure;
  27. ep = PIDVar->e0 - PIDVar->e1;
  28. ei = PIDVar->e0;
  29. ed = PIDVar->e0 - 2*PIDVar->e1 + PIDVar->e2;
  30. out = PIDVar->Kp*ep + PIDVar->Ki*ei + PIDVar->Kd*ed;
  31. if (out > PIDVar->Scope) //输出限幅
  32. {
  33. out = PIDVar->Scope;
  34. }
  35. else if (out < -PIDVar->Scope)
  36. {
  37. out = -PIDVar->Scope;
  38. }
  39. PIDVar->e2 = PIDVar->e1;
  40. PIDVar->e1 = PIDVar->e0;
  41. return out;
  42. }
  43. VAR_TYPE PID_PositionCal(PID_VarType *PIDVar)
  44. {
  45. VAR_TYPE pe, ie, de;
  46. VAR_TYPE out = 0;
  47. PIDVar->e0 = PIDVar->Target - PIDVar->Measure; //计算当前误差
  48. PIDVar->Sum += PIDVar->e0; //误差积分
  49. de = PIDVar->e0 - PIDVar->e1; //误差微分
  50. pe = PIDVar->e0;
  51. ie = PIDVar->Sum;
  52. PIDVar->e1 = PIDVar->e0;
  53. out = pe*(PIDVar->Kp) + ie*(PIDVar->Ki) + de*(PIDVar->Kd);
  54. if (out > PIDVar->Scope) //输出限幅
  55. {
  56. out = PIDVar->Scope;
  57. }
  58. else if (out < -PIDVar->Scope)
  59. {
  60. out = -PIDVar->Scope;
  61. }
  62. return out;
  63. }