脳ティーザー - 移動平均を使用したフィルタリングアルゴリズム

Question

私はMatlabの86400風速（WS）値の1秒のデータセットを持っており、それをフィルタリングする際に助けが必要です。一定のレベルの賢さが必要です。

平均WSを超えた場合：

• 25メートル/秒の30秒の時間間隔
• 3秒の時間間隔

で30メートル/秒で600S時間間隔

28メートル/秒

これらのパラメータのいずれかが満たされている場合、平均WSが300秒の時間間隔で22m/s以下になるまで、WSは「無効」とみなされます。

ここには600秒の要件があります。私は 'データセット'に含まれるデータに対して600〜300秒の移動平均を行います。平均25m/sの最初の出現から22m/s以下の次の出現までの間隔を 'NaN'としてフィルタリングします。フィルタリングの後、別の600秒間の平均を行い、NaNでフラグが立てられた値の間隔はNaNになります。

すなわち

Rolling600avg(:,1) = tsmovavg(dataset(:,2), 's', 600, 1); 

Rolling300avg(:,1) = tsmovavg(dataset(:,2), 's', 300, 1); 

a = find(Rolling600avg(:,2)>25) 

b = find(Rolling300avg(:,2)<22) 

dataset(a:b(a:find(b==1)),2)==NaN; %?? Not sure 

これは、「検索」や、いくつかのインデックスの巧妙な使用を必要とする予定です。誰かが私を助けてくれますか？ 28m/sと30m/sのフィルターは同じ方法に従います。

Answer 1

私があなたの質問に従えば、1つの方法は、forループを使用して、NaNの開始と終了の場所を特定することです。

m = [19 19 19 19 28 28 19 19 28 28 17 17 17 19 29 18 18 29 18 29]; %Example data 
a = find(m>25); 
b = find(m<22); 
m2 = m; 
% Use a loop to isolate segments that should be NaNs; 
for ii = 1:length(a) 
    firstNull = a(ii) 
    lastNull = b(find(b>firstNull,1))-1 % THIS TRIES TO FIND A VALUE IN B GREATER THAN A(II) 
    % IF THERE IS NO SUCH VALUE THEN NANS SHOULD FILL TO THE END OF THE VECTOR 
    if isempty(lastNull), 
     lastNull=length(m); 
    end 
    m2(firstNull:lastNull) = NaN 
end 

これは、tsmovavgが渡されたものと等しい長さのベクトルを戻す場合にのみ機能することに注意してください。そうでなければそれはやりにくく、いくつかの変更が必要になります。

おそらくforループを回避する方法がありますが、これはかなり簡単な解決策です。