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

1 次查看（过去 30 天）
Matthew Worker
Matthew Worker 2023-7-9
编辑： John Kelly 2023-8-2
Hello can you help with the above question
  9 个评论
显示 7更早的评论

请先登录，再进行评论。

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

采纳的回答

Cris LaPierre
Cris LaPierre 2023-7-10
Ran in:
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).';
load acc.mat
yddot = acc;
ydot = cumtrapz(tspan,yddot);
ydot2 = detrend(ydot,10,"SamplePoints",tspan);
y = cumtrapz(tspan,ydot2);
plot(tspan,ydot2,tspan,y)
  2 个评论
Cris LaPierre
Cris LaPierre 2023-7-10
编辑：Cris LaPierre 2023-7-10
Ran in:
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).';
load acc.mat
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
Cris LaPierre 2023-7-11
编辑：John Kelly 2023-8-2
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

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Measurements and Spatial Audio 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by