查看“︁梅林变换”︁的源代码

{{expand English | Mellin transform | time = 2020-4-28}}
在数学中，'''梅林变换'''是一种以幂函数为核的[[积分变换]]，与[[双边拉普拉斯变换]]有密切关联。梅林变换定义式如下：

:<math>\left\{\mathcal{M}f\right\}(s) = \varphi(s)=\int_0^{\infty} x^{s-1} f(x)dx.</math>

而其逆变换为
:<math>\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.</math>

梅林变换有许多应用。出于它与[[狄利克雷级数]]的联系，它也被用以证明[[黎曼ζ函数]]与[[素数计数函数]]有关的的[[黎曼ζ函数#函数方程|函数方程]]；进一步地，它也与解析数论有关，如在[[佩龙公式]]中。

同时，它与[[伽马函数]]密切相关，很多常见函数的梅林变换中都需要用到伽马函数或它衍生出的[[贝塔函数]]；这使得它被运用在[[巴恩斯积分|梅林-巴恩斯积分]]和[[广义超几何函数|超几何函数]]的理论中，衍生出了在计算机代数系统中使用的，可以快速计算大量定积分的[[Meijer_G-函数]]。

== 與其他變換之關係 ==
之所以伽马函数与积分变换的理论联系密切，是因为伽马函数同时是[[指数函数]]的拉普拉斯变换和[[幂函数]]的梅林变换，这也展示了两种积分变换之间的联系。
=== 雙邊拉普拉斯變換 ===
雙邊拉普拉斯變換可以用梅林變換來表示，如下式

:<math> \left\{\mathcal{B} f\right\}(s) = \left\{\mathcal{M} f(-\ln x) \right\}(s)</math>

梅林變換也可以用雙邊拉普拉斯變換來表示，如下式

:<math>\left\{\mathcal{M} f\right\}(s) = \left\{\mathcal{B} f(e^{-x})\right\}(s)</math>

=== 傅立葉變換 ===
傅立葉變換可以用梅林變換來表示，如下式

:<math>\left\{\mathcal{F} f\right\}(-s) = \left\{\mathcal{B} f\right\}(-is) 
= \left\{\mathcal{M} f(-\ln x)\right\}(-is)\ </math>

梅林變換變換也可以用傅立葉來表示，如下式

:<math>\left\{\mathcal{M} f\right\}(s) = \left\{\mathcal{B} f(e^{-x}) \right\}(s) = \left\{\mathcal{F} f(e^{-x})\right\}(-is)\ </math>

== 範例 ==

=== Cahen–Mellin 積分 ===
對於 <math>c>0</math>，<math>\Re(y)>0</math>，且 <math>y^{-s}</math>在主要分支（principal branch）上，我們有

:<math> e^{-y}= \frac{1}{2\pi i} \int_{c-i\infty}^{c+i\infty} \Gamma(s) y^{-s} \; ds </math>

其中 <math>\Gamma(s)</math>為 Γ函數。

=== 數論 ===
假設

:<math>\Re (s+a)<0</math>

我們有

:<math>
f(x)=\begin{cases} 0  & x < 1 \\ x^a & x > 1 \end{cases}
</math>

其中

:<math>
\mathcal M f (s)= - \frac 1 {s+a}
</math>

== 圓柱坐標系下的拉普拉斯算子 ==
在任何維度的圓柱坐標系中，拉普拉斯算子總是會包含下式

:<math>\frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial f}{\partial r} \right)</math>

例如，拉普拉斯算子在二維空間的極坐標表示法

:<math>\nabla^2 f = \frac{1}{r} \frac{\partial}{\partial r} \left(r \frac{\partial f}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 f}{\partial \theta^2}</math>

或是在三維空間的柱坐標表示法

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

而利用梅林變換可以很簡單的處理此項

:<math>\frac{1}{r} \frac{\partial}{\partial r} \left( r \frac{\partial f}{\partial r} \right) = f_{rr} + \frac{f_r}{r}</math>

:<math>\mathcal M \left(r^2 f_{rr} + r f_r,  r \to s \right) = s^2 \mathcal M \left(f,  r \to s \right) = s^2 F</math>

