解析工学インタラクティブ・シミュレータ

数式をビジュアル化して直感的に学ぶ学習プラットフォーム

A

複素周波数解析

複素周波数 $s = \sigma + j\omega$ は、単なる一定の「周波数(振動)」だけでなく、時間の経過とともに減衰または増幅する挙動を包括して表すための重要な概念です。

数理的な意味と日常の例え

物理的なシステム(ブランコや車のサスペンション)を想像してください。

  • 実部 $\sigma$ (減衰因子): 空気抵抗や摩擦に相当します。 $\sigma < 0$ なら時間とともに摩擦で揺れが収まり(安定)、 $\sigma > 0$ ならエネルギーが注入されて揺れが暴走します(不安定)。
  • 虚部 $\omega$ (角周波数): 1秒間にどれだけ速く往復するかという「揺れの細かさ」を表します。

この解析で測定・予測できるもの

  • 過渡応答の収束時間: インパルスや段差(ステップ)による外部衝撃が加わった後、何秒で元の静止状態(安定状態)に収束するかを測定。
  • 減衰振動周波数: 摩擦を伴いながら「実際にシステムが往復する揺れのサイクル数(Hz)」を定量評価。
  • 共振に対する感度(Q値): 外部からの特定の周波数に対して、システムがどれくらい激しく応答(共振)するかを予測。
時間関数: $f(t) = e^{\sigma t} \cos(\omega t)$

インタラクティブ操作

右側の「$s$平面(複素平面)」上をクリック・ドラッグ、または下のスライダーを動かして、極(システムの特性を決める根)の位置を変更してください。応答波形がどう変化するかが分かります。

減衰率 $\sigma$ (実部) -0.20
周波数 $\omega$ (虚部) 4.00

$s$平面 (複素平面)

クリックして極を操作

時間応答波形 $f(t) = e^{\sigma t}\cos(\omega t)$

安定
B

根軌跡を用いた解析

自動制御系において、フィードバック増幅率(ゲイン) $K$ を $0$ から無限大まで変化させたとき、閉ループシステムの安定性を決める特性根(極)が $s$ 平面上でどう移動するかを視覚化したものです。

直感的な理解と「ゲイン」の意味

シャワーのお湯の温度調節を考えてみましょう。

  • ゲイン $K$ が適度(安定): 冷たければ少し温め、速やかにちょうど良い温度に落ち着きます。
  • ゲイン $K$ が大きすぎる(不安定): 「少し冷たい」に対して蛇口を思い切り熱い方に回し、今度は「熱すぎる」と冷たい方に全開で回す動作を繰り返します。これが「制御の発散(ハンチング)」です。

この解析で測定・予測できるもの

  • 安定限界ゲイン ($K_{\text{crit}}$): フィードバックループの「ボリューム(ゲイン)」をどこまで大きくすると、システムが異常振動を起こし制御不能(発散)になるか限界点を正確に予測。
  • 応答のオーバーシュート(行き過ぎ量): 目標値に対して、どれくらい飛び出してしまうかを予測し、適切なダンピング(減衰)比を逆算。
  • ステップ応答時間(応答速度): コントローラによる補正が動作し始めてから、最終定常状態の数パーセント以内に落ち着くまでの「整定時間」を予測。
シミュレーション対象(3次遅れ系):
開ループ伝達関数: $G(s) = \frac{K}{s(s+2)(s+5)}$
特性方程式: $s(s+2)(s+5) + K = 0$
安定限界ゲイン: $K = 70.0$

ゲイン K のコントロール

システムゲイン $K$ 20.00
安定領域です

ゲインが低いため、すべての極が左半面にあります。ステップ応答は最終的に一定値に収束します。

$s$平面上の根軌跡プロット

×: 開ループ極 / ⬤: 閉ループ極

システムのステップ応答

t=0 で入力が1にステップ変化
C

フーリエ級数解析

フーリエ級数とは、「どんな複雑な周期関数(波形)であっても、適切な周波数・振幅を持つ単純な正弦波(サイン波)の足し合わせとして分解・再現できる」という画期的な定理です。

数理的な仕組みと直感の例え

これは「音のイコライザー」「光を七色に分けるプリズム」と同じ原理です。

例えば、バイオリンの音(複雑な波形)は、基本となる低い音の波に、何倍もの細かい高調波(倍音)の波が幾層にも重なることで、あの特有の美しい音色を作っています。これを数式で表したのがフーリエ級数展開です。

この解析で測定・予測できるもの

  • 信号の周波数成分(スペクトル構造): マイクや電気回路から集音・計測した複雑な合成波形に、具体的に「どの周波数(何Hzの成分)」が「どの強度(dB)」で含まれているかを特定・可視化。
  • 高調波歪み率(THD): 入力された音が理想的な「歪みのない純音」からどの程度劣化、変形(歪み)しているかを検出・測定。
  • 不要な電気的・音響的ノイズの位置: 送電線から漏れる電源ノイズ(50Hz / 60Hzなど)や、不快なキーンという高周波不快音のスポットをピンポイントに検出し、排除(フィルタリング)する基準設計に使用。
