記事ソース/フレネル回折からみたレンズの公式 のバックアップ(No.16)
記事ソース/フレネル回折からみたレンズの公式 のバックアップ(No.16)
これはrst2hooktailの記事ソース保存・変換用です（詳細）．
============================================================
フレネル回折からみたレンズの公式
============================================================
天体望遠鏡で大口径のレンズが望まれるのはなぜでしょうか．その理由はレンズの光学伝達関数
から説明できます．ここではフレネル回折の性質を用いた簡単な計算でレンズの公式を導き，さらに
いくつかの簡単な例でフーリエ光学の紹介をします．計算に作用素代数を用いているのが特徴です．
フレネル回折の簡易モデル
============================================================
$z$ 軸方向に進む平面波が $xy$ 平面内に置かれた開口を通過するときのフレネル回折の式
<tex>
u'(x', y') = \frac{A}{i \lambda R} e^{i k R} \int\mspace{-11mu}\int
f(x, y) e^{\frac{i k}{2 R} \{(x' - x)^2 + (y' - y)^2\}} dx dy
</tex>
から出発します[1]．ここで $u'(x', y')$ は点 $(x', y', R) (R > 0)$ における波面（のフェーザ），
$\lambda = \frac{2 \pi}{k}$ は波長， $f(x, y)$ は（開口部で 1，遮蔽部で 0 となる）開口関数です [*]_ ．
.. [*] $z$ 軸に垂直な平面波の波動は $A e^{i (\omega t - k z)}$ のように表されますが，
時間変化を表す $e^{i \omega t}$ は空間内で共通なので，通常は残りの部分（フェーザ）
だけを考えます．
$g(x, y) = e^{i \pi (x^2 + y^2) / \lambda R}$ とおくと，畳み込み積分
<tex>
(f * g)(x', y') = \int\mspace{-11mu}\int f(x, y) g(x' - x, y' - y) dx dy
</tex>
を使って $u'(x', y')$ を
<tex>
u'(x', y') = \frac{A}{i \lambda R} e^{i k R}(f * g)(x', y')
</tex>
と表すことができます．また
$g(x' - x, y' - y) = g(x', y') g(x, y) e^{i 2 \pi (x' x + y' y) / \lambda R}$ ですから
<tex>
u'(x', y') = \frac{A}{i \lambda R} e^{ikR} g(x', y') \int\mspace{-11mu}\int
f(x, y) g(x, y) e^{i 2 \pi (x' x + y' y) / \lambda R} dx dy
</tex>
が成立します．すなわち $u'(x', y')$ は $f(x, y) g(x, y)$ の２次元フーリエ変換を $\lambda R$
倍に拡大した関数に $\frac{A}{i \lambda R} e^{i k R} g(x', y')$ をかけたものに等しくなります．
このことは $f(x, y)$ が文字の形に切り抜かれた開口でも，さらに開口の直前に適当な光学素子を
おいたとして $f(x, y)$ を複素数値の関数 $u(x, y)$ で置換したときでも変わりません． ただし，
フレネル近似が成立するためには $u(x, y)$ の 0 でない部分は $R$ に比べて十分小さいことが
必要です．
人が知覚したり，媒体に記録するのは波面 $u'(x', y')$ そのものでなく，強度分布 $|u'(x', y')|^2$
です．上記の式に含まれる $e^{ikR}$ は $x$, $y$ に依存しない定数なので，フレネル回折の簡易モデル
として，以下では
<tex>
u'(x, y) = \frac{1}{i \lambda R} (u * g)(x, y)
</tex>
を考えましょう（後述の光学素子も同様にモデル化します）．
２次元フーリエ変換の公式
============================================================
点 $(x, y)$ での値が $u(x, y)$ である関数 $u$ を，点 $(x', y')$ での値が
<tex>
u'(x', y') = \int\mspace{-11mu}\int u(x, y) e^{i2 \pi (x x' + y y')} dx dy
</tex>
である関数 $u'$ に変換する写像 $F$ を（２次元の）フーリエ変換といい， $u' = F u$ とかきます．
$u'(x', y') = (F u)(x', y')$, $u'(x, y) = (F u)(x, y)$ です．一般に関数を関数に変換する写像を
作用素といいます．例えば微分作用素はある関数をその導関数に変換します．
$u$ を点 $(x, y)$ での値が $v(x, y) u(x, y)$ である関数や $u(x / a, y / a)$ である関数，
$u(x - a, x - b)$ である関数に変換する写像も作用素です．以下ではこれらの作用素を $\{ v \}$,
$M_a$, $S_{a,b}$ で表します．すなわち，
<tex>
(\{ v \} u)(x, y) = v(x, y) u(x, y)
</tex>
<tex>
(M_a u)(x, y) = u(x / a, y / a)
</tex>
<tex>
(S_{a,b} u)(x, y) = u(x - a, y - b)
</tex>
です．ここで， $v$ は任意の関数， $a$, $b$ は任意の定数で，慣例に従って
<tex>
(a u + b v)(x, y) = a u(x, y) + b v(x, y)
</tex>
と考えます．また，以下でよく使う次の関数もここで定義しておきましょう．
<tex>
\theta _a (x, y) = e^{i \pi a (x^2 + y^2)}
</tex>
<tex>
\phi _{a,b} (x, y) = e^{i2 \pi (a x + b y)}
</tex>
よく知られているように，フーリエ変換について次の公式が成立します．
1. $F(v * u) = (F{v})(F u)$ したがって $v * u = F^{-1} \{ F v \} F u$
2. $F S_{a,b} = \{ \phi _{-a, -b} \} F$, $S_{a,b} F = F \{ \phi _{a, b} \}$
3. $F M_a = a^2 M_{1 / a} F$
4. $F^2 = M_{-1}$ したがって $F^4 u = u$, $F^{-1} u = F M_{-1}$
5. $F \phi _{a,b} = S_{a,b} \delta$, $F \theta _a = i a M_{-1} F \theta _{- 1 / a}$
ここで $\delta$ は
<tex>
\int\mspace{-11mu}\int v(x, y) \delta(x, y) dx dy = v(0, 0) , \hspace{2zw}
v(x, y) = 0 \hspace{1zw} ( (x, y) \neq (0, 0) )
</tex>
という性質をもつ２次元のデルタ関数です．これらの記号を用いると，簡易化したフレネル回折の
作用素 $D_{\lambda d}$ を
<tex>
D_{\lambda d} = F^{-1} \{ \theta _{- \lambda d} \} F
</tex>
で定義でき，さきに積分形式で示したように
<tex>
D_{\lambda d} = \frac{1}{i \lambda d} \{ \theta _{1 / \lambda d} \}
M_{\lambda d} F \{ \theta_{1 / \lambda d} \}
</tex>
が成立します． $D_{\lambda b} D_{\lambda a} = D_{\lambda (a + b)}$ であることは
<tex>
F^{-1} \{ \theta _{- \lambda b} \} F F^{-1} \{ \theta _{- \lambda a} \} F =
F^{-1} \{ \theta _{- \lambda b} \theta _{- \lambda a} \} F =
F^{-1} \{ \theta _{- \lambda (a + b)} \} F
</tex>
で容易に確認できます．
レンズの機能
============================================================
$z$ 軸方向の位相変化を無視すると，開口が十分に広く，焦点距離が $f$ である凸レンズの
機能は作用素 $\{ \theta _{- 1 / \lambda f} \}$ で表され，
<tex>
D_{ \lambda a} \{ \theta _{- 1 / \lambda f} \} D_{ \lambda b} = \{ \frac{1}{i b} \} M_b F
\{ \frac{1}{i a} \theta _{1 / \lambda b} \theta _{- 1 / \lambda f} \theta _{1 / \lambda a} \}
M_a F \{ \theta _a \}
</tex>
より， $\theta _{1 / \lambda b} \theta _{- 1 / \lambda f} \theta _{1 / \lambda a} = \theta _0$
すなわち
<tex>
\frac{1}{a} + \frac{1}{b} = \frac{1}{f}
</tex>
であれば
<tex>
D_{\lambda a} \{ \theta _{- 1 / \lambda f} \} D_{\lambda b} =
\left \{ - \frac{a}{b} \theta _{\lambda b} \right \} M_{- b / a} \{ \theta _{\lambda a} \}
</tex>
<tex>
|( D_{\lambda a} \{ \theta _{- 1 / \lambda f} \} D_{\lambda b} u)(x, y)| =
\left| \frac{a}{b} u \left(- \frac{a}{b} x, - \frac{a}{b} y \right) \right|
</tex>
が成立します．これがレンズの公式です．
レンズの開口関数 $p$ を無視できないとき，レンズの作用素を $\{ p \theta _{ \lambda f} \}$ として解析できます．
<tex>
u' = D_{\lambda b} \{p \theta_{\lambda f} \} D_{\lambda a} u
</tex>
とおくと $1 /a + 1 / b = 1 / f$ のとき
その他の光学素子
============================================================
あとがき
============================================================
フーリエ光学に限らず，情報を変換する物理系や機器の機能を解析するとき，作用素を用いて
モデル化すると分かりやすくなることが少なくありません．このことを例示するのが本資料の
趣旨なので公式の証明は割愛しました．
@@reference: ja.wikipedia.org/wiki/%E3%83%95%E3%83%AC%E3%83%8D%E3%83%AB%E5%9B%9E%E6%8A%98,フレネル回折@@
@@reference: ja.wikipedia.org/wiki/%E3%83%AC%E3%83%B3%E3%82%BA%E3%81%AE%E5%85%AC%E5%BC%8F,レンズの公式@@
@@author: pulsar@@
@@accept: @@
@@category: 波と振動@@
@@id: @@
