AlxisDataBlog-数学・物理学

前々回の式を復習しましょう。

\begin{align} &P_0:\langle 1\rangle_k=\langle 1\rangle_k,\\ &P_1:r^2\Delta\langle 1\rangle_k=\langle X^{-1}\partial+\partial^2\rangle_k,\\ &P_2:(2+r\partial_r)\langle 1\rangle_k=-\langle X\partial\rangle_k,\\ &P_3:(r^2\Delta)^2\langle 1\rangle_k=\mu_2r\partial_r\psi-(\mu_2+\mu_3)\psi+\langle (X^{-1}\partial+\partial^2)^2\rangle_k,\\ &P_4:(r^2\Delta)(2+r\partial_r)\langle 1\rangle_k=\mu_2\psi-\langle (2+X\partial)(X^{-1}\partial+\partial^2)\rangle_k,\\ &P_5:(2+r\partial_r)^2\langle 1\rangle_k=\langle (X\partial)^2\rangle_k,\\ &P_6:(r^2\Delta)(2+r\partial_r)^2\langle 1\rangle_k=\mu_2r\partial_r\psi-(3\mu_2+\mu_3)\psi+\langle (X^{-1}\partial+\partial^2)(X\partial)^2\rangle_k,\\ &P_7:(2+r\partial_r)^3\langle 1\rangle_k=\mu_2\psi-\langle (X\partial)^3\rangle_k \end{align}

今回考えるフレーム内処理では動径微分が禁止されています。これを言葉通りに解釈すると使える方程式は $P_0,P_1,P_3$ のみであり、$\psi,\partial_r\psi$ の両方を計算することはできません。フレーム内処理によって近似的に動径微分を計算することを考えてみましょう。

$P_0,P_1,P_2$ を連立させて

\begin{align} \{c_0+c_1(r^2\Delta)+(2+r\partial_r)\}\langle 1\rangle_k=\langle c_0+c_1(X^{-1}\partial+\partial^2)-X\partial\rangle_k \end{align}

右辺をコーシーの不等式で評価して最小二乗法により $c_0,c_1$ を求めれば

\begin{align} (2+r\partial_r)\langle 1\rangle_k\approx -\{c_0+c_1(r^2\Delta)\}\langle 1\rangle_k \end{align}

のような近似が成立します。数学的には、更なる導関数(例えば $(r^2\Delta)(2+r\partial_r)\langle 1\rangle_k$)に関しては全く主張は得られないのですが、仮に

\begin{align} (2+r\partial_r)\langle 1\rangle_k= -\{c_0+c_1(r^2\Delta)\}\langle 1\rangle_k \end{align}

として議論を進めることにします。すると $P_4,P_5,P_6,P_7$ も使えるようになります。

\begin{align} &(r^2\Delta)(2+r\partial_r)\langle 1\rangle_k=-\{c_0(r^2\Delta)+c_1(r^2\Delta)^2\}\langle 1\rangle_k,\\ &(2+r\partial_r)^2\langle 1\rangle_k=\{c_0^2+2c_1(c_0-1)(r^2\Delta)+c_1^2(r^2\Delta)^2\}\langle 1\rangle_k,\\ &(r^2\Delta)(2+r\partial_r)^2\langle 1\rangle_k=\{c_0^2(r^2\Delta)+2c_1(c_0-1)(r^2\Delta)^2+c_1^2(r^2\Delta)^3\}\langle 1\rangle_k,\\ &(2+r\partial_r)^3\langle 1\rangle_k=-\{c_0^3+c_1(3c_0^2-6c_0+4)(r^2\Delta)+3c_1^2(c_0-2)(r^2\Delta)^2+c_1^3(r^2\Delta)^3\}\langle 1\rangle_k \end{align}

剰余項の生成式を増やすために $P_3,P_4$ を整理しておきましょう。

\begin{align} P_8:\mu_2P_3+(\mu_2+\mu_3)P_4 \end{align}

$P_4,P_8$ を用いて $P_6,P_7$ から環状項を消去します。

\begin{align} &P_9:-P_3+2P_4+P_6,\\ &P_{10}:-P_4+P_7 \end{align}

更に微分してと続けても良いのですが、ラプラシアンの次数に配慮して、この辺りで剰余項の生成を終わりにしましょう。ここからは前々回と同様ですので割愛します。

今回はここまで。次回は動径微分をもう少し良く近似できるようにします。(written by J.Nishigata)

実用的に使えるフレーム内デコンボリューション前半