Из Википедии, бесплатной энциклопедии
Перейти к навигации Перейти к поиску

В математике , то явление Гиббса, обнаружил Генри Wilbraham  ( 1848 ) [1] и вновь открыта Гиббс  ( 1899 ), [2] является своеобразным способом , в котором ряды Фурье из кусочно непрерывно дифференцируемых периодических функций ведет себя в точке А скачок прерывистости . П - й частичной суммы ряда Фурье имеет большие колебания вблизи скачка, что может увеличить максимум частичной суммы выше, чем в самой функции. Перерегулирование не затухает какn увеличивается, но приближается к конечному пределу. [3] Подобное поведение также наблюдалось физиками-экспериментаторами, но считалось, что оно связано с несовершенством измерительного прибора. [4]

Это одна из причин появления артефактов звонка при обработке сигнала .

Описание [ править ]

Функциональная аппроксимация прямоугольной волны с использованием 5 гармоник
Функциональная аппроксимация прямоугольной волны с использованием 25 гармоник
Функциональная аппроксимация прямоугольной волны с использованием 125 гармоник

Феномен Гиббса включает в себя как тот факт, что сумма Фурье выходит за пределы скачка при разрыве скачка , так и то, что это превышение не исчезает по мере добавления дополнительных членов к сумме.

Три изображения справа демонстрируют явление для прямоугольной волны (высоты ), разложение Фурье которой равно

Точнее, это функция f, которая равна между и, и между и для каждого целого числа n ; таким образом, эта прямоугольная волна имеет скачкообразный скачок высоты во всех целых кратных .

Как видно, по мере увеличения числа членов ошибка аппроксимации уменьшается по ширине и энергии, но сходится к фиксированной высоте. Расчет для прямоугольной волны (см. Зигмунд, глава 8.5. Или вычисления в конце этой статьи) дает явную формулу для предела высоты ошибки. Оказывается, ряд Фурье превышает высоту прямоугольной волны на

( OEISA243268 )

или около 9 процентов прыжка. В более общем смысле, в любой точке скачка кусочно-непрерывно дифференцируемой функции со скачком на a , n- й частичный ряд Фурье (для очень большого n ) будет преодолевать этот скачок приблизительно на одном конце и недооценивать его на такую ​​же величину на другом. конец; таким образом, «скачок» в частичном ряду Фурье будет примерно на 18% больше, чем скачок в исходной функции. В месте самого разрыва частичный ряд Фурье будет сходиться к середине скачка (независимо от того, каково фактическое значение исходной функции в этой точке). Количество

( OEISA036792 )

иногда называют константой Уилбрахама – Гиббса .

История [ править ]

Феномен Гиббса был впервые замечен и проанализирован Генри Уилбрахамом в статье 1848 года. [5] Работа привлекала мало внимания до 1914 года, когда она была упомянута в обзоре математического анализа Генриха Буркхардта в энциклопедии Кляйна . [6] В 1898 году Альберт А. Майкельсон разработал устройство, которое могло вычислять и повторно синтезировать ряды Фурье. [7] [8]Широко распространенный миф гласит, что когда в машину вводились коэффициенты Фурье для прямоугольной волны, график колебался на неоднородностях, и что, поскольку это было физическое устройство, подверженное производственным дефектам, Майкельсон был убежден, что перерегулирование было вызвано ошибками. в машине. На самом деле графики, созданные машиной, были недостаточно хороши, чтобы ясно продемонстрировать феномен Гиббса, и Майкельсон, возможно, не заметил этого, поскольку он не упомянул об этом эффекте в своей статье ( Michelson & Stratton 1898 ) о своей машине или своих более поздних письмах. к природе . [1] Вдохновленный некоторым соответствием в природе между Майкельсоном и Лавом о сходимости ряда Фурье квадратной волновой функции, в 1898 г.Дж. Уиллард Гиббс опубликовал короткую заметку, в которой он рассмотрел то, что сегодня назвали бы пилообразной волной, и указал на важное различие между пределом графиков частичных сумм ряда Фурье и графиком функции, которая является предел этих частичных сумм. В своем первом письме Гиббс не заметил явления Гиббса, и предел, который он описал для графиков частичных сумм, был неточным. В 1899 г. он опубликовал поправку, в которой описал выброс в точке разрыва ( Nature : 27 апреля 1899 г., стр. 606). В 1906 году Максим Бохер дал подробный математический анализ этого выброса, придумав термин «феномен Гиббса» [9].и введение этого термина в широкое употребление. [1]

