ボード線図シミュレーター 戻る
制御工学

ボード線図シミュレーター

伝達関数のパラメータを変えてゲイン線図・位相線図をリアルタイム生成。ゲイン余裕・位相余裕・安定判別を直感的に体験できます。

システムタイプ

パラメータ

s
安定
計算結果
ゲイン交差周波数 (rad/s)
位相交差周波数 (rad/s)
位相余裕 (°)
ゲイン余裕 (dB)

ゲイン線図 (Gain Plot)

ゲイン

位相線図 (Phase Plot)

位相
理論・主要公式

$$G(j\omega) = |G(j\omega)| e^{j\angle G(j\omega)}$$

周波数応答:複素数表現。ゲイン \(|G|\)(dB)と位相 \(\angle G\)(deg)を周波数で表示する。

$$|G(j\omega)|_{dB} = 20\log_{10}|G(j\omega)|$$

ゲイン(dB)変換式。0 dBは倍率1倍、20 dBは10倍、−20 dB/decは1次遅れ傾斜。

$$G_M = -|G(j\omega_{pc})|_{dB}, \quad P_M = 180°+\angle G(j\omega_{gc})$$

ゲイン余裕 \(G_M\)(位相交差周波数での)と位相余裕 \(P_M\)(ゲイン交差周波数での安定余裕)。

💬 博士に聞いてみた

🙋
ボード線図ってそもそも何を表してるんですか?なんで対数スケールを使うんですか?
🎓
制御システムに「0.1rad/s の正弦波を入れたらどれだけ増幅・減衰されて、位相がどのくらいずれるか」を周波数ごとに表したものだよ。対数スケールを使う理由は、実用的な周波数レンジが 0.001〜10000 rad/s と広大なのと、直列システムのゲインが掛け算→対数では足し算になって計算が楽になるから。
🙋
位相余裕が少ないと何が起きるんですか?
🎓
ステップ応答で振動が激しくなる。位相余裕が 30° 以下だと大きなオーバーシュートが起きて、0° 以下になると発振(不安定)する。実用的には位相余裕 45〜60°、ゲイン余裕 6dB 以上が設計目標になることが多い。サーボモーター・ロボットアーム・プロセス制御でも同じだよ。
🙋
PIDコントローラを足すと、ゲイン線図の形がどう変わりますか?
🎓
積分要素(I)は低周波のゲインを上げる(-20dB/dec の傾きを加算)ので定常誤差ゼロにできる。微分要素(D)は高周波のゲインを上げて(+20dB/dec)位相を進める。位相が進むと位相余裕が増えて安定化に貢献する。ただし高周波ノイズも増幅されるので Td には上限をつけるのが普通だよ。
🙋
むだ時間(デッドタイム)が入るとなぜ制御が難しくなるんですか?
🎓
むだ時間 e^(-Ls) はゲインには影響しないが位相を -ωL(ラジアン)だけ引き続ける。周波数が上がるほど位相遅れが増えるので、高い周波数帯での位相余裕が激減する。プロセス制御(化学プラント・蒸留塔)では輸送遅れが大きくて有名な難問。スミス予測器などの専用手法が使われるよ。

❓ よくある質問

ボード線図はどんな制御設計に使いますか?

ループゲイン成形(Loop Shaping)と呼ばれる手法で、目標の帯域幅・位相余裕・ゲイン余裕を満たすよう補償器を設計します。PID調整・リードラグ補償・ノッチフィルタ設計など幅広く使われます。

安定判別のためのナイキスト線図との違いは?

ボード線図は伝達関数が安定な系(最小位相系)の安定余裕を直感的に確認するのに優れています。ナイキスト線図は不安定な極を持つ系でも適用できる汎用的な安定判別法です。実用的には多くのケースでボード線図が使われます。

dBとは何ですか?

デシベル(dB)= 20log₁₀(ゲイン比) で電圧・速度などの振幅量を表します。ゲイン2倍≈6dB、10倍=20dB、0.1倍=-20dB。電力比の場合は10log₁₀が使われます。

CAE(FEM)解析との関係は?

CAEで生成した柔軟体(フレキシブル構造)の周波数応答関数(FRF)をボード線図として表示し、制御系設計に組み込む「柔軟体モーダル制御」という手法があります。ロボット・宇宙構造物の制御設計で使われます。

ボード線図シミュレーターとは

伝達関数 \( G(s) = \frac{K \omega_n^2}{s^2 + 2 \zeta \omega_n s + \omega_n^2} \) を基本モデルとし、比例ゲイン \( K \)、固有角周波数 \( \omega_n \)、減衰比 \( \zeta \) をスライダーで調整可能です。この2次系に加え、1次遅れ要素 \( \frac{1}{Ts+1} \) やむだ時間要素 \( e^{-Ls} \) を直列接続でき、複合的な周波数応答を解析します。ゲイン線図は \( 20 \log_{10}|G(j\omega)| \) を、位相線図は \( \angle G(j\omega) \) をプロットし、ゲイン交差周波数 \( \omega_{gc} \) と位相交差周波数 \( \omega_{pc} \) を自動検出します。例えば \( K=10, \zeta=0.3 \) では、ゲイン余裕 \( G_m = -20 \log_{10}|G(j\omega_{pc})| \) が約6dB、位相余裕 \( \phi_m = 180^\circ + \angle G(j\omega_{gc}) \) が約30°となり、安定性の指標を直読できます。これにより、制御系設計におけるパラメータ調整と安定判別の因果関係を直感的に理解できます。

