制御系ステップ応答シミュレーター 戻る
制御工学シミュレーター

制御系ステップ応答・伝達関数シミュレーター

固有角周波数ωnと減衰比ζをスライダーで調整し、2次遅れ系のステップ応答をリアルタイム可視化。過減衰・臨界減衰・不足減衰の3ケースを同時比較し、s平面の極を表示。

システムパラメータ
ωn — 固有角周波数 (rad/s)
rad/s
ζ — 減衰比
K — ゲイン
ライブ数値(追跡中)
0.00
時刻 t [s]
0.000
出力 y(t)
0.0
オーバーシュート %
立上り時間 [s]
整定時間 [s]
リアルタイム・ステップ応答アニメーション
0.00 s
出力 y(t) 目標値(ステップ入力) ±2% 整定帯 頂点 整定時刻
計算結果(最終値)
立上り時間 (s)
オーバーシュート %
整定時間 2% (s)
減衰固有角振動数 ωd
定常ゲイン K
s平面(極配置)
理論・主要公式

$$G(s)=\frac{K\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}$$

$\omega_d = \omega_n\sqrt{1-\zeta^2}$ (不足減衰時)、$M_p=e^{-\zeta\pi/\sqrt{1-\zeta^2}}$、$t_s\approx 4/(\zeta\omega_n)$

検証例:$\zeta=0.5,\ \omega_n=2$ → $M_p\approx16.3\%$, $t_s\approx4.0$ s。

制御系ステップ応答・伝達関数シミュレーターとは

🙋
「減衰比ζ」って何ですか?教科書で見たけど、具体的に応答がどう変わるのかイメージがわかないんです。
🎓
大まかに言うと、振動の収まり具合を決めるパラメータだよ。このシミュレーターで、上のζのスライダーを0.2くらいにしてみて。グラフがガタガタ振動するでしょ?これが「不足減衰」。次にζを1.5にすると、振動せずにゆっくり上がる「過減衰」になる。操作してみると一発で理解できるよ。
🙋
え、そうなんですか!確かに変えました。ζが1の時は一番スムーズに上がってますね。これが「臨界減衰」ってやつですか?実務ではどの値がよく使われるんですか?
🎓
その通り!ζ=1が臨界減衰だ。実務では、ζ=0.7では理論上約4.6%のオーバーシュートで、速さと安定のバランスが良いんだ。ωnのスライダーもいじると、応答の速さがどう変わるか体感できるよ。
🙋
s平面に表示されてる「×」印(極)の位置と、グラフの形が連動して動いてますね。これってどういう関係があるんですか?
🎓
良いところに気づいたね!極の実部(左右の位置)が応答の減衰の速さ、虚部(上下の位置)が振動の速さを決めてるんだ。ζを小さくすると極が上に移動して虚部が大きくなり、グラフの振動が速くなる。ωnを大きくすると極が原点から遠ざかり、全体の応答が速くなる。シミュレーターでパラメータを変えながら極の動きとグラフを見比べるのが、理解の近道だよ。

よくある質問

ブラウザのJavaScriptが有効になっているかご確認ください。また、一部の古いブラウザや拡張機能が原因で動作が遅くなる場合があります。最新のChrome、Edge、Firefoxなどを推奨します。ページをリロードしても解決しない場合は、キャッシュをクリアしてお試しください。
不足減衰(ζ<1)は振動しながら目標値に収束します。臨界減衰(ζ=1)は最も速くオーバーシュートなしで収束します。過減衰(ζ>1)はゆっくりと単調に目標値に達します。凡例の色とs平面の極の位置(複素平面上の虚部の有無)でも判別できます。
すべての極が複素平面の左半面(実部が負)にあれば安定です。実部が負で虚部があると不足減衰(振動的)、実部のみで負なら過減衰または臨界減衰です。実部が正の極があると応答は発散し不安定になります。
あくまで理想的な2次遅れ系のモデルです。実際のシステムには非線形性やノイズ、高次の影響が含まれるため、定性的な傾向把握やパラメータ感度の確認にご利用ください。実機設計には別途詳細なモデリングと検証が必要です。

実世界での応用

産業用ロボット・サーボ制御:アームの位置決め制御に広く利用されます。ζとωnを調整することで、高速動作と振動抑制のトレードオフを設計します。オーバーシュートが大きすぎるとワークを損傷する恐れがあるため、シミュレーターで確認したようなパラメータチューニングが重要です。

自動車のサスペンション設計:車体とタイヤをバネとダンパでつないだモデルは、まさに2次遅れ系です。ζはダンパの減衰力に、ωnはバネの硬さに対応し、乗り心地(振動の収束性)と接地性(応答の速さ)を決めます。

航空機の姿勢制御(オートパイロット):機体のピッチやロールの角度を制御する際の基本モデルとして使われます。乱気流などの外乱に対するロバスト性(安定性)を確保するため、極が左半平面に適切に配置されていることが求められます。

電子回路のフィルタ設計:アクティブフィルタなどの伝達関数は2次系で表現されます。ζ(Q値と関連)によって、バターワース特性(ζ=1/√2)やチェビシェフ特性など、所望の周波数応答を実現するように設計されます。

よくある誤解と注意点

このシミュレーターを使い始めるとき、いくつかつまずきやすいポイントがあるよ。まず「ζを小さくすればするほど応答が速くなる」という誤解。確かにζ=0.2と0.7を比べると、最初の立ち上がりはζ=0.2の方が鋭いけど、振動が収まるまでの「整定時間」はかえって長くなるんだ。速さを求めるなら、ζ=0.7〜1.0の範囲でωnを大きくするのが正解。次にωnの値の解釈。ωn=10 [rad/s]だからと言って、振動周波数が10Hzになるわけじゃない。振動周波数は$\omega_d = \omega_n \sqrt{1-\zeta^2}$で決まるから、例えばζ=0.5なら実質的な振動周波数は約8.7 rad/sになる。あと、実務で一番気をつけるのはモデルと現実の乖離。このシミュレーターは理想的な2次系だけど、実際の機械には摩擦やバックラッシュ、アクチュエータの出力リミットがある。シミュレーションでζ=0.9に設定しても、実際に動かしたら振動が出るなんてことはよくある話だね。

使い方ガイド

  1. 自然角周波数Wnを0.1〜20 rad/sで設定し、応答速度を選びます。
  2. 減衰比Zeta(ζ)を0~2.0の範囲で変更して、過減衰(ζ>1)・臨界減衰(ζ=1)・不足減衰(ζ<1)の応答特性を比較する
  3. 定常ゲインGainを0.1〜5.0で設定し、立上り時間、オーバーシュート、2%整定時間の変化を確認します。
  4. s平面の極配置図を同時観察し、極の位置(実部=減衰の速さ、虚部=振動の速さ)と応答波形の対応を検証する

具体的な計算例

Wn=2.0rad/s、Zeta=0.5、Gain=1.0の不足減衰系を想定。ステップ応答では立上り時間≈0.8秒、オーバーシュート≈16.3%、整定時間2%≈4.0秒となる。減衰固有角振動数ωd=Wn√(1-ζ²)≈1.73rad/sで振動周期T=2π/ωd≈3.6秒。同じWn=2.0で減衰比をZeta=1.0(臨界減衰)に変更すると、振動が消滅し立上り時間≈1.7秒まで延伸するが、オーバーシュート0%を実現。

実務での注意点