13

Mellin変換の基礎的知識

1154
0
$$$$

はじめに

 どうもこんにちは、🐟🍊みかん🍊🐟です。今回は、初めてMathlogを書くついでに、メリン変換に関する重要な結果である"Ramanujan Master Theorem"というものを紹介し、僕が初学の際に躓いた点を解決するにあたって考えたこと書いていこうと思います。この記事は、そこそこ長くなるので、メリン変換を知っている人は前半部分はがっつり読み飛ばしていただいて結構です。

 ちなみに筆者個人としては人名等々は元の言語で記述したい(とはいえ、アルファベットしか書けないので三次方程式解いたことで有名なعمر خیامは普段書いていません)のですが、メリンはひらがなで書くとかわいいので、カタカナで書いている場合が多かったりします。

先に断っておくと、このMathlogは「Ramanujan's master theorem」を完全に理解するという観点からは不適切である可能性があります。というのも、「完全に理解させよう」という意図ではなく、相手が「わかった気になる」くらいの記事を書くことに重きをおき、より多くの人に読みやすいと感じるような記事にしているからです。もし完全に理解したいというのであれば、より発展的な資料を探して読んでみることを強く推奨します。

Mellin変換の触り

 Mathlogを見てみた感じ、メリン変換自体を解説している記事があまりなかったので一応定義しておきます。関数$f(x)$のメリン変換は次の積分で定義されます。

メリン変換

関数$f$のメリン変換は次式で定義される。
$$ \mathcal{M}\{f(x)\}(s)=\int_0^\infty x^{s-1}f(x)dx $$

 初めて見る人では、「こんな変なものを考えて何が面白いのか」と感じる方もいるかもしれません。しかし、数学的につまらないものであれば淘汰されて消えてしまっている可能性が高いので、なにかしら面白みはありそうです。実際今回紹介する定理や、Perronの公式などに応用されています。具体的になにか積分を求めてみましょう。
$$\int_0^\infty\sin(x^2)dx$$
これはFresnel積分という有名な積分で、いくつかの求め方が存在しますが、正弦関数のメリン変換を既知とすると瞬殺できます。(後で頑張って計算します)
\begin{aligned} \int_0^\infty\sin(x^2)dx&=\frac12\int_0^\infty t^{-\frac12}\sin tdt\quad(x^2\mapsto t)\\ &=\frac12\mathcal{M} \{\sin(t)\}\left(\frac{1}{2}\right)\\ &=\frac12\Gamma\left(\frac12\right)\sin\left(\frac\pi4\right)\\ &=\sqrt{\frac\pi8} \end{aligned}
となります。ほとんど何も計算していません。留数計算でも求めることができますが、積分経路が思いつきにくい(八分円を使う)ので、かなり有効な手法だと思います。

Dirichlet積分を求める際にも利用できそうです。次のように計算できます。
\begin{aligned} \int_0^\infty\frac{\sin x}xdx&=\lim_{s\to0^+}\mathcal{M}\{\sin x\}(s)\\ &=\lim_{s\to0^+}\Gamma(s)\sin\left(\frac{s\pi}2\right)\\ &=\lim_{s\to0^+}s\Gamma(s)\frac{\sin\left(\frac{s\pi}2\right)}{\frac{s\pi}2}\frac{\pi}2\\ &=\frac\pi2 \end{aligned}

どうでしょうか。理論的背景はともかくとして、計算上の有用性がありそうな気がしませんか?

Mellin変換の存在条件と逆変換

メリン変換は広義積分で与えられているので、積分値が存在しない場合はメリン変換が定まりません。従って、存在する条件を知っておく必要があります。次の事実が広く知られています。

メリン変換の存在性

局所可積分な関数$f$に対して、$f(x)=O(x^a)\quad(x\to0^+)$$f(x)=O(x^b)\quad(x\to\infty)$が成立するとき、$a<\Re s< b$の範囲で積分$\mathcal Mf(s)$は絶対収束し、その区間内で解析的である。

