フレネル積分
フレネル積分とはなにか？
フレネル積分(Fresnel integrals)とは、次の積分によって表される2つの関数である。
$$ S(x) = \int_{0}^{x}{\sin{(t^2)}dt} $$
$$ C(x) = \int_{0}^{x}{\cos{(t^2)}dt} $$
変数の定数倍の違いであるが、次の定義も使われる。
$$ S(x) = \int_{0}^{x}{\sin{(\frac{\pi}{2}t^2)}dt} $$
$$ C(x) = \int_{0}^{x}{\cos{(\frac{\pi}{2}t^2)}dt} $$
$\sin{(x^2)}$と$\cos{(x^2)}$のグラフ
まず、それぞれの被積分関数である$\sin{(x^2)}$と$\cos{(x^2)}$の様子を確認しておこう。
いずれも偶関数であり、0付近では緩やかな様子だが、両端に近づくにつれて振動が激しくなる (周期が短くなる) 事がわかる。
$S(x)$と$C(x)$のグラフ
次に、フレネル積分の様子をグラフで確認する。
いずれも奇関数であり、有界である。
両端では振動しているが、その振動具合は少しずつ小さくなっており、$x \to \infty$において収束すると考えられる。
$(C(t), S(t))$のグラフ
$x = C(t), y = S(t)$としたときの軌跡もグラフで確認しておこう。
$t \to \pm \infty$において、無限回回転しながら、一点に収束すると考えられる。
この軌跡はオイラー螺旋(Euler spiral)と呼ばれ、またクロソイド(clothoid)もしくは、コルニュ螺旋(Cornu spiral)とも呼ばれる。
ギリシャ神話で糸巻き棒のシンボルで表される女神クロートーに因んでクロソイドという名を付けたのは、チェザロである。
SciPyにおけるフレネル積分
上の図は、SciPyライブラリのscipy.special.fresnelを用いて描画した。
この関数は、冒頭に述べた2番目の定義を実装している。
https://docs.scipy.org/doc/scipy/reference/generated/scipy.special.fresnel.html
上の図は適当な変換をして描画したものである。
フレネル積分の極限
図が示唆するように、フレネル積分は$x \to \infty$において収束する。この極限を求めよう。
$$ \int_{0}^{\infty}{\sin{(t^2)}dt} $$
$$ \int_{0}^{\infty}{\cos{(t^2)}dt} $$
を個別に求める代わりに、
$$ \int_{0}^{\infty}{\{\cos{(t^2) + i \sin{(t^2)}\}}dt} = \int_{0}^{\infty}{e^{it^2}}dt $$
を考える。
教科書的な解法
$R$を正の実数として、次のような積分路を考える。
$$ C_1: 0 \to R $$
$$ C_2: R \to R \exp{[\frac{\pi}{4}i]} $$
$$ C_3: R \exp{[\frac{\pi}{4}i]} \to 0 $$
$$ C = C_1 + C_2 + C_3 $$
ここで、$C_1$は実軸上を直線的に進み、$C_2$は8分の1回転を円周状に進む経路であり、$C_3$はそこから原点に直線的に戻る経路である。
また、各経路に沿った$e^{iz^2}$の複素積分をそれぞれ、$I_1$, $I_2$, $I_3$とする。求めたいのは$I_1$の$R \to \infty$における極限である。
さて、$e^{iz^2}$は正則であり、$C$は複素平面上の単純閉曲線であるから、Cauchyの積分定理により、
$$ I = I_1 + I_2 + I_3 = 0 $$
である。$I_2$および、$I_3$を計算しよう。結論を先に述べると、$I_2$は$0$に収束し、$I_3$はガウス積分となる。
$$ I_2 = \int_{C_2}e^{iz^2}dz = \int_{\theta=0}^{\frac{\pi}{4}} {\exp{[i (R e^{i\theta})^2]} \cdot iRe^{i\theta}d\theta} $$
$$ |I_2| \le \int_{\theta=0}^{\frac{\pi}{4}} {|\exp{[i (R e^{i\theta})^2]} \cdot iRe^{i\theta}|d\theta} = R \int_{\theta=0}^{\frac{\pi}{4}} {|\exp{[i (R e^{i\theta})^2]}|d\theta} = R \int_{\theta=0}^{\frac{\pi}{4}} {\exp{[\Re{(i (R e^{i\theta})^2)}]}d\theta} = R \int_{\theta=0}^{\frac{\pi}{4}} {\exp{[-R^2 \sin{2 \theta}]}d\theta} $$
右辺の積分区間において、$\sin{2\theta}$は上に凸であり、したがって
$$ \sin{2 \theta} \ge \frac{4}{\pi}\theta $$
が成り立つ。そこで、
$$ |I_2| \le R \int_{\theta=0}^{\frac{\pi}{4}} {\exp{[-R^2 \sin{2 \theta}]}d\theta} \le R \int_{\theta=0}^{\frac{\pi}{4}} {\exp{[-R^2 \frac{4}{\pi}\theta]}d\theta} = R \left[ -\frac{\pi}{4R^2} \exp{[-\frac{4R^2}{\pi}\theta]} \right]^{\frac{\pi}{4}}_{0} = \frac{\pi}{4} \frac{1 - \exp{[-R^2]}}{R} $$
となるが、右辺は$R \to \infty$において$0$に収束する。
したがって、$I_2 \to 0$である。
次に、
$$ I_3 = \int_{C_3}e^{iz^2}dz = \int_{r=R}^{0} {\exp{[i (r e^{i\frac{\pi}{4}})^2]} \cdot e^{i\frac{\pi}{4}}dr} = - \int_{r=0}^{R} {\exp{[i (r^2 e^{i\frac{\pi}{2}})]} \cdot (\frac{1}{\sqrt{2}} + \frac{1}{\sqrt{2}}i)dr} = - (\frac{1}{\sqrt{2}} + \frac{1}{\sqrt{2}}i) \int_{r=0}^{R} {\exp{[- r^2]} dr} $$
となり、右辺はガウス積分であるから、$R \to \infty$において
$$ I_3 \to - (\frac{1}{\sqrt{2}} + \frac{1}{\sqrt{2}}i) \frac{\sqrt{\pi}}{2} $$
に収束する。
以上をまとめると、
$$ I = I_1 + I_2 + I_3 = 0 $$
$$ I_2 \to 0 $$
$$ I_3 \to - (\frac{1}{\sqrt{2}} + \frac{1}{\sqrt{2}}i) \frac{\sqrt{\pi}}{2} $$
したがって、
$$ I_1 \to (\frac{1}{\sqrt{2}} + \frac{1}{\sqrt{2}}i) \frac{\sqrt{\pi}}{2} = \int_{0}^{\infty}{\{\cos{(t^2) + i \sin{(t^2)}\}}dt} $$
両辺の実部と虚部を比較して、
$$ \lim_{x \to \infty} S(x) = \int_{0}^{\infty}{\sin{(t^2)}dt} = \frac{1}{2} {\sqrt{\frac{\pi}{2}}} $$
$$ \lim_{x \to \infty} C(x) = \int_{0}^{\infty}{\cos{(t^2)}dt} = \frac{1}{2} {\sqrt{\frac{\pi}{2}}} $$
これはおおよそ、0.626657という値である。
先のオイラー螺旋のグラフを改めて見返すと、実際、最終的にその辺りの値に収束しそうである。
変数変換による方法
上はよく教科書に載っているが、積分路がやや唐突であるのは否めないので、こういった積分路が出てくる事情をもうすこし考えてみよう。
求めたい積分
$$ \int_{0}^{\infty}{e^{it^2}}dt $$
と、既に知っているガウス積分
$$ \int_{0}^{\infty}{e^{-t^2}}dt = \frac{\sqrt{\pi}}{2} $$
は非常によく似ているので、「変数変換によってガウス積分に帰着できそうである」と考える。
そこで次のように変形する。
$$ \int_{0}^{\infty}{e^{it^2}}dt = \int_{0}^{\infty}{\exp{[it^2]}}dt = \int_{0}^{\infty}{\exp{[(-1)(-i)t^2]}}dt = \int_{0}^{\infty}{\exp{[-e^{-\frac{\pi}{2}i}t^2]}}dt = \int_{0}^{\infty}{\exp{[-(e^{-\frac{\pi}{4}i}t)^2]}}dt $$
そこで、変数変換$s = e^{-\frac{\pi}{4}i} t$を考えると、
$$ \int_{0}^{\infty}{e^{it^2}}dt = \int_{0}^{e^{-\frac{\pi}{4}i} \infty}{\exp{[-s^2]}} e^{\frac{\pi}{4}i}ds = e^{\frac{\pi}{4}i} \int_{0}^{e^{-\frac{\pi}{4}i} \infty}{\exp{[-s^2]}}ds $$
となり、複素平面を右下に向かう積分路の広義積分が現れる。
これを求めるために、先ほどと同じように扇形を一周するような積分路でガウス積分を考えてやると、 8分の1回転の部分は0に収束するので結局実軸上でガウス積分をしているのと同じことだとわかる。
そこで、
$$ \int_{0}^{\infty}{e^{it^2}}dt = e^{\frac{\pi}{4}i} \int_{0}^{e^{-\frac{\pi}{4}i} \infty}{\exp{[-s^2]}}ds = e^{\frac{\pi}{4}i} \frac{\sqrt{\pi}}{2} = (\frac{1}{\sqrt{2}} + \frac{1}{\sqrt{2}}i) \frac{\sqrt{\pi}}{2} $$
と積分が求まる。
先に述べた方法で8分の1回転の積分路が出てくる事情はこのようなところにある。
オイラー螺旋の弧長
オイラー螺旋
$$ (x, y) = (C(t), S(t)) $$
の速度ベクトルは、
$$ (\dot{x}, \dot{y}) = (\cos(t^2), \sin(t^2)) $$
である。これの大きさは1であるから、オイラー螺旋は、単位時間あたりに1だけ進むと考えられる。
したがって、時刻$0$からのオイラー螺旋の弧長は
$$ L = t $$
である。
オイラー螺旋の曲率
また、加速度ベクトルは、
$$ (\ddot{x}, \ddot{y}) = (-2t\sin(t^2), 2t\cos(t^2)) $$
なので、曲率半径は、
$$ R = \frac{(\dot{x}^2 + \dot{y}^2)^\frac{3}{2}}{\dot{x}\ddot{y} - \dot{y}\ddot{x}} = \frac{1}{2t} $$
曲率は、
$$ \kappa = \frac{1}{R} = 2t $$
である。
補足
フレネル積分の極限は、少しテクニカルな方法で求めた。 ポイントはガウス積分を見越して8分の1回転の積分路を作るところである。 もう少し別の方法としては、メリン変換(Mellin transform)による方法がある。
オイラー螺旋には、弧長と曲率に比例関係があるという面白い性質がある。 この性質は工学で利用されることがある。 具体的には、直線路からある曲率のカーブに侵入するとき、いきなりカーブに接続してしまうと、曲率の変化率が大きく不快である。 車の運転で言えば、急ハンドルを切らなければならない。 そこで、オイラー螺旋を経由して曲率を少しずつ上げながら(ハンドルを少しずつ回しながら)カーブに接続するという方法が考えられている。 このような「つなぎ」の曲線は緩和曲線(transition curve / spiral easement)と呼ばれている。