Addition-subtraction chain

In applied mathematics, the Kelvin functions ber_ν(x) and bei_ν(x) are the real and imaginary parts, respectively, of

$J_{\nu }(x{e}^{3\pi i/4}),$

where x is real, and J_ν(z), is the ν^th order Bessel function of the first kind. Similarly, the functions Ker_ν(x) and Kei_ν(x) are the real and imaginary parts, respectively, of $K_{\nu }(x{e}^{\pi i/4})$, where $K_{\nu }(z)$ is the ν^th order modified Bessel function of the second kind.

These functions are named after William Thomson, 1st Baron Kelvin.

While the Kelvin functions are defined as the real and imaginary parts of Bessel functions with x taken to be real, the functions can be analytically continued for complex arguments x e^i φ, φ ∈ [0, 2π). With the exception of Ber_n(x) and Bei_n(x) for integral n, the Kelvin functions have a branch point at x = 0.

ber(x)

File:KelvinFunctionBer.png

ber(x) for x between 0 and 10.

File:KelvinFunctionBerNorm.png

$\mathrm{b}\mathrm{e}\mathrm{r}(x)/e^{x/\sqrt{2}}$ for $x$ between 0 and 100.

For integers n, ber_n(x) has the series expansion

${\mathrm{b}\mathrm{e}\mathrm{r}}_n(x)={\left(\frac{x}{2}\right)}^n\sum _{k\ge 0}\frac{\cos \left[(\frac{3n}{4}+\frac{k}{2})\pi \right]}{k!\mathrm{\Gamma }(n+k+1)}{\left(\frac{x^2}{4}\right)}^k$

where Γ(z) is the Gamma function. The special case ber₀(x), commonly denoted as just ber(x), has the series expansion

$\mathrm{b}\mathrm{e}\mathrm{r}(x)=1+\sum _{k\ge 1}\frac{(-1{)}^k(x/2{)}^{4k}}{[(2k)!{]}^2}$

and asymptotic series

$\mathrm{b}\mathrm{e}\mathrm{r}(x)\sim \frac{e^{\frac{x}{\sqrt{2}}}}{\sqrt{2\pi x}}[f_1(x)\cos \alpha +g_1(x)\sin \alpha ]-\frac{\mathrm{k}\mathrm{e}\mathrm{i}(x)}{\pi }$,

where $\alpha =x/\sqrt{2}-\pi /8$, and

$f_1(x)=1+\sum _{k\ge 1}\frac{\cos (k\pi /4)}{k!(8x{)}^k}\prod _{l=1}^k(2l-1{)}^2$

$g_1(x)=\sum _{k\ge 1}\frac{\sin (k\pi /4)}{k!(8x{)}^k}\prod _{l=1}^k(2l-1{)}^2$

bei(x)

File:KelvinFunctionBei.png

bei(x) for $x$ between 0 and 10.

File:KelvinFunctionBeiNorm.png

$\mathrm{b}\mathrm{e}\mathrm{i}(x)/e^{x/\sqrt{2}}$ for $x$ between 0 and 100.

For integers n, bei_n(x) has the series expansion

${\mathrm{b}\mathrm{e}\mathrm{i}}_n(x)={\left(\frac{x}{2}\right)}^n\sum _{k\ge 0}\frac{\sin \left[(\frac{3n}{4}+\frac{k}{2})\pi \right]}{k!\mathrm{\Gamma }(n+k+1)}{\left(\frac{x^2}{4}\right)}^k$

where Γ(z) is the Gamma function. The special case bei₀(x), commonly denoted as just bei(x), has the series expansion

$\mathrm{b}\mathrm{e}\mathrm{i}(x)=\sum _{k\ge 0}\frac{(-1{)}^k(x/2{)}^{4k+2}}{[(2k+1)!{]}^2}$

and asymptotic series

$\mathrm{b}\mathrm{e}\mathrm{i}(x)\sim \frac{e^{\frac{x}{\sqrt{2}}}}{\sqrt{2\pi x}}[f_1(x)\sin \alpha -g_1(x)\cos \alpha ]-\frac{\mathrm{k}\mathrm{e}\mathrm{r}(x)}{\pi }$,

where $\alpha$, $f_1(x)$, and $g_1(x)$ are defined as for ber$(x)$.

ker(x)

For integers n, ker_n(x) has the (complicated) series expansion

