Формула Пуассона и интерполяция

 **Формула суммирования Пуассона** — это фундаментальный закон двойственности математического анализа. Она связывает поведение функции на дискретной сетке с поведением ее непрерывного преобразования Фурье.

Вся теория цифровой обработки сигналов, теоремы отсчетов (интерполяции) и аналитическая теория чисел строятся на одном простом факте: **дискретизация в одном пространстве эквивалентна периодизации в двойственном (частотном) пространстве**.

---

### 1. Формула Пуассона и операция периодизации

Пусть \(f(x)\) — гладкая быстроубывающая функция на вещественной прямой \(\mathbb{R}\), а ее преобразование Фурье определено как:
\[
\widehat{f}(\xi) = \int_{-\infty}^\infty f(x) \, e^{-2\pi i x \xi} \, dx.
\]

#### Классическая формула суммирования Пуассона:
\[
\sum_{n \in \mathbb{Z}} f(n) = \sum_{k \in \mathbb{Z}} \widehat{f}(k).
\]

#### Откуда берется эта формула? (Метод периодизации)
Секрет формулы Пуассона в том, что она представляет собой **разложение периодизованной функции в ряд Фурье**.

1. **Построим периодизацию функции \(f\) с периодом \(T = 1\):**
   Сложим все сдвиги функции на целые числа:
   \[
   F(x) := \sum_{n \in \mathbb{Z}} f(x + n).
   \]
   Функция \(F(x)\) очевидно периодична с периодом \(1\): \(F(x + 1) = F(x)\).

2. **Разложим периодическую функцию \(F(x)\) в ряд Фурье:**
   \[
   F(x) = \sum_{k \in \mathbb{Z}} c_k \, e^{2\pi i k x},
   \]
   где коэффициенты Фурье вычисляются по определению:
   \[
   c_k = \int_0^1 F(x) \, e^{-2\pi i k x} \, dx = \int_0^1 \left( \sum_{n \in \mathbb{Z}} f(x + n) \right) e^{-2\pi i k x} \, dx.
   \]

3. **Трюк разворачивания интеграла (Unfolding trick):**
   Поскольку \(e^{-2\pi i k (x+n)} = e^{-2\pi i k x}\), заменим \(x + n = u\):
   \[
   c_k = \sum_{n \in \mathbb{Z}} \int_n^{n+1} f(u) \, e^{-2\pi i k u} \, du = \int_{-\infty}^\infty f(u) \, e^{-2\pi i k u} \, du = \widehat{f}(k).
   \]
   Коэффициенты ряда Фурье периодической функции \(F(x)\) в точности равны значениям непрерывного преобразования Фурье исходной функции \(\widehat{f}\) в целых точках!

4. **Получение формулы Пуассона:**
   Подставляя \(c_k = \widehat{f}(k)\) обратно в ряд:
   \[
   \sum_{n \in \mathbb{Z}} f(x + n) = \sum_{k \in \mathbb{Z}} \widehat{f}(k) \, e^{2\pi i k x}.
   \]
   Полагая \(x = 0\), мы мгновенно получаем формулу Пуассона:
   \[
   \sum_{n \in \mathbb{Z}} f(n) = \sum_{k \in \mathbb{Z}} \widehat{f}(k).
   \]

---

### 2. Связь с интерполяцией: Теорема Котельникова — Шеннона — Уиттакера

Как формула Пуассона и периодизация порождают формулу восстановления (интерполяции) функции по ее отсчетам?

Рассмотрим масштабированную версию периодизации с шагом \(T > 0\):
\[
\sum_{n \in \mathbb{Z}} f(x + n T) = \frac{1}{T} \sum_{k \in \mathbb{Z}} \widehat{f}\left( \frac{k}{T} \right) e^{2\pi i k x / T}.
\]

Поменяем местами время и частоту: применим формулу периодизации **не к самой функции \(f\), а к ее спектру \(\widehat{f}\)**.

#### Спектральная периодизация:
Периодизуем спектр функции с шагом \(\frac{1}{T}\):
\[
\Phi(\xi) := \sum_{k \in \mathbb{Z}} \widehat{f}\left( \xi + \frac{k}{T} \right) = T \sum_{n \in \mathbb{Z}} f(nT) \, e^{-2\pi i n T \xi}.
\]

Пусть функция \(f(x)\) имеет **финитный спектр** (полоса частот ограничена), то есть:
\[
\operatorname{supp} \widehat{f} \subset [-B, B].
\]

