5

Fourier-Bessel 展開 ・ Gauss Quadrature

219
1
$$\newcommand{BA}[0]{\begin{align*}} \newcommand{BE}[0]{\begin{equation}} \newcommand{bl}[0]{\boldsymbol} \newcommand{BM}[0]{\begin{matrix}} \newcommand{D}[0]{\displaystyle} \newcommand{EA}[0]{\end{align*}} \newcommand{EE}[0]{\end{equation}} \newcommand{EM}[0]{\end{matrix}} \newcommand{h}[0]{\boldsymbol{h}} \newcommand{k}[0]{\boldsymbol{k}} \newcommand{L}[0]{\left} \newcommand{l}[0]{\boldsymbol{l}} \newcommand{m}[0]{\boldsymbol{m}} \newcommand{n}[0]{\boldsymbol{n}} \newcommand{R}[0]{\right} \newcommand{vep}[0]{\varepsilon} $$

内容が薄すぎたので全く関係ない二本立てになりました.

Fourier-Bessel 展開

$\hspace{5pt}$Legendre 多項式などの直行多項式を基底関数とした級数展開は私の最近の興味深いテーマでしたが,興味を無くしたところでふと思いました.Fourier 級数展開(Chebyshev 級数展開)があるなら,Bessel 関数での展開方法もあるのではないかと.調べたら,あるみたいです.『Fourier-Bessel 展開』.以後 FB 展開と略します.
$\hspace{5pt}$ Bessel 関数の直交性と Fourier-Bessel 展開 /東京大学宇宙理論研究室
$\hspace{5pt}$ Fourier-Bessel 展開(Dini 展開) /とつとつとしてろうとせず
$\hspace{5pt}$ Bessel function /Wikipedia
$\hspace{5pt}$また,直行多項式について基本的なことが次の PDF でそこそこ記述されていました.Bessel 関数の直行性も示されています.
$\hspace{5pt}$ 物理数学 II /辻直人(東京大学)

Fourier-Bessel 展開

$\hspace{5pt}$$\nu$次の第一種 Bessel 関数$J_\nu^{}(x)$$n$番目の零点を$\lambda_n^{(\nu)}$とすれば,$\nu>-1$において直行性

$\BA\D \int_0^1 xJ_\nu^{}\L(\lambda_m^{(\nu)}x\R)J_\nu^{}\L(\lambda_n^{(\nu)}x\R)\,dx=\frac{\delta_{m,n}}{2}J_{\nu+1}^{}\L(\lambda_n^{(\nu)}\R)^2 \EA$

が成り立ち,$[0,1]$で定義される任意の関数$f(x)$

$\BA\D f(x)&=\sum_{n=1}^\infty c_n^{}J_{\nu}^{}\L(\lambda_n^{(\nu)}x\R)\\ c_n^{}&=\frac{2}{J_{\nu+1}^{}\L(\lambda_n^{(\nu)}\R)^2}\int_0^1 xf(x)J_{\nu}^{}\L(\lambda_n^{(\nu)}x\R)\,dx \EA$

と展開できます.

$c_n^{}$のことは展開係数と呼ぶことにします.FB 展開係数の計算にあたって Bessel 関数の微分方程式および漸化式が必要になると思いますので,先に列挙してみます.

$\hspace{15pt}$微分漸化式1:$\D 2\frac{d}{dx}J_{\nu}^{}(x)=J_{\nu-1}^{}(x)-J_{\nu+1}^{}(x)$

$\hspace{15pt}$微分漸化式2:$\D 2\nu J_{\nu}^{}(x)=x\L(J_{\nu-1}^{}(x)+J_{\nu+1}^{}(x)\R)$

$\hspace{15pt}$微分漸化式3:$\D \frac{d}{dx}x^{\nu}J_{\nu}^{}(x)=x^{\nu}J_{\nu-1}^{}(x)$

$\hspace{15pt}$微分漸化式4:$\D \frac{d}{dx}x^{-\nu}J_{\nu}^{}(x)=-x^{-\nu}J_{\nu+1}^{}(x)$

$\hspace{15pt}$微分方程式5:$\D \frac{d^2}{dx^2}J_{\nu}^{}(x)+\frac{1}{x}\frac{d}{dx}J_{\nu}^{}(x)-\L(1-\frac{\nu^2}{x^2}\R)J_{\nu}^{}(x)=0$

FB 展開係数の計算

$\hspace{5pt}$先に述べますが,例えば$\sin\pi x$や楕円積分$K(x)$などの場合でも,線形微分方程式を満たすなら漸化式を立てることは可能ですが,多項間漸化式となり,一般項を求めるのは難しいようです.ただ,一部の場合は容易に計算することができそうです.

$f(x)=x^{\nu}$

$\hspace{5pt}\nu>-1$とします.

