查看“︁柱諧函數”︁的源代码

{{unreferenced|time=2014-03-21T16:18:11+00:00}}
{{NoteTA
|G1 = Math
}}

在數學中，'''柱諧函數'''是指在[[柱坐標]]中，[[拉普拉斯方程]]，<math>\nabla^2 V (\rho, \varphi, z) = 0</math> ，的一系列的解。每一個柱諧函數 <math>V_{n,k}(\rho,\varphi,z)</math> 都是三個函數的積：

:<math>V_{n,k}(\rho,\varphi,z)=P_{n,k}(\rho)\Phi_n(\varphi)Z_k(z)\,</math>

其中 <math>(\rho,\varphi,z)</math> 是柱坐標下的坐標（分別為半徑、極角和高度），而 ''n'' 和 ''k'' 則是兩個常數，用以區分不同的柱諧函數。所有的柱諧函數一起，組成一組正交完備的[[基底]]，任何一個拉普拉斯方程的解都可以寫成這些函數的[[線性組合]]。

有時候，柱諧函數也用來指代[[貝塞爾函數]]（柱諧函數最重要的組成部分）。

== 定義 ==
=== 本徵方程的推導 ===
柱坐標下的拉普拉斯方程為：

:<math>\nabla^2 V=\frac{1}{\rho} \frac{\partial}{\partial \rho} \left( \rho \frac{\partial V}{\partial \rho} \right) + \frac{1}{\rho^2} \frac{\partial^2 V}{\partial \varphi^2} + \frac{\partial^2 V}{\partial z^2} =0</math>

使用[[分離變數法]]，設：

:<math>V=P(\rho)\,\Phi(\varphi)\,Z(z)</math>

代入拉普拉斯方程，得到：

:<math>
\frac{\Phi Z}{\rho} \frac{d}{d\rho}(\rho\frac{dP}{d\rho}) + \frac{PZ}{\rho^2}\frac{d^2 \Phi}{d \varphi^2} + P \Phi\frac{d^2 Z}{d z^2} = 0</math>

分離變數後，可以寫成：

:<math>\begin{cases}
\frac{1}{P\rho} \frac{d}{d\rho}(\rho\frac{dP}{d\rho})+\frac{1}{\rho^2}\,(-n^2)+k^2=0  \\
\frac{1}{\Phi}(d^2\Phi/d\varphi^2)=-n^2  \\
\frac{1}{Z}(d^2Z/dz^2)=k^2
\end{cases}</math>，整理得 <math>\begin{cases}
\rho^2 P''+\rho P'+(k^2\rho^2-n^2)P=0  \\
\Phi''+n^2\Phi=0  \\
Z''-k^2Z=0
\end{cases}</math>

=== 本徵方程的求解 ===
這裡，<math>\Phi</math>是一個以<math>2\pi</math>為周期的函數，即滿足周期性邊界條件<math>\Phi(\varphi) = \Phi(\varphi+2\pi)</math>，因此<math>n</math>必須為非負整數。可以解出：
:<math>\Phi_n = \{ \cos(n\varphi) , \sin(n\varphi) \}</math>，<math>n\in\mathbb{N}</math>

或，等價地：
:<math>\Phi_n = \{ e^{in\varphi} , e^{-in\varphi} \}</math>，<math>n\in\mathbb{N}</math>

這裡，花括符表示，兩個解是[[簡併]]的。即對於一個n，方程有兩個線性獨立的解（n=0時除外）。

對於<math>Z</math>的方程，<math>k</math>可以是任意一個複數。對於一個特定的<math>k</math>，方程有兩個線性獨立的解。

若k是一個實數，則：
:<math>Z_k = \{ \cosh(kz) , \sinh(kz) \}</math>

或，等價地：
:<math>Z_k = \{ e^{kz} , e^{-kz} \}</math>


若k是一個純虛數，則：
:<math>Z_k = \{ \cos(|k|z) , \sin(|k|z) \}</math>

或，等價地：
:<math>Z_k = \{ e^{i|k|z} , e^{-i|k|z} \}</math>

對於周期性邊界條件，k取分立值；對於非周期性邊界條件，k取連續值。

而<math>P</math>的方程則是一個[[貝塞爾方程]]，它的解<math>P_{n,k}</math>形式如下。

若<math>k=0</math>，則該方程簡化為一個[[柯西－歐拉方程|歐拉方程]]：
:<math>P_{0,0} = \{1,\ln\rho\} </math>
:<math>P_{n,0} = \{\rho^n,\rho^{-n}\},n\neq0</math>

若<math>k</math>是一個非零實數，則方程的解為[[贝塞尔函数#第一类贝塞尔函数|第一類]]和/或[[贝塞尔函数#第二类贝塞尔函数（诺依曼函数）|第二類貝塞爾函數]]：
:<math>P_{n,k} = \{J_n(k\rho),Y_n(k\rho)\}</math>

若k是一個純虛數，則方程的解為[[贝塞尔函数#修正贝塞尔函数|修正貝塞爾函數]]：
:<math>P_{n,k} = \{I_n(|k|\rho),K_n(|k|\rho)\}</math>

最終，柱諧函數可以表達為以上三個函數的乘積，<math>V_{n,k} = P_{n,k}\Phi_n Z_k</math>。

=== 正交完備性 ===
柱諧函數是正交完備的。正交性是指：
:<math> \int_{0}^{\infty}d\rho \int_{0}^{2\pi}d\varphi \int_{-\infty}^{\infty}dz \left[ V_{n,k}(\rho,\theta,\varphi)V_{n',k'}(\rho,\theta,\varphi)\right] = \frac{1}{C^2_{n,k}} \,\delta_{n,n'}\,\delta_{k,k'} </math>

其中，<math>\delta_{n,n'}</math>和<math>\delta_{k,k'}</math>為[[克羅內克符號]]，<math>C_{n,k}</math>為歸一化系數。

完備性是指，對於柱坐標下的任何一個拉普拉斯方程的解均可以寫成若干個柱諧函數的線性疊加。
:<math>V(\rho,\varphi,z) = \sum_{n,k} A_{n,k}V_{n,k}</math> ，k取分立值
:<math>V(\rho,\varphi,z) = \sum_n \int dk \,A_{n}(k)V_{n,k}</math> ，k取連續值

== 參見 ==
* [[貝塞爾函數]]
* [[球諧函數]]

[[Category:原子物理学|Z]]
[[Category:量子力学|Z]]
[[Category:数学物理|Z]]
[[Category:特殊函數|Z]]