このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。

dsp.MovingRMS

移動平方根平均二乗

説明

dsp.MovingRMS System object™ は、時間の経過に沿って入力信号の移動平方根平均二乗 (RMS) を各チャネルで個別に計算します。このオブジェクトでは、スライディング ウィンドウ法または指数の重み付け法のどちらかを使用して移動 RMS を計算します。スライディング ウィンドウ法では、データ上で指定の長さのウィンドウをサンプルごとに移動させ、ウィンドウ内のデータの RMS を計算します。指数の重み付け法では、オブジェクトはデータ サンプルを二乗し、一連の重み係数を掛け、重み付けしたデータの合計を求めます。次に、オブジェクトは合計の平方根を取ることで RMS を計算します。これらの方法の詳細については、アルゴリズムを参照してください。

入力の移動 RMS を計算するには、次のようにします。

1. dsp.MovingRMS オブジェクトを作成し、そのプロパティを設定します。

2. 関数と同様に、引数を指定してオブジェクトを呼び出します。

System object の機能の詳細については、System object とはを参照してください。

作成

構文

movRMS = dsp.MovingRMS

movRMS = dsp.MovingRMS(Len)

MovRMS = dsp.MovingRMS(Len,Overlap)

movRMS = dsp.MovingRMS(Name,Value)

説明

movRMS = dsp.MovingRMS は、既定のプロパティを使用して、移動 RMS オブジェクト movRMS を返します。

movRMS = dsp.MovingRMS(Len) は WindowLength プロパティを Len に設定します。

MovRMS = dsp.MovingRMS(Len,Overlap) は、WindowLength プロパティを Len に、OverlapLength プロパティを Overlap に設定します。

movRMS = dsp.MovingRMS(Name,Value) は、Name,Value ペアを使用して追加のプロパティを指定します。指定していないプロパティは既定値になります。

例: movRMS = dsp.MovingRMS('Method','Exponential weighting','ForgettingFactor',0.9);

プロパティ

特に指定がない限り、プロパティは "調整不可能" です。つまり、オブジェクトの呼び出し後に値を変更することはできません。オブジェクトは呼び出すとロックされ、ロックを解除するには関数 release を使用します。

プロパティが "調整可能" の場合、その値をいつでも変更できます。

プロパティ値の変更の詳細については、System object を使用した MATLAB でのシステム設計を参照してください。

Method — 移動 RMS メソッド

'Sliding window' (既定値) | 'Exponential weighting'

移動 RMS メソッド。'Sliding window' または 'Exponential weighting' として指定します。

• 'Sliding window' — SpecifyWindowLength によって指定された長さのウィンドウを入力データ上で各チャネルに沿って移動させます。オブジェクトは、ウィンドウが移動する各サンプルについて、ウィンドウ内のデータの RMS を計算します。

• 'Exponential weighting' — オブジェクトはサンプルの二乗と一連の重み係数を乗算します。データが古くなるほど、重み係数の大きさは指数的に小さくなります (0 になることはありません)。RMS を計算するために、このアルゴリズムは重み付けしたデータを合計し、その合計の平方根を取ります。

これらの方法の詳細については、アルゴリズムを参照してください。

SpecifyWindowLength — ウィンドウの長さの指定

true (既定値) | false

ウィンドウの長さを指定するフラグ。スカラーの boolean として指定します。

• true — スライディング ウィンドウの長さが WindowLength プロパティで指定した値と等しくなります。

• false — スライディング ウィンドウの長さが無限になります。このモードでは、RMS は現在のサンプルと以前のすべてのサンプルを使用して計算されます。

依存関係

このプロパティは Method を 'Sliding window' に設定した場合に適用されます。

WindowLength — スライディング ウィンドウの長さ

4 (既定値) | 正のスカラー整数

サンプルのスライディング ウィンドウの長さ。正のスカラー整数として指定します。

依存関係

このプロパティは Method を 'Sliding window'、SpecifyWindowLength を true に設定した場合に適用されます。

データ型: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

OverlapLength — ウィンドウ間のオーバーラップの長さ

WindowLength − 1 (既定値) | 非負の整数

スライディング ウィンドウ間のオーバーラップの長さ。非負の整数として指定します。オーバーラップの長さの値は、[0, WindowLength − 1] の範囲で変化します。指定しない場合、オーバーラップの長さは WindowLength − 1 となります。

依存関係

このプロパティは Method を 'Sliding window'、SpecifyWindowLength を true に設定した場合に適用されます。

データ型: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64

ForgettingFactor — 指数の重み付けの係数

0.9 (既定値) | (0,1] の範囲の正の実数スカラー

指数の重み付けの係数。(0,1] の範囲の正の実数スカラーとして指定します。

忘却係数 0.1 よりも忘却係数 0.9 の方が古いデータに対する重みが大きくなります。忘却係数が 1.0 の場合は無限に記憶します。過去のサンプルに対する重みはすべて同じになります。

