# Convert striaght lines to smooth curves

2 vues (au cours des 30 derniers jours)
Pranav Shende le 22 Août 2022
Commenté : Pranav Shende le 27 Août 2022
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))
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Réponse acceptée

Michael le 22 Août 2022
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 commentairesAfficher 2 commentaires plus anciensMasquer 2 commentaires plus anciens
Michael le 25 Août 2022
@Pranav Shende Did this fix your problem?
Pranav Shende le 26 Août 2022
Modifié(e) : Pranav Shende le 26 Août 2022
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.

Connectez-vous pour commenter.

### Plus de réponses (1)

John D'Errico le 25 Août 2022
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.
T = data(:,1);
Acc = data(:,2);
vel = cumtrapz(T,Acc);
displacement = cumtrapz(T,vel);
plot(vel,displacement,'-')
##### 3 commentairesAfficher 1 commentaire plus ancienMasquer 1 commentaire plus ancien
Michael le 26 Août 2022
@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.
Pranav Shende le 27 Août 2022
@Michael Sure, will try that too. Thanks.

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Spline Postprocessing dans Help Center et File Exchange

R2022a

### Community Treasure Hunt

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

Start Hunting!

Translated by