メリン変換

数学におけるメリン変換(メリンへんかん、: Mellin transform)とは、両側ラプラス変換乗法版と見なされる積分変換である。この変換はディリクレ級数の理論と密接に関連しており、数論漸近展開の理論においてよく用いられる。ラプラス変換フーリエ変換ガンマ関数特殊関数の理論と関係している。

この変換の名はフィンランドの数学者ヒャルマル・メリン(英語版)の名にちなむ。

定義

局所可積分な関数 f のメリン変換は

{ M f } ( s ) = φ ( s ) = 0 x s 1 f ( x ) d x {\displaystyle \left\{{\mathcal {M}}f\right\}(s)=\varphi (s)=\int _{0}^{\infty }x^{s-1}f(x)dx}

により定義される。 任意の小さな正の数 ϵ {\displaystyle \epsilon } に対して、 x + 0 {\displaystyle x\to +0} のとき f ( x ) = O ( x a ϵ ) {\displaystyle f(x)=O(x^{-a-\epsilon })} x + {\displaystyle x\to +\infty } のとき f ( x ) = O ( x b + ϵ ) {\displaystyle f(x)=O(x^{-b+\epsilon })} と評価できるならば、上の積分は絶対収束する。さらに、 { M f } ( s ) {\displaystyle \left\{{\mathcal {M}}f\right\}(s)} a < ( s ) < b {\displaystyle a<\Re (s)<b} で解析的な関数となる。

また、メリン逆変換は

{ M 1 φ } ( x ) = f ( x ) = 1 2 π i c i c + i x s φ ( s ) d s {\displaystyle \left\{{\mathcal {M}}^{-1}\varphi \right\}(x)=f(x)={\frac {1}{2\pi i}}\int _{c-i\infty }^{c+i\infty }x^{-s}\varphi (s)\,ds}

により定義される。記号は、複素平面上の縦軸に沿った線積分を意味している。ここで、 c a < c < b {\displaystyle a<c<b} を満たす任意の実数である。 このような逆が存在するための条件は、メリン逆定理(英語版)で与えられている。

他の変換との関係

両側ラプラス変換は、メリン変換を用いて

{ B f } ( s ) = { M f ( ln x ) } ( s ) {\displaystyle \left\{{\mathcal {B}}f\right\}(s)=\left\{{\mathcal {M}}f(-\ln x)\right\}(s)}

と表すことが出来る。反対に、メリン変換は両側ラプラス変換により

{ M f } ( s ) = { B f ( e x ) } ( s ) {\displaystyle \left\{{\mathcal {M}}f\right\}(s)=\left\{{\mathcal {B}}f(e^{-x})\right\}(s)}

と表される。

メリン変換は、積分核 xs を用いた、乗法的ハール測度 d x x {\displaystyle {\frac {dx}{x}}} についての積分と考えることが出来る。ここで d x x {\displaystyle {\frac {dx}{x}}} は拡張 x a x {\displaystyle x\mapsto ax} について不変であり、したがって d ( a x ) a x = d x x {\displaystyle {\frac {d(ax)}{ax}}={\frac {dx}{x}}} が成り立つ。一方、両側ラプラス変換は加法的ハール測度 d x {\displaystyle dx} についての積分と考えられる。ここで d x {\displaystyle dx} は移動不変であり、したがって d ( x + a ) = d x {\displaystyle d(x+a)=dx} が成り立つ。

同様にフーリエ変換もメリン変換を用いて表すことが出来、またその逆も出来る。もし両側ラプラス変換を上述のように定義するなら、

{ F f } ( s ) = { B f } ( i s ) = { M f ( ln x ) } ( i s ) {\displaystyle \left\{{\mathcal {F}}f\right\}(s)=\left\{{\mathcal {B}}f\right\}(is)=\left\{{\mathcal {M}}f(-\ln x)\right\}(is)}

が成立する。反対に

{ M f } ( s ) = { B f ( e x ) } ( s ) = { F f ( e x ) } ( i s ) {\displaystyle \left\{{\mathcal {M}}f\right\}(s)=\left\{{\mathcal {B}}f(e^{-x})\right\}(s)=\left\{{\mathcal {F}}f(e^{-x})\right\}(-is)}

も成立する。メリン変換はまた、ニュートン級数(英語版)二項変換(英語版)を、ポアソン-メリン-ニュートン・サイクル(英語版)の意味におけるポアソン母関数と結び付ける。

カヘン-メリン積分

c > 0 {\displaystyle c>0} ( y ) > 0 {\displaystyle \Re (y)>0} および主枝(英語版)上の y s {\displaystyle y^{-s}} に対して、

e y = 1 2 π i c i c + i Γ ( s ) y s d s {\displaystyle e^{-y}={\frac {1}{2\pi i}}\int _{c-i\infty }^{c+i\infty }\Gamma (s)y^{-s}\;ds}

