Эффективное формирование динамического окна в MATLAB

Question

Может ли кто-нибудь помочь мне предоставить эффективный способ или помочь мне выполнить код, чтобы сделать те же результаты в минимально возможных шагах. Я буду вам благодарен.

У меня есть оригинальный массив:

A = [1 1 1 4.3 4.5 4 4.3 3 1 0 0 2 6.2 6.3 6 6.2 7.4 8 7.2 2 2 3 3 2];

Выход выглядит как:

A = [1 1 1 4 4 4 4 3 1 0 0 2 6 6 6 6 6 7 7 2 2 3 3 2];

Я применил некоторые ограничения и удалил некоторые значения из массива локальных максимумов, после чего получил несколько новых массивов.

Yposlocfiltered = [6    15    23];
idx             = [4     6     3];
Yneglocfiltered = [2     9    20];
idx_neg         = [1     1     2];

Где я найду idx(локальное значение максимума), я проверим, составляют ли значения позади и впереди окно.

пример

If I will find 4 and 4.5, 4.3 is greater than 4 include it in backward window 
4.3 is greater than 4 include it in forward window.

Мне нужно найти диапазон значений за локальными максимумами и впереди локальных максимумов.

Я пытался написать код, который работает нормально, но проблема в том, что это слишком долго.

Может кто-нибудь предложить мне идею выполнить это действие в минимальные шаги и более быстрыми способами?

Я предоставил код только для положительных локальных максимумов, а для отрицательного кода локальных максимумов. Это просто копия этого.

Код: только для положительных локальных максимумов

 clc
 clear all

 A = [1 0 1 4.3 4.5 5 4.3 3 0 0 0 2 6.2 6.3 7 6.2 7.4 8  7.2 1 2 3 4 2];
 Yposlocfiltered = [ 6    15    23];
 idx = [4     6     3];
 Yneglocfiltered = [2     9    20];
 idx_neg = [1     1     2];

 for b = 1: numel(idx)
 for c = 1:numel(A)
 f = Yposlocfiltered(1,b)-c ;
 g = Yposlocfiltered(1,b)+c ;

 %   if (f > 0 && g <= numel(X))
 if(f > 0)   
 if  (A(Yposlocfiltered(1,b)-c))> idx(1,b)

 else 
   d= f+1;
   z(b)= d;
   backward_window = z;
  break
 end        

 end

 end 

 end 

 for r = 1: numel(idx)
 for s = 1:numel(A)
 u = Yposlocfiltered(1,r)-s ;
 v = Yposlocfiltered(1,r)+s ;

 %   if (f > 0 && g <= numel(X))
 if(v <=numel(A))   
 if  (A(Yposlocfiltered(1,r)+s))> idx(1,r)

 else 
   w= v-1;
   y(r)= w;
   forward_window = y;
 break
 end        

 end

 end 

end 

n=4
for i=1:length(backward_window)
range = backward_window(i): forward_window(i);
p = range 
if n <= numel(p)
        p = range(1:n)
         A( p) = idx(i);
else
%   if (size(range)<= 3)
A( p) = idx(i);
end 
end

Answer 1

С первого взгляда на ваш код, я считаю, что вы можете объединить первые два цикла for в один.

sA = numel(A);
sI = numel(idx);

for i = 1:sI
    f = Yposlocfiltered(i) - [1:sA];
    g = Yposlocfiltered(i) + [1:sA];
    f(f<1) = [];
    g(g>sA) = [];
    backward_window(i) = f(find(A(f) <= idx(i), 1)) + 1;
    forward_window(i) = g(find(A(g) <= idx(i), 1)) - 1;
end

Здесь вы можете использовать find найти элемент массива, соответствующий указанному условию, т.е. g <= numel(X) или же A(f) <= idx(i),

Ваш последний цикл, который изменяет A также может быть интегрирован в тот же цикл, так что вы можете иметь:

sA = numel(A);
sI = numel(idx);

n=4;
for i = 1:sI
    f = Yposlocfiltered(i) - [1:sA];
    g = Yposlocfiltered(i) + [1:sA];
    f(f<1) = [];
    g(g>sA) = [];
    backward_window(i) = f(find(A(f) <= idx(i), 1)) + 1;
    forward_window(i) = g(find(A(g) <= idx(i), 1)) - 1;
    range = backward_window(i) : forward_window(i);
    if n <= numel(range)
        A(range(1:n)) = idx(i);
    else
        A(range) = idx(i);
    end
end