digital integration
digital integration
(OP)
Hello,
I am currently working on a project where I need to find a digital integration for acceleration data. I believe I have found a sufficient integration filter that I have tested in matlab but after the first integration there is a dc offset that makes the second integration to get displacement data impossible. It would be ideal to find another type of filter to remove the dc offset for my project. I have tried basic high pass filters with no luck. Is it possible to do such a thing? and if so could someone advise me on a sufficient high pass filter?
I am currently working on a project where I need to find a digital integration for acceleration data. I believe I have found a sufficient integration filter that I have tested in matlab but after the first integration there is a dc offset that makes the second integration to get displacement data impossible. It would be ideal to find another type of filter to remove the dc offset for my project. I have tried basic high pass filters with no luck. Is it possible to do such a thing? and if so could someone advise me on a sufficient high pass filter?





RE: digital integration
TTFN
FAQ731-376: Eng-Tips.com Forum Policies
RE: digital integration
Dan - Owner

http://www.Hi-TecDesigns.com
RE: digital integration
Tenngineer, I would suggest that you look into methods of numerical integration. From what little you have describe about what you are trying to accomplish, it looks to me as though you are making things too complex by worrying about filtering data. A common form of numerical integration would be Runge Kata (spelling?) and you may consider searching for that term. The book "Math toolkit for real time programming" by Jack Crenshaw is a wonderful resource for a solid how to when it comes to math programming for real time systems.
RE: digital integration
RE: digital integration
Gunnar Englund
www.gke.org
--------------------------------------
100 % recycled posting: Electrons, ideas, finger-tips have been used over and over again...
RE: digital integration
Before you choose an integration routine you need to know if your sample rate meets Nyquist. What is your noise density model? Is your answer within the realm of expectations?
At the risk of copy write infringement I will copy only two lines from his book that I have used in industry.
First you will need both the present acceleration measurement and the past acceleration measurement. Then you produce a predicted velocity. Time now is "I", the past time isi-1, the future is i+1. The time delta is "h", standing for height. But buy the book if you want to know more.
If you want more explanation; Buy the book. Even if you are mildly interested buy the book. It is a great resource for understanding. "Math Toolkit for Real-Time Programming" Jack W. Crenshaw ISBN 1-929629-09-5. I have one from CMP Books, but there is a new publisher according to Jack.
v(i+1) = v(i) + (h/2)*(3*f(i)-f(i-1))
The position is
x(i+1) = x(i) + (h/2)*(v(i+1)+v(i))
Without any other information this is a great thing to do. It is a predictor corrector method, Adam, Bashforth, and Moulton are keywords. Numerical Recipes for C++ recommends PC methods for most problems. There are caveats.
Is your sample interval exactly uniform? If not you could have some jitter problems. So go with a lower order integrator with very small steps. Jack Crenshaw explains how to make a sample rate independent system. I did just that. But it is a lot of work, especially the startup problem. If you do not know what the startup problem is, buy the book!
I am not being facetious. Integration problems are not trivial. I would need to know more to help. I can not help with accelerometers. Can others pipe in?
jsolar