在增加的範圍內查找連續點

Question

我在向量中有一組數據點。例如，

[NaN, NaN, NaN, -1.5363, NaN -1.7664, -1.7475]; 

這些數據來自選擇指定範圍內的3個點（具體爲-0.6和0.6）的代碼。如果列中的三個點不存在於此範圍內，則範圍將逐漸擴大，直到找到三個點。在上面的例子中，範圍增加到-1.8到1.8。然而，我們分析的數據是不穩定的，並且具有隨機的峯值和谷值，導致非鄰接的點被接受到範圍中（元素3被選擇爲有效的，但不是元素4）。

什麼是最好的方式去做這件事？我已經有了一個代碼來遞增地增加範圍以找到三個點，我只需要修改它就不會在任何三點停下來，而是增加範圍直到它找到三個CONTIGUOUS點。如果這樣做是爲了上面的例子，我只是評估斜率去除第三個元素（因爲在3和4之間，斜率是負的）。

謝謝。

Answer 1

假設你的數據在本例中提供的是變量x，您可以使用isnan和findstr像這樣：

x = [NaN, NaN, NaN, -1.5363, NaN -1.7664, -1.7475, 123]; 
~isnan(x) 

ans = 

0  0  0  1  0  1  1  1 

pos = findstr(~isnan(x), [1 1 1]); 

之所以使用findstr像這是我們想找到序列由isnan返回的邏輯數組中的[1 1 1]和findstr將返回此序列出現的輸入數組中位置的索引。

對於您的示例數據，這將返回[]，但如果您將其更改爲我給出的示例中的數據，則它將返回6，並且您可以使用x(pos:pos+2)提取連續區域。如果存在超過3個連續值的情況（如果有4個，它將返回[6 7]）以及存在多個連續區域的情況，您必須小心一點。如果你不需要對這些情況做任何有意義的事情，那麼就使用pos（1）。

如果你要提取的第一相鄰區域，其長度大於或等於3的全部，你可以這樣做：

x = [NaN, NaN, NaN, -1.5363, NaN -1.7664, -1.7475, 123, 456, 789]; 

startPos = []; 
stopPos = []; 

pos = findstr(~isnan(x), [1 1 1]); 
if ~isempty(pos) 
    startPos = pos(1); 
    stopPos = startPos + 2; 

    % Find any cases where we have consecutive numbers in pos 
    if length(pos) > 1 && any(diff(pos) == 1) 
     % We have a contiguous section longer than 3 elements 

     % Find the NaNs 
     nans = find(isnan(x)); 
     % Find the first NaN after pos(1), or the index of the last element 
     stopPos = nans(nans > startPos); 
     if ~isempty(stopPos) 
      stopPos = stopPos(1) - 1; % Don't want the NaN 
     else 
      stopPos = length(x); 
     end 
    end 
end 

x(startPos:stopPos)