が成立する。ここで Γ ( s ) {\displaystyle \Gamma (s)} ガンマ関数である。この積分はカヘン-メリン積分として知られている[1]

数論

数論における重要な応用例として、単関数 f ( x ) = { 0 x < 1 , x a x > 1 {\displaystyle f(x)={\begin{cases}0&x<1,\\x^{a}&x>1\end{cases}}} に対し

M f ( s ) = 1 s + a {\displaystyle {\mathcal {M}}f(s)={\frac {1}{s+a}}}

が成立する、ということが挙げられる。

ゼータ関数

メリン変換を用いることで、リーマンゼータ関数 ζ ( s ) {\displaystyle \zeta (s)} についての公式を得ることができる。 f ( x ) = 1 e x 1 {\displaystyle f(x)={\frac {1}{e^{x}-1}}} としたとき M f ( s ) = 0 x s 1 e x 1 d x = 0 x s 1 e x 1 e x d x = 0 x s 1 n = 1 e n x d x = n = 1 0 x s 1 e n x d x = n = 1 1 n s 0 x s 1 e x d x = n = 1 Γ ( s ) n s = Γ ( s ) ζ ( s ) {\displaystyle {\mathcal {M}}f(s)=\int _{0}^{\infty }{\frac {x^{s-1}}{e^{x}-1}}dx=\int _{0}^{\infty }{\frac {x^{s-1}e^{-x}}{1-e^{-x}}}dx=\int _{0}^{\infty }x^{s-1}\sum _{n=1}^{\infty }e^{-nx}dx=\sum _{n=1}^{\infty }\int _{0}^{\infty }x^{s-1}e^{-nx}dx=\sum _{n=1}^{\infty }{\frac {1}{n^{s}}}\int _{0}^{\infty }x^{s-1}e^{-x}dx=\sum _{n=1}^{\infty }{\frac {\Gamma (s)}{n^{s}}}=\Gamma (s)\zeta (s)} よって

ζ ( s ) = 1 Γ ( s ) 0 x s 1 e x 1 d x {\displaystyle \zeta (s)={\frac {1}{\Gamma (s)}}\int _{0}^{\infty }{\frac {x^{s-1}}{e^{x}-1}}dx}

L2 上のユニタリ作用素として

ヒルベルト空間の研究において、メリン変換は少し異なった方法で定められる。 L 2 ( 0 , ) {\displaystyle L^{2}(0,\infty )} Lp空間を参照されたい)の関数に対して、基本帯(fundamental strip)は常に 1 2 + i R {\displaystyle {\tfrac {1}{2}}+i\mathbb {R} } を含む。そのため、線形作用素 M ~ {\displaystyle {\tilde {\mathcal {M}}}}

M ~ : L 2 ( 0 , ) L 2 ( , ) , { M ~ f } ( s ) := 1 2 π 0 x 1 2 + i s f ( x ) d x {\displaystyle {\tilde {\mathcal {M}}}\colon L^{2}(0,\infty )\to L^{2}(-\infty ,\infty ),\{{\tilde {\mathcal {M}}}f\}(s):={\frac {1}{\sqrt {2\pi }}}\int _{0}^{\infty }x^{-{\frac {1}{2}}+is}f(x)\,dx}

によって定義することが出来る。言い換えると、集合

{ M ~ f } ( s ) := 1 2 π { M f } ( 1 2 i s ) {\displaystyle \{{\tilde {\mathcal {M}}}f\}(s):={\tfrac {1}{\sqrt {2\pi }}}\{{\mathcal {M}}f\}({\tfrac {1}{2}}-is)}

を定義することが出来る。この作用素は通常 M {\displaystyle {\mathcal {M}}} とシンプルに記述され、「メリン変換」と呼ばれる。しかしここでは、上での記述と区別するために M ~ {\displaystyle {\tilde {\mathcal {M}}}} を記号として用いる。このときメリン逆定理(英語版)により、 M ~ {\displaystyle {\tilde {\mathcal {M}}}} は可逆であって、その逆は

M ~ 1 : L 2 ( , ) L 2 ( 0 , ) , { M ~ 1 φ } ( x ) = 1 2 π x 1 2 i s φ ( s ) d s {\displaystyle {\tilde {\mathcal {M}}}^{-1}\colon L^{2}(-\infty ,\infty )\to L^{2}(0,\infty ),\{{\tilde {\mathcal {M}}}^{-1}\varphi \}(x)={\frac {1}{\sqrt {2\pi }}}\int _{-\infty }^{\infty }x^{-{\frac {1}{2}}-is}\varphi (s)\,ds}