После того, как о существовании статьи Генри Уилбрахама стало широко известно, в 1925 году Горацио Скотт Карслоу заметил: «Мы все еще можем называть это свойство ряда Фурье (и некоторых других рядов) феноменом Гиббса; но мы больше не должны утверждать, что это свойство было впервые обнаружено пользователя Гиббса ". [10]

Объяснение [ править ]

Неформально, явление Гиббса отражает трудность , присущую аппроксимирующий разрывную функцию с помощью конечной серии непрерывных волн синуса и косинуса. Важно сделать акцент на слове конечный, потому что даже если каждая частичная сумма ряда Фурье выходит за пределы функции, которую она приближает, предел частичных сумм этого не делает. Значение x, при котором достигается максимальное превышение, смещается все ближе и ближе к разрыву по мере увеличения количества суммированных членов, поэтому, опять же неофициально, после того, как превышение прошло на конкретный x , сходимость при этом значении x возможна.

Нет никакого противоречия в том, что перерегулирование сходится к ненулевой величине, но предел частичных сумм не имеет превышения, потому что местоположение этого перерегулирования смещается. У нас есть поточечная сходимость , но не равномерная . Для кусочной функции C 1 ряд Фурье сходится к функции в каждой точке, кроме скачков. На самих разрывах скачка предел будет сходиться к среднему значению функции по обе стороны от скачка. Это следствие теоремы Дирихле . [11]

Явление Гиббса также тесно связано с принципом, согласно которому уменьшение коэффициентов Фурье функции на бесконечности контролируется гладкостью этой функции; очень гладкие функции будут иметь очень быстро убывающие коэффициенты Фурье (что приведет к быстрой сходимости ряда Фурье), тогда как разрывные функции будут иметь очень медленно убывающие коэффициенты Фурье (что приведет к очень медленной сходимости ряда Фурье). Обратите внимание, например, что коэффициенты Фурье 1, -1/3, 1/5, ... прерывистой прямоугольной волны, описанной выше, затухают только с такой же скоростью, как гармонический ряд , который не является абсолютно сходящимся ; действительно, приведенный выше ряд Фурье оказывается лишь условно сходящимся почти для любого значения  x. Это дает частичное объяснение явления Гиббса, поскольку ряды Фурье с абсолютно сходящимися коэффициентами Фурье будут равномерно сходиться с помощью М-критерия Вейерштрасса и, таким образом, не смогут продемонстрировать вышеуказанное колебательное поведение. Точно так же для разрывной функции невозможно иметь абсолютно сходящиеся коэффициенты Фурье, поскольку функция, таким образом, была бы равномерным пределом непрерывных функций и, следовательно, была бы непрерывной, противоречие. Подробнее об абсолютной сходимости рядов Фурье .

Решения [ править ]

На практике трудности, связанные с явлением Гиббса, можно уменьшить, используя более плавный метод суммирования рядов Фурье, такой как суммирование Фейера или суммирование по Риссу , или используя сигма-аппроксимацию . При использовании непрерывного вейвлет- преобразования феномен вейвлет-Гиббса никогда не превосходит феномен Фурье-Гиббса. [12] Кроме того, при использовании дискретного вейвлет-преобразования с базисными функциями Хаара явление Гиббса вообще не возникает в случае непрерывных данных на скачкообразных скачках [13] и минимально в дискретном случае в точках больших изменений. В вейвлет-анализе это обычно называют феноменом Лонго.. В настройке полиномиальной интерполяции явление Гиббса можно смягчить с помощью алгоритма S-Гиббса. [14]

Формальное математическое описание явления [ править ]

Пусть - кусочно непрерывно дифференцируемая функция, периодическая с некоторым периодом . Предположим, что в какой-то момент левый предел и правый предел функции отличаются ненулевым разрывом :

Для каждого натурального числа N ≥ 1 пусть S N  f будет N- м частным рядом Фурье.

где коэффициенты Фурье задаются обычными формулами

Тогда у нас есть

и

но

В более общем смысле, если это любая последовательность действительных чисел, которая сходится к as , и если пробел a положительный, то

