The Math Forum

Search All of the Math Forum:

Views expressed in these public forums are not endorsed by NCTM or The Math Forum.

Math Forum » Discussions » Software » comp.soft-sys.matlab

Notice: We are no longer accepting new posts, but the forums will continue to be readable.

Topic: numerical integration
Replies: 5   Last Post: Mar 20, 2013 2:43 PM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View   Messages: [ Previous | Next ]
Loren Shure

Posts: 1,067
Registered: 12/7/04
Re: numerical integration
Posted: Mar 20, 2013 8:54 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply

"Nasser M. Abbasi" <> wrote in message
> On 3/18/2013 5:10 AM, JankoZ wrote:
>> Hi I need your help!
>> I have measurement of acceleration in axis X,Y,Z (cca 6000 measurement
>> with spet 0.02sec).
>>I want to get final position with using Numerical Integration in Matlab.
>>Can anybody help me?

> Numerical integration can be hard to get right for this.
> Better to do this in frequency domain using FFT.
> Take the FFT of acceleration data. Since displacement is integration
> 2 times in time domain, then in frequency domain it becomes division
> by -(n*w1)^2 where w1 is the fundamental frequency of the acceleration
> signal (i.e. 2*pi/T) and n is the coefficient (i.e. bin) number.
> i.e. n=0..N-1
> Then you do ifft on the result.
> If you want just velocity, then divide by (I*n*w1) instead and
> do the same.
> This is what I do. It works well. But if you really want to
> do numerical integration, then get the AbraVibe matlab toolbox by
> Anders Brandt, it has a well written timeint() function. Apply
> it 2 times on your data. (for each x,y,z)
> --Nasser

I believe the original data were not evenly sampled, making an FFT (or
filter) solution problematic. Was that just a typo in the data?


Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© The Math Forum at NCTM 1994-2018. All Rights Reserved.