方形波のフーリエ公式:
$f(t) = \frac{4}{\pi} \left( \sin(\omega t) + \frac{1}{3}\sin(3\omega t) + \frac{1}{5}\sin(5\omega t) + \dots \right)$
調和数(合成する波の数) $N$ 4

Nを増やすと、細部のギザギザや急峻な立ち上がりが精密に再現されます。

エピサイクル(回転ベクトルの連鎖)と波形合成アニメーション

左の円の軌跡が、右側に引き伸ばされて波形になります
D

ベッセル関数による解析

円柱(または円形)の形状を持つ構造体における物理現象(波動、熱伝導、電気等)を記述する微分方程式を解く際に現れるのがベッセル関数です。

「丸い太鼓」とベッセル関数の繋がり

直線的な「バイオリンの弦」の揺れ(1次元)はシンプルな $\sin$ 関数で表せます。

しかし、「丸い太鼓の皮(2次元の円形膜)」になると、振動は中心から円周方向へと広がり、かつ周囲は枠にガッチリと固定されている必要があります。この円形に制約された揺れ方を記述するために生まれたのがベッセル関数 $J_m(x)$ です。

この解析で測定・予測できるもの

  • 円筒形状の固有振動数(共鳴点): 丸いドラムの皮、スピーカーの振動板、光ファイバーの導波路、円筒形の排気管などで、どの高周波がエネルギーを損なわずに共振できるか予測可能。
  • 「節線・節円」(揺れないライン)の正確な位置: 振動板上の全く動かない輪状・線状の位置(節)をナノメートル〜ミリメートル単位で予測。これにより、マイクの固定箇所や、構造的な補強ポイントを設計。
  • 熱・電磁場の内部浸透深度: 同軸ケーブルなどのシリンダー内部で、電磁場や熱が中心部から外周に向けてどのように減衰(浸透)して分布しているかを予測。

境界条件: ドラムの最外周(半径 $R$)では、皮が固定されているため変位は必ず $0$ になります。すなわち、振動を維持できるのはベッセル関数の「零点(根)」に対応する特定の振動数のみです。

各パラメータの意味:
m (角方向の節の数): ドラムをケーキのように等分割する「節線(全く動かない直径直線)」の数。
n (半径方向の節の数): ドラムの内部に同心円状に生じる「節円(全く動かない円周)」の数。

第一種ベッセル関数 $J_m(x)$ のプロット

黄色の丸は根(ゼロ点)

ドラムヘッドの振動(3Dメッシュ)

リアルタイム変位表示

4つの解析手法の徹底比較表

紹介した4つの解析手法は、いずれも数学や工学分野における振動現象や動的システムを対象にしていますが、使用されるシチュエーションや適用されるアプローチ、数理的アプローチに明確な違いがあります。

解析手法 物理的・工学的意味 測定・予測・分析できる対象 数理背景 日常的な例え 実社会の応用例
A. 複素周波数解析 システムが持つ固有の周波数と、そのエネルギーが減少(減衰)または増加(発散)する速度を同時に定量評価する。 ・過渡変化時の収束時間
・実際の減衰抵抗(摩擦損失値)
・共鳴周波数に対する感度(Q値)
$s = \sigma + j\omega$
ラプラス変換
特性極の配置
摩擦のあるブランコの揺れが収まる速さ 電気回路の過渡特性設計、地震時のビル構造設計
B. 根軌跡解析 フィードバック制御系の「補正の強さ(ゲイン $K$)」を上げたときに、システム全体がどの瞬間に不安定化するかを軌跡で予測する。 ・限界増幅率(発振ハンチング閾値)
・目標への整定ステップ時間
・行き過ぎ率(オーバーシュート)
$1 + K \cdot G(s) = 0$
閉ループ特性方程式の軌跡
お風呂の温度を慌てて調整して逆に大失敗する様子 ロボットアームの軌道制御、エアコンの自動温度調整設計
C. フーリエ級数解析 どんなに複雑に見える規則正しい「繰り返し波形」でも、様々な波長を持った「きれいなサイン波」の集合体に分解して分析する。 ・波形中の正確な周波数構成(Hz/強度)
・歪みの発生度(全高調波歪み率:THD)
・カットすべきノイズ、共振源の特定
三角関数の直交性
無限和による級数展開
スペクトル
混ざりあった色の光をプリズムでバラバラの光線に分ける 音声ファイルの圧縮(MP3)、画像の周波数解析(JPEG)、オーディオイコライザー
D. ベッセル関数解析 「丸いドラムの皮」や「丸いパイプ内の水流」のように、極座標(円形・球形)の境界条件を持った空間上の波動・熱の伝わり方を解明する。 ・円柱形状での共鳴固有周波数
・振動板の「動かないノード(節)」の位置
・熱や高周波電磁場の同軸内部浸透深さ
ベッセルの微分方程式
第一種ベッセル関数 $J_m(x)$
円筒座標境界値問題
丸いドラムの叩き方で異なる複雑な固有の響き方 光ファイバー内の導波路設計、スピーカーコーンの振動解析、FMラジオの変調設計