Теорема о том, что при подходящих условиях преобразование Фурье свертки двух сигналов является поточечным произведением их преобразований Фурье
В математика , то теорема свертки заявляет, что при подходящих условиях преобразование Фурье из свертка из двух сигналы это точечный продукт их преобразований Фурье. Другими словами, свертка в одном домене (например, область времени ) равно точечному умножению в другой области (например, частотная область ). Версии теоремы о свертке верны для различных Преобразования, связанные с Фурье . Позволять ж { displaystyle f} и грамм { displaystyle g} быть двумя функции с свертка ж ∗ грамм { displaystyle f * g} . (Обратите внимание, что звездочка означает свертку в этом контексте, а не стандартное умножение. В тензорное произведение символ ⊗ { displaystyle otimes} иногда используется вместо этого.)
Если F { displaystyle { mathcal {F}}} обозначает преобразование Фурье оператор , тогда F { ж } { displaystyle { mathcal {F}} {f }} и F { грамм } { Displaystyle { mathcal {F}} {g }} являются преобразованиями Фурье ж { displaystyle f} и грамм { displaystyle g} , соответственно. потом
F { ж ∗ грамм } = F { ж } ⋅ F { грамм } { displaystyle { mathcal {F}} {f * g } = { mathcal {F}} {f } cdot { mathcal {F}} {g }} [1] куда ⋅ { displaystyle cdot} обозначает точечное умножение. Это также работает наоборот:
F { ж ⋅ грамм } = F { ж } ∗ F { грамм } { Displaystyle { mathcal {F}} {е cdot g } = { mathcal {F}} {f } * { mathcal {F}} {g }} Применяя обратное преобразование Фурье F − 1 { Displaystyle { mathcal {F}} ^ {- 1}} , мы можем написать:
ж ∗ грамм = F − 1 { F { ж } ⋅ F { грамм } } { displaystyle f * g = { mathcal {F}} ^ {- 1} { big {} { mathcal {F}} {f } cdot { mathcal {F}} {g }{большой }}} и:
ж ⋅ грамм = F − 1 { F { ж } ∗ F { грамм } } { displaystyle f cdot g = { mathcal {F}} ^ {- 1} { big {} { mathcal {F}} {f } * { mathcal {F}} {g }{большой }}} Приведенные выше соотношения действительны только для формы преобразования Фурье, показанной в Доказательство раздел ниже. Преобразование может быть нормализовано другими способами, и в этом случае постоянные коэффициенты масштабирования (обычно 2 π { displaystyle 2 pi} или же 2 π { displaystyle { sqrt {2 pi}}} ) появится в отношениях выше.
Эта теорема верна и для Преобразование Лапласа , то двустороннее преобразование Лапласа и, при соответствующей модификации, для Преобразование Меллина и Преобразование Хартли (видеть Теорема обращения Меллина ). Его можно продолжить до преобразования Фурье абстрактный гармонический анализ определяется по локально компактные абелевы группы .
Эта формулировка особенно полезна для реализации числовой свертки на компьютер : Стандартный алгоритм свертки имеет квадратичный вычислительная сложность . С помощью теоремы о свертке и быстрое преобразование Фурье , сложность свертки может быть уменьшена с О ( п 2 ) { Displaystyle О влево (п ^ {2} вправо)} к О ( п бревно п ) { Displaystyle О влево (п журнал п вправо)} , с помощью нотация большой O . Это можно использовать для быстрого построения алгоритмы умножения , как в Алгоритм умножения § Методы преобразования Фурье .
Доказательство
Доказательство здесь показано для конкретного нормализация преобразования Фурье. Как упоминалось выше, если преобразование нормализовано иначе, то постоянная коэффициенты масштабирования появится в выводе.
Позволять ж , грамм { displaystyle f, g} принадлежат к Lп -Космос L 1 ( р п ) { Displaystyle L ^ {1} ( mathbb {R} ^ {n})} . Позволять F { displaystyle F} - преобразование Фурье ж { displaystyle f} и грамм { displaystyle G} - преобразование Фурье грамм { displaystyle g} :
F ( ν ) = F { ж } ( ν ) = ∫ р п ж ( Икс ) е − 2 π я Икс ⋅ ν d Икс , грамм ( ν ) = F { грамм } ( ν ) = ∫ р п грамм ( Икс ) е − 2 π я Икс ⋅ ν d Икс , { displaystyle { begin {align} F ( nu) & = { mathcal {F}} {f } ( nu) = int _ { mathbb {R} ^ {n}} f (x ) e ^ {- 2 pi ix cdot nu} , dx, G ( nu) & = { mathcal {F}} {g } ( nu) = int _ { mathbb {R} ^ {n}} g (x) e ^ {- 2 pi ix cdot nu} , dx, end {выровнено}}} где точка между Икс { displaystyle x} и ν { displaystyle nu} указывает на внутренний продукт из р п { Displaystyle mathbb {R} ^ {п}} . Позволять час { displaystyle h} быть свертка из ж { displaystyle f} и грамм { displaystyle g}
час ( z ) = ∫ р п ж ( Икс ) грамм ( z − Икс ) d Икс . { displaystyle h (z) = int _ { mathbb {R} ^ {n}} f (x) g (z-x) , dx.} Также
∬ | ж ( Икс ) грамм ( z − Икс ) | d z d Икс = ∫ ( | ж ( Икс ) | ∫ | грамм ( z − Икс ) | d z ) d Икс = ∫ | ж ( Икс ) | ‖ грамм ‖ 1 d Икс = ‖ ж ‖ 1 ‖ грамм ‖ 1 . { Displaystyle iint | е (х) г (zx) | , dz , dx = int left (| f (x) | int | g (zx) | , dz right) , dx = int | f (x) | , | g | _ {1} , dx = | f | _ {1} | g | _ {1}.} Следовательно Теорема Фубини у нас есть это час ∈ L 1 ( р п ) { Displaystyle ч в L ^ {1} ( mathbb {R} ^ {n})} так что его преобразование Фурье ЧАС { displaystyle H} определяется интегральной формулой
ЧАС ( ν ) = F { час } = ∫ р п час ( z ) е − 2 π я z ⋅ ν d z = ∫ р п ∫ р п ж ( Икс ) грамм ( z − Икс ) d Икс е − 2 π я z ⋅ ν d z . { Displaystyle { begin {align} H ( nu) = { mathcal {F}} {h } & = int _ { mathbb {R} ^ {n}} h (z) e ^ { -2 pi iz cdot nu} , dz & = int _ { mathbb {R} ^ {n}} int _ { mathbb {R} ^ {n}} f (x) g (zx) , dx , e ^ {- 2 pi iz cdot nu} , dz. end {align}}} Обратите внимание, что | ж ( Икс ) грамм ( z − Икс ) е − 2 π я z ⋅ ν | = | ж ( Икс ) грамм ( z − Икс ) | { Displaystyle | е (х) г (г-х) е ^ {- 2 пи из cdot nu} | = | е (х) г (г-х) |} и, следовательно, с помощью приведенных выше аргументов мы можем снова применить теорему Фубини (т.е. поменять порядок интегрирования):
ЧАС ( ν ) = ∫ р п ж ( Икс ) ( ∫ р п грамм ( z − Икс ) е − 2 π я z ⋅ ν d z ) d Икс . { Displaystyle Н ( Nu) = int _ { mathbb {R} ^ {n}} f (x) left ( int _ { mathbb {R} ^ {n}} g (zx) e ^ {-2 pi iz cdot nu} , dz right) , dx.} Подстановка у = z − Икс { Displaystyle у = г-х} дает d у = d z { displaystyle dy = dz} . Следовательно
ЧАС ( ν ) = ∫ р п ж ( Икс ) ( ∫ р п грамм ( у ) е − 2 π я ( у + Икс ) ⋅ ν d у ) d Икс { Displaystyle Н ( Nu) = int _ { mathbb {R} ^ {n}} f (x) left ( int _ { mathbb {R} ^ {n}} g (y) e ^ {-2 pi i (y + x) cdot nu} , dy right) , dx} = ∫ р п ж ( Икс ) е − 2 π я Икс ⋅ ν ( ∫ р п грамм ( у ) е − 2 π я у ⋅ ν d у ) d Икс { displaystyle = int _ { mathbb {R} ^ {n}} f (x) e ^ {- 2 pi ix cdot nu} left ( int _ { mathbb {R} ^ {n }} g (y) e ^ {- 2 pi iy cdot nu} , dy right) , dx} = ∫ р п ж ( Икс ) е − 2 π я Икс ⋅ ν d Икс ∫ р п грамм ( у ) е − 2 π я у ⋅ ν d у . { Displaystyle = int _ { mathbb {R} ^ {n}} е (х) е ^ {- 2 пи ix cdot nu} , dx int _ { mathbb {R} ^ {п }} g (y) e ^ {- 2 pi iy cdot nu} , dy.} Эти два интеграла являются определениями F ( ν ) { Displaystyle F ( nu)} и грамм ( ν ) { Displaystyle G ( nu)} , так:
ЧАС ( ν ) = F ( ν ) ⋅ грамм ( ν ) , { Displaystyle Н ( Nu) = F ( Nu) CDOT G ( Nu),} QED .
Теорема о свертке для обратного преобразования Фурье
Подобный аргумент, как и приведенное выше доказательство, может быть применен к теореме о свертке для обратного преобразования Фурье;
F − 1 { ж ∗ грамм } = F − 1 { ж } ⋅ F − 1 { грамм } { Displaystyle { mathcal {F}} ^ {- 1} {е * г } = { mathcal {F}} ^ {- 1} {е } cdot { mathcal {F}} ^ {-1} {g }} F − 1 { ж ⋅ грамм } = F − 1 { ж } ∗ F − 1 { грамм } { displaystyle { mathcal {F}} ^ {- 1} {f cdot g } = { mathcal {F}} ^ {- 1} {f } * { mathcal {F}} ^ {-1} {g }} так что
ж ∗ грамм = F { F − 1 { ж } ⋅ F − 1 { грамм } } { displaystyle f * g = { mathcal {F}} { big {} { mathcal {F}} ^ {- 1} {f } cdot { mathcal {F}} ^ {- 1 } {g } { big }}} ж ⋅ грамм = F { F − 1 { ж } ∗ F − 1 { грамм } } { displaystyle f cdot g = { mathcal {F}} { big {} { mathcal {F}} ^ {- 1} {f } * { mathcal {F}} ^ {- 1 } {g } { big }}} Теорема свертки для умеренных распределений
Теорема о свертке распространяется на умеренные распределения . Здесь, грамм { displaystyle g} - произвольное умеренное распределение (например, Гребень Дирака )
F { ж ∗ грамм } = F { ж } ⋅ F { грамм } { displaystyle { mathcal {F}} {f * g } = { mathcal {F}} {f } cdot { mathcal {F}} {g }} F { α ⋅ грамм } = F { α } ∗ F { грамм } { Displaystyle { mathcal {F}} { alpha cdot g } = { mathcal {F}} { alpha } * { mathcal {F}} {g }} но ж = F { α } { displaystyle f = F { alpha }} должен быть "быстро убывающим" в сторону − ∞ { displaystyle - infty} и + ∞ { displaystyle + infty} чтобы гарантировать существование как свертки, так и произведения умножения. α = F − 1 { ж } { Displaystyle альфа = F ^ {- 1} {е }} является гладкой "медленно растущей" обычной функцией, она гарантирует существование как произведения умножения, так и произведения свертки.[2] [3] [4]
В частности, каждое умеренное распределение с компактным носителем, такое как Дельта Дирака , "быстро убывает". ограниченные функции , например, функция, которая постоянно 1 { displaystyle 1} гладкие "медленно растущие" обычные функции. Если, например, грамм ≡ III { Displaystyle g Equiv OperatorName {III}} это Гребень Дирака оба уравнения дают Формула суммирования Пуассона и если, кроме того, ж ≡ δ { Displaystyle е экв дельта} дельта Дирака, тогда α ≡ 1 { Displaystyle альфа эквив 1} всегда равно единице, и эти уравнения дают Идентичность гребня Дирака .
Функции дискретных переменных последовательностей
Аналогичный свертка теорема для дискретных последовательностей Икс { displaystyle x} и у { displaystyle y} является:
D Т F Т { Икс ∗ у } = D Т F Т { Икс } ⋅ D Т F Т { у } , { displaystyle scriptstyle { rm {DTFT}} displaystyle {x * y } = scriptstyle { rm {DTFT}} displaystyle {x } cdot scriptstyle { rm {DTFT} } displaystyle {y },} [5] [а] куда DTFT представляет преобразование Фурье с дискретным временем .
Также есть теорема для круговые и периодические свертки :
Икс N ∗ у ≜ ∑ м = − ∞ ∞ Икс N [ м ] ⋅ у [ п − м ] ≡ ∑ м = 0 N − 1 Икс N [ м ] ⋅ у N [ п − м ] , { Displaystyle х _ {_ {N}} * у треугольник сумма _ {м = - infty} ^ { infty} х _ {_ {N}} [м] CDOT у [нм] эквив сумма _ {m = 0} ^ {N-1} x _ {_ {N}} [m] cdot y _ {_ {N}} [нм],} куда Икс N { displaystyle x _ {_ {N}}} и у N { displaystyle y _ {_ {N}}} находятся периодические суммирования последовательностей Икс { displaystyle x} и у { displaystyle y} :
Икс N [ п ] ≜ ∑ м = − ∞ ∞ Икс [ п − м N ] { Displaystyle х _ {_ {N}} [п] треугольник сумма _ {м = - infty} ^ { infty} х [п-мN]} и у N [ п ] ≜ ∑ м = − ∞ ∞ у [ п − м N ] . { Displaystyle y _ {_ {N}} [n] Triangleq sum _ {m = - infty} ^ { infty} y [n-mN].} Теорема:
D F Т { Икс N ∗ у } = D F Т { Икс N } ⋅ D F Т { у N } , { displaystyle scriptstyle { rm {DFT}} displaystyle {x _ {_ {N}} * y } = scriptstyle { rm {DFT}} displaystyle {x _ {_ {N}} } cdot scriptstyle { rm {DFT}} displaystyle {y _ {_ {N}} },} [6] [b] куда DFT представляет собой N-длину Дискретное преобразование Фурье .
И поэтому:
Икс N ∗ у = D F Т − 1 [ D F Т { Икс N } ⋅ D F Т { у N } ] . { displaystyle x _ {_ {N}} * y = scriptstyle { rm {DFT}} ^ {- 1} displaystyle { big [} scriptstyle { rm {DFT}} displaystyle {x_ {_ {N}} } cdot scriptstyle { rm {DFT}} displaystyle {y _ {_ {N}} } { big]}.} За Икс и у последовательности, ненулевая длительность которых меньше или равна N , окончательное упрощение:
Круговая свертка Икс N ∗ у = D F Т − 1 [ D F Т { Икс } ⋅ D F Т { у } ] { displaystyle x _ {_ {N}} * y = scriptstyle { rm {DFT}} ^ {- 1} displaystyle left [ scriptstyle { rm {DFT}} displaystyle {x } cdot scriptstyle { rm {DFT}} displaystyle {y } right]}
При определенных условиях подпоследовательность Икс N ∗ у { displaystyle x _ {_ {N}} * y} эквивалентна линейной (апериодической) свертке Икс { displaystyle x} и у { displaystyle y} , что обычно является желаемым результатом. (видеть Пример ) И когда преобразования эффективно реализованы с помощью Быстрое преобразование Фурье алгоритм, этот расчет намного эффективнее линейной свертки.
Теорема о свертке для коэффициентов ряда Фурье
Существуют две теоремы свертки для Ряд Фурье коэффициенты периодической функции:
Первая теорема о свертке утверждает, что если ж { displaystyle f} и грамм { displaystyle g} находятся в L 1 ( [ − π , π ] ) { Displaystyle L ^ {1} ([- пи, пи])} , коэффициенты ряда Фурье 2π -периодический свертка из ж { displaystyle f} и грамм { displaystyle g} даны: [ ж ∗ 2 π грамм ^ ] ( п ) = 2 π ⋅ ж ^ ( п ) ⋅ грамм ^ ( п ) , { displaystyle [{ widehat {f * _ {2 pi} g}}] (n) = 2 pi cdot { widehat {f}} (n) cdot { widehat {g}} (n ),} [A] куда: [ ж ∗ 2 π грамм ] ( Икс ) ≜ ∫ − π π ж ( ты ) ⋅ грамм [ pv ( Икс − ты ) ] d ты , ( и pv ( Икс ) ≜ аргумент ( е я Икс ) ⏟ основная стоимость ) = ∫ − π π ж ( ты ) ⋅ грамм ( Икс − ты ) d ты , когда грамм ( Икс ) 2 π -периодический. = ∫ 2 π ж ( ты ) ⋅ грамм ( Икс − ты ) d ты , когда обе функции равны 2 π -периодический, а интеграл по любым 2 π интервал. { Displaystyle { begin {align} left [е * _ {2 pi} g right] (x) & треугольникq int _ {- pi} ^ { pi} f (u) cdot g [{ text {pv}} (xu)] , du, && { big (} { text {and}} underbrace {{ text {pv}} (x) треугольникq arg left (e ^ {ix} right)} _ { text {главное значение}} { big)} & = int _ {- pi} ^ { pi} f (u) cdot g (xu ) , du, && scriptstyle { text {when}} g (x) { text {is 2}} pi { text {-periodic.}} & = int _ {2 pi} f (u) cdot g (xu) , du, && scriptstyle { text {когда обе функции 2}} pi { text {-периодические, а интеграл по любому 2}} pi { текст {интервал.}} end {выровненный}}} Вторая теорема о свертке утверждает, что коэффициенты ряда Фурье произведения ж { displaystyle f} и грамм { displaystyle g} даны дискретная свертка из ж ^ { displaystyle { hat {f}}} и грамм ^ { displaystyle { hat {g}}} последовательности: [ ж ⋅ грамм ^ ] ( п ) = [ ж ^ ∗ грамм ^ ] ( п ) . { displaystyle left [{ widehat {f cdot g}} right] (n) = left [{ widehat {f}} * { widehat {g}} right] (n).} Смотрите также
Примечания
^ Масштабный коэффициент всегда равен периоду, 2π в этом случае. Цитирование страниц
Рекомендации
^ McGillem, Clare D .; Купер, Джордж Р. (1984). Анализ непрерывных и дискретных сигналов и систем (2-е изд.). Холт, Райнхарт и Уинстон. п. 118 (3-102). ISBN 0-03-061703-0 . ^ Хорват, Джон (1966). Топологические векторные пространства и распределения . Ридинг, Массачусетс: издательство Addison-Wesley Publishing Company.^ Баррос-Нето, Хосе (1973). Введение в теорию распределений . Нью-Йорк, Нью-Йорк: Деккер. ^ Петерсен, Бент Э. (1983). Введение в преобразование Фурье и псевдодифференциальные операторы . Бостон, Массачусетс: издательство Pitman Publishing. ^ Proakis, John G .; Манолакис, Дмитрий Г. (1996), Цифровая обработка сигналов: принципы, алгоритмы и приложения (3-е изд.), Нью-Джерси: Prentice-Hall International, стр. 297, Bibcode :1996dspp.book ..... P , ISBN 9780133942897 , sAcfAQAAIAAJ ^ Рабинер, Лоуренс Р. ; Золото, Бернард (1975). Теория и применение цифровой обработки сигналов . Энглвуд Клиффс, Нью-Джерси: Prentice-Hall, Inc., стр. 59 (2,163). ISBN 978-0139141010 .Оппенгейм, Алан В. ; Шафер, Рональд В. ; Бак, Джон Р. (1999). Обработка сигналов в дискретном времени (2-е изд.). Верхняя Сэдл-Ривер, штат Нью-Джерси: Prentice Hall. ISBN 0-13-754920-2 . Также доступно на https://d1.amobbs.com/bbs_upload782111/files_24/ourdev_523225.pdf дальнейшее чтение
Кацнельсон, Ицхак (1976), Введение в гармонический анализ , Дувр, ISBN 0-486-63331-4 Ли, Бинг; Бабу, Дж. Йогеш (2019), "Теорема о свертке и асимптотическая эффективность", Аспирантура по статистическому выводу , Нью-Йорк: Springer, стр. 295–327, ISBN 978-1-4939-9759-6 Вайсштейн, Эрик В. .html "Теорема о свертке" . MathWorld .Кратчфилд, Стив (9 октября 2010 г.), "Радость свертки" , Университет Джона Хопкинса , получено 19 ноября, 2010 Дополнительные ресурсы
Для наглядного представления использования теоремы о свертке в обработка сигналов , видеть: