Явление Гиббса: выброс ряда Фурье у разрыва

Явление Гиббса (или эффект Гиббса) - это устойчивый «выброс» частичной суммы ряда Фурье вблизи точки разрыва функции: как бы много гармоник мы ни взяли, около скачка сумма всегда перелетает истинное значение примерно на 9 % величины скачка. Самое неожиданное здесь то, что выброс не уменьшается с ростом числа слагаемых - он лишь сужается, прижимаясь к разрыву, но держит свою высоту. Ниже разберём, откуда берётся этот выброс при разложении прямоугольного сигнала, чему он равен, почему не исчезает и как его величину получить точно через интегральный синус. А пока покрути интерактив - добавляй гармоники сам и следи, как кривая прижимается к прямоугольному сигналу, но «ушко» у разрыва сужается, не опускаясь.

Что такое явление Гиббса

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

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

Эти осцилляции и есть проявление явления Гиббса. Ключевая его особенность: самый высокий «горб» у разрыва не становится ниже, когда мы добавляем гармоники. Он становится только уже и придвигается к точке скачка, но его высота стремится к фиксированной величине, а не к нулю. Именно поэтому говорят, что ряд Фурье сходится к разрывной функции поточечно, но не равномерно в окрестности разрыва.

Откуда берётся выброс: ряд Фурье прямоугольного сигнала

Возьмём канонический пример - нечётный прямоугольный сигнал (меандр), равный $+1$ при $0 < x < \pi$ и $-1$ при $-\pi < x < 0$. В точке $x = 0$ у него скачок величиной $D = 2$ (от $-1$ к $+1$). Его ряд Фурье содержит только нечётные синусы:

$f(x) = \frac{4}{\pi}\sum_{k=0}^{\infty}\frac{\sin\big((2k+1)x\big)}{2k+1} = \frac{4}{\pi}\left(\sin x + \frac{\sin 3x}{3} + \frac{\sin 5x}{5} + \dots\right).$

Частичная сумма $S_n$ - это та же формула, но оборванная на $n$ первых гармониках. Каждая синусоида по отдельности гладкая и ограниченная, но их сумма пытается воспроизвести вертикальную ступеньку. Чем больше слагаемых, тем круче «стенка» суммы у нуля - и тем заметнее, что перед самой стенкой кривая выпрыгивает выше уровня $+1$. Этот первый горб и есть выброс Гиббса.

Удобно представлять каждую гармонику как вектор, вращающийся со своей частотой: первая делает один оборот за период, третья - три, пятая - пять и так далее. Если выстроить эти векторы в цепочку, их общий конец вычерчивает частичную сумму - именно так гармоники «складываются» в прямоугольный сигнал.

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

Чему равен выброс Гиббса

Посчитаем высоту горба в пределе. Максимум частичной суммы $S_n$ около разрыва при $n \to \infty$ стремится не к $1$, а к большему числу:

$\lim_{n\to\infty}\max_x S_n(x) = \frac{2}{\pi}\int_0^{\pi}\frac{\sin t}{t}\,dt = \frac{2}{\pi}\,\mathrm{Si}(\pi) \approx 1{,}17898,$

где $\mathrm{Si}(\pi) \approx 1{,}8519$ - значение интегрального синуса в точке $\pi$. Поскольку истинный уровень сигнала справа от разрыва равен $1$, сумма перелетает его на

$1{,}17898 - 1 = 0{,}17898$

по абсолютной величине. Удобнее измерять выброс в долях полного скачка $D = 2$: тогда перерегулирование составляет

$\frac{0{,}17898}{2} \approx 0{,}0895 = 8{,}95\,\%.$

Отсюда и берётся знаменитое «около 9 %». Эта доля универсальна: она не зависит ни от величины скачка, ни от конкретной функции - для любого разрыва первого рода частичная сумма перелетает его примерно на 8,95 % высоты скачка с каждой стороны.

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

Почему выброс не исчезает с ростом числа гармоник

Интуиция подсказывает: добавим больше гармоник - приближение станет лучше, и выброс уйдёт. Но это не так, и понять почему - самое важное в теме. С ростом $n$ происходят две вещи одновременно. Во-первых, положение первого горба сдвигается к разрыву: пик находится примерно при $x \approx \pi/(2n)$, то есть подбирается к нулю всё ближе. Во-вторых, высота этого горба меняется очень мало и стремится к фиксированному уровню $\tfrac{2}{\pi}\mathrm{Si}(\pi)$.

Получается, что «ушко» не сглаживается, а сжимается по горизонтали, оставаясь почти неизменным по вертикали. В любой фиксированной точке $x \ne 0$ сумма в итоге сходится к правильному значению (поэтому сходимость поточечная), но точка максимального выброса всё время убегает к разрыву, унося с собой свои 9 %. Площадь под горбом стремится к нулю, и в смысле среднеквадратичного приближения ряд сходится - энергия ошибки убывает, что согласуется с неравенством Бесселя для коэффициентов Фурье. Но поточечной близости у самого скачка добиться нельзя: равномерной сходимости здесь нет.

