応用数学 - 特殊関数

概要

• 誤差関数について理解する。
• ガンマ関数について理解する。
• 単位階段関数について理解する。
• デルタ関数について理解する。
• 双曲線関数について理解する。

• 誤差関数

  奇関数である。

  また、 $erfc(t)=1-erf(t)$ は相補誤差関数である。

  
  

• ガンマ関数

  階乗の概念を一般化した関数である。

  自然数nに対して、 $\mathrm{\Gamma }(n+1)=n!$ が成立する。

  
  

• 単位階段関数

  $U(t)=0(t<0)$ または $1(t\ge 0)$

  $U(t)f(t)$ により、 $t<0$ の時 $0$ になる $f(t)$ の変形を作ることができる。

  
  

• デルタ関数

  $\delta (t)=0(t\ne 0)$ または $+\mathrm{\infty }(t=0)$

  $t=0$ において、瞬間的にインパルスが加わった現象を記述できる。

  
  

• 双曲線関数

  $\sinh (ax)=\frac{e^{ax}-e^{-ax}}{2}$

  $\cosh (ax)=\frac{e^{ax}+e^{-ax}}{2}$

  $\tanh (ax)=\frac{e^{ax}-e^{-ax}}{e^{ax}+e^{-ax}}$

特殊関数とは

特殊関数(special functions)とは、特定の分野で何らかの名前や記法が定着している関数のことを指す。
何が特殊関数であるかのはっきりした定義は存在しない。

しばしば特殊関数として扱われるものには、ガンマ関数、ベッセル関数、ゼータ関数等がある。

特殊関数という用語は、初等関数の対義語ではない。
ある関数が初等関数であって同時に特殊関数とされる場合もある。

• ラプラス変換で扱う特殊関数

  ラプラス変換では、以下のような特殊関数を扱う。

  1. 誤差関数(ガウスの誤差関数)
  2. ガンマ関数(オイラーのガンマ関数)
  3. 単位階段関数(ヘビサイドの単位階段関数)
  4. デルタ関数(ディラクの衝撃関数)

誤差関数

定義1:
以下の式で定義された関数erf(t)を誤差関数(error function)という。
$erf(t)=\frac{2}{\sqrt{\pi }}\underset{0}{\overset{t}{\int }}{e}^{-x^2}dx$

定義2:
以下の式で定義された関数erfc(t)を相補誤差関数(complementary error function)という。
$erfc(t)=1-erf(t)$

図. 誤差関数のグラフ

誤差関数は、ガウスの誤差関数(Gauss errorfunction)とも呼ばれる。
ガウスの誤差関数 : $erf\left(\frac{a}{\sigma \sqrt{2}}\right)$

誤差関数には、以下のような特徴がある。

• 誤差関数は奇関数である。
• 測定値が正規分布になっていて、標準偏差が $\sigma$ 、平均が $0$ の場合、1つの測定値の誤差が $-a$ と $a$ の間になる確率は以下で与えられる。

ガンマ関数

定義:
以下の無限積分で定義された関数 $\mathrm{\Gamma }(p)$ をガンマ関数(gamma function)という。
$\mathrm{\Gamma }(p)=\underset{0}{\overset{+\mathrm{\infty }}{\int }}{x}^{p-1}{e}^{-x}dx(p>0)$

$\mathrm{\Gamma }(p)$ を $p$ の初等関数で表すことはできない。
$p$ が正の整数の時や、分母が2である分数の時は簡単に値を求めることができる。

図. ガンマ関数のグラフ

ガンマ関数は、階乗の概念を一般化した特殊関数である。
互いに同値ないくつかの定義が存在する。

1729年、オイラーにより階乗の概念の複素数への一般化として、最初に導入された。

自然数 $n$ に対して、 $\mathrm{\Gamma }(n+1)=n!$ が成立する。
階乗(factorial) : $n!=n(n-1)(n-2)...3\cdot 2\cdot 1(n:\text{ 正 の 整 数 })$