$\begin{array}{rl}{\mathrm{k}\mathrm{e}\mathrm{r}}_n(x)& =\frac{1}{2}{\left(\frac{x}{2}\right)}^{-n}\sum _{k=0}^{n-1}\cos \left[(\frac{3n}{4}+\frac{k}{2})\pi \right]\frac{(n-k-1)!}{k!}{\left(\frac{x^2}{4}\right)}^k-\ln \left(\frac{x}{2}\right){\mathrm{b}\mathrm{e}\mathrm{r}}_n(x)+\frac{\pi }{4}{\mathrm{b}\mathrm{e}\mathrm{i}}_n(x)\\ & +\frac{1}{2}{\left(\frac{x}{2}\right)}^n\sum _{k\ge 0}\cos \left[(\frac{3n}{4}+\frac{k}{2})\pi \right]\frac{\psi (k+1)+\psi (n+k+1)}{k!(n+k)!}{\left(\frac{x^2}{4}\right)}^k\end{array}$

File:KelvinFunctionKer.png

ker(x) for $x$ between 0 and 10.

File:KelvinFunctionKerNorm.png

$\mathrm{k}\mathrm{e}\mathrm{r}(x)e^{x/\sqrt{2}}$ for x between 0 and 100.

where $\psi (z)$ is the Digamma function. The special case ker$0_{}(x)$, commonly denoted as just ker$(x)$, has the series expansion

$\mathrm{k}\mathrm{e}\mathrm{r}(x)=-\ln \left(\frac{x}{2}\right)\mathrm{b}\mathrm{e}\mathrm{r}(x)+\frac{\pi }{4}\mathrm{b}\mathrm{e}\mathrm{i}(x)+\sum _{k\ge 0}(-1{)}^k\frac{\psi (2k+1)}{[(2k)!{]}^2}{\left(\frac{x^2}{4}\right)}^{2k}$

and the asymptotic series

$\mathrm{k}\mathrm{e}\mathrm{r}(x)\sim \sqrt{\frac{\pi }{2x}}{e}^{-\frac{x}{\sqrt{2}}}[f_2(x)\cos \beta +g_2(x)\sin \beta ],$

where $\beta =x/\sqrt{2}+\pi /8$, and

$f_2(x)=1+\sum _{k\ge 1}(-1{)}^k\frac{\cos (k\pi /4)}{k!(8x{)}^k}\prod _{l=1}^k(2l-1{)}^2$

$g_2(x)=\sum _{k\ge 1}(-1{)}^k\frac{\sin (k\pi /4)}{k!(8x{)}^k}\prod _{l=1}^k(2l-1{)}^2.$

kei(x)

For integers n, kei_n(x) has the (complicated) series expansion

${\mathrm{k}\mathrm{e}\mathrm{i}}_n(x)=-\frac{1}{2}{\left(\frac{x}{2}\right)}^{-n}\sum _{k=0}^{n-1}\sin \left[(\frac{3n}{4}+\frac{k}{2})\pi \right]\frac{(n-k-1)!}{k!}{\left(\frac{x^2}{4}\right)}^k-\ln \left(\frac{x}{2}\right){\mathrm{b}\mathrm{e}\mathrm{i}}_n(x)-\frac{\pi }{4}{\mathrm{b}\mathrm{e}\mathrm{r}}_n(x)+\frac{1}{2}{\left(\frac{x}{2}\right)}^n\sum _{k\ge 0}\sin \left[(\frac{3n}{4}+\frac{k}{2})\pi \right]\frac{\psi (k+1)+\psi (n+k+1)}{k!(n+k)!}{\left(\frac{x^2}{4}\right)}^k$

File:KelvinFunctionKei.png

kei(x) for $x$ between 0 and 10.

File:KelvinFunctionKeiNorm.png

$\mathrm{k}\mathrm{e}\mathrm{i}(x)e^{x/\sqrt{2}}$ for $x$ between 0 and 100.

where $\psi (z)$ is the Digamma function. The special case kei$0_{}(x)$, commonly denoted as just kei$(x)$, has the series expansion

$\mathrm{k}\mathrm{e}\mathrm{i}(x)=-\ln \left(\frac{x}{2}\right)\mathrm{b}\mathrm{e}\mathrm{i}(x)-\frac{\pi }{4}\mathrm{b}\mathrm{e}\mathrm{r}(x)+\sum _{k\ge 0}(-1{)}^k\frac{\psi (2k+2)}{[(2k+1)!{]}^2}{\left(\frac{x^2}{4}\right)}^{2k+1}$

and the asymptotic series

$\mathrm{k}\mathrm{e}\mathrm{i}(x)\sim -\sqrt{\frac{\pi }{2x}}{e}^{-\frac{x}{\sqrt{2}}}[f_2(x)\sin \beta +g_2(x)\cos \beta ],$

where $\beta$, $f_2(x)$, and $g_2(x)$ are defined as for ker$(x)$.

See also

• Bessel function

References

• Template:Abramowitz Stegun ref
• Template:Dlmf

External links

• Weisstein, Eric W. "Kelvin Functions." From MathWorld—A Wolfram Web Resource. [1]
• GPL-licensed C/C++ source code for calculating Kelvin functions at codecogs.com: [2]