Как вычислить предел через интегральный синус

Чтобы получить значение $\tfrac{2}{\pi}\mathrm{Si}(\pi)$ строго, частичную сумму удобно записать через ядро Дирихле и перейти к интегралу. Максимум $S_n$ достигается вблизи первого нуля производной, и при больших $n$ сумма у горба ведёт себя как интегральная сумма для функции $\tfrac{\sin t}{t}$:

$\max_x S_n(x) \;\longrightarrow\; \frac{2}{\pi}\int_0^{\pi}\frac{\sin t}{t}\,dt.$

Интеграл $\int_0^{\pi}\frac{\sin t}{t}\,dt$ не берётся в элементарных функциях - это и есть определение интегрального синуса $\mathrm{Si}(\pi)$. Численно $\mathrm{Si}(\pi) = 1{,}851937\ldots$, откуда и получается множитель $1{,}17898$. Подробный вывод через ядро Дирихле опирается на оценку частичной суммы как функции, сходимость которой в среднем гарантирует тот же аппарат ортогональных разложений, что стоит за признаком Дирихле для рядов.

Где встречается на практике и как с ним борются

Явление Гиббса - это не только теоретический курьёз. Оно проявляется везде, где разрывный сигнал восстанавливают из конечного числа частот:

• Обработка изображений. При сжатии (например, JPEG) резкие границы объектов раскладываются по косинусам, и около контрастных краёв появляется «звон» (ringing) - светлые и тёмные каёмки. Это прямой аналог выброса Гиббса.
• Цифровая фильтрация звука. Идеальный фильтр с резким срезом частот даёт переходную характеристику с осцилляциями у фронта - те же горбы.
• Численное решение уравнений. Спектральные методы для задач с разрывами дают паразитные колебания у скачков.

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

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

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

• Считать, что выброс уходит при $n \to \infty$. Не уходит. К нулю стремится ширина горба, а не его высота: перерегулирование остаётся около 9 % скачка при любом конечном $n$.
• Путать поточечную и равномерную сходимость. Ряд Фурье меандра сходится к нему в каждой точке непрерывности, но не равномерно около разрыва - именно из-за выброса.
• Брать долю 9 % от уровня сигнала, а не от скачка. Перерегулирование считается от величины скачка $D$. Для меандра $D = 2$, поэтому абсолютный перелёт $\approx 0{,}179$, а не $0{,}09$.
• Ожидать выброс там, где функция непрерывна. Явление Гиббса возникает только у разрывов первого рода. Для гладкой функции сходимость равномерная, и горбов нет.
• Считать, что в самой точке разрыва сумма «выбрасывается». В точке $x = 0$ ряд сходится к среднему односторонних пределов, то есть к $0$. Выброс наблюдается не в самой точке, а рядом с ней.

FAQ

Чему равен выброс Гиббса в процентах? Около 8,95 % величины скачка с каждой стороны разрыва. Точное значение даётся пределом $\tfrac{2}{\pi}\mathrm{Si}(\pi) \approx 1{,}179$ для максимума суммы при единичном уровне сигнала, что соответствует перелёту на $\approx 0{,}179$ при скачке $D = 2$.

Исчезает ли явление Гиббса при бесконечном числе гармоник? Нет. Высота выброса стремится к фиксированной величине, а к нулю стремится только его ширина. Поэтому при любом конечном числе гармоник около разрыва остаётся горб высотой примерно 9 % скачка.

Как убрать выброс Гиббса? Сменить способ суммирования ряда: суммирование по Фейеру (средние Чезаро), оконные функции или сигма-аппроксимация Ланцоша сглаживают осцилляции и дают равномерное приближение ценой небольшого размытия скачка.

Почему выброс называется именно явлением Гиббса? По имени физика Джозайи Уилларда Гиббса, который в 1899 году объяснил эти осцилляции; ранее их наблюдал Генри Уилбрэхем, поэтому эффект иногда называют явлением Гиббса - Уилбрэхема.

Коротко

Явление Гиббса - устойчивый выброс частичной суммы ряда Фурье у разрыва первого рода: сумма перелетает скачок примерно на 8,95 % его величины. Этот выброс не исчезает с ростом числа гармоник - он лишь сужается к разрыву, сохраняя высоту, потому что ряд сходится поточечно, но не равномерно. Точное значение перелёта даёт предел $\tfrac{2}{\pi}\mathrm{Si}(\pi) \approx 1{,}179$ через интегральный синус. На практике тот же эффект виден как «звон» у краёв на сжатых изображениях и у фронтов фильтрованных сигналов, а подавляют его мягким суммированием ряда - по Фейеру или через оконные функции.