このプロパティは調整可能であるため、その値はオブジェクトがロックされている場合でも変更できます。

調整可能: Yes

依存関係

このプロパティは Method を 'Exponential weighting' に設定した場合に適用されます。

データ型: single | double

使用法

構文

y = movRMS(x)

説明

y = movRMS(x) は、スライディング ウィンドウ法または指数の重み付け法のどちらかを使用して、入力信号 x の移動 RMS を計算します。

入力引数

x — データ入力

ベクトル | 行列

データ入力。ベクトルまたは行列として指定します。x が行列の場合、各列が個別のチャネルとして扱われます。移動 RMS は各チャネルに沿って計算されます。

オブジェクトは可変サイズの入力を受け入れます。オブジェクトがロックされると、各入力チャネルのサイズは変更できますが、チャネルの数は変更できません。

データ型: single | double
複素数のサポート: あり

出力引数

y — 移動 RMS

ベクトル | 行列

入力信号の移動 RMS。ベクトルまたは行列として返されます。

m 行 n 列のサイズの信号をオブジェクトに入力し、Method を 'Sliding window' に、SpecifyWindowLength を true に設定すると、出力は ceil(m/ホップ サイズ) 行 n 列の上限サイズをもちます。ホップ サイズは、ウィンドウの長さからオーバーラップの長さを引いた値です。出力が m 行 n 列のサイズをもつ場合もあります。

このオブジェクトからコードを生成した場合、生成したコードの出力の可変サイズの動作は、入力フレームの長さと入力信号が固定サイズか可変サイズかによって異なります。詳細については、コード生成を参照してください。

データ型: single | double
複素数のサポート: あり

オブジェクト関数

オブジェクト関数を使用するには、System object を最初の入力引数として指定します。たとえば、obj という名前の System object のシステム リソースを解放するには、次の構文を使用します。

release(obj)

すべての System object に共通

step	System object のアルゴリズムの実行
release	リソースを解放し、System object のプロパティ値と入力特性の変更を可能にします。
reset	System object の内部状態のリセット

例

ノイズを含む矩形波信号の移動 RMS

dsp.MovingRMS オブジェクトを使用して、可変振幅を持つノイズを含む矩形波信号の移動 RMS を計算します。

初期化

movrmsWin オブジェクト、movrmsWin_overlap オブジェクト、および movrmsExp オブジェクトを設定します。movrmsWin は、ウィンドウの長さを 20 サンプルとしてスライディング ウィンドウ法を使用します。既定のオーバーラップの長さは 19 サンプルで、これは指定のウィンドウ長より 1 サンプル分小さくなります。movrmsWin_overlap はウィンドウの長さを 20 サンプルとし、オーバーラップの長さは 15 サンプルです。movrmsExp は忘却係数を 0.995 として指数的な重み付け法を使用します。

出力を表示するための時間スコープを作成します。

FrameLength = 10;
Fs = 100;
movrmsWin = dsp.MovingRMS(20);
movrmsWin_overlap = dsp.MovingRMS(20,15);
movrmsExp = dsp.MovingRMS('Method','Exponential weighting',...
    'ForgettingFactor',0.995);
scope  = timescope('SampleRate',[Fs,Fs,Fs/(20-15),Fs],...
    'TimeSpanOverrunAction','Scroll',...
    'TimeSpanSource','Property',...
    'TimeSpan',100,...
    'ShowGrid',true,...
    'YLimits',[-1.0 5.5]);
title = 'Moving RMS';
scope.Title = title;
scope.ChannelNames = {'Original Signal',...
    'Sliding window of 20 samples with default overlap',...
    'Sliding window of 20 samples with an overlap of 15 samples',...
    'Exponential weighting with forgetting factor of 0.995'};

RMS の計算

ノイズを含む矩形波信号を生成します。特定のフレーム数の後、矩形波の振幅を変化させます。スライディング ウィンドウ法および指数の重み付け法をこの信号に対して適用します。時間スコープに出力が表示されます。

count = 1;
Vect = [1/8 1/2 1 2 3 4];
for index = 1:length(Vect)
    V = Vect(index);
    for i = 1:160
        x = V + 0.1 * randn(FrameLength,1);
        y1 = movrmsWin(x);
        y2 = movrmsWin_overlap(x);
        y3 = movrmsExp(x);
        scope(x,y1,y2,y3);
    end
end

アルゴリズム

スライディング ウィンドウ法

スライディング ウィンドウ法では、各入力サンプルの出力は現在のサンプルと以前の Len - 1 個のサンプルの RMS になります。Len はサンプルのウィンドウの長さです。最初の出力サンプルを計算するために、アルゴリズムはホップ サイズの数の入力サンプルを受け取るまで待機します。ホップ サイズは、ウィンドウの長さからオーバーラップの長さを引いた値として定義されます。ウィンドウ内の残りのサンプルはゼロと見なされます。たとえば、ウィンドウの長さが 5 でオーバーラップの長さが 2 の場合、アルゴリズムは最初の出力サンプルを計算するために、入力サンプルを 3 つ受け取るまで待機します。最初の出力の生成後、ホップ サイズの数の入力サンプルごとに後続の出力サンプルを生成します。

