Комплексная форма ряда Фурье: вывод и спектр

Вещественный тригонометрический ряд Фурье удобен для интерпретации, но громоздок в выкладках: каждая гармоника требует двух коэффициентов $a_n$ и $b_n$, а суммирование по положительным $n$ скрывает симметрию спектра. Комплексная форма устраняет оба неудобства: один коэффициент $c_n$ на каждое целое $n$, единый интеграл для вывода и прозрачная запись в виде $\sum_{n=-\infty}^{+\infty} c_n e^{inx}$. Ниже - калькулятор спектра, который покажет структуру коэффициентов в действии.

От вещественной формы к комплексной

Вещественный ряд Фурье функции $f(x)$ с периодом $2\pi$:

$f(x) = \frac{a_0}{2} + \sum_{n=1}^{\infty} \bigl(a_n \cos nx + b_n \sin nx\bigr).$

Заменим косинус и синус через формулы Эйлера:

$\cos nx = \frac{e^{inx} + e^{-inx}}{2}, \qquad \sin nx = \frac{e^{inx} - e^{-inx}}{2i}.$

Подставляя в каждый член суммы:

$a_n \cos nx + b_n \sin nx = \frac{a_n - ib_n}{2}\,e^{inx} + \frac{a_n + ib_n}{2}\,e^{-inx}.$

Обозначим $c_n = \dfrac{a_n - ib_n}{2}$ для $n \geq 1$ и $c_{-n} = \dfrac{a_n + ib_n}{2}$. Тогда весь ряд переписывается как двусторонняя сумма по всем целым $n$:

$\boxed{f(x) = \sum_{n=-\infty}^{+\infty} c_n\, e^{inx}}.$

Для $n = 0$ получаем $c_0 = a_0/2$ - постоянная составляющая.

Формула комплексных коэффициентов

Умножим обе части равенства на $e^{-imx}$ и проинтегрируем по периоду:

$\int_{-\pi}^{\pi} f(x)\, e^{-imx}\, dx = \sum_{n=-\infty}^{+\infty} c_n \int_{-\pi}^{\pi} e^{i(n-m)x}\, dx.$

Интеграл от $e^{i(n-m)x}$ по периоду равен $2\pi$ при $n = m$ и нулю иначе (ортогональность системы $\{e^{inx}\}$). Значит:

$\boxed{c_n = \frac{1}{2\pi} \int_{-\pi}^{\pi} f(x)\, e^{-inx}\, dx}, \qquad n \in \mathbb{Z}.$

Эта единая формула заменяет три отдельных интеграла для $a_0$, $a_n$, $b_n$ в вещественной теории.

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

Ортогональность системы $\{e^{inx}\}$

Ключевой шаг вывода формулы для $c_n$ - ортогональность тригонометрической системы в комплексной форме. Для любых целых $m \neq n$:

$\int_{-\pi}^{\pi} e^{inx} \cdot e^{-imx}\, dx = \int_{-\pi}^{\pi} e^{i(n-m)x}\, dx = \left[\frac{e^{i(n-m)x}}{i(n-m)}\right]_{-\pi}^{\pi} = 0,$

потому что $e^{i(n-m)\pi} = e^{-i(n-m)\pi}$ при целом $n - m$. При $n = m$ интеграл равен $2\pi$. Это означает, что функции $e^{inx}$ попарно ортогональны в пространстве $L^2[-\pi,\pi]$ со скалярным произведением $\langle f, g\rangle = \int_{-\pi}^{\pi} f(x)\overline{g(x)}\, dx$, а их норма $\|e^{inx}\| = \sqrt{2\pi}$. Именно из этой ортогональности «вываливается» формула $c_n = \frac{1}{2\pi}\int f(x) e^{-inx} dx$.

Геометрически: функция $f(x)$ раскладывается по ортонормированному базису $\{\frac{1}{\sqrt{2\pi}}e^{inx}\}_{n\in\mathbb{Z}}$, а коэффициенты - это проекции $f$ на каждую базисную функцию. Теорема Парсеваля - теорема Пифагора в бесконечномерном пространстве.

Связь с вещественными коэффициентами

Из определения $c_n$ немедленно следует:

$a_n = c_n + c_{-n}, \qquad b_n = i(c_n - c_{-n}), \qquad n \geq 1.$

Для вещественной $f(x)$ коэффициенты $a_n$, $b_n$ вещественны, а значит $c_{-n} = \overline{c_n}$ (комплексно-сопряжённые). Спектр «в обе стороны» симметричен по модулю: $|c_n| = |c_{-n}|$. Именно поэтому в формуле двустороннего спектра амплитуда каждой частоты «раздваивается»: физическая амплитуда $A_n = \sqrt{a_n^2 + b_n^2}$ связана с комплексным коэффициентом соотношением $|c_n| = A_n / 2$.

