Get separate index from a column with separate conditions

Question

Jorge Luis Paredes Estacio 2023-7-21

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1998968-get-separate-index-from-a-column-with-separate-conditions

评论： Jorge Luis Paredes Estacio 2023-7-21

Hello, I have the following script with the attached files. I extracting two separate indexes from a column to find specific acceleration values values that matches the "index_initial_time" and "index_final_time" in a for loop. However, the second index ("index_final_time") is not capture properly, which does not allow to estimate the mean of the two acceleration values capture in the next lines. Could you please helpe on this. I would appreciate the help.

clear all
clc
SIGNAL = load('Time_Acceleration.txt');
X_spike_matrix=load('X_spike_matrix.txt');
n_spikes_X=height(X_spike_matrix);
TIME=SIGNAL(:,1);
ACCX=SIGNAL(:,2);
DT=TIME(2,1)-TIME(1,1);
ACCX_spikes_removed=ACCX;
for k=1:n_spikes_X;
    initial_time_ACCX=X_spike_matrix(k,1)-DT; %t-DT time to corresponding acceleration value
    final_time_ACCX=X_spike_matrix(k,1)+DT; %t+DT time to corresponding acceleration value
    index_initial_time=find(TIME==initial_time_ACCX);
    index_final_time=find(TIME==final_time_ACCX);
    initial_ACCX=ACCX(index_initial_time,:); %t-DT acceleration to compute average replacement
    %% line with issue
    final_ACCX=ACCX(index_final_time,:); % t+DT acceleration to compute average replacement 
    %%
    mean_ACCX_replacement=mean([initial_ACCX,final_ACCX]); %estimation of mean neibouring  value to replace peak value
    matrix_ACCX_neighbour_mean(k,1)=mean_ACCX_replacement; %acumulating corrected peak values in a matrix
    ACCX_spikes_removed(TIME==X_spike_matrix(k,1))=mean_ACCX_replacement; %replacement of the peak value with the mean nenighbouring value
end

However, line 20 does not capture the next index value

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Matt J 2023-7-21

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1998968-get-separate-index-from-a-column-with-separate-conditions#answer_1276683

编辑：Matt J 2023-7-21

在 MATLAB Online 中打开

Remember, you are in a world of finite precision computers...

     ftol=@(t) abs(TIME-t)<smallnumber; %compare to within a tolerance
    index_initial_time=find( ftol(initial_time_ACCX) , 1);
    index_final_time=find( ftol(final_time_ACCX)  ,1,'last');