evalfr

指定した周波数でのシステム応答の評価

構文

frsp = evalfr(sys,f)

説明

evalfr は freqresp の簡易バージョンであり、単一の特定周波数に対するラプラス変数値 s または z でのシステム応答を迅速に評価することを意図しています。一連の周波数でのシステム応答を評価するには、freqresp を使用します。周波数応答の振幅と位相のデータをそのプロットと共に取得するには、代わりに bode を使用してください。

例

frsp = evalfr(sys,f) は、連続時間モデルまたは離散時間モデル sys を指定の周波数 f で評価します。

例

すべて折りたたむ

離散時間の伝達関数の評価

ライブスクリプトを開く

次の離散時間伝達関数を作成します。

$H (z) = \frac{z - 1}{z^{2} + z + 1}$

H = tf([1 -1],[1 1 1],-1);

伝達関数を z = 1+j で評価します。

z = 1+j;
evalfr(H,z)

ans = 0.2308 + 0.1538i

与えられた周波数における同定されたモデルの周波数応答の評価

この例では次を使用します。

System Identification Toolbox System Identification Toolbox

ライブスクリプトを開く

次の連続時間の伝達関数モデルを作成します。

$H (s) = \frac{1}{s^{2} + 2 s + 1}$

sys = idtf(1,[1 2 1]);

周波数 0.1 ラジアン/秒で伝達関数を評価します。

w = 0.1;
s = j*w;
evalfr(sys,s)

ans = 0.9705 - 0.1961i

あるいは、freqresp コマンドを使用します。

freqresp(sys,w)

ans = 0.9705 - 0.1961i

MIMO 状態空間モデルの周波数応答

ライブスクリプトを開く

この例では、慣性テンソル J と振幅が 0.2 の減衰力 F で角を中心にして回転しているキューブについて考えます。システムへの入力は駆動トルクであり、角速度が出力です。このキューブの状態空間行列は次のようになります。

$\begin{array}{l} A = - J^{- 1} F, B = J^{- 1}, C = I, D = 0, \\ w h e r e, J = [\begin{array}{ccc} 8 & - 3 & - 3 \\ - 3 & 8 & - 3 \\ - 3 & - 3 & 8 \end{array}] a n d F = [\begin{array}{ccc} 0.2 & 0 & 0 \\ 0 & 0.2 & 0 \\ 0 & 0 & 0.2 \end{array}] \end{array}$

行列 A、B、C および D を指定して、連続時間状態空間モデルを作成します。

J = [8 -3 -3; -3 8 -3; -3 -3 8];
F = 0.2*eye(3);
A = -J\F;
B = inv(J);
C = eye(3);
D = 0;
sys = ss(A,B,C,D);
size(sys)

State-space model with 3 outputs, 3 inputs, and 3 states.

0.2 ラジアン/秒におけるシステムの周波数応答を計算します。sys が連続時間モデルであるため、周波数をラプラス変数 s で表します。

w = 0.2;
s = j*w;
frsp = evalfr(sys,s)

frsp = 3×3 complex

   0.3607 - 0.9672i   0.3197 - 0.5164i   0.3197 - 0.5164i
   0.3197 - 0.5164i   0.3607 - 0.9672i   0.3197 - 0.5164i
   0.3197 - 0.5164i   0.3197 - 0.5164i   0.3607 - 0.9672i

あるいは、freqresp コマンドを使用し、周波数のスカラー値を直接使用して周波数応答を評価することもできます。

H = freqresp(sys,w)

H = 3×3 complex

   0.3607 - 0.9672i   0.3197 - 0.5164i   0.3197 - 0.5164i
   0.3197 - 0.5164i   0.3607 - 0.9672i   0.3197 - 0.5164i
   0.3197 - 0.5164i   0.3197 - 0.5164i   0.3607 - 0.9672i

入力引数

すべて折りたたむ

`sys` — 動的システム
動的システムモデル | モデル配列

動的システム。SISO または MIMO 動的システムモデルか、動的システムモデルの配列として指定します。使用できる動的システムには次のようなものがあります。

ss モデル、tf モデル、zpk モデルなどの LTI モデル。
sparss モデルまたは mechss モデルなどのスパース状態空間モデル。
genss または uss (Robust Control Toolbox) モデルなどの一般化された状態空間モデルまたは不確かさをもつ状態空間モデル。(不確かさをもつモデルを使用するには Robust Control Toolbox™ ソフトウェアが必要です。)
- 調整可能な制御設計ブロックの場合、関数はモデルをその現在の値で評価して周波数応答データを評価します。
- 不確かさをもつ制御設計ブロックの場合、関数はモデルのノミナル値とランダムサンプルで周波数応答を評価します。
idss (System Identification Toolbox) モデルなどの同定された状態空間モデル。(同定されたモデルを使用するには System Identification Toolbox™ ソフトウェアが必要です。)

モデルの完全なリストについては、動的システムモデルを参照してください。

`f` — システム応答を評価する周波数の位置
複素数スカラー

ラプラス変数 s または z として表される、システム応答を評価する周波数の位置。複素数スカラーとして指定します。sys が連続時間モデルであるか離散時間モデルであるかに基づき、それぞれラプラス変数 s または z で周波数を指定します。たとえば、システム sys の周波数応答を周波数値 w rad/s で評価する場合は、以下を使用して f を指定します。