PhaseMargin コマンドは線形システムsys の位相余裕とゲインクロスオーバー周波数を算出します。

このコマンドは通常、マトリクスの形で入力された伝達関数に対応した位相余裕とそのクロスオーバー周波数のリストを、同じくマトリクスの形で出力します。入力が1×1の場合については、1つのリストを出力します。

もしある伝達関数について複数のゲインクロスオーバー周波数が存在する場合は、位相余裕が最小のものが採用されます。

もしある伝達関数についてゲインクロスオーバー周波数が存在しない場合は、リスト[Float(undefined),Float(undefined)]が出力されます。

離散システムの場合は、サンプル時間Tsで伝達関数の z から s への変換 (z = exp(s*Ts))が行われます。この場合、周波数レンジはナイキスト周波数以下に制限されます。

デフォルトでは位相余裕は度[deg]で出力されます。これはomega をゲインクロスオーバー周波数とすると、 180-argument(H(I*omega)) という形で計算されます。

parameters = set( name = complexcons ) または list( name = complexcons )

radians = truefalse

hertz = truefalse

with(DynamicSystems):

sys := TransferFunction (2/(1+s)):

PhaseMargin(sys);

PhaseMargin(sys,radians);

g := PhaseMargin(sys);

gain := g[1];

freq := g[2];

p1 := plot ([[freq, -180], [freq, gain-180]], color=blue, thickness=2):

p2 := PhasePlot (sys):

plots[display](p1, p2);