柱函数-贝塞尔函数

最新推荐文章于 2024-09-09 17:43:12 发布

力语

最新推荐文章于 2024-09-09 17:43:12 发布

阅读量2.3k

点赞数 3

分类专栏：数学物理方法文章标签：线性代数概率论傅立叶分析几何学算法

本文链接：https://blog.csdn.net/woaiwulima/article/details/122190620

版权

本文深入探讨柱函数，特别是Bessel函数，包括其定义、解法、性质和应用。从圆膜振动问题出发，解析Bessel方程的解，详细阐述第一类和第二类柱函数，讨论其特殊性质、零点、递推公式和正交归一化关系。此外，还介绍了其他柱函数如诺依曼函数和变形贝塞尔函数，以及它们的渐近行为。

摘要由CSDN通过智能技术生成

柱函数的提出

考虑固定边界的圆膜振动, 可以归结为定解问题
$\begin{cases}u_{u}=a^{2}\left(u_{x x}+u_{y y}\right) & \left(0 \leqslant x^{2}+y^{2}<l^{2}, t>0\right) \\ \left.u\right|_{x^{2}+y^{2}=l^{2}}=0 & (t \geqslant 0), \\ u(x, y, 0)=\varphi(x, y) & , u_{l}(x, y, 0)=\psi(x, y) \quad \left(0 \leqslant x^{2}+y^{2} \leqslant l^{2}\right),\end{cases}$
在柱坐标下定解问题变为
$\left\{\begin{array}{l}r^{2} R^{\prime \prime}(r)+r R^{\prime}(r)+\left[(k r)^{2}-v^{2}\right] R(r)=0 \\ \Phi^{\prime \prime}(\varphi)+n^{2} \Phi(\varphi)=0, \quad \Phi(\varphi+2 \pi)=\Phi(\varphi) . \\ Z^{\prime \prime}+\mu Z(z)=0\end{array}\right.$
$k^2=\lambda-\mu$ ， $\lambda$ 在分离时间时引入， $\mu$ 在分离 $z$ 变量时引入， $\mu=0$ 表示系统在z方向平移不变， $\lambda=0$ 表示稳定场。

在 $k^2=\lambda-\mu \neq 0$ 时的径向方程为bessel方程。

令 $k r = x$ ， $R (r) = y (x)$ ， $x^{2} y^{\prime \prime}+x y^{\prime}+\left(x^{2}-v^{2}\right) y=0$ ， $v$ 阶Bessel方程

$y^{\prime \prime}+\frac{1}{x} y^{\prime}+\left(1-\frac{v^{2}}{x^{2}}\right) y=0$

Bessel方程的解

$x_0$ 是方程的正则奇点，为了使最后不会得到奇奇怪怪x的小数项，级数展开时都乘以 $x^\rho$ ，则 $y(x)=x^\rho \sum_{k=0}^{\infty }c_k(x-x_0)^k$

$\rho$ 为指标，求 $\rho$ 的方程–指标方程

将 $y=\sum^{\infty}_{k=0}c_kx^{k+\rho}$ 带入Bessel 方程并化简得到

$\left\{\begin{array}{l}\left(\rho^{2}-\nu^{2}\right) c_{0}=0,\\ \left[(\rho+1)^{2}-\nu^{2}\right] c_{1}=0,\\ \left[(\rho+k)^{2}-\nu^{2}\right] c_{k}+c_{k-2}=0, k=2,3,4, \cdots .\end{array}\right.$

由于 $c_0=c_1=0$ 只能得到平凡解，故 $c_0\neq0,c_1\neq0$ ，可以舍去，因此可以得到指标方程为 $\rho^2-v^2=0$ ，所以 $\rho_1=v$ 或 $\rho_2=-v$ 。

当 $\rho_1=v$ 时， $c_0\neq 0$ ， $c_1=0$ ，令 $c_{2 n}=(-1)^{n} \frac{c_{0} \Gamma(\nu+1)}{2^{2 n} n ! \Gamma(\nu+n+1)}$ ， $c_0=\frac{1}{2^v\Gamma(\nu+1)}$ 。

【这样取 $c_0,c_{2n}$ 是为了方便母/子函数表达方便】

则第一个解 $y_1(x)=J_v(x)=\sum_{k=0}^{\infty}(-1)^k\frac{1}{k!\Gamma(v+k+1)}(\frac{x}{2})^{2k+v}$ ， $J_v(x)$ 称作 $v$ 阶Bessel函数。

若 $\rho_2=-v$ ，则称作 $- v$ 阶Bessel 函数， $y_2(x)=J_{-v}(x)=\sum_{k=0}^{\infty}(-1)^k\frac{1}{k!\Gamma(-v+k+1)}(\frac{x}{2})^{2k-v}$

$J_v(x),J_{-v}(x)$ 称为第一类柱函数

①若 $v\neq n$ (整数) ， $J_v(x)$ 与 $J_{-v}(x)$ 是线性无关

$y(x)=a_vJ_v(x)+b_vJ_{-v}(x)$

②若 $v = n$ (整数) ， $J_{-v}(x)=(-1)^nJ_v(x)$

$J_n$ 与 $J_{-n}$ 线性相关，须构建其他与 $J_n(x)$ 线性无关的函数。

Bessel函数的性质

贝塞尔函数的本征值：

$\left\{\begin{array}{l} r R^{\prime \prime}(r)+r R^{\prime}(r)+\left[(k r)^{2}-n^{2}\right] R(r)=0\\ R(r)_{(r\rightarrow0)}有限 \\ \left. \left[\alpha \frac{d R}{d r}+\beta R\right]\right|_{r=a}=0 \end{array}\right.$

$R_m^n(P)=R(k_m^nP)$

$R(r)=J_n(kr)$ ， $k^2=\lambda-\mu$ ，令 $x = k r$

以 $\alpha=0$ 这种第一类齐次边界条件为例

$\left.J_{n}(k r)\right|_{r=a}=J_n(ka)=0$ ，零根 $x_m^n$ 是使 $J_n(ka)=0$ 为 $n$ 阶贝塞尔方程的第 $m$ 个零根， $\left\{\begin{array}{l} k_{m}^{a} a=x_{m}^{n}\\ m=1,2, \cdots \end{array}\right.$ , $k_m^n$ 为 $n$ 阶贝塞尔方程的第 $m$ 个本征值，第 $m$ 个本征函数 $J_n(k_m^nr)=J_n(\frac{x_m^n}{\alpha}r)$