Проверим на примере: для $f(x) = \cos(3x)$ имеем $a_3 = 1$, $b_3 = 0$, все остальные коэффициенты нуль. Тогда $c_3 = (1 - 0)/2 = 1/2$, $c_{-3} = 1/2$, остальные $c_n = 0$. Действительно: $\frac{1}{2}e^{3ix} + \frac{1}{2}e^{-3ix} = \cos(3x)$. Вещественная синусоида - это всегда пара комплексных гармоник с равными амплитудами.

Двусторонний амплитудный спектр

График $|c_n|$ как функция $n$ - это двусторонний амплитудный спектр. Для прямоугольного сигнала (меандра) $f(x) = \operatorname{sign}(x)$:

$c_n = \frac{1}{2\pi}\int_{-\pi}^{\pi} \operatorname{sign}(x)\, e^{-inx}\, dx = \frac{2}{in\pi} \text{ (нечётные } n\text{)}, \quad c_n = 0 \text{ (чётные } n\text{)}.$

Амплитуды: $|c_n| = \dfrac{2}{|n|\pi}$ для нечётных $n$, что порождает характерный «гребень» с убывающими «зубцами» на частотах $\pm 1, \pm 3, \pm 5, \ldots$

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

Теорема Парсеваля для комплексного ряда

Равенство Парсеваля в комплексной форме:

$\frac{1}{2\pi}\int_{-\pi}^{\pi} |f(x)|^2\, dx = \sum_{n=-\infty}^{+\infty} |c_n|^2.$

Это интерпретируется как сохранение энергии при переходе из «временного» представления в спектральное. Для меандра левая часть равна 1 (среднеквадратичное значение сигнала $\pm 1$ равно 1). Значит:

$\sum_{\substack{n=-\infty \\ n \text{ нечётн.}}}^{+\infty} \frac{4}{n^2\pi^2} = 1 \implies \sum_{k=0}^{\infty} \frac{1}{(2k+1)^2} = \frac{\pi^2}{8},$

классический результат, получаемый из теории рядов Фурье бесплатно.

В вещественной форме теорема Парсеваля выглядит иначе:

$\frac{1}{2\pi}\int_{-\pi}^{\pi} |f(x)|^2\, dx = \frac{a_0^2}{4} + \frac{1}{2}\sum_{n=1}^{\infty}(a_n^2 + b_n^2).$

Легко убедиться, что обе записи эквивалентны: $\frac{a_n^2 + b_n^2}{2} = 2|c_n|^2$ (из-за множителя $1/2$ в определении $c_n$). Таким образом, при переходе к комплексной форме «вклад» каждой частоты $n$ равномерно распределяется между $+n$ и $-n$.

Фазовый спектр и полная характеристика гармоники

Помимо амплитудного спектра $|c_n|$, для полного описания сигнала нужен фазовый спектр $\arg(c_n)$. Вместе они дают всю информацию о каждой гармонике. Если записать $c_n = |c_n|\, e^{i\varphi_n}$, то $n$-я «пара» гармоник в разложении:

$c_n e^{inx} + c_{-n} e^{-inx} = 2|c_n| \cos(nx + \varphi_n),$

то есть гармоника с амплитудой $A_n = 2|c_n|$, частотой $n$ и начальной фазой $\varphi_n = \arg(c_n)$. Для вещественной чётной функции (только косинусы) все $c_n$ вещественны и неотрицательны: $\varphi_n = 0$. Для нечётной функции (только синусы) $c_n$ чисто мнимые: $\varphi_n = -\pi/2$. Произвольная функция даёт промежуточные фазы.

Обобщение на произвольный период $2L$

Если функция $f(x)$ имеет период $2L$, формула коэффициентов принимает вид:

$c_n = \frac{1}{2L} \int_{-L}^{L} f(x)\, e^{-in\pi x/L}\, dx, \qquad f(x) = \sum_{n=-\infty}^{+\infty} c_n\, e^{in\pi x/L}.$

Шаг по частоте $\Delta\omega = \pi/L$ становится мельче при увеличении периода. В пределе $L \to \infty$ сумма превращается в интеграл - так рождается преобразование Фурье:

$F(\omega) = \int_{-\infty}^{+\infty} f(x)\, e^{-i\omega x}\, dx, \qquad f(x) = \frac{1}{2\pi}\int_{-\infty}^{+\infty} F(\omega)\, e^{i\omega x}\, d\omega.$

