spa
スペクトル解析を使用した固定の周波数分解能による周波数応答の推定
説明
例
入力引数
出力引数
詳細
アルゴリズム
spa
は、次の手順に従ってブラックマン-テューキー スペクトル解析法を適用します。
共分散と相互共分散を u(t) と y(t) から計算します。
共分散と相互共分散のフーリエ変換を計算します。
ここで、 は幅 (ラグ サイズ) M のハン ウィンドウです。M を指定して推定の周波数分解能を制御できます。2π/M ラジアン/サンプル時間とほぼ等しくなります。
既定では、0 (この値は除く) から π までの等間隔な 128 個の周波数値がこの演算に使用されます。ここで、
w
=[1:128]/128*pi/Ts
であり、Ts
はそのデータ セットのサンプル時間です。ハン ウィンドウの既定のラグ サイズはM = min(length(data)/10,30)
です。既定の周波数の場合は高速フーリエ変換 (FFT) が演算に使用され、ユーザー定義の周波数よりも効率が高くなります。周波数応答関数 と出力ノイズ スペクトル を計算します。
spectrum
は、出力と入力の両方のチャネルのスペクトル行列です。つまり、z = [data.OutputData
, data.InputData]
の場合、spectrum
には z
の行列値のパワー スペクトルがスペクトル データとして格納されます。
ここで、'
は複素共役転置です。
参照
[1] Ljung, Lennart. System Identification: Theory for the User. 2nd ed. Prentice Hall Information and System Sciences Series. Upper Saddle River, NJ: Prentice Hall PTR, 1999.
バージョン履歴
R2006a より前に導入