$\BA\D \int_0^1 xf(x)J_{\nu}^{}\L(\lambda_n^{(\nu)}x\R)\,dx &=\int_0^1 x^{\nu+1}J_{\nu}^{}\L(\lambda_n^{(\nu)}x\R)\,dx\\ &=\frac{1}{\L(\lambda_n^{(\nu)}\R)^{\nu+2}}\int_0^{\lambda_n^{(\nu)}} x^{\nu+1}J_{\nu}^{}\L(x\R)\,dx\\ &=\frac{1}{\L(\lambda_n^{(\nu)}\R)^{\nu+2}}\int_0^{\lambda_n^{(\nu)}} \frac{d}{dx}x^{\nu+1}J_{\nu+1}^{}\L(x\R)\,dx\\ &=\frac{1}{\L(\lambda_n^{(\nu)}\R)^{\nu+2}}\L(\lambda_n^{(\nu)}\R)^{\nu+1}J_{\nu+1}^{}\L(\lambda_n^{(\nu)}\R)\\ &=\frac{J_{\nu+1}^{}\L(\lambda_n^{(\nu)}\R)}{\lambda_n^{(\nu)}} \EA$

よって,

$\BA\D x^{\nu}=\sum_{n=1}^\infty \frac{2J_{\nu}^{}\L(\lambda_n^{(\nu)}x\R)}{\lambda_n^{(\nu)}J_{\nu+1}^{}\L(\lambda_n^{(\nu)}\R)^2} \EA$

となります.ところで,これに似たかたちの Kneser-Sommerfeld expansion というものがあるらしく,次のような式になっています.

$\BA\D \sum_{n=1}^\infty \frac{\L(\lambda_n^{(\nu)}\R)^{2\nu-\alpha}J_\alpha^{}\L(\lambda_n^{(\nu)}x\R)}{\L(\L(\lambda_n^{(\nu)}\R)^2-z^2\R)J_{\nu+1}^{}\L(\lambda_n^{(\nu)}\R)^2}=\frac{\pi}{4J_{\nu}^{}(z)\sin\pi\nu}\L(\frac{\L(\frac{x}{2}\R)^{\alpha-2\nu}J_{\nu}^{}(z)}{\Gamma(1-\nu)\Gamma(1+\alpha-\nu)}{_1}F_2\L[\begin{matrix}1\\1-\nu,1+\alpha-\nu\end{matrix};-\L(\frac{xz}{2}\R)^2\R]-z^{2\nu-\alpha}J_{-\nu}^{}(z)J_{\alpha}^{}(xz)\R) \EA$

詳細は次の文献に記載されています.

$\hspace{5pt}$ Summing Sneddon-Bessel series explicitly /Antonio J. Dur´an, Mario P´erez and Juan L. Varona

FB展開ですが,実際に計算するとなるとのような過程になるのかさっぱりわかりません.

Bessel zeta function

$\hspace{5pt}$Bessel zeta function は次のように定義されます.

$\BA\D \zeta_{\nu}^{}(s)=\sum_{n=1}^\infty \L(\lambda_n^{(\nu)}\R)^{-s} \EA$

Bessel 関数に限らず,特殊関数の零点の和をとったものは様々提案されているようで,次の文献には Bessel 関数のほかに Airy 関数やその q 類似についても記述されています.

$\hspace{5pt}$ Sums of Zeros for Certain Special Functions /Ruiming Zhang

$\zeta_{\nu}^{}(s)$について,$s$が偶数の場合は計算可能らしく,

$\BA\D \zeta_{\nu}^{}(0)=-\frac{1}{2}\L(\nu+\frac{1}{2}\R),\qquad \zeta_{\nu}^{}(2)=-\frac{1}{2^2(\nu+1)},\qquad \zeta_{\nu}^{}(2n)=\frac{1}{n+\nu}\sum_{k=1}^{n-1} \zeta_{\nu}^{}(2k)\zeta_{\nu}^{}(2n-2k)\quad (n=2,3,\cdots) \EA$

となっています.また,$s=-2n$のとき

$\BA\D c_0^{}=c_1^{}=\frac{1}{2}\L(\nu^2-\frac{1}{4}\R) ,\qquad 2c_{n+2}^{}=(n+3)c_{n+1}^{}-\sum_{k=0}^n c_{k}^{}c_{n-k}^{} \EA$

を満たす$c_n^{}$に対して

$\BA\D \zeta_{\nu}^{}(-2n)=\frac{(-1)^nc_{2n-1}}{2} \quad (n=1,2,\cdots) \EA$

となっています.ちなみに Bessel zeta の零点に関する文献は見つかりませんでした.

↓ここから全く別の内容

Gauss quadrature

$\hspace{5pt}$Gauss quadrature(;ガウス求積)とは

