Mellin变换

定义

设 $f(x)$ 是定义在 $(0,+\infty)$ 上的函数，$s=\sigma+i\tau$ 为复参数。若积分在某个竖直条带内收敛，则 $f$ 的 Mellin 变换写作 $\mathcal{M}f(s)=\int_{0}^{\infty}x^{s-1}f(x)\,dx$，其中 $x^{s-1}$ 即变换核，是一族幂函数^[1][2]。

不同的文献在指数写法上并不统一：有的把核写成 $x^{-s}$，即 $\int_{0}^{\infty}F(x)x^{-s}dx$，两者只相差参数 $s$ 的符号替换，性质与结论彼此对应^[3]。中文文献一般把这种变换称为梅林变换，并称 $s=\sigma+i\tau$ 为复数，像函数记作 $M(s)$^[4]。

反演公式为 $f(x)=\frac{1}{2\pi i}\int_{c-i\infty}^{c+i\infty}x^{-s}\mathcal{M}f(s)\,ds$，积分路径是复平面上一条竖直线，$c$ 需落在像函数解析的条带之内^[1][2]。原函数与像函数互为反演的充分条件由 Mellin 反演定理给出^[5]。

原理

变换存在的区域通常是平行于虚轴的一条无限长竖带 $a<\Re s<b$，称为基本条带。在端点 $x\to 0$ 与 $x\to +\infty$ 处，$x^{s-1}$ 对函数分别起放大或压制作用，因此参数实部的上下界由函数在两端的增长阶数决定^[2]。当 $f$ 在每个闭区间上分段连续，且 $\int_{0}^{1}|f(t)|t^{\gamma_1-1}dt$ 与 $\int_{1}^{\infty}|f(t)|t^{\gamma_2-1}dt$ 都有限时，变换在 $\gamma_1<\Re s<\gamma_2$ 内收敛^[6]。

作变量替换 $x=e^{t}$ 后，$\int_{0}^{\infty}x^{s-1}f(x)dx=\int_{-\infty}^{\infty}f(e^{t})e^{st}dt$，即 Mellin 变换化为对数坐标下的双边拉普拉斯变换；把 $s$ 换成 $i\xi$ 则化为傅里叶变换。Mellin 变换因此可看作双边拉普拉斯变换的乘法形式，二者可由换元相互导出^[1][7]。

