インパルス応答のフーリエ変換は、時間領域の応答特性を周波数領域(周波数応答)に変換する操作であり、フィルタ設計・システム解析・音響評価において非常に重要な手法です。
インパルス応答から周波数応答(振幅特性・位相特性)を求めることで、システムが各周波数をどのように処理するかを視覚的に把握できます。
本記事では、インパルス応答とフーリエ変換の関係・周波数応答(H(jω))の意味・振幅特性と位相特性の解釈・スペクトル解析との関係について詳しく解説していきます。
目次
インパルス応答のフーリエ変換と周波数応答の定義
それではまず、インパルス応答のフーリエ変換と周波数応答の定義について解説していきます。
インパルス応答h(t)をフーリエ変換すると、システムの周波数応答H(jω)(または複素伝達関数)が得られます。
フーリエ変換による周波数応答の導出
インパルス応答→周波数応答の変換
H(jω) = F{h(t)} = ∫₋∞^∞ h(t) × e^(-jωt) dt
ω:角周波数(rad/s) = 2πf
H(jω)は複素数:H(jω) = |H(jω)| × e^(jθ(ω))
|H(jω)|:振幅特性(各周波数での利得)
θ(ω):位相特性(各周波数での位相変化)
ラプラス変換との関係
ラプラス変換はs = σ + jωという複素変数を使いますが、虚軸(σ = 0、s = jω)での評価がフーリエ変換に相当します。
システムが安定(全極が左半平面)なとき、ラプラス変換を虚軸上で評価することで周波数応答が得られます。
振幅特性・位相特性の読み方
続いては、周波数応答の振幅特性と位相特性の読み方について確認していきます。
ボード線図(Bode Plot)の構成
振幅特性と位相特性を周波数に対してプロットしたグラフをボード線図(Bode Plot)と呼びます。
横軸に周波数(対数スケール)・上部グラフに振幅(dB)・下部グラフに位相(度)を表示します。
ボード線図から読み取れる主要な情報は以下の通りです。
| 読み取る情報 | 内容 |
|---|---|
| カットオフ周波数 | 振幅が-3 dBになる周波数 |
| ロールオフ率 | 阻止域での減衰の急峻さ(dB/decade) |
| ゲイン余裕 | 位相が-180°のときの振幅余裕(安定性指標) |
| 位相余裕 | 振幅が0 dBのときの位相の-180°からの余裕(安定性指標) |
RCローパスフィルタの周波数応答の例
RCローパスフィルタの周波数応答
H(jω) = 1/(1+jωRC) = 1/(1+j(f/fc))
振幅:|H(jω)| = 1/√(1+(f/fc)²)
位相:θ = -arctan(f/fc)
f = fc での値:|H| = 1/√2 ≈ 0.707(-3 dB)、θ = -45°
f → ∞:|H| → 0(完全遮断)、θ → -90°
インパルス応答のフーリエ変換の実用的な応用
続いては、インパルス応答のフーリエ変換の実用的な応用について確認していきます。
音響インパルス応答の周波数解析
測定した音響インパルス応答h(t)をフーリエ変換することで、その空間(ホール・録音室)の周波数応答特性(どの周波数が強調・減衰されるか)が明らかになります。
この解析により、音響設計の問題点(特定周波数の過度な減衰・共振)を特定し、吸音材・拡散材の配置最適化に活用できます。
FFT(高速フーリエ変換)による実用計算
実際のデジタル信号処理では、測定されたインパルス応答データのFFT(Fast Fourier Transform)によって周波数応答を高速に計算します。
FFTは離散フーリエ変換(DFT)の高速計算アルゴリズムであり、N点のデータをO(N log N)の計算量で処理できます。
まとめ
本記事では、インパルス応答のフーリエ変換・周波数応答H(jω)の意味・振幅特性と位相特性の読み方・ボード線図・実用的な応用について詳しく解説しました。
インパルス応答h(t)をフーリエ変換すると周波数応答H(jω)が得られ、振幅特性|H(jω)|と位相特性θ(ω)に分解してシステムの周波数ごとの特性を解析できます。
ボード線図はカットオフ周波数・ロールオフ率・安定性余裕を視覚的に把握するための標準的なツールです。
FFTを活用したインパルス応答の周波数解析は音響設計・フィルタ設計・制御システム解析の実務において欠かせない技術となるでしょう。