クラマース・クローニッヒの関係式

提供: testwiki
2023年7月26日 (水) 06:58時点におけるimported>Kenlo Nasaharaによる版 (クラマース・クローニッヒの関係式: 本文との整合性のため, 式中の下付き文字のRとIをローマンにした。)
(差分) ← 古い版 | 最新版 (差分) | 新しい版 → (差分)
ナビゲーションに移動 検索に移動

テンプレート:出典の明記 クラマース・クローニッヒの関係式(—かんけいしき、テンプレート:Lang-en-short)とは、線形応答における周波数応答関数の実部虚部ヒルベルト変換で関係づけられていることを示した式である。 1926年ラルフ・クローニッヒ1927年ヘンリク・アンソニー・クラマースによって電磁波分散現象に対して導かれた。

クラマース・クローニッヒの関係式

周波数応答関数H(ω)=HR(ω)+i HI(ω)に対して(ただし、HRHの実部、HIHの虚部である。)

HR(ω)=1π𝒫HI(ω)ωωdωHI(ω)=1π𝒫HR(ω)ωωdω

がクラマース・クローニッヒの関係式である。(𝒫コーシーの主値をとることを表す。)

後述するインパルス応答h(t) が恒に実数であるという条件を付けると、周波数応答関数の実部は偶関数、虚部は奇関数になる。 これを用いて積分範囲を正の部分にするようにクラマース・クローニッヒの関係式を変形すると

HR(ω)=2π𝒫0ωHI(ω)ω2ω2dωHI(ω)=2π𝒫0ωHR(ω)ω2ω2dω

となる。

因果律からの導出

クラマース・クローニッヒの関係式は、刺激よりも前に応答は起こりえないという因果律から導かれる。

線形応答においてはt=0におけるインパルスδ(t)に対する応答h(t)が決まれば、任意の刺激に対する応答が決定される。h(t)を偶関数he(t)と奇関数ho(t)の和の形 テンプレート:Indent に分解すると、因果律よりt<0でh(t)=0なのでho(t)=he(tsgn(t)he(t)=ho(t)·sgn(t)となる。

インパルス応答をフーリエ変換して周波数応答関数を求めると、 テンプレート:Indent となり、偶関数部he(t)のフーリエ変換は周波数応答関数の実部、奇関数部ho(t)のフーリエ変換は周波数応答関数の虚部にあたることが分かる[1]

それぞれに対して積関数のフーリエ変換が畳み込みになることを使えば、クラマース・クローニッヒの関係式が導かれる。ここで、sgn^(ω)は符号関数のフーリエ変換を表す。 テンプレート:Indent

複素関数を用いた導出

クラマース・クローニッヒの関係式の導出における複素積分

またH(ω)複素平面解析接続した複素関数H(z) が、実軸より上側で正則かつ|z|→∞ で一様にH(z)→0 であるとき[2]にはH(ω) がクラマース・クローニッヒの関係式を満たすことを示すことができる。

H(z)/(z-ω) を複素平面上で、以下の4つの区間からなる閉曲線上で複素積分する。

  1. 実軸上の(-R, 0)→(ω - r, 0)
  2. (ω, 0)を中心とする半径r の半円(ω - r, 0)→(ω + r, 0)
  3. 実軸上の(ω + r,0)→(R, 0)
  4. 原点を中心とする半径Rの半円(R, 0)→(-R, 0)

実軸より上側で正則であるという条件から、コーシーの積分定理によりこの閉曲線上の積分は0になる。 ここでR→∞r→0の極限をとると区間4の積分は|z|→∞で一様にH(z)→0の条件より0となる。 区間2の積分はr→0で-iπH(ω)となる。 したがって区間1と3の積分の和はR→∞、r→0の極限で テンプレート:Indent この式の実部と虚部を比較することでクラマース・クローニッヒの関係式が導出される。

応用

クラマース・クローニッヒの関係式を用いることで周波数応答関数の実部か虚部の一方からもう一方を計算で求めることが可能になる。これをクラマース・クローニッヒ解析という。

  • 複素磁化率については、フーリエ変換核磁気共鳴分光法で応用されている。インパルス応答であるFIDn 点のサンプリングデータから周波数応答関数(スペクトル)を離散フーリエ変換で求めた場合、それぞれ独立に実部と虚部のスペクトルデータがn 点得られる。しかし実際にはクラマース・クローニッヒの関係式により実部と虚部は独立ではないので、虚部の自由度を実部に移してスペクトルの分解能を2倍にすることが可能である。(ただしSN比はその分低下する。)このテクニックは時間的制約の大きい二次元NMRなどにおいてデータ点数を二倍にするゼロフィリングなどを用いて実行される。

脚注

テンプレート:Reflist

関連項目


テンプレート:Physics-stub

  1. ここではh(t)は実数であるものとしたが、h(t)が複素数の場合は、 HR=heR+ihoI^,iHI=hoR+iheI^ となり、同様に導出できる。
  2. 実はこの条件は因果律と同値である。実際、t<0でh(t)=0のとき、z=|z|eとして、 テンプレート:Indent であり、逆に、 テンプレート:Indent これが任意の0<θ<πについて|z|→∞で0になるためにはt<0でh(t)=0でなければならない。