Converting acceleration to displacement
4 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Hi,
I have used your code "function dataout = iomega(datain, dt, datain_type, dataout_type)" for conversion of seismic acceleration to displacement. The displacement seems dis-aligned although velocity looks alright. Can you please help.Thanks.
5 comentarios
Dr. Seis
el 1 de Feb. de 2012
Probably should be, the other post has the code that is being referred to in the question. However, the person asking the question here is different from the one who posted the original question in the above link.
Ross Taylor
el 8 de Mzo. de 2013
Hello Elige,
I am a high school student doing a university course in Year 12. I have never used MATLAB till 4 weeks ago. I have chosen a problem and seemed to have bit more than I can handle.
I have accelerometer data which has given acceleration, velocity and the frequency for each data. The logger has recorded this every 30 seconds. Considering I have measured for a few days, I now have 40000 data points. I need this integrated to displacment. But looking at some of these codes, I do not understand how it incorporates the frequency of each wave. For example in the code written below:
function disp_time_data = acc2disp(acc_time_data,dt) % acc_time_data should be acceleration amplitudes evenly spaced in time % dt in units of seconds per sample N1 = length(acc_time_data); N = 2^nextpow2(N1); if N > N1 acc_time_data(N1+1:N) = 0; % pad array with 0's end df = 1 / (N*dt); % frequency increment Nyq = 1 / (2*dt); % Nyquist frequency acc_freq_data = fftshift(fft(acc_time_data)); disp_freq_data = zeros(size(acc_freq_data)); f = -Nyq:df:Nyq-df; % I think this is how fftshift organizes it for i = 1 : N if f(i) ~= 0 disp_freq_data(i) = acc_freq_data(i)/(2*pi*f(i)*sqrt(-1))^2; else disp_freq_data(i) = 0; end end disp_time_data = ifft(ifftshift(disp_freq_data)); disp_time_data = disp_time_data(1:N1); return
If for example, two columns of my data were the following, how can I know the above code is giving displacement values with some accuracy?
Thank you so much for looking into my problem in advance.
Velocity Frequency Acceleration X- mm/s Y-mm/s Z- mm/s X [Hz] Y [Hz] Z [Hz] X [m/s^2] Y [m/s^2] Z [m/s^2] 0.058 0.07 0.064 13.7 102.4 93.1 0.023 0.03 0.022 0.053 0.06 0.056 51.2 93.1 33 0.023 0.029 0.02
Respuestas (1)
Dr. Seis
el 1 de Feb. de 2012
You might need to filter your data. One of the results of going from acceleration to displacement is that you will boost low-frequency noise. I would suggest running a high-pass filter on this type of output from this program and see if that seems to clean things up a bit. Another thing that might help is to pad your data by adding zeros to both the beginning and end (instead of just the beginning as I do in that previous post) and then apply a taper (so that there isn't a harsh transition between the zeros and your data).
0 comentarios
Ver también
Categorías
Más información sobre Axis Labels en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!