и

Если вместо этого зазора является отрицательным, нужно пересадочного предела превосходной с пределом ниже , а также обмена знаками ≤ и ≥, в приведенных выше двух неравенств.

Объяснение обработки сигнала [ править ]

Функция sinc , импульсная характеристика идеального фильтра нижних частот . Масштабирование сужает функцию и, соответственно, увеличивает величину (которая здесь не показана), но не уменьшает величину недостижения, которая является интегралом хвоста.

С точки зрения обработки сигналов явление Гиббса - это переходная характеристика фильтра нижних частот , а колебания называются « звонкими» или « звенящими» артефактами . Усечение преобразования Фурье сигнала на действительной линии или ряда Фурье периодического сигнала (эквивалентно сигналу на круге) соответствует фильтрации более высоких частот идеальным ( кирпичная стена ) низкочастотным / высокочастотным фильтром. вырезать фильтр. Это можно представить как свертку исходного сигнала с импульсной характеристикой фильтра (также известной как ядро ), которая является функцией sinc.. Таким образом, явление Гиббса можно рассматривать как результат свертки ступенчатой ​​функции Хевисайда (если периодичность не требуется) или прямоугольной волны (если периодическая) с функцией sinc: колебания в функции sinc вызывают рябь на выходе.

Синус интеграла , демонстрируя явление Гиббса для ступенчатой функции на вещественной прямой.

В случае свертки со ступенчатой ​​функцией Хевисайда результирующая функция является в точности интегралом функции sinc, интегралом синуса ; для прямоугольной волны описание не так просто. Для ступенчатой ​​функции величина недорега, таким образом, является в точности интегралом от (левого) хвоста, интегрируемого с первым отрицательным нулем: для нормализованного sinc периода единичной выборки это значение, соответственно, имеет такую ​​же величину: интеграл правого хвоста, или, что то же самое, разность между интегралом от отрицательной бесконечности до первого положительного нуля, минус 1 (значение без выхода за пределы).

Перерегулирование и недорегулирование можно понять так: ядра обычно нормализуются, чтобы иметь интеграл 1, поэтому они приводят к отображению постоянных функций в постоянные функции - в противном случае они имеют усиление . Значение свертки в точке представляет собой линейную комбинацию входного сигнала с коэффициентами (весами) значений ядра. Если ядро ​​неотрицательно, например, для ядра Гаусса , тогда значение отфильтрованного сигнала будет выпуклой комбинациейвходных значений (коэффициенты (ядро) интегрируются до 1 и являются неотрицательными) и, таким образом, будут находиться между минимумом и максимумом входного сигнала - он не будет недооценивать или выходить за пределы допустимого диапазона. Если, с другой стороны, ядро ​​принимает отрицательные значения, такие как функция sinc, тогда значение отфильтрованного сигнала вместо этого будет аффинной комбинацией входных значений и может выходить за пределы минимума и максимума входного сигнала. , что приводит к недостижению и перерегулированию, как в феномене Гиббса.

Более длительное расширение - резка на более высокой частоте - соответствует в частотной области расширению кирпичной стены, что во временной области соответствует сужению функции sinc и увеличению ее высоты в тот же раз, при этом интегралы между соответствующими точками остаются неизменными. . Это общая особенность преобразования Фурье: расширение в одной области соответствует сужению и увеличению высоты в другой. Это приводит к тому, что колебания в sinc становятся более узкими и высокими, а в отфильтрованной функции (после свертки) получают более узкие колебания и, следовательно, меньшую площадь, но не уменьшают величину:обрезание на любой конечной частоте приводит к sinc-функции, даже узкой, с теми же хвостовыми интегралами. Это объясняет постоянство перерегулирования и недорегулирования.

  • Колебания можно интерпретировать как свертку с синк.

  • Более высокое значение отсечки делает sinc более узким, но более высоким, с такими же величинами хвостовых интегралов, что приводит к более высокочастотным колебаниям, но величина которых не исчезает.

Таким образом, особенности явления Гиббса интерпретируются следующим образом:

  • недолет происходит из-за импульсной характеристики, имеющей отрицательный хвостовой интеграл, что возможно, поскольку функция принимает отрицательные значения;
  • перерегулирование компенсирует это за счет симметрии (общий интеграл не изменяется при фильтрации);
  • Устойчивость колебаний обусловлена ​​тем, что увеличение отсечки сужает импульсный отклик, но не уменьшает его интеграл - колебания, таким образом, движутся к разрыву, но не уменьшаются по величине.