Комплексный ряд Фурье, таким образом, - это дискретный прообраз преобразования Фурье.

Практика: вычислить $c_n$ для пилообразной функции

Пусть $f(x) = x$ на $(-\pi, \pi)$ с нечётным продолжением. Тогда $a_n = 0$, а вещественные коэффициенты $b_n = 2(-1)^{n+1}/n$. Переходим к комплексным:

$c_n = \frac{0 - i\cdot 2(-1)^{n+1}/n}{2} = \frac{(-1)^{n+1}}{in} = \frac{i(-1)^n}{n}, \qquad n \neq 0.$

Для $n = 0$: $c_0 = a_0/2 = 0$ (нечётная функция). Проверка: $c_{-n} = \overline{c_n}$ - выполняется, ведь $f(x)$ вещественна.

Можно вычислить $c_n$ непосредственно из интегральной формулы:

$c_n = \frac{1}{2\pi}\int_{-\pi}^{\pi} x\, e^{-inx}\, dx.$

Интегрируя по частям (при $n \neq 0$):

$c_n = \frac{1}{2\pi}\left[\frac{x\, e^{-inx}}{-in}\right]_{-\pi}^{\pi} + \frac{1}{2\pi}\int_{-\pi}^{\pi} \frac{e^{-inx}}{in}\, dx = \frac{1}{2\pi}\cdot\frac{\pi\, e^{-in\pi} - (-\pi)\, e^{in\pi}}{-in}.$

Поскольку $e^{\pm in\pi} = (-1)^n$, числитель равен $\pi(-1)^n + \pi(-1)^n = 2\pi(-1)^n$, откуда $c_n = \frac{i(-1)^n}{n}$ - совпадает с результатом через $b_n$.

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

• Забывают делитель $2\pi$ в формуле $c_n = \frac{1}{2\pi}\int\ldots$ - путаница с нормировкой на $[-\pi,\pi]$ и на $[0, 2\pi]$. Всегда уточняйте, какой промежуток интегрирования.
• Неправильный знак в показателе: в $c_n$ стоит $e^{-inx}$, а в разложении $e^{+inx}$. Перепутать знак - и получить комплексно-сопряжённые коэффициенты.
• Считают $|c_n|$ вдвое больше реальной амплитуды: физическая амплитуда $n$-й гармоники равна $A_n = 2|c_n|$, потому что в сигнал входят $c_n e^{inx}$ и $c_{-n} e^{-inx}$ совместно.
• Применяют комплексный ряд к функции с разрывами без оговорок: в точках разрыва ряд сходится к среднему $\frac{f(x^+) + f(x^-)}{2}$ (условие Дирихле), не к одному из значений функции.
• Путают двусторонний спектр с односторонним: в некоторых задачах по электротехнике спектр рисуют только для $n \geq 0$, а амплитуды на $n > 0$ удваивают. Уточняйте конвенцию.

FAQ

Зачем вводить комплексные коэффициенты, если функция вещественна? Комплексная форма компактнее и удобнее в выкладках: вместо двух формул для $a_n$ и $b_n$ одна для $c_n$. Условие вещественности $f(x)$ автоматически выражается соотношением $c_{-n} = \overline{c_n}$, которое легко контролируется.

Как связан комплексный ряд Фурье с преобразованием Фурье? Комплексный ряд - предельный случай при периоде $2L \to \infty$. Дискретный набор частот $n\pi/L$ переходит в непрерывную переменную $\omega$, а сумма по $n$ - в интеграл. Весь аппарат (формулы, теорема Парсеваля, свёртка) сохраняется в обоих случаях.

Почему в электротехнике пишут $e^{j\omega t}$, а в математике $e^{inx}$? Это одно и то же с разными обозначениями: $j = i$ (мнимая единица), $\omega$ играет роль $n$, а $t$ - роль $x$. В физике и технике мнимую единицу называют $j$, чтобы не путать с током $i$. Формулы идентичны.

Коротко

Комплексная форма ряда Фурье $f(x) = \sum_{n \in \mathbb{Z}} c_n e^{inx}$ - это алгебраически компактное переписывание вещественного ряда через формулу Эйлера. Единый коэффициент $c_n = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x) e^{-inx} dx$ охватывает все гармоники разом; связь $c_n = (a_n - ib_n)/2$ позволяет переходить между двумя формами без потери информации. Двусторонний спектр симметричен для вещественных функций, а теорема Парсеваля гарантирует сохранение «энергии» при спектральном разложении.