在 MATLAB 中以高效的方式动态 window 形成

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（局部最大值）我会检查后面和前面的值是否更大使 window。

示例

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

在 MATLAB 中以高效的方式动态 window 形成

Dynamic window forming in efficient way in MATLAB

sorting

matlab

signal-processing