Motor control

Dependencies:   mbed QEI HIDScope biquadFilter MODSERIAL FastPWM

Committer:
BasB
Date:
Mon Oct 14 09:06:08 2019 +0000
Revision:
9:08a7a8e59a6a
Parent:
8:c6c94d55b088
Child:
10:b871d1b05787
Voordat we heel veel aan gaan passen. Voordat we putty er bij gingen betrekken

Who changed what in which revision?

UserRevisionLine numberNew contents of line
freek100 0:e4858e2df9c7 1 #include "mbed.h"
freek100 3:07fedd2e252c 2 #include "HIDScope.h"
freek100 0:e4858e2df9c7 3 #include "QEI.h"
freek100 0:e4858e2df9c7 4 #include "MODSERIAL.h"
freek100 3:07fedd2e252c 5 #include "BiQuad.h"
freek100 0:e4858e2df9c7 6 #include "FastPWM.h"
freek100 0:e4858e2df9c7 7
freek100 6:1c0b6e55e900 8 // Button and potmeter1 control
freek100 0:e4858e2df9c7 9 InterruptIn button1(D11);
freek100 0:e4858e2df9c7 10 InterruptIn button2(D10);
freek100 6:1c0b6e55e900 11 AnalogIn potmeter1(A0);
freek100 6:1c0b6e55e900 12 AnalogIn potmeter2(A1);
freek100 8:c6c94d55b088 13 AnalogIn potmeter3(A2);
BasB 9:08a7a8e59a6a 14 AnalogIn potmeter4(A3);
freek100 0:e4858e2df9c7 15 // Encoder
freek100 0:e4858e2df9c7 16 DigitalIn encA(D13);
freek100 0:e4858e2df9c7 17 DigitalIn encB(D12);
freek100 0:e4858e2df9c7 18 QEI encoder(D13,D12,NC,64,QEI::X4_ENCODING);
freek100 2:d7286c36595f 19 float Ts = 0.01;
freek100 1:08e8cc33fcae 20 float angle;
freek100 1:08e8cc33fcae 21 float omega;
freek100 1:08e8cc33fcae 22
freek100 0:e4858e2df9c7 23
freek100 0:e4858e2df9c7 24 // Motor
freek100 0:e4858e2df9c7 25 DigitalOut motor2Direction(D4);
freek100 0:e4858e2df9c7 26 FastPWM motor2Power(D5);
freek100 0:e4858e2df9c7 27 DigitalOut motor1Direction(D7);
freek100 0:e4858e2df9c7 28 FastPWM motor1Power(D6);
freek100 0:e4858e2df9c7 29
freek100 5:17aa878564d0 30 volatile int motor1Toggle = 1;
freek100 5:17aa878564d0 31
freek100 0:e4858e2df9c7 32 //Motorcontrol
freek100 0:e4858e2df9c7 33 bool motordir;
freek100 0:e4858e2df9c7 34 double motorpwm;
BasB 9:08a7a8e59a6a 35 float u1= 0;
freek100 1:08e8cc33fcae 36 double u2;
freek100 0:e4858e2df9c7 37 double potValue;
freek100 1:08e8cc33fcae 38 double pi2= 6.283185;
freek100 2:d7286c36595f 39 float e; //e = error
freek100 8:c6c94d55b088 40 float Kp;
freek100 8:c6c94d55b088 41 float Ki;
freek100 7:20a802dfe664 42 float Kd;
freek100 2:d7286c36595f 43 float u_k;
freek100 2:d7286c36595f 44 float u_i;
freek100 2:d7286c36595f 45 float u_d;
freek100 3:07fedd2e252c 46
freek100 4:e7d50c6a7c53 47 //Hidscope
freek100 4:e7d50c6a7c53 48 HIDScope scope(3); //Going to send 3 channels of data. To access data go to 'http:/localhost:18082/' after starting HIDScope application.
freek100 0:e4858e2df9c7 49 // PC connection
freek100 0:e4858e2df9c7 50 MODSERIAL pc(USBTX, USBRX);
freek100 0:e4858e2df9c7 51
freek100 0:e4858e2df9c7 52 // Intializing tickers
freek100 0:e4858e2df9c7 53 Ticker motorTicker;
freek100 0:e4858e2df9c7 54 Ticker controlTicker;
freek100 0:e4858e2df9c7 55 Ticker directionTicker;
freek100 1:08e8cc33fcae 56 Ticker encoderTicker;
freek100 4:e7d50c6a7c53 57 Ticker scopeTicker;
freek100 0:e4858e2df9c7 58
freek100 0:e4858e2df9c7 59 const float PWM_period = 1e-6;
freek100 0:e4858e2df9c7 60
freek100 0:e4858e2df9c7 61 volatile int counts; // Encoder counts
freek100 0:e4858e2df9c7 62 volatile int countsPrev = 0;
freek100 0:e4858e2df9c7 63 volatile int deltaCounts;
freek100 0:e4858e2df9c7 64
freek100 0:e4858e2df9c7 65 float factorin = 6.23185/64; // Convert encoder counts to angle in rad
freek100 0:e4858e2df9c7 66 float gearratio = 131.25; // Gear ratio of gearbox
freek100 0:e4858e2df9c7 67
freek100 2:d7286c36595f 68
freek100 2:d7286c36595f 69 float PID_controller(float e){
freek100 2:d7286c36595f 70 static float error_integral=0;
freek100 3:07fedd2e252c 71 static float e_prev=e;
freek100 3:07fedd2e252c 72 static BiQuad LowPassFilter(0.0640,0.1279,0.0640,-1.1683,0.4241);
freek100 2:d7286c36595f 73
freek100 2:d7286c36595f 74 //Proportional part:
BasB 9:08a7a8e59a6a 75 //Kp=2*potmeter1.read();
BasB 9:08a7a8e59a6a 76 Kp=101;
freek100 2:d7286c36595f 77 u_k=Kp*e;
freek100 2:d7286c36595f 78
freek100 2:d7286c36595f 79 //Integral part
BasB 9:08a7a8e59a6a 80 Ki=10;
freek100 2:d7286c36595f 81 error_integral=error_integral+e*Ts;
freek100 2:d7286c36595f 82 u_i=Ki*error_integral;
freek100 2:d7286c36595f 83
freek100 3:07fedd2e252c 84 //Derivative part
BasB 9:08a7a8e59a6a 85 //Kd=2*potmeter2.read();
BasB 9:08a7a8e59a6a 86 Kd=10;
freek100 3:07fedd2e252c 87 float error_derivative =(e-e_prev)/Ts;
freek100 3:07fedd2e252c 88 float filtered_error_derivative = LowPassFilter.step(error_derivative);
freek100 3:07fedd2e252c 89 u_d=Kd*filtered_error_derivative;
freek100 3:07fedd2e252c 90 e_prev=e;
freek100 3:07fedd2e252c 91
freek100 2:d7286c36595f 92 // Sum and return
freek100 3:07fedd2e252c 93 return u_k+u_i+u_d;
freek100 2:d7286c36595f 94 }
freek100 2:d7286c36595f 95
freek100 4:e7d50c6a7c53 96
freek100 4:e7d50c6a7c53 97 void readEncoder()
freek100 4:e7d50c6a7c53 98 {
freek100 4:e7d50c6a7c53 99 counts = encoder.getPulses();
freek100 4:e7d50c6a7c53 100 deltaCounts = counts - countsPrev;
freek100 4:e7d50c6a7c53 101
freek100 4:e7d50c6a7c53 102 countsPrev = counts;
freek100 4:e7d50c6a7c53 103 }
BasB 9:08a7a8e59a6a 104
BasB 9:08a7a8e59a6a 105 void togglehoek(){
BasB 9:08a7a8e59a6a 106 u1= 0.25*pi2+u1;
BasB 9:08a7a8e59a6a 107 }
BasB 9:08a7a8e59a6a 108
freek100 0:e4858e2df9c7 109 void motorControl()
freek100 0:e4858e2df9c7 110 {
BasB 9:08a7a8e59a6a 111 button1.fall(&togglehoek);
freek100 1:08e8cc33fcae 112 angle = counts * factorin / gearratio; // Angle of motor shaft in rad
freek100 2:d7286c36595f 113 omega = deltaCounts / Ts * factorin / gearratio; // Angular velocity of motor shaft in rad/s
freek100 6:1c0b6e55e900 114 potValue= potmeter1.read();
BasB 9:08a7a8e59a6a 115 //u1= (potValue*2*pi2)-pi2;
freek100 2:d7286c36595f 116 e=u1-angle;
freek100 2:d7286c36595f 117
freek100 2:d7286c36595f 118 u2=PID_controller(e);
freek100 2:d7286c36595f 119
freek100 2:d7286c36595f 120 motorpwm= abs(u2);
freek100 1:08e8cc33fcae 121 if (u2<0){
freek100 0:e4858e2df9c7 122 motordir= 0;}
freek100 0:e4858e2df9c7 123 else {
freek100 0:e4858e2df9c7 124 motordir= 1;}
freek100 5:17aa878564d0 125 motor1Power.pulsewidth(motorpwm * PWM_period*motor1Toggle );
freek100 0:e4858e2df9c7 126 motor1Direction= motordir;
freek100 0:e4858e2df9c7 127 }
freek100 0:e4858e2df9c7 128
freek100 4:e7d50c6a7c53 129 void Plotje()
freek100 1:08e8cc33fcae 130 {
freek100 4:e7d50c6a7c53 131 scope.set(0,u1); //gewenste hoek
freek100 4:e7d50c6a7c53 132 scope.set(1,angle); //Gemeten hoek
freek100 4:e7d50c6a7c53 133 scope.set(2,e); //verschil in gewenste en gemeten hoek
freek100 1:08e8cc33fcae 134
freek100 4:e7d50c6a7c53 135 scope.send(); //send what's in scope memory to PC
freek100 1:08e8cc33fcae 136 }
freek100 0:e4858e2df9c7 137
freek100 5:17aa878564d0 138 void toggleMotor()
freek100 5:17aa878564d0 139 {
freek100 5:17aa878564d0 140 motor1Toggle = !motor1Toggle;
freek100 5:17aa878564d0 141 }
freek100 4:e7d50c6a7c53 142
freek100 0:e4858e2df9c7 143 int main()
freek100 0:e4858e2df9c7 144 {
freek100 0:e4858e2df9c7 145 pc.baud(115200);
freek100 0:e4858e2df9c7 146 pc.printf("\r\nStarting...\r\n\r\n");
freek100 0:e4858e2df9c7 147
freek100 0:e4858e2df9c7 148 motor1Power.period(PWM_period);
freek100 0:e4858e2df9c7 149 motorTicker.attach(motorControl, 0.01);
freek100 4:e7d50c6a7c53 150 scopeTicker.attach(Plotje, 0.01);
freek100 2:d7286c36595f 151 encoderTicker.attach(readEncoder, Ts);
freek100 5:17aa878564d0 152
freek100 5:17aa878564d0 153 button2.fall(&toggleMotor);
BasB 9:08a7a8e59a6a 154
freek100 0:e4858e2df9c7 155 while (true) {
freek100 0:e4858e2df9c7 156
freek100 7:20a802dfe664 157 //pc.printf("Potmeter: %d \r\n", potValue,);
freek100 7:20a802dfe664 158 //pc.printf("Counts: %i DeltaCounts: %i\r\n", counts, deltaCounts);
freek100 1:08e8cc33fcae 159 pc.printf("Angle: %f Omega: %f\r\n", angle, omega);
freek100 2:d7286c36595f 160 pc.printf("U1: %f Error: %f \r\n",u1, e);
freek100 7:20a802dfe664 161 pc.printf("Kp: %f Kd: %f Ki: %f \r\n", Kp, Kd, Ki);
freek100 6:1c0b6e55e900 162
freek100 0:e4858e2df9c7 163 wait(0.5);
freek100 0:e4858e2df9c7 164 }
freek100 0:e4858e2df9c7 165 }