• ガンマ関数の性質
  
  • $\mathrm{\Gamma }(p+1)=p\mathrm{\Gamma }(p)(p>0)$

    $p=n-1$ とおけば、以下が成立する。

    $\mathrm{\Gamma }(n)=(n-1)\mathrm{\Gamma }(n-1)$

  • $\mathrm{\Gamma }(1)=1$
  • $\mathrm{\Gamma }(n)=(n-1)!(n=2,3,4,\cdots )$
  • $\mathrm{\Gamma }\left(\frac{1}{2}\right)=\sqrt{\pi }$
  • $\mathrm{\Gamma }\left(\frac{n}{2}\right)=\left(\frac{\sqrt{\pi }}{2^m}\right)(2m-1)(2m-3)\cdots 3\cdot 1(m,n:\text{ 正 整 数 },n=2m+1:\text{ 奇 数 })$

式 $\mathrm{\Gamma }(p+1)=p\mathrm{\Gamma }(p)$ を使用すると、 $p<0$の時でも $\mathrm{\Gamma }(p)$ を求めることができる。

図. 定義域を拡張したガンマ関数のグラフ

例題:
以下のガンマ関数の値を求めよ。
(1) $\mathrm{\Gamma }(3)$
(2) $\mathrm{\Gamma }\left(\frac{3}{2}\right)$

解答:
(1) $\mathrm{\Gamma }(3)=(3-1)!=2$

(2) 
$\begin{array}{rl}\mathrm{\Gamma }\left(\frac{3}{2}\right)& =(\frac{3}{2}-1)\mathrm{\Gamma }(\frac{3}{2}-1)\\ & =\frac{1}{2}\mathrm{\Gamma }\left(\frac{1}{2}\right)\\ & =\frac{\sqrt{\pi }}{2}\end{array}$

単位階段関数