ウィンドウの長さを指定しない場合、アルゴリズムはウィンドウの長さを無限と見なします。このモードでは、出力はチャネルの現在のサンプルと以前のすべてのサンプルの移動 RMS になります。

スライディング ウィンドウ法を使用してストリーミング入力データの移動 RMS を計算する例を考えます。アルゴリズムはウィンドウの長さ 4 とオーバーラップの長さ 3 を使用します。受け取った各入力サンプルで、長さ 4 のウィンドウがデータに沿って移動します。

指数の重み付け法

指数の重み付け法では、次の式を使用して再帰的に移動 RMS が計算されます。

$$\begin{array}{l}{w}_{N,\lambda }=\lambda w_{N-1,\lambda }+1\\ x\_rm{s}_{N,\lambda }=\sqrt{\left(1-\frac{1}{w_{N,\lambda }}\right)x\_rm{s}_{N-1,\lambda }+\left(\frac{1}{w_{N,\lambda }}\right)x^2{}_N}\end{array}$$

• $$x\_rm{s}_{N,\lambda }$$ — 現在のサンプルの移動 RMS

• $$x^2{}_N$$ — 現在の入力データ サンプルの二乗

• $$x\_rm{s}_{N-1,\lambda }$$ — 前のサンプルの移動 RMS

• λ — 忘却係数

• $$w_{N,\lambda }$$ — 現在のデータ サンプルに適用される重み係数

• $$\left(1-\frac{1}{w_{N,\lambda }}\right)x\_rm{s}_{N-1,\lambda }$$ — RMS に対する前のデータの影響

アルゴリズムは、最初のサンプル (N = 1) に対して$$w_{N,\lambda }$$ = 1 を選択します。次のサンプルについては、再帰方程式に従って更新された重み係数を使用して RMS を計算します。データが古くなるほど、重み係数の大きさは指数的に小さくなります (0 になることはありません)。つまり、現在の RMS に対する影響は古いデータよりも新しいデータの方が大きくなります。

重み係数の変化率は忘却係数の値で決まります。忘却係数 0.1 よりも忘却係数 0.9 の方が古いデータに対する重みが大きくなります。忘却係数が 1.0 の場合は無限に記憶します。以前のサンプルに対する重みはすべて同じになります。

以下は、指数の重み付け法を使用した移動 RMS の計算例です。忘却係数は 0.9 です。

参照

[1] Bodenham, Dean. “Adaptive Filtering and Change Detection for Streaming Data.” PH.D. Thesis. Imperial College, London, 2012.

拡張機能

C/C++ コード生成

MATLAB® Coder™ を使用して C および C++ コードを生成します。

使用上の注意および制限:

• MATLAB コード生成における System object (MATLAB Coder)を参照してください。

• Method を 'Sliding window' に、SpecifyWindowLength を true に設定してこのオブジェクトからコードを生成した場合の出力の可変サイズの動作は、入力フレームの長さと入力信号が固定サイズか可変サイズかによって異なります。

  詳細については、次の表を参照してください。

  入力信号	入力サイズ	出力信号
  固定サイズ	m 行 n 列	入力フレーム長がホップ サイズの倍数である場合、出力信号は (m/ホップ サイズ) 行 n 列の固定サイズとなります。 入力フレーム長がホップ サイズの倍数でない場合、出力信号は可変サイズとなり、ceil(m/ホップ サイズ) 行 n 列の上限をもちます。
  可変サイズ	m 行 n 列	出力は可変サイズの信号です。 出力は、ceil(m/ホップ サイズ) 行 n 列の上限サイズをもちます。

バージョン履歴

R2016b で導入

R2022b: 生成したコードの出力信号の可変サイズの動作の変更

R2022b 以降、Method プロパティを 'Sliding window' に、SpecifyWindowLength プロパティを true に設定してこのオブジェクトからコードを生成し、フレーム長がホップ サイズ (ウィンドウの長さからオーバーラップの長さを引いた値) の倍数となる信号を入力した場合、このオブジェクトは、生成したコードに固定サイズの出力信号を生成します。詳細については、コード生成を参照してください。

参考

オブジェクト

• dsp.MovingAverage | dsp.MovingMaximum | dsp.MovingMinimum | dsp.MovingStandardDeviation | dsp.MovingVariance | dsp.MedianFilter

ブロック

• Moving RMS | RMS | Moving Average | Moving Maximum | Moving Minimum | Moving Standard Deviation | Moving Variance | Median Filter

トピック

• What Are Moving Statistics?
• ストリーミング信号の統計
• スライディング ウィンドウ法と指数の重み付け法
• Energy Detection in the Time Domain