Преобразование Меллина: определение, полоса и примеры

Преобразование Меллина переводит функцию, заданную на положительной полуоси, в функцию комплексной переменной $s$. Оно естественно появляется там, где важна не сдвиговая, а масштабная симметрия: асимптотика интегралов, аналитическая теория чисел, дзета-функция, теория вероятностей произведений случайных величин. Формально это близкий родственник преобразований Фурье и Лапласа, но «работает» он не со сдвигом аргумента, а с растяжением. Разберём строгое определение, где живёт образ (фундаментальная полоса), как считается классический пример с экспонентой и как устроено обратное преобразование.

Ниже - интерактивный калькулятор: подвигайте параметр и посмотрите, как образ экспоненты совпадает с гамма-функцией, а полоса сходимости задаёт область определения.

Определение преобразования Меллина

Преобразование Меллина функции $f(x)$, заданной при $x > 0$, определяется как интеграл

$$\mathcal{M}\{f\}(s) = F(s) = \int_0^{\infty} x^{s-1} f(x)\,dx,$$

где $s = \sigma + i t$ - комплексная переменная. Ядро $x^{s-1}$ - это степенная функция, и именно она делает преобразование чувствительным к масштабу: при замене $x \to \lambda x$ интеграл лишь домножается на $\lambda^{-s}$. Это масштабная (мультипликативная) симметрия - в отличие от Фурье, где базовая симметрия сдвиговая (аддитивная).

Чтобы интеграл сходился, на функцию накладывают условия степенного поведения на обоих концах: $f(x) \sim x^{-a}$ при $x \to 0^{+}$ и $f(x) \sim x^{-b}$ при $x \to \infty$, причём $a < b$. Тогда образ определён в вертикальной полосе $a < \Re(s) < b$ комплексной плоскости.

Фундаментальная полоса сходимости

Главная особенность Меллина - образ существует не на всей плоскости, а в вертикальной полосе $\alpha < \Re(s) < \beta$, которую называют фундаментальной полосой. Её границы задаются поведением функции на нуле и на бесконечности.

Рассчитать для своей задачи

Разберём на интеграле $\int_0^{\infty} x^{s-1} f(x)\,dx$. Разобьём его на $\int_0^1$ и $\int_1^{\infty}$:

• На отрезке $(0,1)$ сходимость определяет поведение в нуле. Если $f(x) = O(x^{-\alpha})$ при $x \to 0$, то $\int_0^1 x^{s-1} f(x)\,dx$ сходится при $\Re(s) > \alpha$ - это левая граница полосы.
• На $(1, \infty)$ сходимость определяет поведение на бесконечности. Если $f(x) = O(x^{-\beta})$ при $x \to \infty$, то $\int_1^{\infty}$ сходится при $\Re(s) < \beta$ - правая граница.

Если $\alpha < \beta$, полоса непуста, и в ней $F(s)$ - аналитическая функция. Вне полосы интеграл расходится, но образ часто продолжается аналитически - именно так гамма-функция продолжается с полуплоскости $\Re(s) > 0$ на всю плоскость с полюсами в неположительных целых точках.

Образ экспоненты: гамма-функция

Канонический пример, на котором стоит проверять все формулы, - преобразование Меллина функции $f(x) = e^{-x}$:

$$\mathcal{M}\{e^{-x}\}(s) = \int_0^{\infty} x^{s-1} e^{-x}\,dx = \Gamma(s).$$

Это в точности интегральное определение гамма-функции. Здесь $e^{-x} \to 1$ при $x \to 0$ (то есть $\alpha = 0$) и убывает быстрее любой степени при $x \to \infty$ (то есть $\beta = +\infty$). Значит, фундаментальная полоса - это правая полуплоскость $\Re(s) > 0$, и именно там интеграл сходится.

Зная это, легко получить целый ворох частных значений: $\Gamma(1) = 1$, $\Gamma(n) = (n-1)!$ для натуральных $n$, $\Gamma(1/2) = \sqrt{\pi}$. Функциональное уравнение $\Gamma(s+1) = s\,\Gamma(s)$ и позволяет аналитически продолжить образ за пределы полосы сходимости.

Свойства: масштаб, степень, производная

Преобразование Меллина обладает набором правил, которые делают его удобным для асимптотики. Пусть $F(s) = \mathcal{M}\{f\}(s)$. Тогда:

