フレネル積分

提供: miniwiki
移動先:案内検索
ファイル:Fresnel Integrals (Unnormalised).svg
S(x)C(x)C(x) の最大値は約 0.977451424。t2 の代わりに πt2/2 を使うと、図は水平および垂直方向に縮小される(下図)

フレネル積分(フレネルせきぶん、: Fresnel integrals)とは、オーギュスタン・ジャン・フレネルの名を冠した2つの超越関数 S(x) と C(x) であり、光学で使われている。近接場フレネル回折現象を説明する際に現れ、以下のような積分で定義される。

[math]S(x)=\int_0^x \sin(t^2)\,dt,\quad C(x)=\int_0^x \cos(t^2)\,dt[/math]

S(x) と C(x) をパラメトリック方程式として描画したものがクロソイド曲線である。

定義

ファイル:Fresnel Integrals (Normalised).svg
正規化したフレネル積分 S(x)C(x)。三角関数の引数を上図での t2 ではなく、πt2/2 にしている。

フレネル積分は、全ての x について収束する次の冪級数展開式で表せる。

[math]S(x)=\int_0^x \sin(t^2)\,dt=\sum_{n=0}^{\infin}(-1)^n\frac{x^{4n+3}}{(4n+3)(2n+1)!}[/math]
[math]C(x)=\int_0^x \cos(t^2)\,dt=\sum_{n=0}^{\infin}(-1)^n\frac{x^{4n+1}}{(4n+1)(2n)!}[/math]

Abramowitz and Stegun, (eqs 7.3.1 – 7.3.2) などの書籍では、S(x) と C(x) を定義する積分の冪指数として [math]\frac{\pi}{2}t^2[/math] を使っている。同じ関数を得るには、積分式に [math]\sqrt{\frac{2}{\pi}}[/math] をかけ、同じ係数で引数 x を割ればよい。

オイラーの螺旋

ファイル:Cornu Spiral.svg
オイラーの螺旋 (xy) = (C(t), S(t))。t が正および負の無限大に近づくと、曲線は2つの穴の中心に収束していく。

オイラーの螺旋は、コルニュ螺旋またはクロソイドとも呼ばれ、S(t) と C(t) をパラメトリックにプロットした曲線である。コルニュ螺旋はマリー・アルフレッド・コルニュが回折の計算用にノモグラムとして考案したものである。

ここで、

[math]C'(t)^2 + S'(t)^2 = \sin^2(t^2) + \cos^2(t^2) = 1 \,[/math]

であるから、この曲線の接ベクトル単位長で、t は原点 (0,0) からの曲線に沿った弧長である。したがって、どちらの方向の曲線も無限の長さがある。

この曲線は、任意の点の曲率が原点からの曲線に沿った長さに比例するという特徴がある。このため、高速道路や鉄道で緩和曲線として用いられる。この曲線上で乗り物が一定速度で走行すると、角加速度が一定のレートとなる。クロソイド曲線(の一部)は例えばローラーコースターのループの形状にも使われている。

属性

  • C(x) と S(x) は x奇関数である。
  • CS整関数である。
  • 上述の冪級数展開式を使うと、フレネル積分は定義域を複素数に拡張でき、複素数値の解析関数となる。フレネル積分は誤差関数を使って以下のように表現できる。
    • [math]S(x)=\frac{\sqrt{\pi}}{4} \left( \sqrt{i}\,\operatorname{erf}(\sqrt{i}\,x) + \sqrt{-i}\,\operatorname{erf}(\sqrt{-i}\,x) \right)[/math]
    • [math]C(x)=\frac{\sqrt{\pi}}{4} \left( \sqrt{-i}\,\operatorname{erf}(\sqrt{i}\,x) + \sqrt{i}\,\operatorname{erf}(\sqrt{-i}\,x) \right)[/math]
  • C(x) と S(x) を定義している積分式は、特別な場合を除いては、初等関数を使って閉形式で評価することができない。x が無限大に漸近したときのこれらの関数の極限は次のようになることが知られている。
    • [math]\int_{0}^{\infty} \cos t^2\,dt = \int_{0}^{\infty} \sin t^2\,dt = \frac{\sqrt{2\pi}}{4} = \sqrt{\frac{\pi}{8}}[/math]

評価

ファイル:Fresnel Integral Contour.svg
扇形の輪郭線を使い、フレネル積分の極限を計算する。

引数が無限大に漸近したときの CS の極限は、複素解析の手法で求められる。それには、正のx-軸、半直線 y = x, x ≥ 0、原点を中心とした半径 R の円で囲まれた複素平面での扇形の領域の境界線に沿って、次の関数の扇形積分を使う。

[math]\oint \mathrm{d}z ~ e^{iz^2} = 0[/math]

この周回積分が0になるのは領域内にがないためである。R が無限大の極限を考えると、[math]\gamma_2[/math](円弧部分)上の積分は0になり、[math]\gamma_3[/math]上の積分はガウス積分から

[math] \int_{\gamma_3} \mathrm{d}z ~ e^{iz^2} = \int_0^R \mathrm{d}t ~ \frac{1+i}{\sqrt{2}} e^{-t^2} \to \sqrt{\frac {\pi}{8}} (1+i) \quad (R \to \infty) [/math]

となる。よって、[math]\gamma_1[/math](実軸)上の積分の実部と虚部を取ることで、[math]C(\infty) = S(\infty) = \sqrt{\pi/8}[/math]が求められる。

関連項目

参考文献