Effacer les filtres
Effacer les filtres

Trying to change plot style

2 vues (au cours des 30 derniers jours)
Nathaniel Porter
Nathaniel Porter le 19 Déc 2021
Commenté : Star Strider le 19 Déc 2021
%Plot of entire patient one time series
load ('glucose.mat')
glucose_mgdl = glucose * 18;
plot(date1+time,glucose_mgdl),
xlabel('Days'), ylabel('Glucose')
title('Glucose readings vs days/time')
M = mean(glucose_mgdl)
M = 183.4936
S = std(glucose_mgdl)
S = 87.6448
Trying to display correlogram as the picture seen below what do I have to add to my code:
%Plot of correlogram
[r,lags] = xcorr(glucose_mgdl,'coeff');
figure
plot(lags, r, '.')
grid
xlabel('Lags')
ylabel('Croorelation Coefficient')
title('Autocorrelation')

Réponse acceptée

Star Strider
Star Strider le 19 Déc 2021
Trying to display correlogram as the picture seen below what do I have to add to my code:
The plot in the documentation is a stem plot. I plotted that originally, however the data are so densely packed that it simply appears as a filled area below the curve (even when using '.' as the marker).
One option would be to select a subset of values and plot them —
%Plot of entire patient one time series
T1 = readtable('https://www.mathworks.com/matlabcentral/answers/uploaded_files/838945/glucose.csv');
glucose = T1.glucose;
date_time = T1.date + T1.time;
glucose_mgdl = glucose * 18;
plot(date_time,glucose_mgdl),
xlabel('Days'), ylabel('Glucose')
title('Glucose readings vs days/time')
M = mean(glucose_mgdl)
M = 183.4936
S = std(glucose_mgdl)
S = 87.6448
%Plot of correlogram
[r,lags] = xcorr(glucose_mgdl,'coeff');
sampstep = 10; % Sampling Step
rs = r(1:sampstep:end);
lagss = lags(1:sampstep:end);
figure
stem(lagss, rs, '.')
grid
xlabel('Lags')
ylabel('Crorelation Coefficient')
title('Autocorrelation')
Change the value of ‘sampstep’ to get different stem plot results.
.
  2 commentaires
Nathaniel Porter
Nathaniel Porter le 19 Déc 2021
Oh okay thank you
Star Strider
Star Strider le 19 Déc 2021
As always, my pleasure!
.

Connectez-vous pour commenter.

Plus de réponses (0)

Produits


Version

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by