フレネル積分

フレネル積分（フレネルせきぶん、英: Fresnel integrals）とは、オーギュスタン・ジャン・フレネルの名を冠した2つの超越関数 S(x) と C(x) であり、光学で使われている。近接場のフレネル回折現象を説明する際に現れ、以下のような積分で定義される。 S ( x ) = ∫ 0 x sin ⁡ ( t 2 ) d t , C ( x ) = ∫ 0 x cos ⁡ ( t 2 ) d t {\displaystyle S(x)=\int _{0}^{x}\sin(t^{2})\,dt,\quad C(x)=\int _{0}^{x}\cos(t^{2})\,dt} S(x) と C(x) をパラメトリック方程式として描画したものがクロソイド曲線である。 フレネル積分は、全ての x について収束する次の冪級数展開式で表せる。 S ( x ) = ∫ 0 x sin ⁡ ( t 2 ) d t = ∑ n = 0 ∞ ( − 1 ) n x 4 n + 3 ( 4 n + 3 ) ( 2 n + 1 ) ! {\displaystyle S(x)=\int _{0}^{x}\sin(t^{2})\,dt=\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{4n+3}}{(4n+3)(2n+1)!}}} C ( x ) = ∫ 0 x cos ⁡ ( t 2 ) d t = ∑ n = 0 ∞ ( − 1 ) n x 4 n + 1 ( 4 n + 1 ) ( 2 n ) ! {\displaystyle C(x)=\int _{0}^{x}\cos(t^{2})\,dt=\sum _{n=0}^{\infty }(-1)^{n}{\frac {x^{4n+1}}{(4n+1)(2n)!}}} Abramowitz and Stegun, (eqs 7.3.1 – 7.3.2) などの書籍では、S(x) と C(x) を定義する積分の冪指数として π 2 t 2 {\displaystyle {\frac {\pi }{2}}t^{2}} を使っている。同じ関数を得るには、積分式に 2 π {\displaystyle {\sqrt {\frac {2}{\pi }}}} をかけ、同じ係数で引数 x を割ればよい。 オイラーの螺旋は、コルニュ螺旋またはクロソイドとも呼ばれ、S(t) と C(t) をパラメトリックにプロットした曲線である。コルニュ螺旋はマリー・アルフレッド・コルニュが回折の計算用にノモグラムとして考案したものである。 ここで、 C ′ ( t ) 2 + S ′ ( t ) 2 = sin 2 ⁡ ( t 2 ) + cos 2 ⁡ ( t 2 ) = 1 {\displaystyle C'(t)^{2}+S'(t)^{2}=\sin ^{2}(t^{2})+\cos ^{2}(t^{2})=1\,} であるから、この曲線の接ベクトルは単位長で、t は原点 (0,0) からの曲線に沿った弧長である。したがって、どちらの方向の曲線も無限の長さがある。 この曲線は、任意の点の曲率が原点からの曲線に沿った長さに比例するという特徴がある。このため、高速道路や鉄道で緩和曲線として用いられる。この曲線上で乗り物が一定速度で走行すると、角加速度が一定のレートとなる。クロソイド曲線（の一部）は例えばローラーコースターのループの形状にも使われている。 C(x) と S(x) は x の奇関数である。 C と S は整関数である。 上述の冪級数展開式を使うと、フレネル積分は定義域を複素数に拡張でき、複素数値の解析関数となる。フレネル積分は誤差関数を使って以下のように表現できる。 S ( x ) = π 4 ( i erf ⁡ ( i x ) + − i erf ⁡ ( − i x ) ) {\displaystyle S(x)={\frac {\sqrt {\pi }}{4}}\left({\sqrt {i}}\,\operatorname {erf} ({\sqrt {i}}\,x)+{\sqrt {-i}}\,\operatorname {erf} ({\sqrt {-i}}\,x)\right)} C ( x ) = π 4 ( − i erf ⁡ ( i x ) + i erf ⁡ ( − i x ) ) {\displaystyle C(x)={\frac {\sqrt {\pi }}{4}}\left({\sqrt {-i}}\,\operatorname {erf} ({\sqrt {i}}\,x)+{\sqrt {i}}\,\operatorname {erf} ({\sqrt {-i}}\,x)\right)} C(x) と S(x) を定義している積分式は、特別な場合を除いては、初等関数を使って閉形式で評価することができない。x が無限大に漸近したときのこれらの関数の極限は次のようになることが知られている。 ∫ 0 ∞ cos ⁡ t 2 d t = ∫ 0 ∞ sin ⁡ t 2 d t = 2 π 4 = π 8 {\displaystyle \int _{0}^{\infty }\cos t^{2}\,dt=\int _{0}^{\infty }\sin t^{2}\,dt={\frac {\sqrt {2\pi }}{4}}={\sqrt {\frac {\pi }{8}}}} 引数が無限大に漸近したときの C と S の極限は、複素解析の手法で求められる。それには、正のx-軸、半直線 y = x, x ≥ 0、原点を中心とした半径 R の円で囲まれた複素平面での扇形の領域の境界線に沿って、次の関数の扇形積分を使う。 ∮ ⁡ d z   e i z 2 = 0 {\displaystyle \oint \mathrm {d} z~e^{iz^{2}}=0} この周回積分が0になるのは領域内に極がないためである。R が無限大の極限を考えると、 γ 2 {\displaystyle \gamma _{2}} （円弧部分）上の積分は0になり、 γ 3 {\displaystyle \gamma _{3}} 上の積分はガウス積分から ∫ γ 3 d z   e i z 2 = ∫ 0 R d t   1 + i 2 e − t 2 → π 8 ( 1 + i ) ( R → ∞ ) {\displaystyle \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 )} となる。よって、 γ 1 {\displaystyle \gamma _{1}} （実軸）上の積分の実部と虚部を取ることで、 C ( ∞ ) = S ( ∞ ) = π / 8 {\displaystyle C(\infty )=S(\infty )={\sqrt {\pi /8}}} が求められる。 オーギュスタン・ジャン・フレネル クロソイド曲線 回折レンズ Weisstein, Eric W. "Fresnel Integrals". MathWorld（英語）. CS1 maint: Multiple names: authors list Weisstein, Eric W. "Cornu Spiral". MathWorld（英語）. CS1 maint: Multiple names: authors list R. Nave, The Cornu spiral, Hyperphysics (2002) （t2 の代わりに πt2/2 を使っている） Milton Abramowitz and Irene A. Stegun, eds. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. New York: Dover, 1972. (See Chapter 7)