梅林变换

来自testwiki
跳转到导航 跳转到搜索

Template:Expand English 在数学中,梅林变换是一种以幂函数为核的积分变换,与双边拉普拉斯变换有密切关联。梅林变换定义式如下:

{f}(s)=φ(s)=0xs1f(x)dx.

而其逆变换为

{1φ}(x)=f(x)=12πicic+ixsφ(s)ds.

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

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

與其他變換之關係

之所以伽马函数与积分变换的理论联系密切,是因为伽马函数同时是指数函数的拉普拉斯变换和幂函数的梅林变换,这也展示了两种积分变换之间的联系。

雙邊拉普拉斯變換

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

{f}(s)={f(lnx)}(s)

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

{f}(s)={f(ex)}(s)

傅立葉變換

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

{f}(s)={f}(is)={f(lnx)}(is) 

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

{f}(s)={f(ex)}(s)={f(ex)}(is) 

範例

Cahen–Mellin 積分

對於 c>0(y)>0,且 ys在主要分支(principal branch)上,我們有

ey=12πicic+iΓ(s)ysds

其中 Γ(s)為 Γ函數。

數論

假設

(s+a)<0

我們有

f(x)={0x<1xax>1

其中

f(s)=1s+a

圓柱坐標系下的拉普拉斯算子

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

1rr(rfr)

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

2f=1rr(rfr)+1r22fθ2

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

2f=1rr(rfr)+1r22fφ2+2fz2

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

1rr(rfr)=frr+frr
(r2frr+rfr,rs)=s2(f,rs)=s2F

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

r2frr+rfr+fθθ=0

或是

1rr(rfr)+1r22fθ2=0

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

Fθθ+s2F=0

通解為

F(s,θ)=C1(s)cos(sθ)+C2(s)sin(sθ)

若給定邊界條件

f(r,θ0)=a(r),f(r,θ0)=b(r)

其梅林變換為

F(s,θ0)=A(s),F(s,θ0)=B(s)

則通解可以寫成

F(s,θ)=A(s)sin(s(θ0θ))sin(2θ0s)+B(s)sin(s(θ0+θ))sin(2θ0s)

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

1(sin(sφ)sin(2θ0s);sr)=12θ0rmsin(mφ)1+2rmcos(mφ)+r2m

其中

m=π2θ0

可以得到

f(r,θ)=rmcos(mθ)2θ00{a(x)x2m+2rmxmsin(mθ)+r2m+b(x)x2m2rmxmsin(mθ)+r2m}xm1dx

应用

因为具有Template:Link-en,梅林变换广泛应用于计算机科学的算法分析。[1]对于纯虚输入,原函数函数的梅林变换与对其进行尺度伸缩后函数的梅林变换幅度相同。尺度不变性类似于傅里叶变换的时移不变性,即原函数与对其进行时移的函数的傅里叶变换幅度相同。这一性质对图像识别非常有用:当物体与摄像机的距离发生变化时,图像尺度会发生变化。

在量子力学(特别是量子场论)中,由于动量和位置之间存在傅里叶变化的关系,傅里叶空间被广泛应用。2011年,A. Liam Fitzpatrick、Jared Kaplan、Template:Link-enTemplate:Link-en和Balt C. van Rees证明了梅林空间在AdS/CFT对偶中具有类似的效用。[2][3][4]

梅林变换表

下表展示了部分函数的梅林变化结果,在Template:HarvtxtTemplate:Harvtxt中可以找到。

梅林变化表
函数 f(x) 梅林变换

f~(s)={f}(s)

收敛域 注释
ex Γ(s) 0<s<
ex1 Γ(s) 1<s<0
ex1+x Γ(s) 2<s<1 一般来说, Γ(s)exn=0N1(1)nn!xn, for N<s<N+1的梅林变换。[5]
ex2 12Γ(12s) 0<s<
erfc(x) Γ(12(1+s))πs 0<s<
e(lnx)2 πe14s2 <s<
δ(xa) as1 <s< a>0,δ(x)狄拉克函数
u(1x)={1if0<x<10if1<x< 1s 0<s< u(x)单位阶跃函数
u(x1)={0if0<x<11if1<x< 1s <s<0
u(1x)xa={xaif0<x<10if1<x< 1s+a a<s<
u(x1)xa={0if0<x<1xaif1<x< 1s+a <s<a
u(1x)xalnx={xalnxif0<x<10if1<x< 1(s+a)2 a<s<
u(x1)xalnx={0if0<x<1xalnxif1<x< 1(s+a)2 <s<a
11+x πsin(πs) 0<s<1
11x πtan(πs) 0<s<1
11+x2 π2sin(12πs) 0<s<2
ln(1+x) πssin(πs) 1<s<0
sin(x) sin(12πs)Γ(s) 1<s<1
cos(x) cos(12πs)Γ(s) 0<s<1
eix eiπs/2Γ(s) 0<s<1
J0(x) 2s1πsin(πs/2)[Γ(s/2)]2 0<s<32 J0(x)第一类贝塞尔函数
Y0(x) 2s1πcos(πs/2)[Γ(s/2)]2 0<s<32 Y0(x)第二类贝塞尔函数
K0(x) 2s2[Γ(s/2)]2 0<s< K0(x)第二类修正贝塞尔函数

参考文献

Template:Reflist

  1. 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.
  2. Template:Cite journal
  3. Template:Cite journal
  4. A. Liam Fitzpatrick. "AdS/CFT and the Holographic S-Matrix" Template:Wayback, video lecture.
  5. Jacqueline Bertrand, Pierre Bertrand, Jean-Philippe Ovarlez. The Mellin Transform. The Transforms and Applications Handbook, 1995, 978-1420066524. ffhal-03152634f