基本的MATLAB变量操作

Question

我有一个dwc = [3001 x 2 double] ，这或多或少是一个sin函数，我有一个for循环， for在dwc(:,2)查找顶部值。 可以说dwc(531,2)中有一个最大值，这是最好的方法，或者是什么很容易采用dwc(531,1)和dwc(531,2)并使M = [num_of_top_points x 2 double] ？

对于以下循环，我该怎么办？

j = 0;
for i = 2:size(dwcL01,1)-1
    if dwcL01(i,2) > dwcL01(i-1,2) && dwcL01(i,2) > dwcL01(i+1,2)
        j = j+1;
        ?? = dwcL01(i,:);
    end
end

Answer 1

这是您完成循环的方式

j = 0;
M = [];
for i = 2:size(dwcL01,1)-1
    if dwcL01(i,2) > dwcL01(i-1,2) && dwcL01(i,2) > dwcL01(i+1,2)
        j = j+1;
        M(j, :) = dwcL01(i, :);
    end
end

但是您可以通过向量化更有效地执行此操作

%//Some example data
x = -4*pi:0.5:4*pi;
y = cos(x);
dwcL01 = [x(:), y(:)]; %// (:) just makes it a column

%// Finding the peaks using diff and sign. Note that I add the first element to the beginning as diff reduces the size by one so this prevents offsetting
F = diff(sign(diff([dwcL01(1,2);dwcL01(:,2)]))) < 0;
M = [dwcL01(F,:)', dwcL01(F,:)'];    

plot(x, y, M(:,1), M(:,2), '*r')

首先，我们发现每个元素连续元素对的区别。 现在，当符号改变时，这意味着我们达到了最大值或最小值。 如果符号变化为负，则梯度从正变为负，这是最大值。 因此，我使用diff(sign())来找到符号变化的点，然后使用> 0来创建一个逻辑矩阵，在所有期望的最大值处都为false。 然后，我使用逻辑索引来提取最大值。

Answer 2

您可以将其附加到矩阵中（我们称其为dwcL01_max ）-这不是最快的方法，因为矩阵大小会改变每个循环，但它可以工作：

dwcL01_max = [dwcL01_max dwcL01(i,:)];

另一个选择是使用内置的findpeaks （来自signal proc工具箱）

[~, dwcL01_peaks] = findpeaks(dwcL01(:,2));
dwcL01_max = dwcL01(dwcL01_peaks, :);