Bessel函数
2022/9/6 23:25:47
本文主要是介绍Bessel函数,对大家解决编程问题具有一定的参考价值,需要的程序猿们随着小编来一起学习吧!
在这篇文章中,我们将会罗列Bessel函数的一些基本性质。
A. Definition and Basic Properties
We define the Bessel function $J_{\nu}$ of order $\nu$ by its Poisson representation formula $$J_{\nu}(t) = \frac{(t/2)^{\nu}}{\Gamma(\nu + 1/2)\Gamma(1/2)}\int_{-1}^1e^{its}(1 - s^2)^{\nu}\frac{{\rm d}s}{\sqrt{1 - s^2}},$$ where ${\rm Re}\ \nu > -1/2$ and $t \geq 0$. Note that in this case $J_{\nu}(t)$ is a real number.
Here are some basic properties of Bessel functions.
- (Recurrence formula)$$\frac{{\rm d}}{{\rm d}t}(t^{-\nu}J_{\nu}(t)) = -t^{-\nu}J_{\nu + 1}(t), \qquad {\rm Re} \ \nu > -\frac{1}{2}; $$
- (Companion recurrence formula) $$\frac{{\rm d}}{{\rm d}t}(t^{\nu}J_{\nu}(t)) = t^{\nu}J_{\nu - 1}(t), \qquad {\rm Re} \ \nu > \frac{1}{2}; $$
- $J_{\nu}(t)$ satisfies the Bessel equation $$t^2f''(t) + tf'(t) + (t^2 - \nu^2)f(t) = 0;$$
- If $\nu \in \mathbb{Z}_+$, then we have the following identity, which was taken by Bessel as the definition of $J_{\nu}$ for integer $\nu$: $$J_{\nu}(t) = \frac{1}{2\pi}\int_0^{2\pi}e^{it\sin\theta}e^{-i\nu\theta} \,{\rm d}\theta = \frac{1}{2\pi}\int_0^{2\pi}\cos(t\sin\theta - \nu\theta) \,{\rm d}\theta;$$
- For ${\rm Re}\ \nu > -1/2$ we have the following identity: $$J_{\nu}(t) = \frac{(t/2)^{\nu}}{\Gamma(1/2)}\sum\limits_{j = 0}^{\infty}(-1)^j\frac{\Gamma(j + 1/2)}{\Gamma(j + \nu + 1)}\frac{t^{2j}}{(2j)!} = \sum\limits_{j = 0}^{\infty}\frac{(-1)^j}{j!}\frac{(t/2)^{2j + \nu}}{\Gamma(j + \nu + 1)}; $$
- For ${\rm Re}\ \nu > 1/2$ we have the following identity: $$\frac{{\rm d}}{{\rm d}t}(J_{\nu}(t)) = \frac{1}{2}(J_{\nu - 1}(t) - J_{\nu + 1}(t)).$$
Remark. The property 1, 2, 3 and 4 still hold when $t \in \mathbb{C}$.
Proposition. Let ${\rm Re} \ \mu > -1/2, {\rm Re} \ \nu > -1$ and $t > 0$. Then the following identity is valid: \begin{equation}\label{1}\int_0^1J_{\mu}(ts)s^{\mu + 1}(1 - s^2)^{\nu} \,{\rm d}s = \frac{\Gamma(\nu + 1)2^{\nu}}{t^{\nu + 1}}J_{\mu + \nu + 1}(t).\end{equation}
B. The Fourier Transform of Surface Measure on $\mathbb{S}^{n - 1}$
Let ${\rm d}\sigma$ denote the surface on $\mathbb{S}^{n - 1}$ for $n \geq 2$. Then the following is true: $$\widehat{{\rm d}\sigma}(\xi) = \int_{\mathbb{S}^{n - 1}}e^{-2\pi i\xi \cdot \theta} \,{\rm d}\theta = \frac{2\pi}{|\xi|^{\frac{n - 2}{2}}}J_{\frac{n - 2}{2}}(2\pi|\xi|).$$ Using polar coordinates, we can also obtain the Fourier transform of a radial function on $\mathbb{R}^n$: $$\widehat{f}(\xi) = \frac{2\pi}{|\xi|^{\frac{n - 2}{2}}}\int_0^{\infty}J_{\frac{n - 2}{2}}(2\pi r|\xi|)r^{\frac{n}{2}} \,{\rm d}r,$$ where $f(x) = f_0(|x|)$.
Example. Consider the radial function $f(x) = \chi_{B_1(0)}(x)$ on $\mathbb{R}^n$. It follows that $$(\chi_{B_1(0)})^{\wedge}(\xi) = \frac{2\pi}{|\xi|^{\frac{n - 2}{2}}}\int_0^1J_{\frac{n - 2}{2}}(2\pi|\xi|r)r^{\frac{n}{2}} \,{\rm d}r = \frac{J_{\frac{n}{2}}(2\pi|\xi|)}{|\xi|^{\frac{n}{2}}},$$ where we use identity \eqref{1}. More generally, for ${\rm Re}\ \lambda > -1$, let \begin{equation*}m_{\lambda}(\xi) = \begin{cases}(1 - |\xi|^2)^{\lambda} \quad &|\xi| \leq 1, \\ 0 \quad &|\xi| > 1. \end{cases}\end{equation*} Then $$m_{\lambda}^{\vee}(x) = \frac{\Gamma(\lambda + 1)}{\pi^{\lambda}}\frac{J_{\frac{n}{2} + \lambda}(2\pi|x|)}{|x|^{\frac{n}{2} + \lambda}}.$$
C. Asymtotics of Bessel Funtions
Let ${\rm Re}\ \nu > -1/2$. We have the following results:
\begin{equation*}J_{\nu}(t) = \begin{cases}\frac{t^{\nu}}{2^{\nu}\Gamma(\nu + 1)} + O(t^{{\rm Re}\ \nu + 1}) \qquad &t \rightarrow 0^+, \\ \sqrt{\frac{2}{t\pi}}\cos\left(t - \frac{\nu\pi}{2} - \frac{\pi}{4}\right) + O(t^{-\frac{3}{2}}) \qquad &t \rightarrow \infty. \\ \end{cases}\end{equation*} In particular, for fixed $\nu$, we have $J_{\nu}(t) = O(t^{-\frac{1}{2}})$ as $t \rightarrow \infty$.
Ref. Grafakos, L. Classical Fourier Analysis.
苗长兴, 现代调和分析及其应用讲义.
这篇关于Bessel函数的文章就介绍到这儿,希望我们推荐的文章对大家有所帮助,也希望大家多多支持为之网!
- 2024-05-01为什么公共事业机构会偏爱 TiDB :TiDB 数据库在某省妇幼健康管理系统的应用
- 2024-04-26敏捷开发:想要快速交付就必须舍弃产品质量?
- 2024-04-26静态代码分析的这些好处,我竟然都不知道?
- 2024-04-26你在测试金字塔的哪一层?(下)
- 2024-04-26快刀斩乱麻,DevOps让代码评审也自动起来
- 2024-04-262024年最好用的10款ER图神器!
- 2024-04-2203-为啥大模型LLM还没能完全替代你?
- 2024-04-2101-大语言模型发展
- 2024-04-17基于SpringWeb MultipartFile文件上传、下载功能
- 2024-04-14个人开发者,Spring Boot 项目如何部署