Convert striaght lines to smooth curves

3 次查看（过去 30 天）
Pranav Shende
Pranav Shende 2022-8-22
评论： Pranav Shende 2022-8-27
Hi, I have acceleration and time data and want to create a phase diagram (displacement vs velocity).
My current plot looks like below fig 1, with straight lines connecting each point. Is it possible to change these to smooth curves as fig 2?
fig 1
fig 2
My current code is as below.
v = diff(Acc1).*diff(X_Value); %Acc1 = acceleration data and X_Value = Time
s = diff(v).*diff(X_Value(1:end-1));
figure
plot(s,v(1:end-1))

请先登录，再回答此问题。

采纳的回答

Michael
Michael 2022-8-22
It looks like your sample rate was too low. I'd recollect the data at a higher sample rate.
You could try doing a cubic spline interpolation on your acceleration data to estimate intermediate time points and then rerun your code.
dt = mean(diff(X_value));
t_hiRes = 0:dt/100:max(X_value);
hiresAccel1 = interp1(X_value, Acc1, t_hirRes,'cubic');
%Same for the rest
v = diff(Acc1).*diff(X_Value); %Acc1 = acceleration data and X_Value = Time
s = diff(v).*diff(X_Value(1:end-1));
figure
plot(s,v(1:end-1))
  4 个评论
显示 2更早的评论
Pranav Shende
Pranav Shende 2022-8-26
编辑：Pranav Shende 2022-8-26
Thanks for your help Micheal, I understand using 'diff' comes out confusing to find velocity from acceleration. I have tried using trapz and cumtrapz for my data. The problem I faced is all the values (velocty and displacement) increases/decreases with time. I believe it should follow the same trend as accleration if plotted against time as my sensors are attached to the rig fixed on a table. Ideally as below figure. I think as you mentioned the sampling frequency is too low for my data. I will be using a higher sampling frequency and try again. Or else I am missing something fundamental and need to go back to basics.

请先登录，再进行评论。

更多回答（1 个）

John D'Errico
John D'Errico 2022-8-25
Ran in:
Why in the name of god and little green apples are you differentiating acceleration? You need to INTEGRATE acceleration, once to get velocity, then again to get displacement.
Differentiating acceleration results in noisy garbage, as it probably should. You differentiated the accelerations TWICE.
data = xlsread('starved_1000_001.lvm.xlsx');
T = data(:,1);
Acc = data(:,2);
vel = cumtrapz(T,Acc);
displacement = cumtrapz(T,vel);
plot(vel,displacement,'-')
  3 个评论
显示 1更早的评论
Michael
Michael 2022-8-26
@Pranav Shende This might be due a DC bias in your accel measurement. This can happen if you have a bit of gravity measured in the axis of the accel. Try doing
a = a_raw - mean(a_raw)
to eliminate the bias. Also, be sure to accept the answer.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Descriptive Statistics 的更多信息

标签

产品


版本

R2022a

Community Treasure Hunt

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

Start Hunting!

Translated by