100赫兹で1秒間サンプリングされる正弦波信号を生成します。ノイズをシミュレーションするため,高周波数の正弦波を追加します。

fs = 100;t = 0:1 / fs: 1;x =罪(2 *π* t * 3) + 0.25 *罪(2 *π* t * 40);

10次のメディアンフィルターを使用して,信号を滑らかにします。結果をプロットします。

y = medfilt1 (x, 10);情节(t t, x,, y)传奇(“原始”，“过滤”)传说(“boxoff”）

異なる周波数の正弦波で構成される2チャネル信号を生成します。ランダムな位置にスパイクを配置します。南を使用し、欠損サンプルをランダムに追加します。再現可能な結果が必要な場合は、乱数発生器をリセットします。信号をプロットします。

rng (“默认”n = 59;x =罪(π。/ 10[15]‘* (1:n) +π/ 3)';randi spk = (2 * n, 9,1);(spk) = x (spk) * 2;x (randi (2 * n, 6,1)) =南;情节(x)

medfilt1を既定の設定で使用して,信号をフィルター処理します。フィルター処理された信号をプロットします。既定では、フィルターは欠損サンプルを含むすべてのセグメントの中央値に NaN を代入します。

y = medfilt1 (x);情节(y)

元の信号を転置します。関数が行に沿って動作するように指定して、再度フィルター処理を行います。中央値を計算する際には、欠損サンプルを排除します。2 番目の引数を空のままにした場合、medfilt1は既定のフィルター次数3を使用します。

y = medfilt1 (x ' [] [], 2,“omitnan”）;情节(y ')

関数は,南のみを含むセグメントに値を代入することはできません。この問題を解決するには,セグメント長を増やします。さらに,この変更によって,外れ値がより徹底的に削除されるようになります。

y = medfilt1 (x 4“omitnan”）;情节(y)

既定のゼロパディングは,関数がエッジにおける信号値を過小に推定する原因になります。縮小ウィンドウを使用して両端における中央値を計算することで,この影響が低減されます。

y = medfilt1 (x 4“omitnan”，“截断”）;情节(y)