これは広義積分の収束条件から簡単に証明できるので証明は与えないことにします。「必要十分じゃないじゃん」と思われる方もいらっしゃるかもしれませんが、十分条件が分かっているため間違って使うことはないので大丈夫だと思います。

さて、メリン変換には逆変換が存在します。次のように与えられます。

メリン逆定理

関数$\varphi$$a<\Re s< b$の範囲で解析的であり、$\Im s\to\pm\infty$において$\varphi(s)$$0$に一様収束するとき、$a< c< b$を満たす実数$c$に対して、
$$ f(x)=\{\mathcal M^{-1}\varphi\}=\frac1{2\pi i}\int_{c-i\infty}^{c+i\infty}x^{-s}\varphi(s)ds $$
が成り立つ。

定理の条件が満たされているとき、$\varphi(s)=\mathcal Mf(s)$とおくと、Diracのdelta関数を用いて
\begin{aligned} \frac1{2\pi i}\int_{c-i\infty}^{c+i\infty}x^{-s}\varphi(s)ds&=\frac1{2\pi i}\int_{c-i\infty}^{c+i\infty}\left(\int_0^\infty f(u)u^{s-1}du\right)x^{-s}ds\\ &=\frac1{2\pi i}\int_0^\infty\int_{-\infty}^\infty\left(\frac ux\right)^{c+i\sigma}id\sigma f(u)\frac{du}{u}\quad(s\mapsto c+i\sigma, u\mapsto xe^t)\\ &=\frac1{2\pi}\int_{-\infty}^\infty\int_{-\infty}^\infty e^{t(c+i\sigma)}d\sigma f(xe^t)e^{ct}dt\\ &=\int_{-\infty}^\infty\left(\frac1{2\pi}\int_{-\infty}^\infty e^{it\sigma}d\sigma\right)f(xe^t)e^{ct}dt\\ &=\int_{-\infty}^\infty\delta(t)f(xe^t)e^{ct}dt\\ &=f(x) \end{aligned}
となるので、成立する。

この証明は個人的にかなり好きです。

Ramanujan's Master Theorem

ようやくメインディッシュです。メリン変換の定義式を見ても、結局のところ元の関数に対して定まるメリン変換というのは(実解析を考えている限り)そこまで明らかなものではなく、例えば正弦関数のメリン変換

$$ \mathcal{M}\{\sin x\}(s)=\Gamma(s)\sin\left(\frac{\pi s}2\right)\quad(0<\Re s<1) $$

にしてもかなり大変そうです。Laplace変換の有名な等式

$$ \int_0^\infty \mathcal{L}(f)(x)g(x)dx=\int_0^\infty f(x)\mathcal{L}(g)(x)dx $$

を用いれば

\begin{aligned} \mathcal{M}\{\sin x\}(s) &=\int_0^\infty\mathcal{L}^{-1}(\sin)(x)\mathcal{L}\left(\frac1{x^{1-s}}\right)(x)dx\\ &=\frac1{2\Gamma(1-s)}\int_0^\infty\frac{u^{-\frac s2-\frac12}}{1+u}du\quad(x^2\mapsto u)\\ &=\frac1{2\Gamma(1-s)}\operatorname{B}\left(\frac12-\frac s2,\frac12+\frac s2\right)\\ &=\frac1{2\Gamma(1-s)}\Gamma\left(\frac12-\frac s2\right)\Gamma\left(\frac12+\frac s2\right)\\ &=\frac{1}{2\Gamma(1-s)}\frac\pi{\sin\left(\frac\pi2+\frac{\pi s}2\right)}\\ &=\frac{\Gamma(s)\sin(\pi s)}{2\pi}\frac\pi{\cos\left(\frac{\pi s}2\right)}\\ &=\Gamma(s)\sin\left(\frac{\pi s}2\right) \end{aligned}

と計算できますが、Laplace変換を知っていないとなかなか思いつけない変形ですし、結構面倒ではありませんか?面倒ですよね?そう、面倒なんですよ。

そこで、お役立ちの定理、Ramanujan's Master Theoremが出てくるわけですね。まず紹介をします。

Ramanujan's Master Theorem