よくある質問

シミュレーター上でゲイン線図が0dBを横切る周波数(ゲイン交差周波数)と、位相線図が-180°を横切る周波数(位相交差周波数)が自動検出されます。そのときのゲイン余裕(位相交差周波数でのゲインの0dBからの差)と位相余裕(ゲイン交差周波数での位相の-180°からの差)が表示されるので、数値を確認して安定性を判断してください。
むだ時間要素 e^{-Ls} のゲインは周波数によらず一定(0dB)ですが、位相は周波数に比例して遅れます。そのため、他の要素と組み合わせると高周波域で位相が急激に回転し、ゲイン交差周波数付近で位相余裕が減少します。これにより安定性が低下し、ゲイン線図自体は振動しませんが、位相線図が周期的に変化する様子が観察できます。
Kを大きくするとゲイン線図が全体的に上昇し、ゲイン交差周波数が高周波側に移動します。これにより位相余裕が減少し、システムが不安定になりやすくなります。逆にKを小さくすると位相余裕が増加し安定性は向上しますが、定常偏差が大きくなるトレードオフがあるため、適切な値をスライダーで探索してください。
減衰比ζが小さいと、固有角周波数ω_n付近で共振が発生し、ゲイン線図に鋭いピーク(共振ピーク)が現れます。これはシステムの極が虚軸に近づき、応答が持続振動的になるためです。ζ=0では理論上無限大のピークとなり不安定になります。シミュレーターでζを0.1程度まで下げると、その様子をリアルタイムで確認できます。

実世界での応用

産業での実際の使用例
自動車業界では、電動パワーステアリング(EPS)のモーター制御系設計に本シミュレーターが活用されています。伝達関数のゲイン余裕と位相余裕をリアルタイムで確認しながら、ハンドル操作に対する応答の安定性を調整。例えば、トヨタや日産の車両開発現場で、EPSの共振周波数を回避するフィルタ設計に利用され、実車試験前のパラメータ最適化を効率化しています。また、産業用ロボットのサーボモーター制御でも、負荷変動に対する位相遅れを視覚的に補償する用途で導入されています。

研究・教育での活用
大学の制御工学実験では、学生がボード線図の理論を直感的に理解するための教材として使用。例えば、PID制御の比例ゲインを変更した際のゲイン余裕の変化を即座に観察でき、安定判別(ナイキスト安定判別との関連)を体感的に学べます。研究分野では、ドローン姿勢制御系の位相補償設計において、シミュレーター上で位相進み補償器の効果を試行錯誤し、理論値と実機応答の差異を検証するケースが増えています。

CAE解析との連携や実務での位置付け
本ツールは、MATLAB/SimulinkやANSYSなどのCAEツールで得られた伝達関数モデルを即座にインポートし、安定性評価を行うブリッジ的役割を担います。例えば、構造解析で得られた機械共振の伝達関数を本シミュレーターに読み込み、制御系設計者がゲイン余裕を6dB以上確保するようフィルタを調整。実務では、試作前の制御パラメータ決定プロセスに組み込まれ、ハードウェア試験の回数を削減する標準ワークフローの一部として定着しています。

よくある誤解と注意点

「ゲインが0dBを下回っていれば安定」と思いがちですが、実際はゲイン余裕と位相余裕の両方を確認する必要があります。ゲイン線図が0dBを下回っていても、位相が-180度に達する周波数でのゲイン値が負であれば安定ですが、位相余裕が小さすぎると、現実のシステムではノイズや非線形性により不安定化するリスクがあります。

「位相余裕が正なら必ず安定」と誤解されがちですが、実際には非最小位相系やむだ時間を含むシステムでは、ゲイン線図と位相線図の対応が複雑になるため注意が必要です。特に、ゲインが0dBを横切る周波数が複数ある場合、位相余裕の定義が曖昧になり、安定判別を誤る可能性があります。

「シミュレーター上で安定だから実機でも安定」と思いがちですが、実際の制御系にはセンサノイズ、アクチュエータの飽和、温度ドリフトなどシミュレーションで考慮されていない要素が存在します。ボード線図は線形近似であり、実機の非線形特性や高周波域の寄生要素が安定性に影響を与える点に注意が必要です。

使い方ガイド

  1. 比例ゲインK、時定数τ、二次系ゲインK2、固有角周波数ωnをスライダーまたは数値入力で設定する
  2. 伝達関数G(s)=K/(τs+1)×K2/(s²+2ζωns+ωn²)のボード線図をリアルタイム生成
  3. ゲイン線図(dB)と位相線図(度)から安定判別、ゲイン余裕、位相余裕を自動計算
  4. ナイキスト線図との連動で制御系の安定性を視覚的に確認

具体的な計算例

モータ速度制御系でK=5、τ=0.1秒、K2=2、ωn=10rad/sを設定。100Hzクロスオーバー周波数でゲイン余裕GM=6dB、位相余裕PM=45度を確認。鋼製歯車減速機(慣性モーメント0.05kg·m²)の場合、開ループゲイン20dBで過度応答時間0.4秒以下を達成。逆に位相余裕が20度未満ではリミットサイクル発生リスク増加。

実務での注意点