Zbigniew Druzbacki / Mbed 2 deprecated SPI_Slave

SPI slave program to enable communication between the FPGA and the STM32L432 board.

Dependencies:   mbed

IMUs.cpp@11:366f1186c121, 2019-03-19 (annotated)

Committer:
Zbyszek
Date:
Tue Mar 19 01:26:11 2019 +0000
Revision:
11:366f1186c121
Parent:
7:0e9af5986488
Child:
13:c7e8e277f884
Several IMUs can now be read using STM32L432.

Who changed what in which revision?

UserRevisionLine numberNew contents of line
Zbyszek 6:0ebecfecadc9 1 #include "Structures.h"
Zbyszek 7:0e9af5986488 2 #include "Quaternions.h"
Zbyszek 6:0ebecfecadc9 3 #include "IMUs.h"
Zbyszek 6:0ebecfecadc9 4 #include "mbed.h"
Zbyszek 6:0ebecfecadc9 5
Zbyszek 6:0ebecfecadc9 6 //IMU Class Constructor
Zbyszek 6:0ebecfecadc9 7 IMU::IMU(char IMU_ID, double OffsetAX, double OffsetAY, double OffsetAZ, double OffsetGX, double OffsetGY, double OffsetGZ, char SSFA, char SSFG) {
Zbyszek 6:0ebecfecadc9 8
Zbyszek 6:0ebecfecadc9 9 AccelerometerOffset.x = OffsetAX;
Zbyszek 6:0ebecfecadc9 10 AccelerometerOffset.y = OffsetAY;
Zbyszek 6:0ebecfecadc9 11 AccelerometerOffset.z = OffsetAZ;
Zbyszek 6:0ebecfecadc9 12
Zbyszek 6:0ebecfecadc9 13 GyroscopeOffset.x = OffsetGX;
Zbyszek 6:0ebecfecadc9 14 GyroscopeOffset.y = OffsetGY;
Zbyszek 6:0ebecfecadc9 15 GyroscopeOffset.z = OffsetGZ;
Zbyszek 6:0ebecfecadc9 16
Zbyszek 6:0ebecfecadc9 17
Zbyszek 6:0ebecfecadc9 18 IMU_Identifier = IMU_ID;
Zbyszek 6:0ebecfecadc9 19
Zbyszek 6:0ebecfecadc9 20 switch(SSFA) {
Zbyszek 6:0ebecfecadc9 21 case 0:
Zbyszek 6:0ebecfecadc9 22 accelSSF = 0.00006103515625f;
Zbyszek 6:0ebecfecadc9 23 break;
Zbyszek 6:0ebecfecadc9 24 case 1:
Zbyszek 6:0ebecfecadc9 25 accelSSF = 0.0001220703125f;
Zbyszek 6:0ebecfecadc9 26 break;
Zbyszek 6:0ebecfecadc9 27 case 2:
Zbyszek 6:0ebecfecadc9 28 accelSSF = 0.000244140625f;
Zbyszek 6:0ebecfecadc9 29 break;
Zbyszek 6:0ebecfecadc9 30 case 3:
Zbyszek 6:0ebecfecadc9 31 accelSSF = 0.00048828125f;
Zbyszek 6:0ebecfecadc9 32 break;
Zbyszek 6:0ebecfecadc9 33 default:
Zbyszek 6:0ebecfecadc9 34 break;
Zbyszek 6:0ebecfecadc9 35 }
Zbyszek 6:0ebecfecadc9 36
Zbyszek 6:0ebecfecadc9 37 switch(SSFG) {
Zbyszek 6:0ebecfecadc9 38 case 0:
Zbyszek 6:0ebecfecadc9 39 gyroSSF = 0.00763358778625954198473282442748f;
Zbyszek 6:0ebecfecadc9 40 break;
Zbyszek 6:0ebecfecadc9 41 case 1:
Zbyszek 6:0ebecfecadc9 42 gyroSSF = 0.01526717557251908396946564885496f;
Zbyszek 6:0ebecfecadc9 43 break;
Zbyszek 6:0ebecfecadc9 44 case 2:
Zbyszek 6:0ebecfecadc9 45 gyroSSF = 0.03048780487804878048780487804878f;
Zbyszek 6:0ebecfecadc9 46 break;
Zbyszek 6:0ebecfecadc9 47 case 3:
Zbyszek 6:0ebecfecadc9 48 gyroSSF = 0.06097560975609756097560975609756f;
Zbyszek 6:0ebecfecadc9 49 break;
Zbyszek 6:0ebecfecadc9 50 default:
Zbyszek 6:0ebecfecadc9 51 break;
Zbyszek 6:0ebecfecadc9 52 }
Zbyszek 6:0ebecfecadc9 53
Zbyszek 7:0e9af5986488 54 this->t.start();
Zbyszek 7:0e9af5986488 55
Zbyszek 6:0ebecfecadc9 56 }
Zbyszek 6:0ebecfecadc9 57
Zbyszek 6:0ebecfecadc9 58 //void IMU::CFAngle() {
Zbyszek 6:0ebecfecadc9 59
Zbyszek 6:0ebecfecadc9 60
Zbyszek 6:0ebecfecadc9 61
Zbyszek 6:0ebecfecadc9 62 //}
Zbyszek 6:0ebecfecadc9 63
Zbyszek 6:0ebecfecadc9 64 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 6:0ebecfecadc9 65 IMU_Data IMU::concatenateData(int16_t SamplesPieces[12]) {
Zbyszek 6:0ebecfecadc9 66 //Local-Variables------Local-Variables------Local-Variables------Local-Variables
Zbyszek 6:0ebecfecadc9 67
Zbyszek 6:0ebecfecadc9 68 int16_t MSB = 0; //Store Most Significant Byte of data piece in this variable for processing
Zbyszek 6:0ebecfecadc9 69 int16_t LSB = 0; //Store Least Significant Byte of data piece in this variable for processing
Zbyszek 6:0ebecfecadc9 70 char arrPointer = 0; //Array Pointer
Zbyszek 6:0ebecfecadc9 71
Zbyszek 6:0ebecfecadc9 72 IMU_Data LocalD;
Zbyszek 6:0ebecfecadc9 73
Zbyszek 6:0ebecfecadc9 74 //Procedure--------------Procedure------------Procedure-------------Procedure---
Zbyszek 6:0ebecfecadc9 75 for(char x = 0; x <= 5; x++) {
Zbyszek 6:0ebecfecadc9 76 MSB = SamplesPieces[arrPointer]; //Odd data pieces are MSBs
Zbyszek 6:0ebecfecadc9 77 MSB &= ~(255<<8); //Mask the MSB bits as they are not part of data
Zbyszek 6:0ebecfecadc9 78 MSB = MSB << 8; //Shift the Value as its the MSB of the data piece
Zbyszek 6:0ebecfecadc9 79 arrPointer++; //Increment array pointer
Zbyszek 6:0ebecfecadc9 80 LSB = SamplesPieces[arrPointer]; //Even data pieces are LSBs
Zbyszek 6:0ebecfecadc9 81 LSB &= ~(255 << 8); //Mask the MSB bits as they are not part of data
Zbyszek 6:0ebecfecadc9 82 arrPointer++; //Increment array pointer
Zbyszek 6:0ebecfecadc9 83
Zbyszek 6:0ebecfecadc9 84 switch(x) {
Zbyszek 6:0ebecfecadc9 85 case 0:
Zbyszek 6:0ebecfecadc9 86 LocalD.Ax = (MSB + LSB) + AccelerometerOffset.x; //Combine Accelerometer x-axis data
Zbyszek 6:0ebecfecadc9 87 break;
Zbyszek 6:0ebecfecadc9 88 case 1:
Zbyszek 6:0ebecfecadc9 89 LocalD.Ay = (MSB + LSB) + AccelerometerOffset.y; //Combine Accelerometer y-axis data
Zbyszek 6:0ebecfecadc9 90 break;
Zbyszek 6:0ebecfecadc9 91 case 2:
Zbyszek 6:0ebecfecadc9 92 LocalD.Az = (MSB + LSB) + AccelerometerOffset.z; //Combine Accelerometer z-axis data
Zbyszek 6:0ebecfecadc9 93 break;
Zbyszek 6:0ebecfecadc9 94 case 3:
Zbyszek 6:0ebecfecadc9 95 LocalD.Gx = (MSB + LSB) + GyroscopeOffset.x; //Combine Gyroscope x-axis data
Zbyszek 6:0ebecfecadc9 96 break;
Zbyszek 6:0ebecfecadc9 97 case 4:
Zbyszek 6:0ebecfecadc9 98 LocalD.Gy = (MSB + LSB) + GyroscopeOffset.y; //Combine Gyroscope y-axis data
Zbyszek 6:0ebecfecadc9 99 break;
Zbyszek 6:0ebecfecadc9 100 case 5:
Zbyszek 6:0ebecfecadc9 101 LocalD.Gz = (MSB + LSB) + GyroscopeOffset.z; //Combine Gyroscope z-axis data
Zbyszek 6:0ebecfecadc9 102 break;
Zbyszek 6:0ebecfecadc9 103 default:
Zbyszek 6:0ebecfecadc9 104 break;
Zbyszek 6:0ebecfecadc9 105 }//switch(x)
Zbyszek 6:0ebecfecadc9 106 }//for(char x = 0; x <= 5; x++)
Zbyszek 11:366f1186c121 107
Zbyszek 11:366f1186c121 108 //printf("Accel X: %+2f, Accel Y: %+2f, Accel Z: %+2f, Gyro X: %+2f, Gyro Y: %+2f, Gyro Z: %+2f\n\r", LocalD.Ax, LocalD.Ay, LocalD.Az, LocalD.Gx, LocalD.Gy, LocalD.Gz);
Zbyszek 6:0ebecfecadc9 109
Zbyszek 6:0ebecfecadc9 110 return LocalD;
Zbyszek 6:0ebecfecadc9 111 }
Zbyszek 6:0ebecfecadc9 112 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 6:0ebecfecadc9 113
Zbyszek 6:0ebecfecadc9 114
Zbyszek 6:0ebecfecadc9 115
Zbyszek 6:0ebecfecadc9 116
Zbyszek 6:0ebecfecadc9 117
Zbyszek 6:0ebecfecadc9 118 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 6:0ebecfecadc9 119 IMU_Data IMU::SSFmultiply(IMU_Data RawData) {
Zbyszek 6:0ebecfecadc9 120
Zbyszek 6:0ebecfecadc9 121 IMU_Data localD;
Zbyszek 6:0ebecfecadc9 122
Zbyszek 6:0ebecfecadc9 123
Zbyszek 6:0ebecfecadc9 124 localD.Ax = RawData.Ax * accelSSF;
Zbyszek 6:0ebecfecadc9 125 localD.Ay = RawData.Ay * accelSSF;
Zbyszek 6:0ebecfecadc9 126 localD.Az = RawData.Az * accelSSF;
Zbyszek 6:0ebecfecadc9 127
Zbyszek 6:0ebecfecadc9 128 localD.Gx = RawData.Gx * gyroSSF;
Zbyszek 6:0ebecfecadc9 129 localD.Gy = RawData.Gy * gyroSSF;
Zbyszek 6:0ebecfecadc9 130 localD.Gz = RawData.Gz * gyroSSF;
Zbyszek 6:0ebecfecadc9 131
Zbyszek 7:0e9af5986488 132 //printf("Accel X: %+2f, Accel Y: %+2f, Accel Z: %+2f, Gyro X: %+2f, Gyro Y: %+2f, Gyro Z: %+2f\n\r", LocalD.Ax, LocalD.Ay, LocalD.Az, LocalD.Gx, LocalD.Gy, LocalD.Gz);
Zbyszek 6:0ebecfecadc9 133 return localD;
Zbyszek 6:0ebecfecadc9 134 }
Zbyszek 7:0e9af5986488 135 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 136
Zbyszek 7:0e9af5986488 137
Zbyszek 7:0e9af5986488 138
Zbyszek 7:0e9af5986488 139
Zbyszek 7:0e9af5986488 140
Zbyszek 7:0e9af5986488 141 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 142 IMU_Data IMU::getAngles(IMU_Data SSF_Format_Values, float dt) {
Zbyszek 7:0e9af5986488 143
Zbyszek 7:0e9af5986488 144 IMU_Data newValues, AngleData;
Zbyszek 7:0e9af5986488 145
Zbyszek 7:0e9af5986488 146 newValues = SSF_Format_Values;
Zbyszek 7:0e9af5986488 147
Zbyszek 7:0e9af5986488 148 AngleData.Ax = (atan2f(newValues.Ay, newValues.Az) * 57.29577951f);
Zbyszek 7:0e9af5986488 149 AngleData.Ay = (atan2f((-newValues.Ax), sqrt(pow(newValues.Ay, 2) + pow(newValues.Az, 2) )) * 57.29577951f);
Zbyszek 7:0e9af5986488 150 AngleData.Az = 0; //Cannot calculate angle for this one as it remains constant all the time
Zbyszek 7:0e9af5986488 151
Zbyszek 7:0e9af5986488 152 AngleData.Gx = (newValues.Gx * dt);
Zbyszek 7:0e9af5986488 153 AngleData.Gy = (newValues.Gy * dt);
Zbyszek 7:0e9af5986488 154 AngleData.Gz = (newValues.Gz * dt);
Zbyszek 7:0e9af5986488 155
Zbyszek 7:0e9af5986488 156 //printf("Accel X: %+2f, Accel Y: %+2f, Accel Z: %+2f, Gyro X: %+2f, Gyro Y: %+2f, Gyro Z: %+2f dt: %2f\n\r", AngleData.Ax, AngleData.Ay, AngleData.Az, AngleData.Gx, AngleData.Gy, AngleData.Gz, dt);
Zbyszek 7:0e9af5986488 157 return AngleData;
Zbyszek 7:0e9af5986488 158
Zbyszek 7:0e9af5986488 159 }
Zbyszek 7:0e9af5986488 160 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 161
Zbyszek 7:0e9af5986488 162
Zbyszek 7:0e9af5986488 163
Zbyszek 7:0e9af5986488 164
Zbyszek 7:0e9af5986488 165
Zbyszek 7:0e9af5986488 166 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 167 vector IMU::getAccelAngles(vector SSF_Accel) {
Zbyszek 7:0e9af5986488 168
Zbyszek 7:0e9af5986488 169 vector newValues, AngleData;
Zbyszek 7:0e9af5986488 170
Zbyszek 7:0e9af5986488 171 newValues = SSF_Accel;
Zbyszek 7:0e9af5986488 172
Zbyszek 7:0e9af5986488 173 AngleData.x = (atan2f(newValues.y, newValues.z) * 57.29577951f);
Zbyszek 7:0e9af5986488 174 AngleData.y = (atan2f((-newValues.x), sqrt(pow(newValues.y, 2) + pow(newValues.z, 2) )) * 57.29577951f);
Zbyszek 7:0e9af5986488 175 //AngleData.z = 0; //Cannot calculate angle for this one as it remains constant all the time
Zbyszek 7:0e9af5986488 176
Zbyszek 7:0e9af5986488 177 //printf("Accel X: %+2f, Accel Y: %+2f, Accel Z: %+2f, Gyro X: %+2f, Gyro Y: %+2f, Gyro Z: %+2f dt: %2f\n\r", AngleData.Ax, AngleData.Ay, AngleData.Az, AngleData.Gx, AngleData.Gy, AngleData.Gz, dt);
Zbyszek 7:0e9af5986488 178 return AngleData;
Zbyszek 7:0e9af5986488 179
Zbyszek 7:0e9af5986488 180 }
Zbyszek 7:0e9af5986488 181 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 182
Zbyszek 7:0e9af5986488 183
Zbyszek 7:0e9af5986488 184
Zbyszek 7:0e9af5986488 185
Zbyszek 7:0e9af5986488 186
Zbyszek 7:0e9af5986488 187
Zbyszek 7:0e9af5986488 188
Zbyszek 7:0e9af5986488 189
Zbyszek 7:0e9af5986488 190 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 191 IMU_Data IMU::CalculateAngles(int16_t SamplesPieces[12]){
Zbyszek 7:0e9af5986488 192
Zbyszek 7:0e9af5986488 193 IMU_Data ConcatenatedValues, SSFValues, SimpleAngles;
Zbyszek 7:0e9af5986488 194 static IMU_Data PrevSimpleAngles;
Zbyszek 7:0e9af5986488 195 float deltaT;
Zbyszek 7:0e9af5986488 196
Zbyszek 7:0e9af5986488 197 ConcatenatedValues = this->concatenateData(SamplesPieces);
Zbyszek 7:0e9af5986488 198 SSFValues = this->SSFmultiply(ConcatenatedValues);
Zbyszek 7:0e9af5986488 199
Zbyszek 7:0e9af5986488 200 this->t.stop();
Zbyszek 7:0e9af5986488 201 deltaT = this->t.read();
Zbyszek 7:0e9af5986488 202 SimpleAngles = this->getAngles(SSFValues, deltaT);
Zbyszek 7:0e9af5986488 203
Zbyszek 7:0e9af5986488 204 PrevSimpleAngles.Gx = PrevSimpleAngles.Gx + SimpleAngles.Gx;
Zbyszek 7:0e9af5986488 205 PrevSimpleAngles.Gy = PrevSimpleAngles.Gy + SimpleAngles.Gy;
Zbyszek 7:0e9af5986488 206 PrevSimpleAngles.Gz = PrevSimpleAngles.Gz + SimpleAngles.Gz;
Zbyszek 7:0e9af5986488 207
Zbyszek 7:0e9af5986488 208 this->t.reset();
Zbyszek 7:0e9af5986488 209 this->t.start();
Zbyszek 7:0e9af5986488 210
Zbyszek 7:0e9af5986488 211 //printf("Accel X: %+2f, Accel Y: %+2f, Accel Z: %+2f, Gyro X: %+2f, Gyro Y: %+2f, Gyro Z: %+2f\n\r", PrevSimpleAngles.Ax, PrevSimpleAngles.Ay, PrevSimpleAngles.Az, PrevSimpleAngles.Gx, PrevSimpleAngles.Gy, PrevSimpleAngles.Gz);
Zbyszek 7:0e9af5986488 212 return PrevSimpleAngles;
Zbyszek 7:0e9af5986488 213 }
Zbyszek 7:0e9af5986488 214 //----------------------------------------------------------------------------------------------------------------------------------------------------------
Zbyszek 7:0e9af5986488 215
Zbyszek 7:0e9af5986488 216
Zbyszek 7:0e9af5986488 217 vector IMU::CalculateCFAngles(int16_t SamplesPieces[12]) {
Zbyszek 7:0e9af5986488 218
Zbyszek 7:0e9af5986488 219 IMU_Data ConcatenatedValues,NewAngle, SSFValues;
Zbyszek 7:0e9af5986488 220 float deltaT;
Zbyszek 7:0e9af5986488 221
Zbyszek 7:0e9af5986488 222
Zbyszek 7:0e9af5986488 223 this->t.stop();
Zbyszek 7:0e9af5986488 224 deltaT = this->t.read();
Zbyszek 7:0e9af5986488 225 ConcatenatedValues = this->concatenateData(SamplesPieces);
Zbyszek 7:0e9af5986488 226 SSFValues = this->SSFmultiply(ConcatenatedValues);
Zbyszek 7:0e9af5986488 227 NewAngle = this->getAngles(SSFValues, deltaT);
Zbyszek 7:0e9af5986488 228
Zbyszek 11:366f1186c121 229 this->CFAngle.x = 0.98f*(CFAngle.x + NewAngle.Gx) + 0.02f*NewAngle.Ax;
Zbyszek 11:366f1186c121 230 this->CFAngle.y = 0.98f*(CFAngle.y + NewAngle.Gy) + 0.02f*NewAngle.Ay;
Zbyszek 11:366f1186c121 231 this->CFAngle.z = 0.98f*(CFAngle.z + NewAngle.Gz) + 0.02f*NewAngle.Az;
Zbyszek 7:0e9af5986488 232
Zbyszek 11:366f1186c121 233 t.reset();
Zbyszek 11:366f1186c121 234 t.start();
Zbyszek 7:0e9af5986488 235
Zbyszek 7:0e9af5986488 236 //printf("Accel X: %+2f, Accel Y: %+2f, Accel Z: %+2f, dt: %+2f\n\r", CFAngle.x, CFAngle.y, CFAngle.z, deltaT);
Zbyszek 7:0e9af5986488 237 printf("%+2f,%+2f,%+2f\n\r", CFAngle.x, CFAngle.y, CFAngle.z);
Zbyszek 7:0e9af5986488 238 return CFAngle;
Zbyszek 7:0e9af5986488 239
Zbyszek 7:0e9af5986488 240 }
Zbyszek 7:0e9af5986488 241
Zbyszek 7:0e9af5986488 242
Zbyszek 7:0e9af5986488 243
Zbyszek 7:0e9af5986488 244 vector IMU::CalculateQCFAngles(int16_t SamplesPieces[12]) {
Zbyszek 7:0e9af5986488 245
Zbyszek 7:0e9af5986488 246 IMU_Data ConcatenatedValues,NewAngle, SSFValues;
Zbyszek 7:0e9af5986488 247 float deltaT;
Zbyszek 7:0e9af5986488 248 static Quaternion CF = {.w = 1.0f, .x = 0.0001f, .y = 0.0001f, .z = 0.0001f};
Zbyszek 7:0e9af5986488 249 vector Accel_Angles, Eangles, accelSSFVals, SSFAngularRate, QuaternionAngles;
Zbyszek 7:0e9af5986488 250
Zbyszek 7:0e9af5986488 251
Zbyszek 7:0e9af5986488 252 this->t.stop();
Zbyszek 7:0e9af5986488 253 deltaT = this->t.read();
Zbyszek 7:0e9af5986488 254 ConcatenatedValues = this->concatenateData(SamplesPieces);
Zbyszek 7:0e9af5986488 255 SSFValues = this->SSFmultiply(ConcatenatedValues);
Zbyszek 7:0e9af5986488 256
Zbyszek 7:0e9af5986488 257 //convert to radians per second
Zbyszek 7:0e9af5986488 258 SSFValues.Gx *= 0.0174533f;
Zbyszek 7:0e9af5986488 259 SSFValues.Gy *= 0.0174533f;
Zbyszek 7:0e9af5986488 260 SSFValues.Gz *= 0.0174533f;
Zbyszek 7:0e9af5986488 261
Zbyszek 7:0e9af5986488 262 //Angular rates used for Quaternion gyro integration
Zbyszek 7:0e9af5986488 263 SSFAngularRate.x = SSFValues.Gx;
Zbyszek 7:0e9af5986488 264 SSFAngularRate.y = SSFValues.Gy;
Zbyszek 7:0e9af5986488 265 SSFAngularRate.z = SSFValues.Gz;
Zbyszek 7:0e9af5986488 266
Zbyszek 7:0e9af5986488 267 CF = updateQuaternion(CF, SSFAngularRate, deltaT);
Zbyszek 7:0e9af5986488 268 CF = normaliseQuaternion(CF);
Zbyszek 7:0e9af5986488 269 Eangles = eulerA(CF);
Zbyszek 7:0e9af5986488 270
Zbyszek 7:0e9af5986488 271
Zbyszek 7:0e9af5986488 272 //Get accel ssf values into vector format
Zbyszek 7:0e9af5986488 273 accelSSFVals.x = SSFValues.Ax;
Zbyszek 7:0e9af5986488 274 accelSSFVals.y = SSFValues.Ay;
Zbyszek 7:0e9af5986488 275 accelSSFVals.z = SSFValues.Az;
Zbyszek 7:0e9af5986488 276 Accel_Angles = getAccelAngles(accelSSFVals); //Convert Accel data to angle
Zbyszek 7:0e9af5986488 277
Zbyszek 7:0e9af5986488 278
Zbyszek 7:0e9af5986488 279 CF.x = 0.98f*(Eangles.x) + 0.02f*Accel_Angles.x;
Zbyszek 7:0e9af5986488 280 CF.y = 0.98f*(Eangles.y) + 0.02f*Accel_Angles.y;
Zbyszek 7:0e9af5986488 281 CF.z = 0.98f*(Eangles.z) + 0.02f*Accel_Angles.z;
Zbyszek 7:0e9af5986488 282
Zbyszek 7:0e9af5986488 283 QuaternionAngles.x = CF.x;
Zbyszek 7:0e9af5986488 284 QuaternionAngles.y = CF.y;
Zbyszek 7:0e9af5986488 285 QuaternionAngles.z = CF.z;
Zbyszek 7:0e9af5986488 286
Zbyszek 7:0e9af5986488 287 //Convert back to radians
Zbyszek 7:0e9af5986488 288 CF.x *= 0.0174533f;
Zbyszek 7:0e9af5986488 289 CF.y *= 0.0174533f;
Zbyszek 7:0e9af5986488 290 CF.z *= 0.0174533f;
Zbyszek 7:0e9af5986488 291
Zbyszek 7:0e9af5986488 292 CF = euler2Quaternion(CF); //Conver CF angle back to quaternion for next iteration
Zbyszek 7:0e9af5986488 293 CF = normaliseQuaternion(CF); //Normalise the quaternion
Zbyszek 7:0e9af5986488 294
Zbyszek 7:0e9af5986488 295 this->t.reset();
Zbyszek 7:0e9af5986488 296 this->t.start();
Zbyszek 7:0e9af5986488 297
Zbyszek 7:0e9af5986488 298 printf("Quat Angle X: %+2f, Quat Angle Y: %+2f, Quat Angle Z: %+2f, dt: %+2f\n\r", QuaternionAngles.x, QuaternionAngles.y, QuaternionAngles.z, deltaT);
Zbyszek 7:0e9af5986488 299 return QuaternionAngles;
Zbyszek 7:0e9af5986488 300
Zbyszek 7:0e9af5986488 301 }