核心性质是尺度律：若 $g(x)=f(ax)$，$a>0$，则 $\mathcal{M}g(s)=a^{-s}\mathcal{M}f(s)$。也就是说，原变量的伸缩在像域变成乘以幂因子，这与傅里叶变换把平移变成相位因子相平行。此外还有 $\mathcal{M}[x^{\alpha}f](s)=\mathcal{M}f(s+\alpha)$，以及在边界项消失时 $\mathcal{M}[xf'(x)](s)=-s\,\mathcal{M}f(s)$，欧拉微分算子 $x\frac{d}{dx}$ 由此化为乘法^[8]。

对 Mellin 卷积 $(f*g)(y)=\int_{0}^{\infty}f(x)g\left(\frac{y}{x}\right)\frac{dx}{x}$，变换把卷积变成像函数的普通乘积，$\mathcal{M}[f*g]=\mathcal{M}f\cdot\mathcal{M}g$，这使含卷积的积分方程可以在像域中按代数方式处理^[9][2]。伽马函数是最基本的例子：$\Gamma(s)=\mathcal{M}[e^{-x}](s)$，收敛域为 $0<\Re s<\infty$^[1]。

发展历程

变换及其反演式最早出现在黎曼研究素数分布的论文中，其后 Cahen 给出了更明确的表述；把这些关系系统化、并指出其若干应用的是芬兰数学家 Hjalmar Mellin，变换因此以其姓氏命名^[10]。Mellin 本人从无穷乘积与伽马函数的关系入手，进而研究超几何函数与黎曼 ζ 函数，并在复平面上使用围道积分与留数方法^[11]。

关于发表年份，文献记载并不完全一致：一种说法是 Mellin 于 1897 年在芬兰科学学会的刊物上引入这一变换^[12]。另有研究指出，定义变换与反演的那一对关系实际出现在他更早的论文中，两式相隔数页，而他本人并未写过专门讨论最一般情形的文章^[11]。

1917 年，Hardy 在引用 Mellin 论文的《On Mellin's inversion formula》中，把这对关系的归属明确归于 Mellin，此后这一定名被普遍沿用^[11]。1921 年，Fowler 把它引入理论物理问题；1936 年，Meijer 将 Mellin 关于超几何函数的工作推广为 Meijer G 函数，后者本身即以 Gamma 函数之比的逆 Mellin 变换隐式定义^[11][1]。到 20 世纪 30 年代中期，这一变换开始进入分析学教科书^[11]。

20 世纪 70 年代，Slater 定理得到完整证明并用于积分计算，Marichev 等人将其完善为可操作的算法，成为符号积分软件处理大量定积分的基础之一^[13]。1990 年代，Flajolet 与 Sedgewick 等人系统地把 Mellin 变换用于算法平均情况分析，Knuth 与 De Bruijn 较早将其引入这一领域^[14][15]。

应用

在解析数论中，Mellin 变换与狄利克雷级数天然对应，可用于证明黎曼 ζ 函数与素数计数函数有关的函数方程，佩龙公式也是这一联系的体现^[1]。像函数解析延拓后的极点位置与留数，能够刻画出原函数在 $x\to 0$ 或 $x\to +\infty$ 时的渐近主项，这是渐近展开理论中的常用手段^[2]。

在概率论中，Mellin 变换用于研究随机变量之积的分布：若 $X$ 与 $Y$ 独立，则乘积的 Mellin 变换等于两者变换之积，这一性质使其成为处理乘法型随机模型的有力工具^[5]。

在计算机科学中，由于变换模在纯虚参数下对尺度伸缩不变，Mellin 变换被广泛用于算法分析，可处理数字查找树、四叉树等结构的平均代价的高阶渐近估计^[1][15]。尺度不变性类似于傅里叶变换模的时移不变性，因而在图像识别中也有用处：物体与相机距离变化时图像尺度随之改变，而变换模保持不变^[1]。

在特殊函数与积分计算方面，许多常见函数的 Mellin 变换都含伽马函数或贝塔函数，由此发展出 Mellin–Barnes 积分，并在计算机代数系统中通过 Meijer G 函数快速计算大量定积分^[1]。求解极坐标或柱坐标下的拉普拉斯方程时，Mellin 变换可把径向部分化为常系数方程^[1]。2011 年，Fitzpatrick 等人在 AdS/CFT 对偶的研究中指出 Mellin 空间具有类似于傅里叶空间的作用^[1]。

局限

Mellin 变换只作用于正半轴上的函数，定义域因此天然受限；复平面上双侧都有定义的函数需要拆分为正负部分分别处理，才能定义其 Mellin 变换^[5]。

积分在两端都需收敛，仅当参数落在基本条带内变换才存在，条带为空时变换无定义；这与拉普拉斯变换、傅里叶变换在适用条件上的限制类似^[2][14]。

反演要求像函数在竖直线 $\Re s=c$ 上绝对可积，并通过沿竖直线的复积分还原原函数，计算通常依赖留数定理和围道变形，对像函数的解析延拓与极点分布有较高要求^[2][6]。

此外，这一变换长期处于拉普拉斯变换与傅里叶变换的次要位置，专门文献相对分散；历史文献中关于 Mellin 本人发表年月的记载也多有出入，引用时需核对原始刊物^[11]。

参见

• 拉普拉斯变换 —— 双边拉普拉斯变换经换元即可与 Mellin 变换相互转化。

• 傅里叶变换 —— Mellin 变换在对数坐标下即为一种傅里叶变换。

• 伽马函数 —— 指数函数 $e^{-x}$ 的 Mellin 变换恰为伽马函数。

• 黎曼ζ函数 —— 其函数方程可借助 Mellin 变换证明。

• 狄利克雷级数 —— 与 Mellin 变换存在直接的对应关系。

• Meijer G函数 —— 由 Gamma 函数之比的逆 Mellin 变换隐式定义。