Пример прямоугольной волны [ править ]

Анимация аддитивного синтеза прямоугольной волны с увеличивающимся числом гармоник. Феномен Гиббса особенно заметен при большом количестве гармоник.

Без ограничения общности мы можем предположить случай прямоугольной волны, в котором период L равен , разрыв находится в нуле, а скачок равен . Для простоты рассмотрим случай, когда N четно (случай нечетного N очень похож). Тогда у нас есть

Подставляя , получаем

как заявлено выше. Далее мы вычисляем

Если ввести нормированную синк функцию , мы можем переписать это как

Но выражение в квадратных скобках - это приближение суммы Римана к интегралу (точнее, приближение по правилу средней точки с интервалом ). Поскольку функция sinc непрерывна, это приближение сходится к фактическому интегралу как . Таким образом, мы имеем

что было заявлено в предыдущем разделе. Аналогичное вычисление показывает

Последствия [ править ]

При обработке сигналов явление Гиббса нежелательно, потому что оно вызывает артефакты, а именно отсечение из-за перерегулирования и недорегулирования, а также артефакты звона от колебаний. В случае фильтрации нижних частот их можно уменьшить или устранить с помощью различных фильтров нижних частот.

В МРТ феномен Гиббса вызывает артефакты в присутствии соседних областей с заметно различающейся интенсивностью сигнала. Это наиболее часто встречается при МРТ позвоночника, где феномен Гиббса может имитировать появление сирингомиелии .

Феномен Гиббса проявляется в виде артефакта перекрестного рисунка в дискретном преобразовании Фурье изображения [15], где у большинства изображений (например, микрофотографий или фотографий) есть резкий разрыв между границами вверху / внизу и слева / справа изображения. Когда периодические граничные условия накладываются в преобразовании Фурье, этот скачкообразный разрыв представлен континуумом частот вдоль осей в обратном пространстве (т. Е. Перекрестная картина интенсивности в преобразовании Фурье).

См. Также [ править ]

  • σ-приближение, которое корректирует суммирование Фурье, чтобы исключить явление Гиббса, которое в противном случае возникло бы на разрывах
  • Феномен Пинского
  • Феномен Рунге (аналогичное явление в полиномиальных приближениях)
  • Интеграл синуса
  • Полосы Маха

