固有角周波数ωnと減衰比ζをスライダーで調整し、2次遅れ系のステップ応答をリアルタイム可視化。過減衰・臨界減衰・不足減衰の3ケースを同時比較し、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}$ （不足減衰時）
このシミュレーターの核心は、2次遅れ系の伝達関数です。入力（目標値）と出力（応答）の動的関係をラプラス変換したs領域で表しています。
$$G(s)=\frac{Y(s)}{U(s)}=\frac{K\omega_n^2}{s^2+2\zeta\omega_n s+\omega_n^2}$$$G(s)$: 伝達関数
$Y(s)$: 出力のラプラス変換
$U(s)$: 入力（ステップ入力）のラプラス変換
$K$: 定常ゲイン（最終的に到達する値の比率）
$\omega_n$: 固有角周波数 [rad/s]（系の「固さ」や「速さ」を決める）
$\zeta$: 減衰比（振動の収まり具合を決める無次元数）
$s$: ラプラス演算子
伝達関数の分母多項式（特性方程式）の根を「極」と呼び、s平面上の位置が応答の性質を直接決定します。これはシミュレーターのs平面プロットに対応しています。
$$s^2+2\zeta\omega_n s+\omega_n^2 = 0 \quad \Rightarrow \quad s = -\zeta\omega_n \pm j\omega_n\sqrt{1-\zeta^2}$$極の実部 $-\zeta\omega_n$: これが負（左半平面）であることが安定条件。絶対値が大きいほど減衰が速い。
極の虚部 $\pm \omega_n\sqrt{1-\zeta^2}$: $\zeta<1$（不足減衰）の時に存在し、振動の角周波数 $\omega_d$ を表す。虚部の絶対値が大きいほど振動が速い。
この式から、$\zeta$と$\omega_n$を変えることが、極の位置＝システムの挙動を直接変える操作であることがわかります。
産業用ロボット・サーボ制御：アームの位置決め制御に広く利用されます。ζとω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に設定しても、実際に動かしたら振動が出るなんてことはよくある話だね。
2次遅れ系の理解は、制御工学の枠を超えて、さまざまな物理現象の解析に直結するんだ。例えば機械工学の振動解析。バネ（剛性k）とダンパ（減衰c）と質量（m）から成る単振動系は、まさにこの2次系そのもの。ここでは$\omega_n = \sqrt{k/m}$, $\zeta = c / (2\sqrt{mk})$ と対応する。同じ式が電気工学のRLC直列回路にも現れるよ。インダクタL、抵抗R、キャパシタCの回路の電流応答は、Lを質量m、Rをダンパc、1/Cをバネkと見なせば、全く同じ形の微分方程式で記述できる。さらに構造物の耐震設計でも、建物の固有周期と減衰定数はこの概念で考えられる。このシミュレーターで「極の位置」をいじる操作は、機械ならダンパを交換し、電気なら抵抗値を変え、建築なら制震ダンパーを調整する作業に、それぞれ対応しているんだ。
このツールで遊んで感覚が掴めたら、次のステップに進もう。まずは数学的な背景の深堀りがおすすめ。ラプラス変換でなぜ代数方程式になるのか、インパルス応答や周波数応答（ボード線図）との関係はどうなっているのか、を理解すると視野が広がる。次に「1次遅れ系」や「むだ時間要素」との組み合わせを学ぶと実務に近づく。多くの実システムは、2次系だけじゃなくてこれらの要素が直列になった高次系だからね。最終的にはフィードバック制御系の設計を目標にしよう。このシミュレーターでいじっているのは「制御対象」の部分だけ。ここにPID制御器を追加して閉ループ系を組むと、目標値の変更や外乱に対する応答をどう設計するか、という本当の制御設計の世界に入れる。その時、このツールで培った「極の位置と応答の関係」の直感が、非常に強力な武器になるんだ。