$$ f(z)=\sum_{k=0}^\infty{\varphi(k)}(-z)^k $$
と展開されるとき、そのメリン変換は
$$ \mathcal{M}\{f(x)\}(s)=\frac{\pi}{\sin{\pi s}}\varphi(-s) $$
で与えられる。

これを見てどう思いますか?僕は、「そんなわけないのでは?????」と思いました。考えてみましょう。$\varphi(s)=\frac1{s\Gamma(s)}$としてみます。このとき、$f(z)=e^{-z}$ですから積分は確かにガンマ関数に一致します。

しかしながらです。$\varphi(s)=\frac{\cos(2\pi s)}{s\Gamma(s)}$としても$f(z)$は変わりませんが、そのメリン変換の結果が変わってしまいます。不合理です。従ってこの命題は成立しないと思いました。

しかし、証明を見ないとわからないことがあると思い、証明を見てみることにしました。証明の大まかな流れとしては、

(1)最後のメリン変換の式をメリン逆変換する
(2)留数定理を使う

という流れでした。ここで、よくよく考えてみると$\frac{\pi}{\sin{\pi s}}\varphi(-s)$というのはいつでもメリン逆変換できるわけではないですよね。それを考えて厳密に考え直すと定理の内容を修正する必要があります。要するに、定理の主張だけ何故だか収束のことをあまり気にしていないんですね。なので、追加するべき定理の条件としては、

(1)$\frac{\pi}{\sin{\pi s}}\varphi(-s)$がメリン逆変換可能な単葉関数

とでもなるでしょうか。このようにすると、不安はある程度解消されそうな感覚はしませんか?例えば先の余弦関数を掛ける例だと、$\varphi(s)=\cos(2\pi s)\Gamma(s)$$\Im s\to\infty$で0に一様収束してくれないんですよね。実際に確かめてみると、次の主張が従うことを示すことができます。(証明が間違っていたら申し訳ございません)

二つの関数$\psi, \varphi$の通常型母関数が収束し、$^\forall k\in\mathbb{Z}_{\ge0}[\psi(k)=\varphi(k)]$が成り立ちそれぞれがメリン逆変換可能であれば、$\psi=\varphi$となる。

証明はチャレンジしてみてください。(僕だったら自分で証明しないときっと納得しなかったので)ともかく、上記の命題を認める限り先ほど出した「余弦関数を掛けてやったら違うの作れるじゃん」という問題は完全に消し去れます。やったね!!!

ちなみに、精密な議論ではない直感的な証明を記すと、次のようにすることができます。(一般的に正しくない議論がいくつか含まれています)

Ramanujan's master theorem

$$ \varphi(s)=\sum_{a\ge0}\lambda(a)a^s $$
という形で書けるものとすると、
\begin{aligned} \int_0^\infty t^{s-1}f(t)dt &=\int_0^\infty t^{s-1}\sum_{k=0}^\infty\varphi(k)(-t)^kdt\\ &=\int_0^\infty t^{s-1}\sum_{k=0}^\infty\sum_{a\ge0}\lambda(a)a^k(-t)^kdt \\ &=\sum_{a\ge0}\lambda(a)\int_0^\infty t^{s-1}\sum_{k=0}^\infty(-at)^kdt\\ &=\sum_{a\ge0}\lambda(a)\int_0^\infty\frac{t^{s-1}}{1+at}dt\\ &=\sum_{a\ge0}\lambda(a)a^{-s}\int_0^\infty\frac{u^{s-1}}{1+u}du &(u\mapsto at)\\ &=\frac{\pi}{\sin(\pi s)}\varphi(-s) \end{aligned}

だいぶ直感的な導出ですね。

よく使う等式

いくつかメリン変換を求めるうえでよく使う等式を示しておきたいと思います。なお、ここは補足に近いです。また、$f(x)$はメリン変換可能な解析関数としておきます。

$a>0$のとき、
$$ \mathcal{M}\{f(ax)\}(s)=a^{-s}\mathcal{M}\{f(x)\}(s) $$

積分の表示において$ax\mapsto t$とすることで得られる。

