Распределение Накагами

Распределение Накагами — это непрерывное распределение вероятностей, используемое в теории связи и статистике для моделирования замираний сигнала в каналах с многолучевым распространением радиоволн. Оно было предложено японским инженером и математиком Минору Накагами в 1940-х годах и обобщает ряд других распределений, таких как распределение Рэлея и распределение Райса. В отличие от последних, распределение Накагами не требует предположения о равномерном распределении фаз отражённых лучей, что делает его более гибким для описания реальных условий распространения сигнала.

История

Распределение Накагами было впервые опубликовано Минору Накагами в 1940 году в японском научном журнале. В 1958 году работа была переведена на английский язык и получила широкое признание в западной научной литературе. Накагами разработал это распределение для описания замираний сигнала в коротковолновых радиоканалах, где многолучевое распространение приводит к значительным флуктуациям амплитуды принимаемого сигнала. Впоследствии распределение стало стандартным инструментом в теории мобильной связи, спутниковой связи и радиолокации.

Определение

Распределение Накагами описывает амплитуду (огибающую) сигнала. Его плотность вероятности задаётся формулой:

\[ f(r) = \frac{2 m^m r^{2m-1}}{\Gamma(m) \Omega^m} \exp\left(-\frac{m r^2}{\Omega}\right), \quad r \ge 0 \]

где:

• \( r \) — амплитуда сигнала,
• \( m \) — параметр формы (параметр Накагами), \( m \ge 0.5 \),
• \( \Omega = \mathbb{E}[r^2] \) — средняя мощность сигнала,
• \( \Gamma(m) \) — гамма-функция.

Функция распределения (кумулятивная функция) имеет вид:

\[ F(r) = \frac{\gamma\left(m, \frac{m r^2}{\Omega}\right)}{\Gamma(m)} \]

где \( \gamma(\cdot, \cdot) \) — неполная гамма-функция.

Параметры

Параметр формы \( m \)

Параметр \( m \) определяет степень замираний в канале:

• \( m = 0.5 \) — одностороннее нормальное распределение (максимальная степень замираний).
• \( m = 1 \) — распределение Рэлея (стандартный случай для многолучевого распространения без прямой видимости).
• \( m > 1 \) — распределение, приближающееся к распределению Райса; чем больше \( m \), тем слабее замирания и тем ближе распределение к гауссовскому. При \( m \to \infty \) распределение стремится к дельта-функции (отсутствие замираний).

Параметр мощности \( \Omega \)

Параметр \( \Omega \) равен средней мощности сигнала и может быть оценён по выборке как выборочное среднее квадратов амплитуд. В практических приложениях \( \Omega \) часто нормируют к единице.

Свойства

Моменты

Момент \( k \)-го порядка для распределения Накагами вычисляется по формуле:

\[ \mathbb{E}[r^k] = \frac{\Gamma\left(m + \frac{k}{2}\right)}{\Gamma(m)} \left(\frac{\Omega}{m}\right)^{k/2} \]

В частности:

• Среднее значение: \( \mathbb{E}[r] = \frac{\Gamma\left(m + 0.5\right)}{\Gamma(m)} \sqrt{\frac{\Omega}{m}} \)
• Дисперсия: \( \text{Var}[r] = \Omega - \left(\mathbb{E}[r]\right)^2 \)

Медиана

Медиана распределения не выражается в элементарных функциях, но может быть найдена численно из уравнения \( F(r_{0.5}) = 0.5 \).

Мода

Мода (наиболее вероятное значение) равна:

\[ r_{\text{mode}} = \sqrt{\frac{\Omega (2m - 1)}{2m}}, \quad m \ge 0.5 \]

При \( m = 0.5 \) мода находится в нуле.

Связь с другими распределениями

• Распределение Рэлея — частный случай при \( m = 1 \).
• Распределение Райса — аппроксимируется распределением Накагами при \( m > 1 \); точное преобразование: \( m = \frac{(1 + K)^2}{1 + 2K} \), где \( K \) — фактор Райса.
• Гамма-распределение — квадрат амплитуды \( r^2 \) имеет гамма-распределение с параметром формы \( m \) и масштабом \( \Omega/m \).
• Нормальное распределение — при \( m \to \infty \) распределение Накагами стремится к нормальному с дисперсией \( \Omega/(2m) \).

Применение

Теория связи

Распределение Накагами широко используется для моделирования замираний в беспроводных каналах связи, особенно в условиях:

• Городской застройки (многолучевое распространение без прямой видимости).
• Спутниковой связи (сложные условия распространения через ионосферу).
• Подводной акустической связи.

Параметр \( m \) позволяет гибко описывать различные степени замираний — от очень сильных (\( m = 0.5 \)) до слабых (\( m > 10 \)).

Радиолокация

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

Статистика и эконометрика

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

Оценка параметров

Для оценки параметров \( m \) и \( \Omega \) по выборке амплитуд \( r_1, r_2, \dots, r_n \) используется метод моментов или метод максимального правдоподобия.

Метод моментов

Оценки по методу моментов:

\[ \hat{\Omega} = \frac{1}{n} \sum_{i=1}^n r_i^2 \] \[ \hat{m} = \frac{\hat{\Omega}^2}{\frac{1}{n} \sum_{i=1}^n (r_i^2 - \hat{\Omega})^2} \]

Метод максимального правдоподобия

Оценка максимального правдоподобия для \( m \) требует численного решения уравнения:

\[ \ln m - \psi(m) = \ln\left(\frac{1}{n} \sum_{i=1}^n r_i^2\right) - \frac{1}{n} \sum_{i=1}^n \ln r_i^2 \]

где \( \psi(m) \) — дигамма-функция. Оценка для \( \Omega \) остаётся той же, что и в методе моментов.

Критика и ограничения

Основное ограничение распределения Накагами — оно не описывает корреляцию между последовательными отсчётами сигнала (временные замирания). Для моделирования коррелированных замираний используются более сложные модели, такие как процесс Дженкинса или модель Кларка. Кроме того, распределение Накагами не учитывает эффект Доплера, что ограничивает его применение в каналах с подвижными объектами.

Интересные факты

• Распределение Накагами иногда называют «m-распределением» из-за параметра формы \( m \).
• В японской научной литературе распределение часто обозначается как «распределение Накагами — Райса», хотя Райс предложил своё распределение независимо в 1944 году.
• Распределение Накагами является частным случаем обобщённого гамма-распределения.

Источники

• Nakagami M. The m-distribution — A general formula of intensity distribution of rapid fading // Statistical Methods in Radio Wave Propagation. — Pergamon Press, 1960. — С. 3–36.
• Proakis J. G. Digital Communications. — 5th ed. — McGraw-Hill, 2007. — Глава 2.
• Simon M. K., Alouini M.-S. Digital Communication over Fading Channels. — 2nd ed. — Wiley, 2005. — Глава 2.
• Ахметов Р. Г., Крейнделин В. Б. Статистическая теория связи. — М.: Радио и связь, 1985. — Глава 4.