Complementary filter and relative orientation with MPU9250

Dupont cables

MPU9250

Arduino Nano

Mini Breadboard

Arduino IDE

Main article: <a href="https://www.hibit.dev/posts/67/complementary-filter-and-relative-orientation-with-mpu9250?utm_source=arduino&utm_medium=social&utm_campaign=project_hub">Complementary filter and relative orientation with MPU9250</a>MPU-9250 is one of the most advanced combined accelerometer, gyroscope and compass small size sensors currently available. It replaces the popular MPU-9150 lowering the power consumption, improving gyro noise and compass full scale range performance. It has many advanced features, including low pass filtering, motion detection and even a programmable specialized processor.Internally it includes the MPU-6500, which contains a 3-axis gyroscope plus a 3-axis accelerometer, and the AK8963, the market leading 3-axis digital compass. The MPU-9250 uses 16-bit analog-to-digital converters (ADCs) for digitizing all 9 axes.The MPU-9250 is a high performance accelerometer, gyroscope and magnetometer option for use in applications such as gesture recognition, self-balancing robots, cell phones, fitness monitoring and similar applications where detection of movement direction and magnitude along with rotation is desired. In <a href="https://www.hibit.dev/posts/36/reading-mpu9250-sensors-with-arduino">reading MPU9250 sensors with Arduino </a>article we've introduced the module and how to interact with it, now we will see how to use extracted data and get relative orientation.<h3 >Wiring schema</h3>Connect the power supply pins and SDA, SCL pins according to the connection diagram shown below.<figure><img src=https://projects.arduinocontent.cc/87bff5fd-3caf-4259-b1c6-ff5049b82b60.png /><figcaption>Arduino Nano wiring with MPU9250</figcaption></figure>Note: the module operates with 3.3 volts, although some versions have a regulator that allows to connect it to 5V. Make sure to check your module specification before using it.<h3 >Relative orientation</h3>Relative orientation is the recovery of the position and orientation of one object relative to another. Depending on the direction there are three types of angular rate measurements:<ul ><li >Yaw: the horizontal rotation on a flat surface when seen the object from above.</li></ul><ul ><li >Pitch: vertical rotation as seen the object from front.</li></ul><ul ><li >Roll: horizontal rotation when seen the object from front.</li></ul>A graphical representation for each axis using a plane:<figure><img src=https://projects.arduinocontent.cc/630e30c2-d301-41f2-bd3b-0bebde0fe858.png /><figcaption> </figcaption></figure>Angles can not be directly extracted from MPU-9250 module but with correct formulas we can obtain them.<h3 >Accelerometer angles</h3>Commonly, g is used as a unit for acceleration, relative to standard gravity (1g = 9.80665m/s2). Acceleration for perfectly aligned IMU will be 1g for Z axis. Terrestrial gravity and trigonometry allow to use the accelerometer readings and know the angle of inclination. The formulas to calculate the angles are:<blockquote >pitch = atan2(accelerometer_x, sqrt(accelerometer_y^2 + accelerometer_z^2))</blockquote><blockquote >roll = atan2(accelerometer_y, sqrt(accelerometer_x^2 + accelerometer_z^2))</blockquote>You can not obtain yaw angles from accelerometer because yaw motions occur in a plane and so no matter what X and Y readings you have.<h3 >Gyroscope angles</h3>The gyroscope measures the angular velocity, i.e., the number of degrees rotated in a unit of time. Having initial angle of the IMU, angular velocity and the time passed between samples, we can know the angle change. It's very similar to linear speed, that multiplied by time give us the total distance passed. The formulas to calculate the angles are:<blockquote >pitch = pitch + gyroscope_x * elapsed_time</blockquote><blockquote >roll = roll + gyroscope_y * elapsed_time</blockquote><blockquote >yaw = yaw + gyroscope_z * elapsed_time</blockquote>Where elapsed time is the time interval between samples. The Z axis is usually ignored, since we don't have yaw angle from accelerometer and we can not rely only on gyroscope sensor.<h3 >Noise and errors</h3>In the ideal word, sensors data with applied formulas would provide us precise and exact angles. The reality is different as some factors affect the sensors output.Typically, when you move around with an accelerometer it experiences movement accelerations. Additional acceleration values may affect the orientation accuracy. Another accelerometer related problem is the noise: unwanted disturbance in an electrical signal. The accelerometer is able to measure any angle, however its readings are noisy and have a certain margin of error even with a low pass filter.On the other hand, gyroscopes are subject to bias instabilities, in which the initial zero reading of the gyroscope will cause drift over time due to integration of inherent imperfections and noise within the device.So, what can we do? There are different algorithms to solve this problems. The one we are going to use is known as complementary filter. Idea behind complementary filter is to take slow moving signals from accelerometer and fast moving signals from a gyroscope and combine them. It is ideal to implement with Arduino: easy to use, low cost of processing and good precision.<h3 >Complementary filters</h3>The complementary filter can be thought of as a union of two different filters: a high-pass filter for the gyroscope and a low-pass filter for the accelerometer. The first lets only pass the values ​​above a certain limit, unlike the low-pass filter, which only allows those below.Accelerometer gives a good indicator of orientation in static conditions and gyroscope gives a good indicator of tilt in dynamic conditions. The formula resulting from combining the two filters is:<blockquote >angle = (1 - α) * (angle + gyroscope * dt) + α * accelerometer</blockquote>A common value for α is 0.98, which means that 98% of the weight lays on the gyroscope measurements.<blockquote >pitch = 0.98 * (pitch + gyroscope_x * dt) + 0.02* accelerometer_x</blockquote><blockquote >roll = 0.98 * (roll + gyroscope_y * dt) + 0.02* accelerometer_y</blockquote>As the data changes very rapidly we can sample for some amount of time and take the average for more precise results.<h3 >Install Arduino library for I2C</h3>We will use I2C protocol to interact with registers and read/write data. The I2C library provides a very simple interface for that purpose and can be reused in other projects that use I2C protocol. It can be downloaded from our <a href="https://github.com/hibit-dev/mpu9250/raw/master/lib/I2C.zip">official repository.</a>To import a library, open the Arduino IDE, go to Sketch > Include Library > Add.ZIP Library and select the library file downloaded from our GitHub repository.<figure><img src=https://projects.arduinocontent.cc/616277f7-f520-4af9-8350-56752f5a770c.png /><figcaption> </figcaption></figure>Then you can simply use include statement:<pre ><code >#include "I2C.h" </code></pre>It will include the library with predefined functions to interact with registers.<h3 >Arduino code</h3>We are going to reuse most of the code generated from the raw sensors reading. To be more specific, everything related with data reading and normalization, together with defined structs.<blockquote >The sampling theorem states that a real signal which is band-limited to f HZ can be reconstructed without error from samples taken uniformly at a rate R > 2f. So, it's important to choose an output rate that will be at least two times the DLPF, otherwise the data will be corrupted (aliased).</blockquote>A new structure has been defined, called angle, to represent object position. Make sure to follow the previous step and import I2C library in order to use it with the include statement:<pre ><code >#include "Wire.h" #include "I2C.h" #define MPU9250_IMU_ADDRESS 0x68 #define MPU9250_MAG_ADDRESS 0x0C #define GYRO_FULL_SCALE_250_DPS 0x00 #define GYRO_FULL_SCALE_500_DPS 0x08 #define GYRO_FULL_SCALE_1000_DPS 0x10 #define GYRO_FULL_SCALE_2000_DPS 0x18 #define ACC_FULL_SCALE_2G 0x00 #define ACC_FULL_SCALE_4G 0x08 #define ACC_FULL_SCALE_8G 0x10 #define ACC_FULL_SCALE_16G 0x18 #define TEMPERATURE_OFFSET 21 // As defined in documentation #define INTERVAL_MS_PRINT 1000 #define G 9.80665 struct gyroscope_raw { int16_t x, y, z; } gyroscope; struct accelerometer_raw { int16_t x, y, z; } accelerometer; struct magnetometer_raw { int16_t x, y, z; struct { int8_t x, y, z; } adjustment; } magnetometer; struct temperature_raw { int16_t value; } temperature; struct { struct { float x, y, z; } accelerometer, gyroscope, magnetometer; float temperature; } normalized; struct angle { float x, y, z = 0; }; angle position; unsigned long lastPrintMillis = 0; unsigned long lastSampleMicros = 0; void setup() { Wire.begin(); Serial.begin(115200); // Data Output Rate = 1000 / (1 + SRD) // Output rate must be at least 2x DLPF rate I2CwriteByte(MPU9250_IMU_ADDRESS, 25, 0x01); // Set the SRD to 1 // Set the DLPF to 184HZ by default I2CwriteByte(MPU9250_IMU_ADDRESS, 26, 0x01); // Configure gyroscope and accelerometer scale I2CwriteByte(MPU9250_IMU_ADDRESS, 27, GYRO_FULL_SCALE_1000_DPS); I2CwriteByte(MPU9250_IMU_ADDRESS, 28, ACC_FULL_SCALE_2G); // Enable interrupt pin for raw data I2CwriteByte(MPU9250_IMU_ADDRESS, 56, 0x01); } void loop() { unsigned long currentMillis = millis(); readSample(); if (currentMillis - lastPrintMillis > INTERVAL_MS_PRINT) { Serial.print("TEMP:\\\\t"); Serial.print(normalized.temperature, 2); Serial.print("\\\\xC2\\\\xB0"); //Print degree symbol Serial.print("C"); Serial.println(); Serial.print("Pitch:\\\\t"); Serial.print(getPitch()); Serial.print("\\\\xC2\\\\xB0"); //Print degree symbol Serial.println(); Serial.print("Roll:\\\\t"); Serial.print(getRoll()); Serial.print("\\\\xC2\\\\xB0"); //Print degree symbol Serial.println(); Serial.println(); lastPrintMillis = currentMillis; } } bool readSample() { if (isImuReady() == false) { return false; } unsigned long sampleMicros = (lastSampleMicros > 0) ? micros() - lastSampleMicros : 0; lastSampleMicros = micros(); readRawImu(); normalize(gyroscope); normalize(accelerometer); normalize(temperature); angle accelerometer = calculateAccelerometerAngles(); angle gyroscope = calculateGyroscopeAngles(sampleMicros); detectPitch(gyroscope, accelerometer); detectRoll(gyroscope, accelerometer); return true; } </code></pre>Note: the snippet is part of Arduino project located in our GitHub <a href="https://github.com/hibit-dev/mpu9250/tree/master/src/complementary_filter">repository </a>with the code separated in different logical files.<h3 >Testing</h3>The serial monitor will print latest available sensors data every INTERVAL_MS_PRINT milliseconds (in the example above, once per second) and it should be similar to:<figure><img src=https://projects.arduinocontent.cc/320bd60d-1d75-4b21-a46d-c09aedb8473e.png /><figcaption> </figcaption></figure>Usually the orientation is drawn on the physical module so you can easily detect pitch and roll.

Complementary filter and relative orientation with MPU9250

Getting relative orientation using MPU9250 and Arduino Nano