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

sectorplot

R2020b

周波数の関数としてのセクター インデックスの計算またはプロット

構文

sectorplot(H,Q)

sectorplot(H,Q,w)

sectorplot(H1,H2,...,HN,Q)

sectorplot(H1,H2,...,HN,Q,w)

sectorplot(H1,LineSpec1,...,HN,LineSpecN,Q)

sectorplot(H1,LineSpec1,...,HN,LineSpecN,Q,w)

[index,wout] = sectorplot(H,Q)

index = sectorplot(H,Q,w)

説明

sectorplot(H,Q) は動的システム H と指定されたセクター行列 Q の相対セクター インデックスをプロットします。これらのインデックスはセクター境界がどの程度満たされているか (1 未満のインデックス) または違反されているか (指定された周波数で 1 より大きいインデックス) によって測定します (セクター インデックスの意味の詳細についてはセクター境界とセクター インデックスについてを参照してください)。sectorplot は H のダイナミクスに基づいて自動的に周波数の範囲と点の数を選択します。

次の式を正の部分と負の部分への対称行列 Q の直交分解にします。

$$Q=W_1{W}_1^T-W_2{W}_2^T,W_1^T{W}_2=0.$$

セクター インデックスのプロットは、$$W_2^{T}H$$ に固有の安定した逆行列がある場合にのみ意味があります。この場合、セクター インデックスは次の式の特異値です。

$$\left(W_1^{T}H\left(j\omega \right)\right){\left(W_2^{T}H\left(j\omega \right)\right)}^{-1}.$$

sectorplot(H,Q,w) は w で指定した周波数のセクター インデックスをプロットします。

• w が形式 {wmin,wmax} の cell 配列の場合、sectorplot はセクター インデックスを wmin から wmax の範囲の周波数でプロットします。

• w が周波数のベクトルの場合、sectorplot はそれぞれの指定された周波数でセクター インデックスをプロットします。

sectorplot(H1,H2,...,HN,Q) と sectorplot(H1,H2,...,HN,Q,w) は同じプロット上で複数の動的システム H1,H2,...,HN のセクター インデックスをプロットします。

sectorplot(H1,LineSpec1,...,HN,LineSpecN,Q) と sectorplot(H1,LineSpec1,...,HN,LineSpecN,Q,w) は、プロット内の各システムに色、ライン スタイル、マーカーを指定します。

[index,wout] = sectorplot(H,Q) はベクトル wout の各周波数でセクター インデックスを返します。出力 index は行列です。値 index(:,k) はセクター インデックスを周波数 w(k) で降順で返します。この構文はプロットを描画しません。

index = sectorplot(H,Q,w) は w で指定される周波数でセクター インデックスを返します。

例

セクター インデックスと周波数のプロット

セクター インデックスをプロットして、次のように定義されるセクター内に $$G\left(s\right)=\left(s+2\right)/\left(s+1\right)$$ の I/O 軌跡が収まる周波数を可視化します。

$$S=\left\{\left(y,u\right):0.1u^2<uy<10u^2\right\}.$$

U/Y 領域では、このセクターは次の図の影付きの領域になります。

このセクターの Q 行列は次のようになります。

a = 0.1;  
b = 10; 
Q = [1 -(a+b)/2 ; -(a+b)/2 a*b];

軌跡 $$y\left(t\right)=Gu\left(t\right)$$ は、すべての T > 0 について次が成り立つとき、セクター S 内にあります。

$$0.1{\int }_0^{T}u{\left(t\right)}^2<{\int }_0^{T}u\left(t\right)y\left(t\right)dt<10{\int }_0^{T}u{\left(t\right)}^{2}dt.$$

周波数領域では、同じ条件を以下のように表すことができます。

$${\left(\begin{array}{c}G\left(j\omega \right)\\ 1\end{array}\right)}^{H}Q\left(\begin{array}{c}G\left(j\omega \right)\\ 1\end{array}\right)<0.$$

任意の周波数において G がこの条件を満たしているか違反しているかを確認するには、H = [G;1] のセクター インデックスをプロットします。

G = tf([1 2],[1 1]); 
sectorplot([G;1],Q)

プロットは、すべての周波数でセクター インデックスが 1 未満であることを示しています。したがって、G(s) の軌跡はすべての周波数において指定したセクター Q の範囲内に収まります。

MIMO システムによるセクター プロット

特定のセクターについて、2 出力 2 入力システムのセクター プロットを調べます。