$$\begin{aligned} \mathcal{M}\{f(\lambda x)\}(s) &= \lambda^{-s} F(s), \\ \mathcal{M}\{x^{a} f(x)\}(s) &= F(s + a), \\ \mathcal{M}\{f(x^{p})\}(s) &= \tfrac{1}{p}\, F\!\left(\tfrac{s}{p}\right), \\ \mathcal{M}\{x f'(x)\}(s) &= -s\, F(s). \end{aligned}$$

Первое свойство - ключевое: масштаб аргумента превращается в простой множитель $\lambda^{-s}$. Сравните с Фурье, где сдвиг аргумента даёт множитель $e^{-i\omega \tau}$. Именно поэтому Меллин - «преобразование масштаба»: задачи, инвариантные относительно растяжения, в его образе становятся алгебраическими.

Рассчитать для своей задачи

Связь с Фурье и Лапласом

Преобразование Меллина - не самостоятельный «чужак», а замена переменной в уже знакомых преобразованиях. Подстановка $x = e^{-u}$ переводит интеграл Меллина в двусторонний интеграл по $u \in (-\infty, \infty)$:

$$F(s) = \int_{-\infty}^{\infty} e^{-s u}\, f(e^{-u})\,du.$$

При $s = \sigma + i t$ это в точности двустороннее преобразование Лапласа функции $g(u) = f(e^{-u})$, взятое в точке $s$. А если зафиксировать $\sigma$ и менять $t$, получится преобразование Фурье функции $e^{-\sigma u} g(u)$. Поэтому абсцисса $\sigma = \Re(s)$ внутри фундаментальной полосы играет ту же роль, что и параметр сходимости в преобразовании Лапласа: её выбирают так, чтобы интеграл сходился.

Эта связь объясняет, почему многие свойства Меллина дословно повторяют свойства Фурье и Лапласа - они получаются друг из друга логарифмической заменой переменной.

Обратное преобразование Меллина

Восстановление функции по её образу выполняется интегралом по вертикальной прямой внутри фундаментальной полосы:

$$f(x) = \frac{1}{2\pi i} \int_{c - i\infty}^{c + i\infty} x^{-s}\, F(s)\,ds,$$

где $c$ - любое число из интервала $\alpha < c < \beta$. Это контурный интеграл по прямой $\Re(s) = c$, параллельной мнимой оси. Результат не зависит от конкретного $c$, пока прямая лежит внутри полосы сходимости.

На практике этот интеграл считают через вычеты: контур замыкают полуокружностью и собирают вклады полюсов $F(s)$. Каждый полюс даёт слагаемое асимптотического разложения $f(x)$. Например, полюсы гамма-функции в точках $s = 0, -1, -2, \dots$ дают ряд Тейлора $e^{-x} = \sum_{n \ge 0} \frac{(-1)^n}{n!} x^n$ - это и есть метод Меллина–Барнса для асимптотики.

Где применяется

Преобразование Меллина - рабочий инструмент сразу нескольких областей:

• Аналитическая теория чисел. Дзета-функция Римана $\zeta(s)$ выражается через Меллин-образ, а суммы по простым числам сводятся к контурным интегралам.
• Асимптотика интегралов и сумм. Метод Меллина–Барнса даёт полные асимптотические разложения через полюсы образа.
• Теория вероятностей. Меллин-образ плотности - это аналог производящей функции для произведений независимых случайных величин: образ произведения равен произведению образов.
• Цифровая обработка и распознавание. Лог-полярное преобразование плюс Фурье даёт инвариантность к масштабу - основа Фурье–Меллин-дескрипторов.

Частые ошибки

• Забыть про полосу сходимости. Формула $\mathcal{M}\{e^{-x}\}(s) = \Gamma(s)$ верна как интеграл только при $\Re(s) > 0$. За пределами полосы равенство держится лишь как аналитическое продолжение, а сам интеграл расходится.
• Перепутать показатель ядра. Ядро Меллина - $x^{s-1}$, а не $x^{s}$. Лишняя или потерянная единица сдвигает образ на единицу по $s$ и ломает все формулы.
• Считать обратный интеграл по «неправильной» прямой. Абсцисса $c$ должна лежать внутри фундаментальной полосы; прямая вне полосы даёт неверный результат или расходимость.
• Игнорировать поведение на обоих концах. Левую границу полосы задаёт нуль, правую - бесконечность. Проверять надо обе асимптотики, иначе полоса определена неверно.
• Путать масштаб и сдвиг. Меллин реагирует на растяжение $x \to \lambda x$, а не на сдвиг. Перенос свойств Фурье «как есть» приводит к ошибкам в множителях.

FAQ

Чем преобразование Меллина отличается от Фурье? Фурье разлагает функцию по гармоникам и инвариантно к сдвигу аргумента: сдвиг даёт фазовый множитель. Меллин разлагает по степеням и инвариантно к масштабу: растяжение $x \to \lambda x$ даёт множитель $\lambda^{-s}$. Связаны они логарифмической заменой $x = e^{-u}$, после которой Меллин превращается в Фурье (или Лаплас) функции $f(e^{-u})$.

Почему образ существует только в полосе, а не на всей плоскости? Сходимость интеграла на нуле требует $\Re(s) > \alpha$, а на бесконечности - $\Re(s) < \beta$. Эти два условия вместе вырезают вертикальную полосу $\alpha < \Re(s) < \beta$. Вне неё интеграл расходится, но аналитическое продолжение часто расширяет область определения образа.

Как преобразование Меллина связано с гамма-функцией? Гамма-функция - это в точности Меллин-образ экспоненты: $\Gamma(s) = \int_0^{\infty} x^{s-1} e^{-x}\,dx$. Поэтому гамма-функция играет роль базового «словаря»: через неё и её сдвиги выражаются образы степенных и экспоненциальных функций.

Коротко

Преобразование Меллина $F(s) = \int_0^{\infty} x^{s-1} f(x)\,dx$ переводит функцию с положительной полуоси в функцию комплексного $s$ и устроено вокруг масштабной симметрии. Образ живёт в фундаментальной полосе $\alpha < \Re(s) < \beta$, границы которой задают асимптотики на нуле и бесконечности. Канонический пример - образ экспоненты, равный гамма-функции $\Gamma(s)$ при $\Re(s) > 0$. Логарифмическая замена связывает Меллин с Фурье и Лапласом, а обратная формула - контурный интеграл по вертикальной прямой внутри полосы, который удобно считать через вычеты.