定義:
以下で定義された関数U(t)を単位階段関数(unit step function)という。
$U(t)=\{\begin{array}{ll}0& (t<0)\\ 1& (t\ge 0)\end{array}$

単位階段関数は、ヘビサイドの階段関数(Heaviside step function)とも呼ばれる。
ただし、 $t=0$ の場合の扱いが異なるもの(例えば、 $t=0$ の値を未定義とする)をそう呼ぶ場合もある。

図. 単位階段関数のグラフ

単位階段関数 $U(t)$ は、時刻 $t$ の関数として用いられることが多い。
$t=0$ においてグラフを(別の)グラフ $f(t)$ に変えたい時、積 $U(t)f(t)$ を作れば、その関数を表現することができる。

さらに、時刻 $t=a$ において突然変化するような現象を表したいときは、 $U(t)f(t)$ をt軸へa方向だけ移動して $U(t-a)f(t-a)$ とすればよい。

図. $U(t)f(t)$ : $t=0$ から別のグラフ $f(t)$ のグラフになる関数

例題:
以下の関数のグラフを書け。
$y=U(t)\sin t$

解答:
$y=U(t)\sin t$ のグラフは以下のようになる。
ある関数に単位階段関数 $U(t)$ を乗算する $t<0$ の時、値が $0$ になる。

図. $y=U(t)\sin t$ のグラフ

デルタ関数

定義:
以下の式で表される関数δ(t)を、デルタ関数(delta function)という。
$\delta (t)=\{\begin{array}{ll}0& (t\ne 0)\\ +\mathrm{\infty }& (t=0)\end{array}$

$\underset{-\mathrm{\infty }}{\overset{+\mathrm{\infty }}{\int }}\delta (t)dt=1$

図. デルタ関数のグラフ

デルタ関数は、ディラクのデルタ関数(Dirac deltafunction)、あるいは、単位衝撃関数(unit impulsefunction)とも呼ばれる。

デルタ関数は、1930年頃ディラク(P. Dirac)によって最初に導入され、(離散的な)クロネッカーデルタ(Kronecker delta)の連続版として導入された。

デルタ関数はt = 0において、瞬間的に衝撃(インパルス)が加わった現象を記述することができる。
例えば、t = a, 2a, 3a, ... (a≧0)において衝撃が加わっている現象は、次式と下図のように表すことができる。
$\sum _{n=0}^{\mathrm{\infty }}\delta (t-na)(n:\text{ 自 然 数 })$

図. $t=a,2a,3a,...(a\ge 0)$ において衝撃が加わっている現象

• デルタ関数の性質
  以下において、aは定数とする。
  • $\delta (x-a)=\frac{d}{dx}U(x-a)$
  • $\underset{-\mathrm{\infty }}{\overset{+\mathrm{\infty }}{\int }}f(x)\delta (x-a)=f(a)$
  • $x\delta (x-a)=a\delta (x-a)$
  • $f(x)\delta (x-a)=f(a)\delta (x-a)$

双曲線関数

双曲線関数の概要

以下に示す2つの指数関数を組み合わせることによって作られる関数を考える。

• $y=e^{ax},y=e^{-ax}(a\ne 0)$

以下に示す3つを双曲線関数(hyperbolic function)という。

• $\sinh (ax)=\frac{e^{ax}-e^{-ax}}{2}$
• $\cosh (ax)=\frac{e^{ax}+e^{-ax}}{2}$
• $\tanh (ax)=\frac{e^{ax}-e^{-ax}}{e^{ax}+e^{-ax}}$

双曲線関数 $\sinh (ax),\cosh (ax),\tanh (ax)$ は、それぞれ三角関数 $\sin (ax),\cos (ax),\tan (ax)$ とよく似た性質をもつ。

特に、 $\cosh (ax)$ は、懸垂曲線(カテナリー, catenary)とも呼ばれる。
懸垂曲線は柔軟な線状のものが垂れ下がった状態を表す。(懸垂曲線は放物線と似ているが異なる)
例 : 送電線が垂れ下がった曲線。

指数関数を偶関数の部分と奇関数の部分に分けた時、次式のようになり、偶関数部分が cosh(x) で、奇関数部分が sinh(x) である。

• $e^x=\cosh (x)+\sinh (x)$
• $e^{-x}=\cosh (x)-\sinh (x)$

座標上の点( $\cosh (x),\sinh (x)$ )は、双曲線 $x^2-y^2=1$ 上の点である。
そのため以下が成立する。
${\cosh }^2(x)-{\sinh }^2(x)=1$

図. 双曲線関数のグラフ

双曲線関数の性質

(1) ${\cosh }^2(ax)-{\sinh }^2(ax)=1$
(2) $\tanh (ax)=\frac{\sinh (ax)}{\cosh (ax)}$
(3) $\sinh (x\pm y)=\sinh (x)\cosh (y)\pm \cosh (x)\sinh (y)$
(4) $\cosh (x\pm y)=\cosh (x)\cosh (y)\pm \sinh (x)\sinh (y)$
(5) $\tanh (x\pm y)=\frac{\tanh (x)\pm \tanh (y)}{1\pm \tanh (x)\tanh (y)}$

双曲線関数の微分:
(1) $\frac{d}{dx}\sinh (ax)=a\cosh (ax)$
(2) $\frac{d}{dx}\cosh (ax)=a\sinh (ax)$

双曲線関数の積分:
(1) $\int \sinh (ax)dx=\frac{\cosh (ax)}{a}+C$
(2) $\int \cosh (ax)dx=\frac{\sinh (ax)}{a}+C$

双曲線関数の例題

例題:
双曲線関数の定義に基づいて以下が成立することを示せ。
$\frac{d}{dx}\sinh (ax)=a\cosh (ax)$

解答:
$\begin{array}{rl}\frac{d}{dx}\sinh (ax)& =\frac{d}{dx}\frac{e^{ax}-e^{-ax}}{2}\\ & =\frac{a{e}^{ax}-(-a)e^{-ax}}{2}\\ & =\frac{a(e^{ax}+e^{-ax})}{2}\\ & =a\cosh (ax)\end{array}$