と得られることが分かる。さらにこの作用素は等長であること、すなわち M ~ f L 2 ( , ) = f L 2 ( 0 , ) {\displaystyle \|{\tilde {\mathcal {M}}}f\|_{L^{2}(-\infty ,\infty )}=\|f\|_{L^{2}(0,\infty )}} がすべての f L 2 ( 0 , ) {\displaystyle f\in L^{2}(0,\infty )} に対して成立することが分かる(この性質のために係数 1 / 2 π {\displaystyle 1/{\sqrt {2\pi }}} が用いられている)。したがって、 M ~ {\displaystyle {\tilde {\mathcal {M}}}} ユニタリ作用素である。

確率論において

確率論におけるメリン変換は、確率変数の積の分布の研究によく用いられる[2]X を確率変数とし、X+ = max{X,0} をその正の部分、X − = max{−X,0} をその負の部分としたとき、X のメリン変換は

M X ( s ) = 0 x s d F X + ( x ) + γ 0 x s d F X ( x ) , {\displaystyle {\mathcal {M}}_{X}(s)=\int _{0}^{\infty }x^{s}dF_{X^{+}}(x)+\gamma \int _{0}^{\infty }x^{s}dF_{X^{-}}(x),}

として定義される[3]。ここで γ は、γ2 = 1 を満たすもの(formal indeterminate)である。この変換は、複素帯領域 D = {s: a ≤ Re(s) ≤ b}(ただしa ≤ 0 ≤ b)内のすべての s に対して存在する[3]

確率変数 X のメリン変換 M X ( i t ) {\displaystyle \scriptstyle {\mathcal {M}}_{X}(it)} は、その分布関数 FX を一意に定める[3]。確率論におけるメリン変換が持つ重要な性質として、次が挙げられる: X および Y を二つの独立な確率変数としたとき、それらの積のメリン変換は、それぞれのメリン変換の積と等しい[4]。すなわち、

M X Y ( s ) = M X ( s ) M Y ( s ) {\displaystyle {\mathcal {M}}_{XY}(s)={\mathcal {M}}_{X}(s){\mathcal {M}}_{Y}(s)}

が成立する。

応用

メリン変換は、そのスケール不変性のため、計算機科学の分野で広く用いられている。あるスケール変換を施された関数のメリン変換の絶対値は、もとの関数の絶対値と等しい。このスケール不変性は、フーリエ変換のシフト不変性とも同様である。時間に関してシフトされた関数のフーリエ変換の絶対値は、もとの関数のそれと等しい。

この性質は、画像認識を行う際に役に立つ。物体の画像は、その物体がカメラに近づいたり離れたりするだけで簡単にスケールが変わってしまうからである。

その他の例

関連項目

注釈

  1. ^ Hardy, G. H.; Littlewood, J. E. (1916). “Contributions to the Theory of the Riemann Zeta-Function and the Theory of the Distribution of Primes”. Acta Mathematica 41 (1): 119–196. doi:10.1007/BF02422942.  (See notes therein for further references to Cahen's and Mellin's work, including Cahen's thesis.)
  2. ^ Galambos & Simonelli (2004, p. 15)
  3. ^ a b c Galambos & Simonelli (2004, p. 16)
  4. ^ Galambos & Simonelli (2004, p. 23)

参考文献

  • Galambos, Janos; Simonelli, Italo (2004). Products of random variables: applications to problems of physics and to arithmetical functions. Marcel Dekker, Inc.. ISBN 0-8247-5402-6 
  • Paris, R. B.; Kaminski, D. (2001). Asymptotics and Mellin-Barnes Integrals. Cambridge University Press 
  • Polyanin, A. D.; Manzhirov, A. V. (1998). Handbook of Integral Equations. Boca Raton: CRC Press. ISBN 0-8493-2876-4 
  • Flajolet, P.; Gourdon, X.; Dumas, P. (1995). “Mellin transforms and asymptotics: Harmonic sums”. Theoretical Computer Science 144 (1-2): 3–58. 
  • Tables of Integral Transforms at EqWorld: The World of Mathematical Equations.
  • Weisstein, Eric W. "Mellin Transform". mathworld.wolfram.com (英語).

外部リンク

  • Philippe Flajolet, Xavier Gourdon, Philippe Dumas, Mellin Transforms and Asymptotics: Harmonic sums.
  • Antonio Gonzáles, Marko Riedel Celebrando un clásico, newsgroup es.ciencia.matematicas
  • Juan Sacerdoti, Funciones Eulerianas (in Spanish).
  • Mellin Transform Methods, Digital Library of Mathematical Functions, 2011-08-29, National Institute of Standards and Technology
  • Antonio De Sena and Davide Rocchesso, A FAST MELLIN TRANSFORM WITH APPLICATIONS IN DAFX