[ガウス求積 - Wikipedia](https://www.bing.com/ck/a?!&&p=7e345525f9b84443JmltdHM9MTcwMTIxNjAwMCZpZ3VpZD0zNmRlM2Y0OS01N2VkLTY4MTUtMDZlMy0yYzlmNTYwNTY5ODgmaW5zaWQ9NTE3Nw&ptn=3&ver=2&hsh=3&fclid=36de3f49-57ed-6815-06e3-2c9f56056988&psq=Gaussian+Quadrature+%e3%81%a8%e3%81%af&u=a1aHR0cHM6Ly9qYS53aWtpcGVkaWEub3JnL3dpa2kvJUUzJTgyJUFDJUUzJTgyJUE2JUUzJTgyJUI5JUU2JUIxJTgyJUU3JUE5JThE&ntb=1) ガウス求積 - Wikipedia

らしいです.英語表記では,Gaussian quadratureも見るのですが,termとしてはどちらを用いればよいのでしょうか.

$\omega(x)$を重み関数にもつ$n$次の直交多項式$p_n^{}(x)$に対して
$p_n^{}(x)$$n$個の零点をノード(分点)として選ぶと,次数が$2n-1$以下の任意の多項式について正確な積分値を与える$n$個の重み$w_1^{},w_2^{},\cdots,w_n^{}$を選ぶことができる.さらに,それらのノードには重複がなくすべて開区間$(a,b)$にある.

らしいです.

$\bigcirc$$w_i^{}$の計算方法
$\hspace{5pt}$重み$w_i^{}$は次のように表されます.

$\BA\D w_i^{}=\frac{A_n^{}}{A_{n-1}^{}}\frac{1}{p_n'(x_i^{})p_{n-1}^{}(x_i^{})}\int_a^b p_{n-1}^{}(x)^2\omega(x)\,dx \EA$

ただし,$A_n^{}$$p_n^{}(x)$における$x^n$の係数です.
長々と何か書いてあるみたいですが,結果だけ抜き出すと,

$\hspace{0pt}$Gauss quadrature

$\hspace{5pt}$$2n-1$次以下の多項式$f(x)$に対して,重み関数$\omega(x)$の直交多項式$p_n^{}(x)=A_n^{}(x-x_1^{})(x-x_2^{})\cdots(x-x_n^{})$とし,

$\BA\D w_i^{}=\frac{A_n^{}}{A_{n-1}^{}}\frac{1}{p_n'(x_i^{})p_{n-1}^{}(x_i^{})}\int_a^b p_{n-1}^{}(x)^2\omega(x)\,dx \EA$

$(a,b)$は各$p_n^{}(x)$ごとに決まる区間

とすれば,

$\BA\D \int_a^b \omega(x)f(x)\,dx=\sum_{k=1}^n w_k^{}f(x_k^{}). \EA$

ということだと思います.間違っていたら申し訳ないです.
以下にいくつかの直交多項式と$w_i$をまとめてみました.

$p_n^{}(x)$$\omega(x)$$(a,b)$$w_i^{}$
$\textrm{Legendre~Polynomial}:P_n^{}(x)$$1$$(-1,1)$$\D\frac{2}{\D (1-x_i^2)P_n'(x_i^{})^2}$
$\textrm{Chebyshev~Polynomial}:T_n^{}(x)$$\D\frac{1}{\sqrt{1-x^2}}$$(-1,1)$$\D\frac{\pi}{n}$
$\textrm{Chebyshev~Polynomial}:U_n^{}(x)$$\sqrt{1-x^2}$$(-1,1)$$\D\frac{\pi}{n+1}\sin^2\frac{\pi i}{n+1}$
$\textrm{Gegenbauer~Polynomial}:C_n^{(\alpha)}(x)$$(1-x^2)^{-\frac{1}{2}+\alpha}$$(-1,1)$$\D \frac{\pi\Gamma(2\alpha)}{\Gamma(\alpha)^2}\frac{(2\alpha)_n^{}}{n!}\frac{\D 2^{2-2\alpha}(1-x_i^2)}{\D n(2n-1+\alpha)C_{n-1}^{(\alpha)}(x_i^{})^2}$
$\textrm{Laguerre~Polynomial}:L_n^{}(x)$$e^{-x}$$(0,\infty)$$\D\frac{x_i^{}}{\D (n+1)^2L_{n+1}^{}(x_i^{})^2}$
$\textrm{Hermite~Polynomial}:H_n^{}(x)$$e^{-x^2}$$(-\infty,\infty)$$\D\frac{\D 2^{n-1}n!\sqrt{\pi}}{\D n^2H_{n-1}^{}(x_i^{})^2}$

また,$\D\gamma_n^{}=\int_a^b p_{n}^{}(x)^2\omega(x)\,dx$とおくとき

$\BA\D \sum_{k=0}^n \frac{p_k^{}(x)^2}{\gamma_k^{}}=\frac{A_n^{}}{A_{n+1}^{}\gamma_n^{}}\L(p_n^{}(x)p_{n+1}'(x)-p_n'(x)p_{n+1}^{}(x)\R) \EA$

が成り立つらしいです.
$f(x)$$2n$次以上のときは,積分と級数の間に誤差が生じ,$(a,b)$上のどこかの数$\xi$を用いて

$\BA\D E_n^{}=\int_a^b \omega(x)f(x)\,dx-\sum_{k=1}^n w_k^{}f(x_k^{})=\frac{\gamma_n^{}}{A_n^2}\frac{f^{(2n)}(\xi)}{(2n)!} \EA$

と書けるらしいです.

投稿日:20231129
更新日:20231130

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中