舉例來說，二維拉普拉斯方程的極坐標表示法具有以下形式

:<math> r^2 f_{rr} +  r f_r + f_{\theta \theta} = 0</math>

或是

:<math>\frac{1}{r} \frac{\partial}{\partial r} \left(r \frac{\partial f}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 f}{\partial \theta^2} = 0</math>

利用梅林變換，可以轉換成一個簡諧振子的形式

:<math> F_{\theta \theta} + s^2 F = 0</math>

通解為

:<math> F (s, \theta) = C_1(s) \cos (s\theta) + C_2(s) \sin (s \theta)</math>

若給定邊界條件

:<math> f(r,-\theta_0) = a(r), \quad  f(r,\theta_0) = b(r) </math>

其梅林變換為

:<math> F(s,-\theta_0) = A(s), \quad  F(s,\theta_0) = B(s) </math>

則通解可以寫成

:<math> F (s, \theta) = A(s) \frac {\sin(s (\theta_0 - \theta))}{\sin (2 \theta_0 s)}+ B(s) \frac {\sin(s (\theta_0 + \theta))}{\sin (2 \theta_0 s)}</math>

最後利用逆變換以及卷積定理

:<math>\mathcal M^{-1} \left( \frac {\sin (s \varphi)}{\sin (2 \theta_0 s)}; s \to r \right) = \frac 1 {2 \theta_0} \frac{r^m \sin (m \varphi)}{1+2r^m \cos(m \varphi) + r^{2m}}</math>

其中

:<math>m= \frac \pi {2 \theta_0}</math>

可以得到

:<math> f(r, \theta) = \frac{r^m \cos (m \theta)}{2 \theta_0} \int_0^\infty \left \{ \frac{a(x)}{x^{2m} + 2r^m x^m \sin(m \theta) + r^{2m}} + \frac{b(x)}{x^{2m} - 2r^m x^m \sin(m \theta) + r^{2m}} \right \} x^{m-1} \, dx </math>
:

== 应用 ==
因为具有{{Link-en|尺度不变性|Scale_invariance}}，梅林变换广泛应用于计算机科学的算法分析。<ref>Philippe Flajolet and Robert Sedgewick. The Average Case Analysis of Algorithms: Mellin Transform Asymptotics. Research Report 2956. 93 pages. Institut National de Recherche en Informatique et en Automatique (INRIA), 1996.</ref>对于纯虚输入，原函数函数的梅林变换与对其进行尺度伸缩后函数的梅林变换幅度相同。尺度不变性类似于傅里叶变换的时移不变性，即原函数与对其进行时移的函数的傅里叶变换幅度相同。这一性质对[[圖型識別|图像识别]]非常有用：当物体与摄像机的距离发生变化时，图像尺度会发生变化。

