I am very new to DSP so please forgive any ignorance within this post.
I am working on a device which has a 6-axis IMU measuring acceleration and rotation. The device itself is attached to a physical bracket and is calibrated on a desk. During calibration, the system waits for a state of motionless, which is defined by a change in 1.5 degrees in any axes, before proceeding.
Because of the bracket, I've noticed it is susceptible to mechanical vibration. For example, when resting on a desk and people are walking by, the device can sometimes report movement in 2-5 degrees in some axes.
I've thought about relaxing that threshold from 1.5 deg to something higher, but thinking about it this past weekend I realized I don't want to do that. If the system is calibrating, it shouldn't be experiencing any notable motion during this step, or that won't truly be 'calibrated'. That would lead to errors for anything relying on that calibrated state.
I've thought about running some accelerometer readings (as in the lsb's) through an FFT to see if there are specific frequencies that appear when people walk by that I can filter out. I talked myself out of this since I think the vibrations itself would appear random due to differences in body weight, height, and stride of a human being.
The state of motionless mentioned above is determined by the difference between the current orientation and previous orientation for each axis. So, it's a old/new angle comparison really. The algorithm used is a median filter.
1) The new angle is the output of a median filter with window size of three. A new reading is rolled into the median filter at 100Hz (10ms), then the filter is sorted and the middle value is taken. The actual IMU samplings occur much quicker, but for this step (every 10ms) we using take the latest sample. I believe this is essentially down-sampling?
2) The old angle is simply the previous 'new' angle.
3) If new - old is greater than 1.5 degrees, the system considers itself to be under motion (rightfully so in my opinion).
Is there any way to account or compensate for this mechanical vibration in software? I am pretty lost as to where to go forward from this point. Any help would be appreciated.
EDIT I should probably mention the IMU is mounted rather well. The issue lies with the overall device, not the mounting of the IMU to the device PCB.
As requested, I've gotten around to collect some data surrounding my problem. Sorry this took so long, other priorities were shifted.
I collected two sets of data:
1) Natural vibration of the system during calibration
- pastebin.com/bEN0Fc0r (Raw) - pastebin.com/kyi5eET3 (Median Filter)
2) Vibration observed while people were walking nearby
- pastebin.com/XvahfKP3 (Raw) - pastebin.com/kJDPWrg0 (Median Filter)
Raw accelerometer data was sampled at 100Hz for about 12 seconds. I also saved the median-filter result (sampled at same frequency and period). I wrote a python script (dependent on scipy & numpy) to process both sets of logs. What's shown below is as follows:
1) Raw 3-axis accelerometer data (in g's)
2) Median-filter output of 3-axis accelerometer data (in g's)
3) FFT of raw 3-axis accelerometer data (in g's)
4) New/Old Angle Deltas (in degrees)
- Any measurement above 1.5 degrees would be considered to be in motion.
The first set of figures was collected when the physical bracket was set on the desk and allowed to naturally vibrate. You can see that it eventually settles but it does take several seconds. This was without anyone moving nearby.
The second set of figures was collected while myself and 2 other people walked around nearby.
* In the people-walking-nearby set, the same frequency buckets are shown but with a much smaller amplitude. Unfortunately it's still picked up by the system as shown in the deltas in the next diagram.
I'd like to add that once a state of motionless is achieved (<1.5deg in all axes), the system waits one full second before calibrating. If a state of motion is entered (>1.5 deg) during this period, the system must reenter a motionless state, and then another full second waiting period occurs. In the diagrams shown above, the calibration takes place right around the 12s mark.