freqresp
グリッド上の周波数応答
構文
[H,wout] = freqresp(sys)
H = freqresp(sys,w)
H = freqresp(sys,w,units)
[H,wout,covH] = freqresp(idsys,...)
説明
[
は、周波数 H
,wout
] = freqresp(sys
)wout
における動的システム モデル sys
の周波数応答を返します。freqresp
コマンドは、sys
のダイナミクスに基づいて周波数を自動的に判断します。
は、ベクトル H
= freqresp(sys
,w
)w
で指定された実周波数グリッドの周波数応答を返します。
は、H
= freqresp(sys
,w
,units
)units
を使用して、w
の周波数単位を明示的に指定します。
[
は、同定されたモデル H
,wout
,covH
] = freqresp(idsys
,...)idsys
の周波数応答の共分散 covH
も返します。
入力引数
|
任意の動的システム モデルまたはモデル配列。 |
|
周波数応答を評価する実周波数のベクトル。周波数を |
|
入力周波数ベクトル
既定値: |
|
任意の同定されたモデル。 |
出力引数
|
周波数応答値を含む配列。
|
|
|
|
応答 |
例
詳細
アルゴリズム
伝達関数モデルまたは零点-極-ゲイン モデルの場合、freqresp
は指定された周波数点で分子 (s) と分母 (s) を評価します。連続時間の状態空間モデル (A, B, C, D) の周波数応答は次のようになります。
効率を高めるため、A は上ヘッセンベルグ形式に低減され、ヘッセンベルグ構造を利用しつつ、線形方程式 (jω − A)X = B を各周波数点で解きます。ヘッセンベルグ形式に低減することで、効率と信頼性の両方の妥協を取ることができます。この技法の詳細は、[1] を参照してください。
代替方法
evalfr
を使用して、個々の周波数または少数の周波数の周波数応答を評価します。freqresp
は中~大規模の周波数ベクトル用に最適化されています。
参考文献
[1] Laub, A.J., "Efficient Multivariable Frequency Response Computations," IEEE® Transactions on Automatic Control, AC-26 (1981), pp. 407-408.