$\nu\ne0$に対して
$$ \mathcal{M}\{f(x^\nu)\}(s)=\frac1{\lvert\nu\rvert}\mathcal{M}\{f(x)\}\left(\frac s\nu\right) $$

積分の表示において$x\nu\mapsto t$とすることで得られる。

$$ \mathcal{M}\{x^\nu f(x)\}(s)=\mathcal{M}\{f(x)\}\left(s+\nu\right) $$

定義から明らかです。

$$ \mathcal{M}\{\ln f(x)\}(s)=\frac d{ds}\mathcal{M}\{f(x)\}\left(s\right) $$

関係式
$$ \frac d{ds}x^s=\ln sx^s $$
を用いて、(偏)微分と積分を交換して得ます。

次のものは英語版のWikipediaにも載っていなかったので需要はあると思います。

$$ \mathcal{M}\{f(-x)\}(s)=\cos(\pi s)\mathcal{M}\{f(x)\}\left(s\right) $$

$f(x)$は解析的なので
$$ f(x)=\sum_{n=0}^\infty \varphi(n)x^n $$
としてみると、
\begin{aligned} f(-x)&=\sum_{n=0}^\infty \varphi(n)(-x)^n\\ f(x)&=\sum_{n=0}^\infty \cos(\pi n)\varphi(n)(-x)^n \end{aligned}
となるので、Ramanujan's master theroemを用いて得る。

よく使うのはこれくらいです。興味があれば英語版WikipediaのMellin変換の公式を見ておくといいと思います。

具体例

具体例に関して記述しようとは思ったのですが、おそらく僕がここに書くよりも、こちらの Wikipedia(Ramanujan's master theorem) やそのExternal linksの欄に掲載されているpdfを見ていただくほうがより多くの具体例に触れられると思います。

一応二つの例を示しておきます。最初に示した正弦関数のメリン変換と、超幾何関数のメリン変換について触れておきます。

正弦関数のほうに関しては簡単に計算できて、
\begin{aligned} \mathcal{M}\{\sin(x)\}(s) &=\mathcal{M}\left\{x\sum_{n=0}^\infty\frac1{\Gamma(2n+2)}(-x)^{2n}\right\}(s)\\ &=\frac12\mathcal{M}\left\{\sum_{n=0}^\infty\frac1{\Gamma(2n+2)}(-x)^n\right\}\left(\frac{s+1}2\right)\\ &=\frac12\frac\pi{\sin( \frac{s+1}2\pi)}\frac1{\Gamma(1-s)}\\ &=\Gamma(s)\sin\left(\frac{\pi s}2\right) \end{aligned}
となります。 超幾何関数(Wikipedia) に関しては、よく使う等式とPochhammer記号のガンマ関数表記
$$ (x)_n=\frac{\Gamma(x+n)}{\Gamma(x)} $$
を考えると、$\Gamma(\boldsymbol{a})=\prod_k\Gamma(a_k)$と略記すれば即座に以下の結果を得ます。

超幾何関数のメリン変換

次の等式が成り立つ。
$$ \int_0^\infty {_rF_s}\left[\begin{array}{c} \boldsymbol c\\\boldsymbol d\end{array};-x\right]x^{s-1}dx=\Gamma(s)\frac{\Gamma(\boldsymbol{c}-s)\Gamma(\boldsymbol{d})}{\Gamma(\boldsymbol{c})\Gamma(\boldsymbol{d}-s)} $$

多くの関数は超幾何級数で表すことができるので、超幾何級数表示が判明している関数のメリン変換は(命題5から命題9の結果も併せて考えれば)瞬殺できることになります。逆に、メリン変換の結果がガンマ関数で書けているなら、元の関数の超幾何級数表示もわかりますね。

おわりに

初めて書いたmathlogで慣れないことも多くありましたが、「わかった気になる」をお届けするという観点からはうまくできていると嬉しいです。最後まで読んでいただきありがとうございました!

投稿日:2023226

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

怠惰なB1らしいです.どうしてB1になってしまったのであろうか

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中