rng(4,'twister');
H = rss(3,4,2); 
Q = [-5.12   2.16  -2.04   2.17
      2.16  -1.22  -0.28  -1.11
     -2.04  -0.28  -3.35   0.00
      2.17  -1.11   0.00   0.18];
sectorplot(H,Q)

H は 2 行 2 列なので、セクター プロットには 2 本のラインがあります。セクター インデックスの最大値は約 0.5 rad/s より低いところと、3 rad/s 前後の狭い帯域内で 1 を超えています。したがって、H は Q で表されるセクター境界を満たしていません。

入力引数

H — 解析するモデル

動的システムのモデル | モデル配列

セクター境界に対して解析するモデル。tf モデル、ss モデル、または genss モデルなどの動的システム モデルとして指定します。H は連続または離散のどちらにすることもできます。H が調整可能なブロックまたは不確かさをもつブロックを含む一般化モデルの場合、sectorplot は H の現在のノミナル値を解析します。

線形システム G のすべての I/O 軌跡 (u(t),y(t) が特定のセクターにあるかどうかを解析するには、H = [G;I] を使用します。ここで I = eyes(nu) のとき、nu は G の入力の数です。

H がモデル配列の場合、sectorplot は同じプロット内の配列にあるすべてのモデルのセクター インデックスをプロットします。出力引数を使用してセクター インデックスのデータを取得する場合、H は単一モデルでなければなりません。

Q — セクターの形状

行列 | LTI モデル

セクターの形状。次のように指定します。

• 定数のセクター形状の場合は行列。Q は 1 辺が ny の対称正方行列です。ここで、ny は H の出力数です。

• 周波数依存のセクター形状の場合は LTI モデル。Q は Q(s)’ = Q(–s) を満たします。つまり、Q(s) は各周波数でエルミート行列と評価されます。

行列 Q は適切に定義された円錐セクターを表すために不定でなければなりません。不定行列には、正と負両方の固有値があります。

詳細については、セクター境界とセクター インデックスについてを参照してください。

w — 周波数

{wmin,wmax} | ベクトル

インデックスを計算およびプロットする周波数。cell 配列 {wmin,wmax} または周波数値のベクトルとして指定します。

• w が形式 {wmin,wmax} の cell 配列の場合、関数は wmin から wmax の範囲の周波数でインデックスを計算します。

• w が周波数のベクトルの場合、関数は指定された各周波数でインデックスを計算します。たとえば、logspace を使用すると、対数的に等間隔な周波数値の行ベクトルを生成できます。

周波数はラジアン/TimeUnit 単位で指定します。ここで TimeUnit はモデルの TimeUnit プロパティです。

LineSpec — ライン スタイル、マーカー、色

文字ベクトル | string

ライン スタイル、マーカー、色。1、2、または 3 文字の string または文字ベクトルとして指定します。文字が表示される順序は任意です。3 つの特性 (ライン スタイル、マーカーおよび色) をすべて指定する必要はありません。たとえば、ライン スタイルを省略してマーカーを指定した場合、プロットはラインなしでマーカーのみを表示します。この引数の設定の詳細については、関数 plot の入力引数 LineSpec を参照してください。

例: 'r--' は赤い破線を指定します。

例: '*b' は青いアスタリスク マーカーを指定します。

例: 'y' は黄色いラインを指定します。

出力引数

index — セクター インデックス

行列

周波数の関数としてのセクター インデックス。行列として返されます。index には、これらを指定した場合は w、指定しない場合は wout の周波数で計算されたセクター インデックスが含まれます。index には w または wout にある値と同数の列と、H にある入力と同数の行が含まれます。このため、値 index(:,k) はセクター インデックスを周波数 w(k) で降順で返します。

たとえば、G が 3 入力 3 出力システム、Q が適切なセクター行列、w が 1 行 30 列の周波数のベクトルの場合、次の構文は 3 行 30 列の行列 index を返します。

H = [G;eyes(3)]
index = sectorplot(H,Q,w);

エントリ index(:,k) には H の 3 つのセクター インデックスが周波数 w(k) で降順で含まれます。

詳細については、セクター境界とセクター インデックスについてを参照してください。

wout — 周波数

ベクトル

インデックスが計算される周波数。ベクトルとして返されます。関数は、周波数範囲と点数をモデルのダイナミクスに基づいて自動的に選択します。

バージョン履歴

R2016a で導入

参考

トピック

• セクター境界とセクター インデックスについて