マシュー方程式 $\ddot{x}+ (\delta + \varepsilon\cos 2t)x = 0$ の安定領域チャートと時刻歴応答をリアルタイム計算。δとεを操作して不安定領域を体験しましょう。
赤：不安定領域 / グレー：安定領域 / 黄色十字：現在の(δ,ε)
機械・構造物の設計：回転軸や歯車など、回転に伴って剛性が周期的に変化する部品では、パラメトリック共振による予期せぬ大振幅振動（ねじれ振動など）が破損の原因となります。設計段階でδ-εパラメータを安全領域に収めることが重要です。
マイクロ・ナノスケールのデバイス：MEMS（微小電気機械システム）の共振子など、高周波で動作するデバイスでは、意図的にパラメトリック励振を用いて感度を高めたり、逆にその発生を抑制する制御が研究されています。
土木構造物：風や水流によって橋梁の空力特性（実効的な剛性）が周期的に変調され、パラメトリック励起が起こる可能性があります。タコマナローズ橋の悲劇的な崩壊には、このメカニズムも関与していたと考えられています。
量子光学：光共振器内の電場や、イオントラップ中の荷電粒子の運動も、パラメータが時間変調する系として記述され、パラメトリック共振の概念が応用されています。
このシミュレーターを使い始めるとき、いくつかつまずきやすいポイントがあるよ。まず、「εが大きいほど必ず不安定になる」と思いがちだけど、そう単純じゃないんだ。δ=1のくさびの先端付近では、εが小さいと確かに不安定だけど、εをもっと大きくしていくと、逆に安定領域が現れる「島」が存在する場合がある。例えば、δ=1、γ=0でεを0から2までゆっくり動かしてみて。不安定→安定→不安定と、交互に変わるのが見えるはず。これは、変調が激しすぎるとかえってエネルギー注入のタイミングが合わなくなるからなんだ。
次に、シミュレーション上の「発散」と現実の「破損」はイコールではないということ。このツールは線形モデルだから、発散は理論上無限大まで成長することを示す。でも実際の構造物には非線形性（例えば、バネが限界以上に伸びないとか）が必ずあり、振幅はある有限値で落ち着くか、別の破壊モードに至る。シミュレーションで赤い領域に入ったから即NGと判断せず、「ここは注意深く非線形応答を調べる必要がある危険ゾーン」と捉えよう。
最後にパラメータ設定のコツ。「γ（減衰）を0に近い値でシミュレーションするのは、現実的ではないことが多い」ことを覚えておいて。多くの機械構造物には必ず何らかの減衰がある。研究で原理を理解するならγ=0もありだけど、実務を想定するなら、材料や結合部から見積もった現実的な減衰係数（例えば鋼構造物で0.01〜0.05程度）を設定してからチャートを見る癖をつけよう。そうしないと、必要以上に保守的（重く、コスト高）な設計になってしまうからね。
梁の横揺れパラメトリック共振（f_0=2Hz、m=50kg）でδ=0.8、ε=0.6、減衰0.1、初期変位2mmを設定した場合、フロケ指数λが0.15/sとなり不安定領域に入る。10秒の応答計算では振幅が2mmから約45mmまで指数増幅され、マシュー方程式d²x/dt²+2ζω₀(dx/dt)+(ω₀²(1+εcos(Ωt)))x=0の予測と一致する。減衰を0.25に増加させると同一パラメータで安定化し、振幅は20mm以下に抑制される。