ナイキスト線図・安定余裕とは
🙋
ナイキスト線図って、複素平面上にぐるぐる描くあの曲線ですよね？これを見てどうやって「安定」か判断するんですか？
🎓
大まかに言うと、そのぐるぐるした軌跡が「臨界点」と呼ばれる(-1, 0)をどう回るかで決まるんだ。このシミュレーターで、上のゲインKのスライダーを大きくしていくと、軌跡がだんだん(-1, 0)に近づいて、ついには囲み始める様子がリアルタイムで見えるよ。囲むと不安定になるんだ。
🙋
え、囲むとダメなんですか？「ゲイン余裕」や「位相余裕」って表示されてますけど、これとどう関係あるんですか？
🎓
いいところに気づいたね。軌跡がギリギリ囲まない状態でも、(-1, 0)に近すぎると実用上は危険なんだ。例えばロボットの制御で振動が起きたりする。そこで「あとどれだけ余裕があるか」を数値化したのがゲイン余裕(GM)と位相余裕(PM)だ。シミュレーターで極`a`や零点`z`を変えると、この余裕がどう変わるか、すぐに確かめられるよ。
🙋
なるほど！でも、右側に同時に出てるボード線図と、ナイキスト線図ってどう対応してるんですか？見方が大きく異なる気がするんですが…。
🎓
実は表裏一体なんだ。ボード線図でゲインが0dBになる周波数（ゲイン交差周波数）での位相を見れば、それが位相余裕(PM)だ。逆に、位相が-180度になる周波数でのゲインを見れば、それがゲイン余裕(GM)だ。シミュレーターでパラメータをいじりながら、両方の図で同じ値が強調表示されるのを確認してみて。現場のエンジニアは、周波数特性の直感的理解にボード線図、安定性の厳密な判定や余裕の視覚化にナイキスト線図を使い分けてるんだ。
物理モデルと主要な数式
このシミュレーターで扱う標準形の開ループ伝達関数 $L(s)$ です。パラメータK, a, b, zは画面上のスライダーで変更できます。
$$L(s) = G(s)H(s) = K \frac{(s/z + 1)}{(s/a + 1)(s/b + 1)}$$
$K$: ゲイン（定数）, $z$: 零点の値, $a, b$: 極の値。ここで $s = j\omega$ (虚数単位×角周波数) を代入することで、周波数応答 $L(j\omega)$ が得られ、ナイキスト線図が描けます。
ナイキストの安定判別則。閉ループ系の安定性を、開ループの情報だけから判定できる強力な定理です。
$$Z = N + P$$
$P$: 開ループ伝達関数 $L(s)$ の右半平面（実部が正）にある極の数。
$N$: ナイキスト軌跡 $L(j\omega)$ ($\omega: -\infty \to +\infty$) が点 $(-1, 0)$ を反時計回りに囲む回数。
$Z$: 求めたい閉ループ系の不安定極の数。$Z=0$ なら閉ループ系は安定です。
よくある質問
ゲイン余裕は位相が-180°のときのゲイン[dB]の絶対値、位相余裕はゲインが0dBのときの位相と-180°の差です。ボード線図上で該当周波数の値を確認するか、画面右の自動計算結果を参照してください。余裕が正なら安定、負なら不安定です。
スライダーをドラッグ後、マウスを離すか、数値入力欄でEnterキーを押すとリアルタイム更新されます。ブラウザのJavaScriptが有効か確認し、ページをリロードしても改善しない場合は、キャッシュをクリアして再度アクセスしてください。
ゲインKを小さくするとゲイン余裕が増え安定化しやすくなります。また、零点zを大きく（左半平面に移動）すると位相進み効果で位相余裕が改善します。極a,bを小さくすると応答が遅くなりますが安定性は向上します。各パラメータを試しながらナイキスト線図が(-1,0)を左に見るように調整してください。
Pは開ループ伝達関数L(s)の右半平面（実部>0）にある極の数です。本シミュレーターでは極a,bが正の値だと右半平面極となります。例えばa=2（正）ならP=1です。画面上の「極の値」が正か負かを確認し、正の個数を数えてください。
実世界での応用
航空機・ドローンの姿勢制御：飛行中の機体は常に外乱（風など）を受けます。ナイキスト線図で安定余裕（GM, PM）を十分に確保することで、この外乱に対しても振動することなく安定した飛行を維持できる制御系を設計します。余裕が小さいとパイロットの操作に対して危険な振動（ピオー現象）が発生する可能性があります。
自動車のクルーズコントロール：設定した速度を維持するためのエンジンやブレーキの制御です。坂道など負荷が変動しても速度が大きくぶれないようにするため、制御系の安定性と応答性が重要です。ナイキスト線図を用いて、追従性能と安定性のバランスの取れたコントローラゲインを決定します。
ロボットアームの位置決め制御：工場の生産ラインで精密な部品を組み立てるロボットは、指定された位置に素早く、かつオーバーシュート（行き過ぎ）なく移動する必要があります。位相余裕(PM)はこの過渡応答の振動特性と強く関連しており、設計の重要な指標となります。
電力系統の安定度解析：大規模な送電網では、発電機同士の同期を保つことがシステム全体の安定性に直結します。ナイキスト判別法は、このような高次元で複雑なシステムの安定性を周波数領域で評価するための古典的かつ有効な手法として今も用いられています。
よくある誤解と注意点
まず、「ナイキスト線図が(-1,0)を囲まなければ絶対に安定」とは限らないという点に注意だ。これはナイキストの安定判別則を適用する大前提として、開ループ系が安定（P=0）である場合の話だからだ。例えば、倒立振子のような元々不安定なプラント（P>0）を制御する場合、この前提が崩れる。シミュレーターで極`a`や`b`を負の値（右半平面に極がある状態）に設定して試してみると、話が変わってくるのがわかるよ。
次に、安定余裕の数値だけを盲信しないこと。例えば、ゲイン余裕(GM)が10dB、位相余裕(PM)が45度あれば教科書的には十分そうに見える。しかし、ナイキスト軌跡のカーブの形状が鋭く尖っていたり、(-1,0)のすぐ近くをかすめていくような形だと、モデルのわずかな誤差や外乱で簡単に不安定領域に飛び込んでしまう。このツールでKを上げ下げしながら、余裕の「質」、つまり軌跡の(-1,0)への近づき方も目視で確認する習慣をつけよう。
最後に、シミュレーターの伝達関数モデルと現実のギャップを理解しておくこと。このツールで扱っているのはごくシンプルな2次系モデルだ。実際の制御対象には、むだ時間やより高次の振動モード、非線形性が必ず含まれる。例えば、むだ時間 $e^{-Ls}$ が加わると位相がさらに遅れ、ナイキスト軌跡はらせん状に無限に回り込む。ツールで学んだ基本原理を、より現実に近い複雑なモデルへどう拡張していくかを常に考えてほしい。