在量子力学（特别是量子场论）中，由于动量和位置之间存在傅里叶变化的关系，傅里叶空间被广泛应用。2011年，A. Liam Fitzpatrick、Jared Kaplan、{{Link-en|João Penedones|João_Penedones}}、{{Link-en|Suvrat Raju|Suvrat_Raju}}和Balt C. van Rees证明了梅林空间在[[AdS/CFT对偶]]中具有类似的效用。<ref>{{Cite journal |last=Fitzpatrick |first=A. Liam |last2=Kaplan |first2=Jared |last3=Penedones |first3=Joao |last4=Raju |first4=Suvrat |last5=van Rees |first5=Balt C. |date=2011-11 |title=A Natural Language for AdS/CFT Correlators |url=http://arxiv.org/abs/1107.1499 |journal=Journal of High Energy Physics |volume=2011 |issue=11 |doi=10.1007/JHEP11(2011)095 |issn=1029-8479 |access-date=2015-11-27 |archive-date=2022-12-25 |archive-url=https://web.archive.org/web/20221225110343/https://arxiv.org/abs/1107.1499 |dead-url=no }}</ref><ref>{{Cite journal |last=Fitzpatrick |first=A. Liam |last2=Kaplan |first2=Jared |date=2012-10 |title=Unitarity and the Holographic S-Matrix |url=http://arxiv.org/abs/1112.4845 |journal=Journal of High Energy Physics |volume=2012 |issue=10 |doi=10.1007/JHEP10(2012)032 |issn=1029-8479 |access-date=2015-11-27 |archive-date=2023-11-16 |archive-url=https://web.archive.org/web/20231116215229/http://arxiv.org/abs/1112.4845 |dead-url=no }}</ref><ref>A. Liam Fitzpatrick.  [http://online.kitp.ucsb.edu/online/qgravity15/fitzpatrick/ "AdS/CFT and the Holographic S-Matrix"] {{Wayback|url=http://online.kitp.ucsb.edu/online/qgravity15/fitzpatrick/ |date=20180929091045 }}, video lecture.</ref>

== 梅林变换表 ==
下表展示了部分函数的梅林变化结果，在{{harvtxt|Bracewell|2000}}与{{harvtxt|Erdélyi|1954}}中可以找到。
{| class="wikitable"
|+梅林变化表
!函数 <math>f(x)</math>
!梅林变换
<math>\tilde{f}(s) = \mathcal{M}\{f\}(s)</math>
!收敛域
!注释
|-
|<math> e^{-x} </math>
|<math> \Gamma(s) </math>
|<math> 0 < \Re s < \infty </math>
|
|-
|<math> e^{-x}-1 </math>
|<math> \Gamma(s) </math>
|<math> -1 < \Re s < 0 </math>
|
|-
|<math> e^{-x}-1 +x </math>
|<math> \Gamma(s) </math>
|<math> -2 < \Re s < -1 </math>
|一般来说， <math>\Gamma(s)</math>是<math>e^{-x}-\sum_{n=0}^{N-1} \frac{(-1)^n}{n !} x^n, \text{ for } -N< \Re s <-N+1</math>的梅林变换。<ref>Jacqueline Bertrand, Pierre Bertrand, Jean-Philippe Ovarlez. The Mellin Transform. The Transforms and Applications Handbook, 1995, 978-1420066524. ffhal-03152634f</ref> 
|-
|<math> e^{-x^2} </math>
|<math> \tfrac{1}{2}\Gamma(\tfrac{1}{2}s) </math>
|<math> 0 < \Re s < \infty </math>
|
|-
|<math> \mathrm{erfc}(x) </math>
|<math> \frac{\Gamma(\tfrac{1}{2}(1+s))}{\sqrt{\pi}\;s} </math>
|<math> 0 < \Re s < \infty </math>
|
|-
|<math> e^{-(\ln x)^2} </math>
|<math> \sqrt{\pi} \, e^{\tfrac{1}{4}s^2} </math>
|<math> -\infty < \Re s < \infty </math>
|
|-
|<math> \delta(x-a) </math>
|<math> a^{s-1} </math>
|<math> -\infty < \Re s < \infty </math>
|<math> a>0, \; \delta(x)</math> 是[[狄拉克δ函数|狄拉克函数]]。
|-
|<math> u(1-x) = \left\{ \begin{aligned} &1 &&\;\text{if}\; 0<x<1 &\\ &0 &&\;\text{if}\; 1<x<\infty &\end{aligned} \right. </math>
|<math> \frac{1}{s} </math>
|<math> 0 < \Re s < \infty </math>
|<math>u(x)</math> 是[[单位阶跃函数]]。
|-
|<math> -u(x-1) = \left\{ \begin{aligned} &0 &&\;\text{if}\; 0<x<1 &\\ &-1 &&\;\text{if}\; 1<x<\infty &\end{aligned} \right. </math>
|<math> \frac{1}{s} </math>
|<math> -\infty < \Re s < 0 </math>
|
|-
|<math> u(1-x)\,x^a = \left\{ \begin{aligned} &x^a &&\;\text{if}\; 0<x<1 &\\ &0 &&\;\text{if}\; 1<x<\infty &\end{aligned} \right. </math>
|<math> \frac{1}{s+a} </math>
|<math> -\Re a < \Re s < \infty </math>
|
|-
|<math> -u(x-1)\,x^a = \left\{ \begin{aligned} &0 &&\;\text{if}\; 0<x<1 &\\ &-x^a &&\;\text{if}\; 1<x<\infty &\end{aligned} \right. </math>
|<math> \frac{1}{s+a} </math>
|<math> -\infty < \Re s < -\Re a </math>
|
|-
|<math> u(1-x)\,x^a \ln x = \left\{ \begin{aligned} &x^a \ln x &&\;\text{if}\; 0<x<1 &\\ &0 &&\;\text{if}\; 1<x<\infty &\end{aligned} \right. </math>
|<math> \frac{1}{(s+a)^2} </math>
|<math> -\Re a < \Re s < \infty </math>
|
|-
|<math> -u(x-1)\,x^a \ln x = \left\{ \begin{aligned} &0 &&\;\text{if}\; 0<x<1 &\\ &-x^a \ln x &&\;\text{if}\; 1<x<\infty &\end{aligned} \right. </math>
|<math> \frac{1}{(s+a)^2} </math>
|<math> -\infty < \Re s < -\Re a </math>
|
|-
|<math> \frac{1}{1+x} </math>
|<math> \frac{\pi}{\sin(\pi s)} </math>
|<math> 0 < \Re s < 1 </math>
|
|-
|<math> \frac{1}{1-x} </math>
|<math> \frac{\pi}{\tan(\pi s)} </math>
|<math> 0 < \Re s < 1 </math>
|
|-
|<math> \frac{1}{1+x^2} </math>
|<math> \frac{\pi}{2\sin(\tfrac{1}{2}\pi s)} </math>
|<math> 0 < \Re s < 2 </math>
|
|-
|<math> \ln(1+x) </math>
|<math> \frac{\pi}{s\,\sin(\pi s)} </math>
|<math> -1 < \Re s < 0 </math>
|
|-
|<math> \sin(x) </math>
|<math> \sin(\tfrac{1}{2}\pi s) \, \Gamma(s) </math>
|<math> -1 < \Re s < 1 </math>
|
|-
|<math> \cos(x) </math>
|<math> \cos(\tfrac{1}{2}\pi s) \, \Gamma(s) </math>
|<math> 0 < \Re s < 1 </math>
|
|-
|<math> e^{ix} </math>
|<math> e^{i\pi s/2} \, \Gamma(s) </math>
|<math> 0 < \Re s < 1 </math>
|
|-
|<math> J_0(x) </math>
|<math> \frac{2^{s-1}}{\pi} \, \sin(\pi s/2) \, \left[\Gamma(s/2)\right]^2 </math>
|<math> 0 < \Re s < \tfrac{3}{2} </math>
|<math> J_0(x) </math>是[[贝塞尔函数#第一类贝塞尔函数|第一类贝塞尔函数]]。
|-
|<math> Y_0(x) </math>
|<math> -\frac{2^{s-1}}{\pi} \, \cos(\pi s/2) \, \left[\Gamma(s/2)\right]^2 </math>
|<math> 0 < \Re s < \tfrac{3}{2} </math>
|<math> Y_0(x) </math>是[[贝塞尔函数#第二类贝塞尔函数（诺依曼函数）|第二类贝塞尔函数]]。
|-
|<math> K_0(x) </math>
|<math> 2^{s-2} \, \left[\Gamma(s/2)\right]^2 </math>
|<math> 0 < \Re s < \infty </math>
|<math> K_0(x) </math>是[[贝塞尔函数#修正贝塞尔函数|第二类修正贝塞尔函数]]。
|}

== 参考文献 ==
*{{cite book
  | last1 = Galambos  | first1 = Janos
  | last2 = Simonelli | first2 = Italo
  | year = 2004
  | title = Products of random variables: applications to problems of physics and to arithmetical functions
  | url = https://archive.org/details/productsofrandom00gala  | publisher = Marcel Dekker, Inc.
  | isbn = 0-8247-5402-6
  | ref = harv }}
{{reflist}}

[[Category:积分变换]]
[[Category:复分析]]