Displacement from acceleration measurements

3 次查看(过去 30 天)
I have the attached acceleretion measurements and I have used the following code to get the displacement. I suspected to get the displacement curve having multiple cycles, why I got only one cycle (It must contain multiple cycles)?
Any suggestions?
data = dlmread("test4.txt"); x = data(:, 1);
acc1 = detrend(x,'linear');
Fs = 1000 ; % Sampling frequency Ts = 1/Fs ; % Sampling period L = length (acc1) ; % Length of signal t = 0:Ts:(Ts*L)-Ts; dt = mean(diff(t)); % Average dt %fs = 1/dt; % Frequency [Hz] or sampling rate
% some additionnal high pass filtering %N = 4; %fc = 0.05; % Hz %[B,A] = butter(N,2*fc/Fs,'high'); %%acc2 = filter(B,A,acc); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
velocity1 = cumtrapz(dt,acc1); velocity1 = detrend(velocity1,'linear'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% disp1 = cumtrapz(dt,velocity1);
figure(1) subplot(311),plot(t,acc1); subplot(312),plot(t,velocity1); subplot(313),plot(t,disp1);

回答(1 个)

Star Strider
Star Strider 2023-3-25
The slopes are due to a constant that is being integrated can be removed by subtracting the mean of the variable to be integrated (thus removing any trends due to the integration), and then calling cumtrapz to do the integration.
No further signal processing is necessary —
A = readmatrix('https://www.mathworks.com/matlabcentral/answers/uploaded_files/1334739/test4.txt');
A = fillmissing(A,'linear');
Fs = 1000;
L = size(A,1);
t = linspace(0, L-1, L).'/Fs;
V = cumtrapz(t, A-mean(A));
D = cumtrapz(t, V-mean(V));
figure
tiledlayout(size(A,2),3)
for k = 1:size(A,2)
nexttile
plot(t, A(:,k))
grid
ylim([-0.5 2.5])
ylabel('Acceleration')
xlabel('Time')
title("Column "+k)
nexttile
plot(t, V(:,k))
grid
ylim([-0.02 0.02])
ylabel('Velocity')
xlabel('Time')
nexttile
plot(t, D(:,k))
grid
ylim([-10 30]*1E-4)
ylabel('Displacement')
xlabel('Time')
end
I used ylim to force all the y-axis ranges to be the same, relative to each integration. This may make them a bit easier to interpret.
.

类别

Help CenterFile Exchange 中查找有关 Biological and Health Sciences 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by