#### Условие Найквиста — Котельникова (отсутствие наложения):
Если мы выберем шаг дискретизации \(T\) достаточно малым, так чтобы период частотной сетки был больше ширины спектра:
\[
\frac{1}{T} \ge 2B \iff T \le \frac{1}{2B},
\]
то сдвинутые копии спектра \(\widehat{f}\left(\xi + \frac{k}{T}\right)\) **не перекрываются**.

В пределах основного частотного диапазона \(\xi \in [-B, B]\) все слагаемые при \(k \neq 0\) равны нулю, и периодизованный спектр совпадает с истинным:
\[
\widehat{f}(\xi) = \mathbf{1}_{[-B, B]}(\xi) \cdot \Phi(\xi) = \mathbf{1}_{[-B, B]}(\xi) \cdot T \sum_{n \in \mathbb{Z}} f(nT) \, e^{-2\pi i n T \xi}.
\]

#### Восстановление функции (Интерполяция):
Применяя обратное преобразование Фурье к обеим частям:
\[
f(x) = \int_{-B}^B \widehat{f}(\xi) \, e^{2\pi i x \xi} \, d\xi = \sum_{n \in \mathbb{Z}} f(nT) \left( T \int_{-B}^B e^{2\pi i (x - nT) \xi} \, d\xi \right).
\]

Интеграл от экспоненты по отрезку \([-B, B]\) при \(T = \frac{1}{2B}\) дает знаменитое кардинальное ядро \(\operatorname{sinc}\):
\[
T \int_{-B}^B e^{2\pi i (x - nT) \xi} \, d\xi = \frac{\sin \pi (x/T - n)}{\pi (x/T - n)} = \operatorname{sinc}\left( \frac{x - nT}{T} \right).
\]

Мы получили классическую формулу кардинальной интерполяции:
\[
f(x) = \sum_{n \in \mathbb{Z}} f(nT) \, \operatorname{sinc}\left( \frac{x - nT}{T} \right).
\]

> **Интерполяционный вывод:** Теорема Котельникова — Шеннона — это прямое следствие того, что коэффициенты ряда Фурье периодизованного спектра равны точечным отсчетам функции во времени \(f(nT)\).

---

### 3. Наложение спектров (Aliasing) на языке периодизации

Что происходит, если шаг интерполяции \(T\) выбран слишком большим (\(T > \frac{1}{2B}\))?

В формуле Пуассона для спектра:
\[
\sum_{k \in \mathbb{Z}} \widehat{f}\left( \xi + \frac{k}{T} \right)
\]
соседние копии \(\widehat{f}\) начинают накладываться друг на друга (явление **элайзинга / aliasing**). Высокие частоты маскируются под низкие, спектр искажается, и восстановить функцию \(f(x)\) через простую \(\operatorname{sinc}\)-интерполяцию становится невозможно без ошибок.

---

### 4. Фундаментальный принцип двойственности

Формула Пуассона утверждает универсальный закон гармонического анализа:

\[
\begin{array}{ccc}
\textbf{Прямое пространство (Время / Пространство)} & & \textbf{Двойственное пространство (Частоты)} \\
\hline
\text{Периодизация с шагом } T & \longleftrightarrow & \text{Дискретизация с шагом } 1/T \\
\text{Дискретизация с шагом } T & \longleftrightarrow & \text{Периодизация с шагом } 1/T \\
\text{Интерполяция по сетке } nT & \longleftrightarrow & \text{Фильтрация (срезка) периода спектра}
\end{array}
\]

---

### 5. Многомерные обобщения: Решетки и упаковка сфер

На многомерном евклидовом пространстве \(\mathbb{R}^d\) формула Пуассона формулируется для произвольной дискретной решетки \(\Lambda \subset \mathbb{R}^d\) и ее **дуальной решетки** \(\Lambda^*\):
\[
\sum_{x \in \Lambda} f(x) = \frac{1}{\operatorname{vol}(\mathbb{R}^d / \Lambda)} \sum_{y \in \Lambda^*} \widehat{f}(y).
\]

Именно эта многомерная связь между периодизацией на решетках, интерполяцией и преобразованием Фурье лежит в основе современных прорывов:
* **Интерполяция Фурье Радченко — Вязовской (2019):** Использование модулярных форм для построения ядер интерполяции, обобщающих \(\operatorname{sinc}\)-функции на случай функций с неограниченным спектром;
* **Плотнейшая упаковка сфер в размерностях 8 и 24:** Решение задачи через формулу Пуассона для решеток \(E_8\) и Лича.