ゼルニケ多項式

ゼルニケ多項式のうち、初めの21個を示す。

ゼルニケ多項式 (ゼルニケたこうしき、英語: Zernike polynomials）とは、単位円上で定義された直交多項式である。

とくに光学において軸対称な光学収差を回折理論に基づいて解析的に取り扱う際に用いられる。。

呼称は、位相差顕微鏡の発明によって1953年にノーベル物理学賞を受賞した光物理学者フリッツ・ゼルニケに由来する。

定義

ゼルニケ多項式 $Z_{n}^{m}(\rho ,\varphi )$ は、

$Z_{n}^{m}(\rho ,\varphi )={\begin{cases}R_{n}^{m}(\rho )\,\cos(m\,\varphi )&m\geq 0\\R_{n}^{|m|}(\rho )\,\sin(|m|\,\varphi )&m<0\end{cases}}$

により定義される。ここで、 n は非負整数、 m は n ≧ |m| なる整数であり、 ρ は 動径 (0 ≦ ρ ≦ 1)、 φ は偏角である。ゼルニケ多項式は $|Z_{n}^{m}(\rho ,\varphi )|\leq 1$ の範囲の値を取る。ここで、動径多項式 $R_{n}^{m}(\rho )$ は、 n − m が偶数の場合、

$R_{n}^{m}(\rho )=\sum _{k=0}^{\tfrac {n-m}{2}}{\frac {(-1)^{k}\,(n-k)!}{k!\left({\tfrac {n+m}{2}}-k\right)!\left({\tfrac {n-m}{2}}-k\right)!}}\;\rho ^{n-2\,k}$

また奇数の場合0として定義される。

他の定義

動径多項式は、二項係数を用いて、

$R_{n}^{m}(\rho )=\sum _{k=0}^{\tfrac {n-m}{2}}(-1)^{k}{\binom {n-k}{k}}{\binom {n-2k}{{\tfrac {n-m}{2}}-k}}\rho ^{n-2k}$.

と書き表すことができ、これより多項式の係数はすべて整数であることが示される。

ガウスの超幾何関数を用いて表現することもできる。この表現は、漸化式や微分方程式の導出の他、本多項式がヤコビの多項式の一部であることを示すのに有用である。

${\begin{aligned}R_{n}^{m}(\rho )&={\binom {n}{\tfrac {n+m}{2}}}\rho ^{n}\ {}_{2}F_{1}\left(-{\tfrac {n+m}{2}},-{\tfrac {n-m}{2}};-n;\rho ^{-2}\right)\\&=(-1)^{\tfrac {n-m}{2}}{\binom {\tfrac {n+m}{2}}{m}}\rho ^{m}\ {}_{2}F_{1}\left(1+{\tfrac {n+m}{2}},-{\tfrac {n-m}{2}};1+m;\rho ^{2}\right)\end{aligned}}$

動径多項式 $R_{n}^{m}(\rho )$ に含まれる項 $\rho ^{n-2k}$ は、バーンスタイン基底関数を用いて展開できる。 n が偶数の場合は $b_{s,n/2}(\rho ^{2})$ 、奇数の場合は $b_{s,(n-1)/2}(\rho ^{2})$ と $\rho$ の積で展開される。ここで、 s は $\lfloor n/2\rfloor -k\leq s\leq \lfloor n/2\rfloor$ の範囲をとる。これより、動径多項式は有限次のバーンスタイン関数として表される。

$R_{n}^{m}(\rho )={\frac {1}{\binom {\lfloor n/2\rfloor }{\lfloor m/2\rfloor }}}\rho ^{n\mod 2}\sum _{s=\lfloor m/2\rfloor }^{\lfloor n/2\rfloor }(-1)^{\lfloor n/2\rfloor -s}{\binom {s}{\lfloor m/2\rfloor }}{\binom {(n+m)/2}{s+\lceil m/2\rceil }}b_{s,\lfloor n/2\rfloor }(\rho ^{2}).$

Nollによる記法

2つの指数 n,m を並べて、1つの指数　j に統合する方法として、Nollにより提案されたのは、

$j={\frac {n(n+1)}{2}}+|m|+\left\{{\begin{array}{ll}0,&m>0\land n\in \{0,1\}{\pmod {4}};\\0,&m<0\land n\in \{2,3\}{\pmod {4}};\\1,&m\geq 0\land n\in \{2,3\}{\pmod {4}};\\1,&m\leq 0\land n\in \{0,1\}{\pmod {4}}.\end{array}}\right.$

とするものである。初めの20項を下表に示す。

n,m	0,0	1,1	1,−1	2,0	2,−2	2,2	3,−1	3,1	3,−3	3,3	4,0	4,2	4,−2	4,4	4,−4	5,1	5,−1	5,3	5,−3	5,5
j	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20

OSA/ANSIによる記法

OSA / ANSI標準ゼルニケ多項式は、以下のように定義される。

$j={\frac {n(n+2)+m}{2}}$

初めの20項を下表に示す。

n,m	0,0	1,-1	1,1	2,-2	2,0	2,2	3,-3	3,-1	3,1	3,3	4,-4	4,-2	4,0	4,2	4,4	5,-5	5,-3	5,-1	5,1	5,3
j	0	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19

Fringeによる記法

フリンジゼルニケ多項式は、以下のように定義される。

$j=\left(1+{\frac {n+|m|}{2}}\right)^{2}-2|m|+{\frac {1-\operatorname {sgn} m}{2}}$

この定義は光学設計ソフトウェアや光学検査で多く用いられる。初めの20項を下表に示す。

n,m	0,0	1,1	1,−1	2,0	2,2	2,-2	3,1	3,-1	4,0	3,3	3,-3	4,2	4,−2	5,1	5,−1	6,0	4,4	4,-4	5,3	5,-3
j	1	2	3	4	5	6	7	8	9	10	11	12	13	14	15	16	17	18	19	20

性質

直交性

動径部分は以下の直交関係を満たす。

$\int _{0}^{1}\rho {\sqrt {2n+2}}R_{n}^{m}(\rho )\,{\sqrt {2n'+2}}R_{n'}^{m}(\rho )\,d\rho =\delta _{n,n'}$

偏角部分については、初等的な計算により、

$\int _{0}^{2\pi }\cos(m\varphi )\cos(m'\varphi )\,d\varphi =\epsilon _{m}\pi \delta _{|m|,|m'|},$

$\int _{0}^{2\pi }\sin(m\varphi )\sin(m'\varphi )\,d\varphi =(-1)^{m+m'}\pi \delta _{|m|,|m'|};\quad m\neq 0,$

$\int _{0}^{2\pi }\cos(m\varphi )\sin(m'\varphi )\,d\varphi =0$

ここで $\epsilon _{m}$ は $m=0$ のとき2、$m\neq 0$ のとき1と定義される。これらより、単位円状でのゼルニケ多項式の直交性

$\int Z_{n}^{m}(\rho ,\varphi )Z_{n'}^{m'}(\rho ,\varphi )\,d^{2}r={\frac {\epsilon _{m}\pi }{2n+2}}\delta _{n,n'}\delta _{m,m'}$

が導かれる。ここで $d^{2}r=\rho \,d\rho \,d\varphi$ であり $n-m$ と $n'-m'$ はいずれも偶数と仮定している。

対称性

x軸に関する軸対称性より、

$Z_{n}^{m}(\rho ,\varphi )=(-1)^{m}Z_{n}^{m}(\rho ,-\varphi )$

原点に関する点対称性より、

$Z_{n}^{m}(\rho ,\varphi )=(-1)^{m}Z_{n}^{m}(\rho ,\varphi +\pi )$

ここで、 $n-m$ は偶数であると仮定しているので、 $(-1)^{m}$ は $(-1)^{n}$ と書き換えることができる。動径多項式は、n,mに応じて偶関数または奇関数である。

$R_{n}^{m}(\rho )=(-1)^{n}R_{n}^{m}(-\rho )=(-1)^{m}R_{n}^{m}(-\rho )$

三角関数の周期性より、原点を中心とした m 回回転対称性が生じる。

$Z_{n}^{m}\left(\rho ,\varphi +{\tfrac {2\pi k}{m}}\right)=Z_{n}^{m}(\rho ,\varphi ),\qquad k=0,\pm 1,\pm 2,\dots$

漸化式

動径多項式は、以下の漸化式を満たす。

${\begin{aligned}R_{n}^{m}(\rho )+R_{n-2}^{m}(\rho )=\rho \left[R_{n-1}^{\left|m-1\right|}(\rho )+R_{n-1}^{m+1}(\rho )\right]{\text{ .}}\end{aligned}}$

動径多項式の定義より、$R_{m}^{m}(\rho )=\rho ^{m},\quad R_{m+2}^{m}(\rho )=((m+2)\rho ^{2}-(m+1))\rho ^{m}$ である。これと、以下の三項間漸化式 により、すべての $R_{n}^{m}(\rho )$ を計算することができる。

$R_{n}^{m}(\rho )={\frac {2(n-1)(2n(n-2)\rho ^{2}-m^{2}-n(n-2))R_{n-2}^{m}(\rho )-n(n+m-2)(n-m-2)R_{n-4}^{m}(\rho )}{(n+m)(n-m)(n-2)}}{\text{ .}}$

この式より、動径多項式の導関数を、2つの動径多項式から計算することができる。

${\frac {\operatorname {d} }{\operatorname {d} \!\rho }}R_{n}^{m}(\rho )={\frac {(2nm(\rho ^{2}-1)+(n-m)(m+n(2\rho ^{2}-1)))R_{n}^{m}(\rho )-(n+m)(n-m)R_{n-2}^{m}(\rho )}{2n\rho (\rho ^{2}-1)}}{\text{ .}}$

例

動径多項式

動径多項式は以下のような式となる。

$R_{0}^{0}(\rho )=1\,$

$R_{1}^{1}(\rho )=\rho \,$

$R_{2}^{0}(\rho )=2\rho ^{2}-1\,$

$R_{2}^{2}(\rho )=\rho ^{2}\,$

$R_{3}^{1}(\rho )=3\rho ^{3}-2\rho \,$

$R_{3}^{3}(\rho )=\rho ^{3}\,$

$R_{4}^{0}(\rho )=6\rho ^{4}-6\rho ^{2}+1\,$

$R_{4}^{2}(\rho )=4\rho ^{4}-3\rho ^{2}\,$

$R_{4}^{4}(\rho )=\rho ^{4}\,$

$R_{5}^{1}(\rho )=10\rho ^{5}-12\rho ^{3}+3\rho \,$

$R_{5}^{3}(\rho )=5\rho ^{5}-4\rho ^{3}\,$

$R_{5}^{5}(\rho )=\rho ^{5}\,$

$R_{6}^{0}(\rho )=20\rho ^{6}-30\rho ^{4}+12\rho ^{2}-1\,$

$R_{6}^{2}(\rho )=15\rho ^{6}-20\rho ^{4}+6\rho ^{2}\,$

$R_{6}^{4}(\rho )=6\rho ^{6}-5\rho ^{4}\,$

$R_{6}^{6}(\rho )=\rho ^{6}\,$

ゼルニケ多項式

ゼルニケ多項式は以下のような式となる。なお、各式は

$\int _{0}^{2\pi }\int _{0}^{1}Z_{j}^{2}\,\rho \,d\rho \,d\theta =\pi$

を満たすよう規格化されている。

	OSA/ANSI 指数 ($j$)	Noll 指数 ($j$)	Fringe 指数 ($j$)	動径 次数 ($n$)	偏角 次数 ($m$)	$Z_{j}$	名称
$Z_{0}^{0}$	00	01	1	0	00	$1$	Piston
$Z_{1}^{-1}$	01	03	3	1	−1	$2\rho \sin \theta$	Tilt (Y-Tilt, Vertical tilt)
$Z_{1}^{1}$	02	02	2	1	+1	$2\rho \cos \theta$	Tip (X-Tilt, Horizontal tilt)
$Z_{2}^{-2}$	03	05	6	2	−2	${\sqrt {6}}\rho ^{2}\sin 2\theta$	Oblique astigmatism
$Z_{2}^{0}$	04	04	4	2	00	${\sqrt {3}}(2\rho ^{2}-1)$	Defocus
$Z_{2}^{2}$	05	06	5	2	+2	${\sqrt {6}}\rho ^{2}\cos 2\theta$	Vertical astigmatism
$Z_{3}^{-3}$	06	09	11	3	−3	${\sqrt {8}}\rho ^{3}\sin 3\theta$	Vertical trefoil
$Z_{3}^{-1}$	07	07	8	3	−1	${\sqrt {8}}(3\rho ^{3}-2\rho )\sin \theta$	Vertical coma
$Z_{3}^{1}$	08	08	7	3	+1	${\sqrt {8}}(3\rho ^{3}-2\rho )\cos \theta$	Horizontal coma
$Z_{3}^{3}$	09	10	10	3	+3	${\sqrt {8}}\rho ^{3}\cos 3\theta$	Oblique trefoil
$Z_{4}^{-4}$	10	15	18	4	−4	${\sqrt {10}}\rho ^{4}\sin 4\theta$	Oblique quadrafoil
$Z_{4}^{-2}$	11	13	13	4	−2	${\sqrt {10}}(4\rho ^{4}-3\rho ^{2})\sin 2\theta$	Oblique secondary astigmatism
$Z_{4}^{0}$	12	11	9	4	00	${\sqrt {5}}(6\rho ^{4}-6\rho ^{2}+1)$	Primary spherical
$Z_{4}^{2}$	13	12	12	4	+2	${\sqrt {10}}(4\rho ^{4}-3\rho ^{2})\cos 2\theta$	Vertical secondary astigmatism
$Z_{4}^{4}$	14	14	17	4	+4	${\sqrt {10}}\rho ^{4}\cos 4\theta$	Vertical quadrafoil