After brute force trial and error, I think this is working. Can someone advise on the units?
A 2hz sine wave with an amplitude of 2:
y=2*sin(2 * 2 * pi * t)
is producing a velocity wave with amplitude 0.1625.
Here is the updated code for anyone curious to try it out.
clear
home...
So I have the frequency axis working ok after taking an FFT, without using fftshift. But, I'm still stuck. I think the I need to use the entire frequency vector (freq3) in the omega arithmetic : G=Y./(2*pi*freq3*1i);
Not just the first half. In the code below I tried to make the 2nd half...
Thanks Pete. You're right. I'm too used to working with the data after using "fftshift". From the Matlab website:
I will investigate more tomorrow when I have the code in front of me, but I have a good feeling about the code working once I define the freq vector correctly for when fftshift...
If you don't have time to read the whole blog, figure 12 and figure 13 summarize the difference in results between Omega Arithmetic Method and filter->integrate-> filter again.
It looks to me like there isn't much difference, especially if you were to manually remove the DC offset at the end...
Thanks for the reply. The reasons I would like to use the frequency domain are outlined here:
http://blog.prosig.com/2010/12/16/methods-of-conversion-between-acceleration-velocity-and-displacement/
So far, I've used the FFT-> modify the freq domain signal-> IFFT back to time domain...
No. The application is engine vibration on a v6 engine. My company is the supplier of an accessory drive. I've taken a lot of accelerometer data "for fun" that's not related to our parts, and I'd like to do something with it.
They didn't teach anything close to this where I went to school..
Hello, this is my first post.
I've searched and seen some similar quesitons in the past, but thus far I'm still stuck, so I'll ask.
I'm trying to convert acceleration results into velocity using
"Omega Arithmetic". The Matlab code I wrote is something like this...
