Can my integration code for an accelerometer data work for a Square wave?

1 view (last 30 days)
I have a code written to integrate an accelerometer data to displacement. this code integrates the signal.
The code includes a butterworth high pass filter and the "cumtrapz" integration code.
Applying the code to a square wave, the wave becomes distorted. Does it mean my code is not good enough or the square wave will always be distorted by the filter.
Thanks in advance.

Accepted Answer

Dr. Seis
Dr. Seis on 11 May 2012
Instead of trying to convert using your code, try mine located here:
Say your acceleration data is called "square_acceleration_data", then (once you save my iomega.m file in a location Matlab can find it, like your Matlab home directory) you would produce your displacements like this:
square_displacement_data = iomega(square_acceleration_data,dt,3,1);
where "dt" is equal to your time increment (i.e., 1/Fs or 1 over your sampling frequency).
Example:
  2 Comments
kunle
kunle on 14 May 2012
Thanks Elige and Walter for your responses. Elige, thanks for you iomega.m program but it only works for my square wave. When I apply it to the acceleration data, it gives a 'U' shaped curve which shouldn't be.
Dr. Seis
Dr. Seis on 14 May 2012
Can you post a picture of the acceleration data with the resulting displacement data? It might be necessary to run a high-pass filter on the acceleration data to avoid introducing low-frequency noise into the displacement data.

Sign in to comment.

More Answers (1)

Walter Roberson
Walter Roberson on 11 May 2012
It could mean both, but Yes, a square wave will always be distorted by any frequency filter.
square waves (and all other kinds of waves with straight edges... other than the constant signal of course) require infinite signal bandwidth to represent properly. Remove any frequency from that bandwidth (by way of filtering) and the square wave will be distorted. You often get "overshoot" on the leading and trailing edges as well.

Categories

Find more on Historical Contests in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!