Заметки [ править ]

  1. ^ a b c Хьюитт, Эдвин; Хьюитт, Роберт Э. (1979). «Феномен Гиббса-Уилбрахама: эпизод в анализе Фурье». Архив истории точных наук . 21 (2): 129–160. DOI : 10.1007 / BF00330404 . S2CID  119355426 .Доступно в Интернете по адресу: Национальный университет Цзяо Дун: Open Course Ware: Hewitt & Hewitt, 1979. Архивировано 4 марта 2016 г. в Wayback Machine.
  2. ^ Эндрю Dimarogonas (1996). Вибрация для инженеров . ISBN 978-0-13-462938-4.
  3. ^ HS Carslaw (1930). «Глава IX» . Введение в теорию рядов и интегралов Фурье (Третье изд.). Нью-Йорк: Dover Publications Inc.
  4. ^ Vretblad 2000 Раздел 4.7.
  5. ^ Уилбрахам, Генри (1848) «Об одной периодической функции», Кембриджский и Дублинский математический журнал , 3  : 198–201.
  6. ^ Encyklopädie der Mathematischen Wissenschaften mit Einschluss ihrer Anwendungen (PDF) . Том II Т. 1 H 1. Висбаден: Vieweg + Teubner Verlag. 1914. с. 1049 . Проверено 14 сентября 2016 года .
  7. ^ Хаммак, Билл; Кранц, Стив; Карпентер, Брюс (2014-10-29). Гармонический анализатор Альберта Майкельсона: Визуальный тур по машине девятнадцатого века, которая выполняет анализ Фурье . Артикулируйте шумовые книги. ISBN 9780983966173. Проверено 14 сентября 2016 года .
  8. ^ Вольфрам, Стивен (2002). Новый вид науки . Wolfram Media, Inc. стр. 899 . ISBN 978-1-57955-008-0.
  9. ^ Bocher, Maxime (апрель 1906) "Введение в теорию рядов Фурье" , Анналы Mathethematics , второй серии, 7 (3): 81-152. Феномен Гиббса обсуждается на страницах 123–132; Роль Гиббса упоминается на странице 129.
  10. ^ Carslaw, HS (1 октября 1925). «Историческая справка о феномене Гиббса в рядах и интегралах Фурье» . Бюллетень Американского математического общества . 31 (8): 420–424. DOI : 10.1090 / s0002-9904-1925-04081-1 . ISSN 0002-9904 . Проверено 14 сентября 2016 года . 
  11. М. Пинский (2002). Введение в анализ Фурье и всплески . Соединенные Штаты Америки: Брукс / Коул. п. 27 .
  12. ^ Расмуссен, Хенрик О. "Феномен Вейвлета Гиббса". В " Вейвлетах, фракталах и преобразованиях Фурье ", ред. М. Фардж и др ., Clarendon Press, Oxford, 1993.
  13. ^ Келли, Сьюзен Э. «Феномен Гиббса для всплесков». Прикладной и вычислительный гармонический анализ 3, 1995. «Архивная копия» (PDF) . Архивировано из оригинального (PDF) 09.09.2013 . Проверено 31 марта 2012 . CS1 maint: archived copy as title (link)
  14. ^ Де Марчи, Стефано; Маркетти, Франческо; Перраккионе, Эмма; Поджиали, Давиде (2020). «Полиномиальная интерполяция через сопоставленные базы без повторной выборки». J. Comput. Прил. Математика . 364 : 112347. дои : 10.1016 / j.cam.2019.112347 . ISSN 0377-0427 . 
  15. ^ Р. Hovden, Y. Jiang, HL Xin, LF Kourkoutis (2015). "Периодическое уменьшение артефактов в преобразованиях Фурье изображений полного поля с атомным разрешением" . Микроскопия и микроанализ . 21 (2): 436–441. DOI : 10.1017 / S1431927614014639 . PMID 25597865 . CS1 maint: multiple names: authors list (link)

Ссылки [ править ]

  • Гиббса, Дж Willard (1898), "Серия Фурье" , Природа , 59 (1522): 200, DOI : 10.1038 / 059200b0 , ISSN  0028-0836 , S2CID  4004787
  • Гиббса, Дж Willard (1899), "Серия Фурье", Природа , 59 (1539): 606, DOI : 10.1038 / 059606a0 , ISSN  0028-0836 , S2CID  13420929
  • Михельсон, AA; Страттон, SW (1898), «Новый гармонический анализатор», Philosophical Magazine , 5 (45): 85–91
  • Антони Зигмунд , Тригонометрическая серия , Дуврские публикации, 1955.
  • Уилбрахам, Генри (1848), "Об одной периодической функции" , Cambridge and Dublin Mathematical Journal , 3 : 198–201
  • Пол Дж. Нахин , Сказочная формула доктора Эйлера, Princeton University Press, 2006. Ch. 4, разд. 4.
  • Вретблад, Андерс (2000), Анализ Фурье и его приложения , Тексты для выпускников по математике, 223 , Нью-Йорк: Springer Publishing , стр. 93, ISBN 978-0-387-00836-3
  • Вольфрам, Стивен (2002), Новый вид науки , Wolfram Media, Inc., стр. 1152 , ISBN 1-57955-008-8

Внешние ссылки [ править ]

  • СМИ, связанные с феноменом Гиббса, на Викискладе?
  • "Феномен Гиббса" , Математическая энциклопедия , EMS Press , 2001 [1994]
  • Вайсштейн, Эрик В., « Феномен Гиббса ». Материал из MathWorld - веб-ресурса Wolfram.
  • Прандони, Паоло, " Феномен Гиббса ".
  • Радаэлли-Санчес, Рикардо и Ричард Баранюк, « Феномен Гиббса ». Проект Connexions. (Лицензия Creative Commons Attribution)
  • Горацио С. Карслав: Введение в теорию рядов и интегралов Фурье.pdf (Introductiontot00unkngoog.pdf) на archive.org
  • Python реализация алгоритма S-Гиббса смягчающего Явление Гиббса https://github.com/pog87/FakeNodes .