# please, how can i make the plot on the x axis from 15 to 25 hertz to rest on the x axis. that is to make them start from zero. i have used detrend but, it is not working

2 vues (au cours des 30 derniers jours)
Matthew Worker le 9 Juil 2023
Modifié(e) : John Kelly le 2 Août 2023
Hello can you help with the above question
##### 9 commentairesAfficher 7 commentaires plus anciensMasquer 7 commentaires plus anciens
Cris LaPierre le 18 Juil 2023
We are several steps removed from the original question. I suggest asking a new question in the forum.
Matthew Worker le 24 Juil 2023
some of the data i dropped are from an ongoing research, which i should not have dropped

Connectez-vous pour commenter.

### Réponse acceptée

Cris LaPierre le 10 Juil 2023
This is likely due to the resolution of your data. The signal is not able to capture everything, so over time, the error accumulates. There are likely techniques you can use when collecting the data to minimize this. However, since you already have the signal, I think you should try higher-order detrending. Here I picked 10 based solely on how the ouput looks
fs= 1000;
T=7;
dt= 1/fs;
tspan= (0:dt:T).';
yddot = acc;
ydot = cumtrapz(tspan,yddot);
ydot2 = detrend(ydot,10,"SamplePoints",tspan);
y = cumtrapz(tspan,ydot2);
plot(tspan,ydot2,tspan,y)
##### 2 commentairesAfficher AucuneMasquer Aucune
Cris LaPierre le 10 Juil 2023
Modifié(e) : Cris LaPierre le 10 Juil 2023
Another approach might be to smooth your signal and then subtract the smoothed signal from the raw signal.
fs= 1000;
T=7;
dt= 1/fs;
tspan= (0:dt:T).';
yddot = acc;
ydot_raw = cumtrapz(tspan,yddot);
smoothedData = smoothdata(ydot_raw,"lowess","SmoothingFactor",0.125);
ydot = ydot_raw - smoothedData;
y = cumtrapz(tspan,ydot);
plot(tspan,ydot,tspan,y)
Cris LaPierre le 11 Juil 2023
Modifié(e) : John Kelly le 2 Août 2023
Post that was deleted
Thank you very much for your response. what of the other question i asked? did you figured it out also?
______________________________________________________________________
Response
It seems you already have the code written for that, but you are only capturing the final result from your nested for loops.
You can capture it easily enough. You just need to make Energy a matrix instead of a vector.
Change
for c = 1.5:1:5.5
...
Energy(ii) = (trapz(tspan,c*(zdot).^2));
end
to something like
dcoeffs = 1.5:1:5.5;
for jj = 1:length(dcoeffs)
c = dcoeffs(jj);
...
Energy(ii,jj) = (trapz(tspan,c*(zdot).^2));
end

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Programming dans Help Center et File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by