rooteig

固有ベクトル法を使用した周波数とパワー成分

ページ内をすべて折りたたむ

構文

[w,pow] = rooteig(x,p)
[w,pow] = rooteig(___,'corr')
[f,pow] = rooteig(___,fs)

説明

[w,pow] = rooteig(x,p) では、入力信号 x の周波数成分が推定され、ラジアン/サンプル単位の周波数のベクトル w と、対応する信号強度のベクトル pow が返されます。入力引数 p を使用して、信号部分空間の次元を指定できます。

p の 2 番目のエントリにある追加のしきい値パラメーターを使用すると、ノイズ部分空間と信号部分空間の割り当てを柔軟にコントロールできるようになります。

例

[w,pow] = rooteig(___,'corr') では、入力引数 x が信号データの行列ではなく、相関行列として解釈されます。この構文では、x は正方行列でなければなりません。また、すべての固有値は非負でなければなりません。この構文には、前の構文の入力引数を含めることができます。

メモ

'corr' は p の後の任意の位置に配置できます。

[f,pow] = rooteig(___,fs) では、Hz 単位で計算された周波数ベクトル f が返されます。サンプリング周波数 fs は Hz 単位で与えられます。fs を空ベクトル [] として指定すると、サンプリング周波数の既定値である 1 Hz が使用されます。

例

すべて折りたたむ

ライブ スクリプトを開く

3 つの複素指数成分とノイズから構成される信号の周波数成分を求めます。修正共分散法を使用して、固有ベクトル法で使用される相関行列を計算します。再現可能な結果が必要な場合は、乱数発生器をリセットします。

rng default
n = 0:99;   
s = exp(1i*pi/2*n)+2*exp(1i*pi/4*n)+exp(1i*pi/3*n)+randn(1,100);

X = corrmtx(s,12,'mod'); 
[W,P] = rooteig(X,3)
W = 3×1

    0.7883
    1.5674
    1.0429


P = 3×1

    4.1748
    1.0572
    1.2419

入力引数

すべて折りたたむ

ベクトルまたは行列として指定される入力信号。x がベクトルの場合は、信号の 1 つの観測値として扱われます。x が行列の場合、x の各行は信号の個別の観測値を表します。たとえば、x'*x が相関行列の推定となるような、配列処理の場合のように各行がセンサー配列の 1 つの出力である場合などです。

入力データ x が複素数の場合、pow と w は同じ長さです。実数値入力データ x に対して、対応する出力ベクトル pow の長さは 0.5*length(w) です。

メモ

関数 corrmtx の出力を使用して、このような配列 x を作成できます。

複素数のサポート: あり

部分空間の次元。実数の正の整数または 2 要素ベクトルとして指定します。p が実数の正の整数である場合は、部分空間の次元として扱われます。p が 2 要素ベクトルの場合は、p の 2 番目の要素は、信号の相関行列での推定固有値の最小値、λminを乗算したしきい値を表しています。しきい値 λmin*p(2) より小さい固有値は、ノイズ部分空間に割り当てられます。この場合、p(1) では、信号部分空間の最大次元が指定されます。p の 2 番目のエントリにある追加のしきい値パラメーターを使用すると、ノイズ部分空間と信号部分空間の割り当てを柔軟にコントロールできるようになります。

サンプル レート。正のスカラーで指定します。サンプル レート fs を Hz 単位で指定できます。fs を空ベクトル [] として指定すると、サンプル レートの既定値である 1 Hz が使用されます。

出力引数

すべて折りたたむ

ラジアン/サンプル単位の出力周波数。ベクトルとして返されます。ベクトル w の長さは、信号部分空間の計算された次元です。

信号強度。ベクトルとして返されます。

Hz 単位の出力周波数。ベクトルとして返されます。サンプリング周波数 fs は Hz 単位で与えられます。fs を空ベクトル [] として指定すると、サンプリング周波数の既定値である 1 Hz が使用されます。

アルゴリズム

関数 rooteig で使用する固有ベクトル法は、関数 peig で使用するものと同じです。このアルゴリズムでは、信号の周波数成分を推定するために、信号の相関行列の固有空間解析が行われます。

関数 peig と rooteig の違いは、以下のとおりです。

  • 関数 peig では、すべての周波数サンプルで疑似スペクトルが返されます。

  • 関数 rooteig は、離散周波数スペクトルの推定を対応する信号のパワー推定と共に返します。

関数 rooteig は、加法性ホワイト ガウス ノイズに正弦波が組み込まれているような信号の周波数の推定にもっとも有効です。

拡張機能

すべて展開する

バージョン履歴

R2006a より前に導入

参考

| | |