Преобразование Гильберта - Hilbert transform

В математика И в обработка сигналов, то Преобразование Гильберта это конкретный линейный оператор который принимает функцию, ты(т) реальной переменной и производит другую функцию реальной переменной ЧАС(ты)(т). Этот линейный оператор задается формулой свертка с функцией :

то несобственный интеграл понимается в основная стоимость смысл. Преобразование Гильберта имеет особенно простое представление в виде частотная область: Придает сдвиг фазы −90 ° на каждую фурье-компоненту функции. Например, преобразование Гильберта , куда ω> 0, является .

Преобразование Гильберта важно при обработке сигналов, поскольку оно позволяет получить аналитическое представление сигнала с действительным знаком ты(т). В частности, преобразование Гильберта ты это его гармоническое сопряжение v, функция действительной переменной т так что сложный -значная функция ты + я в допускает расширение комплекса верхняя полуплоскость удовлетворение Уравнения Коши – Римана. Преобразование Гильберта было впервые введено Дэвид Гильберт в этой настройке, чтобы решить частный случай Проблема Римана – Гильберта для аналитических функций.

Вступление

Преобразование Гильберта ты можно рассматривать как свертка из ты(т) с функцией час(т) = 1/ π т, известный как Ядро Коши. Потому что1т не является интегрируемый через т = 0 , интеграл, определяющий свертку, не всегда сходится. Вместо этого преобразование Гильберта определяется с использованием Главное значение Коши (здесь обозначается p.v.). Явно преобразование Гильберта функции (или сигнала) ты(т) дан кем-то

при условии, что этот интеграл существует как главное значение. Это как раз свертка ты с умеренное распределение p.v. 1/ π т (из-за Шварц (1950); видеть Панди (1996), Глава 3)). В качестве альтернативы, изменяя переменные, интеграл главного значения может быть записан явно (Зигмунд 1968, §XVI.1) в качестве

Когда преобразование Гильберта применяется дважды подряд к функции ты, результат отрицательный ты:

при условии, что интегралы, определяющие обе итерации, сходятся в подходящем смысле. В частности, обратное преобразование имеет вид ЧАС. Этот факт легче всего увидеть, рассмотрев влияние преобразования Гильберта на преобразование Фурье ты(т) (видеть Связь с преобразованием Фурье, ниже).

Для аналитическая функция в верхняя полуплоскость, преобразование Гильберта описывает связь между действительной и мнимой частью граничных значений. То есть, если ж (z) аналитична в верхней полукомплексной плоскости {z : ℐм z > 0 } , и ты(т) = ℛе ж ( т + 0·я ) , тогда м ж ( т + 0·я ) = ЧАС(ты)(т) с точностью до аддитивной константы, если существует это преобразование Гильберта.

Обозначение

В обработка сигналов преобразование Гильберта ты(т) обычно обозначается как (например., Брэндвуд 2003, п. 87). Однако в математике это обозначение уже широко используется для обозначения преобразования Фурье ты(т) (например., Штайн и Вайс, 1971 г. ). Иногда преобразование Гильберта можно обозначить как . Кроме того, многие источники определяют преобразование Гильберта как отрицательное по сравнению с определенным здесь (например, Bracewell 2000, п. 359).

История

Преобразование Гильберта возникло в 1905 году в работе Гильберта над проблемой, поставленной Риманом относительно аналитических функций (Кресс (1989); Бицадзе (2001)), который стал известен как Проблема Римана – Гильберта. Работа Гильберта в основном была связана с преобразованием Гильберта для функций, определенных на окружности (Хведелидзе 2001; Гильберт 1953 ). Некоторые из его более ранних работ, связанных с дискретным преобразованием Гильберта, восходят к лекциям, которые он читал. Гёттинген. Результаты были позже опубликованы Германом Вейлем в его диссертации (Харди, Литтлвуд и Поля 1952, §9.1). Шур улучшил результаты Гильберта о дискретном преобразовании Гильберта и распространил их на интегральный случай (Харди, Литтлвуд и Поля 1952, §9.2). Эти результаты были ограничены пространствами L2 и 2. В 1928 г. Марсель Рис доказал, что преобразование Гильберта можно определить для ты в L п(ℝ) за 1 < п < ∞ , что преобразование Гильберта является ограниченный оператор на L п(ℝ) за 1 < п < ∞ , и что аналогичные результаты верны для преобразования Гильберта на окружности, а также для дискретного преобразования Гильберта (Рис 1928 ). Преобразование Гильберта было вдохновляющим примером для Антони Зигмунд и Альберто Кальдерон во время изучения сингулярные интегралы (Кальдерон и Зигмунд 1952 ). Их исследования сыграли фундаментальную роль в современном гармоническом анализе. Различные обобщения преобразования Гильберта, такие как билинейные и трилинейные преобразования Гильберта, все еще являются активными областями исследований сегодня.

Связь с преобразованием Фурье

Преобразование Гильберта - это оператор множителя (Duoandikoetxea 2000, Глава 3). Множитель ЧАС является σЧАС(ω) = −я sgn (ω) , куда sgn это сигнум функция. Следовательно:

куда обозначает преобразование Фурье. С sgn (Икс) = sgn (2πИкс) , следует, что этот результат применим к трем общим определениям .

К Формула Эйлера,

Следовательно, ЧАС(ты)(т) имеет эффект сдвига фазы отрицательная частота компоненты ты(т) на + 90 ° (π2 радиан) и фазу положительных частотных составляющих на −90 °. И я·ЧАС(ты)(т) имеет эффект восстановления положительных частотных составляющих при сдвиге отрицательных частотных составляющих еще на + 90 °, что приводит к их отрицанию (то есть умножению на -1).

Когда преобразование Гильберта применяется дважды, фаза отрицательной и положительной частотных составляющих ты(т) сдвинуты соответственно на + 180 ° и -180 °, что эквивалентно. Сигнал отменен; т.е. ЧАС(ЧАС(ты)) = −ты , потому что

Таблица избранных преобразований Гильберта

В следующей таблице частота параметр реально.

Сигнал
Преобразование Гильберта[fn 1]

[fn 2]
[fn 2]


(видеть Функция Доусона )
Функция Sinc
Прямоугольная функция
Дельта-функция Дирака
Характеристика Функция
Примечания
  1. ^ Некоторые авторы (например, Брейсвелл) используют наш -ЧАС как их определение прямого преобразования. Как следствие, правая колонка этой таблицы будет инвертирована.
  2. ^ а б Преобразование Гильберта функций sin и cos можно определить, взяв главное значение интеграла на бесконечности. Это определение согласуется с результатом определения преобразования Гильберта по распределению.

Доступна обширная таблица преобразований Гильберта (Король 2009b Отметим, что преобразование Гильберта константы равно нулю.

Область определения

Ни в коем случае не очевидно, что преобразование Гильберта вообще правильно определено, поскольку определяющий его несобственный интеграл должен сходиться в подходящем смысле. Однако преобразование Гильберта хорошо определено для широкого класса функций, а именно для функций из L п(ℝ) за 1 < п < ∞ .

Точнее, если ты в L п(ℝ) за 1 < п < ∞ , то предел, определяющий несобственный интеграл

существует для почти каждый т . Функция предела также находится в L п(ℝ) и фактически также является пределом в среднем несобственного интеграла. То есть,

в качестве ε → 0 в L п норме, а также поточечно почти всюду Теорема Титчмарша (Титчмарш 1948, Глава 5).

В случае L = 1 преобразование Гильберта по-прежнему сходится поточечно почти всюду, но само может не быть интегрируемым даже локально (Титчмарш 1948, §5.14). В частности, сходимости в среднем в этом случае вообще не происходит. Преобразование Гильберта L1 функция сходится, однако, в L1-слабым, а преобразование Гильберта является ограниченным оператором из L1 к L1, Вт (Штайн и Вайс, 1971 г., Лемма V.2.8). (В частности, поскольку преобразование Гильберта также является оператором множителя на L2, Интерполяция Марцинкевича и аргумент двойственности предоставляет альтернативное доказательство того, что ЧАС ограничен L п.)

Характеристики

Ограниченность

Если 1 < п < ∞ , то преобразование Гильберта на L п(ℝ) это ограниченный линейный оператор, что означает, что существует постоянная Cп такой, что

для всех тыL п(ℝ) . Эта теорема связана с Рисс (1928 г., VII); смотрите также Титчмарш (1948, Теорема 101).

Лучшая константа дан кем-то

Этот результат обусловлен (Пихориды 1972 ); смотрите также Графакос (2004 г., Замечание 4.1.8). Простой способ найти лучшее за степень двойки обусловлена ​​так называемой идентичностью Котлара, которая для всех действительных f. Те же самые лучшие константы верны для периодического преобразования Гильберта.

Из ограниченности преобразования Гильберта следует L п(ℝ) сходимость симметричного оператора частичной суммы

к ж в L п(ℝ)см. например (Duoandikoetxea 2000, п. 59).

Антисамосопряженность

Преобразование Гильберта является анти-самосопряженный оператор относительно двойственности пары между L п(ℝ) и двойное пространство Lq(ℝ), куда п и q находятся Конъюгаты Гёльдера и 1 < п, q < ∞ . Символично,

за тыL п(ℝ) и vLq(ℝ)(Титчмарш 1948, Теорема 102).

Обратное преобразование

Преобразование Гильберта - это антиинволюция (Титчмарш 1948, п. 120), означающий, что

при условии, что каждое преобразование четко определено. С ЧАС сохраняет пространство L п(ℝ), из этого, в частности, следует обратимость преобразования Гильберта на L п(ℝ), и это

Сложная структура

Потому что ЧАС2 = −Ⅰ   (“  " это оператор идентификации ) на реальном Банахово пространство из настоящий-значные функции в L п(ℝ)преобразование Гильберта определяет линейная сложная структура на этом банаховом пространстве. В частности, когда п = 2 ,   преобразование Гильберта дает гильбертово пространство действительных функций в L2(ℝ) структура сложный Гильбертово пространство.

Комплекс) собственные состояния преобразования Гильберта допускают представления в виде голоморфные функции в верхней и нижней полуплоскостях в Харди космос ЧАС2 посредством Теорема Пэли – Винера..

Дифференциация

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

Повторяя эту личность,

Это строго верно, как указано при условии ты и его первый k производные принадлежат L п(ℝ) (Панди 1996, §3.3). Это легко проверить в частотной области, где дифференцирование превращается в умножение на ω.

Свертки

Преобразование Гильберта формально можно реализовать как свертка с умеренное распределение (Duistermaat & Kolk 2010, п. 211)

Таким образом, формально

Тем не мение, априори это может быть определено только для ты распределение компактная опора. С этим можно работать несколько строго, поскольку функции с компактным носителем (которые являются распределениями a fortiori) находятся плотный в L п . В качестве альтернативы можно использовать тот факт, что час(т) это производная по распределению функции бревно остроумие

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

Это строго верно, если ты и v являются распределениями с компактным носителем, поскольку в этом случае

Таким образом, переходя к соответствующему пределу, это также верно, если тыL п и vLq при условии, что

из теоремы в силу Титчмарш (1948, Теорема 104).

Инвариантность

Преобразование Гильберта обладает следующими свойствами инвариантности на L2(ℝ).

  • Он общается с переводчиками. То есть коммутирует с операторами Та ж (Икс) = ж (Икс + а) для всех а в .
  • Он коммутирует с положительными расширениями. То есть коммутирует с операторами Mλ ж (Икс) = ж (λ x) для всех λ > 0 .
  • Это антикоммуты с отражением R f (Икс) = ж (-X) .

С точностью до мультипликативной константы преобразование Гильберта является единственным ограниченным оператором на L2 с этими свойствами (Штейн 1970, §III.1).

На самом деле существует более широкий набор операторов, которые коммутируют с преобразованием Гильберта. Группа SL (2, ℝ) действует унитарными операторами Uграмм на пространстве L2(ℝ) по формуле

Этот унитарное представительство является примером представление основной серии из SL (2, ℝ). В этом случае оно приводимо, разбиваясь как ортогональная сумма двух инвариантных подпространств, Харди космос ЧАС2(ℝ) и его сопряженный. Это пространства L2 граничные значения голоморфных функций на верхней и нижней полуплоскостях. ЧАС2(ℝ) и его сопряжение состоит именно из тех L2 функции с преобразованием Фурье, исчезающим на отрицательной и положительной частях вещественной оси соответственно. Поскольку преобразование Гильберта равно ЧАС = −я (2п − Ⅰ) , с п ортогональная проекция из L2(ℝ) на ЧАС2(ℝ), и то оператор идентификации, следует, что ЧАС2(ℝ) а его ортогональные - собственные подпространства ЧАС для собственных значений ±я. Другими словами, ЧАС коммутирует с операторами Uграмм. Ограничения операторов Uграмм к ЧАС2(ℝ) и его сопряженные дают неприводимые представления SL (2, ℝ) - так называемой предел представлений дискретной серии.[1]

Расширение области определения

Преобразование Гильберта распределений

Кроме того, можно распространить преобразование Гильберта на некоторые пространства распределения (Панди 1996, Глава 3). Поскольку преобразование Гильберта коммутирует с дифференцированием и является ограниченным оператором на L п, ЧАС ограничивает непрерывное преобразование обратный предел из Соболевские пространства:

Тогда преобразование Гильберта можно определить на двойственном пространстве к , обозначенный , состоящий из L п раздачи. Это достигается за счет спаривания дуальности:
За , определять:

Можно определить преобразование Гильберта на пространстве умеренные распределения а также подходом из-за Гельфанд и Шилов (1968), но требуется гораздо больше осторожности из-за особенности в интеграле.

Преобразование Гильберта ограниченных функций

Преобразование Гильберта можно определить для функций из L(ℝ) также, но это требует некоторых изменений и оговорок. При правильном понимании преобразование Гильберта отображает L(ℝ) к Банахово пространство из ограниченное среднее колебание (BMO) классы.

Наивно интерпретируемое преобразование Гильберта ограниченной функции явно неверно определено. Например, с ты = sgn (Икс), интеграл, определяющий ЧАС(ты) расходится почти везде, чтобы ±∞. Чтобы облегчить такие трудности, преобразование Гильберта L функция поэтому определяется следующим упорядоченный форма интеграла

где как указано выше час(Икс) = 1/π Икс и

Модифицированное преобразование ЧАС согласуется с исходным преобразованием функций компактного носителя из общего результата Кальдерон и Зигмунд (1952); видеть Фефферман (1971). Кроме того, полученный интеграл поточечно сходится почти всюду и относительно нормы BMO к функции ограниченного среднего колебания.

А глубокий результат из Фефферман (1971) и Фефферман и Штейн (1972) состоит в том, что функция имеет ограниченное среднее колебание тогда и только тогда, когда она имеет вид ж + ЧАС(грамм) для некоторых ж, граммL(ℝ) .

Сопряженные функции

Преобразование Гильберта можно понять в терминах пары функций ж (Икс) и грамм(Икс) так что функция

граничное значение голоморфная функция F(z) в верхней полуплоскости (Титчмарш 1948, Глава V). В этих условиях, если ж и грамм достаточно интегрируемы, то одно является преобразованием Гильберта другого.

Предположим, что жL п(ℝ) . Тогда по теории Интеграл Пуассона, ж допускает единственное гармоническое продолжение в верхнюю полуплоскость, и это продолжение дается формулой

что является сверткой ж с Ядро Пуассона

Кроме того, существует уникальная гармоническая функция v определенная в верхней полуплоскости такая, что F(z) = ты(z) + я в(z) голоморфен и

Эта гармоническая функция получается из ж взяв свертку с сопряженное ядро ​​Пуассона

Таким образом

Действительно, действительная и мнимая части ядра Коши равны

так что F = ты + я в голоморфен Интегральная формула Коши.

Функция v получен из ты таким образом называется гармоническое сопряжение из ты. Граничный (не касательный) предел v(Икс,у) в качестве у → 0 преобразование Гильберта ж. Таким образом, кратко,

Теорема Титчмарша

Теорема Титчмарша (названная в честь E.C. Титчмарш включивший его в свою работу 1937 г.) уточняет связь между граничными значениями голоморфных функций в верхней полуплоскости и преобразованием Гильберта (Титчмарш 1948, Теорема 95). Он дает необходимые и достаточные условия для комплекснозначного интегрируемый с квадратом функция F(Икс) на действительной прямой как граничное значение функции в Харди космос ЧАС2(U) голоморфных функций в верхней полуплоскости U.

Теорема утверждает, что следующие условия для комплекснозначной квадратично интегрируемой функции F : ℝ → ℂ эквивалентны:

  • F(Икс) это предел как zИкс голоморфной функции F(z) в верхней полуплоскости такая, что
  • Реальная и мнимая части F(Икс) являются преобразованиями Гильберта друг друга.
  • В преобразование Фурье исчезает для Икс < 0 .

Более слабый результат верен для функций класса L п за п > 1 (Титчмарш 1948, Теорема 103). В частности, если F(z) голоморфная функция такая, что

для всех у, то существует комплексная функция F(Икс) в L п(ℝ) такой, что F(Икс + я у) → F(Икс) в L п норма как у → 0 (а также удерживая точечно почти всюду ). Более того,

куда ж является вещественной функцией в L п(ℝ) и грамм - преобразование Гильберта (класса L п) из ж.

Это не так в случае п = 1 . Фактически, преобразование Гильберта L1 функция ж не обязательно сходиться в среднем к другому L1 функция. Тем не менее, (Титчмарш 1948, Теорема 105)преобразование Гильберта ж сходится почти всюду к конечной функции грамм такой, что

Этот результат прямо аналогичен результату Андрей Колмогоров для функций Харди на диске (Дюрен 1970, Теорема 4.2). Хотя этот результат обычно называют теоремой Титчмарша, результат объединяет многие работы других авторов, в том числе Харди, Пэли и Винера (см. Теорема Пэли – Винера. ), а также работы Рисса, Хилле и Тамаркина (см. раздел 4.22 в Король (2009a) ).

Проблема Римана – Гильберта

Одна форма Проблема Римана – Гильберта стремится идентифицировать пары функций F+ и F такой, что F+ является голоморфный в верхней полуплоскости и F голоморфна в нижней полуплоскости, так что при Икс вдоль действительной оси,

куда ж (Икс) некоторая заданная действительная функция от Икс ∈ ℝ . Левую часть этого уравнения можно понимать либо как разность пределов F± из соответствующих полуплоскостей или как гиперфункция распределение. Две функции этого вида являются решением проблемы Римана – Гильберта.

Формально, если F± решить проблему Римана – Гильберта

то преобразование Гильберта ж (Икс) дан кем-то

(Панди 1996, Глава 2).

Преобразование Гильберта на окружности

Для периодической функции ж определяется круговое преобразование Гильберта:

Круговое преобразование Гильберта используется при описании пространства Харди и при изучении сопряженной функции в рядах Фурье. Ядро,

известен как Ядро Гильберта поскольку именно в этой форме преобразование Гильберта было первоначально изучено (Хведелидзе 2001 ).

Ядро Гильберта (для кругового преобразования Гильберта) можно получить, сделав ядро ​​Коши1Икс периодический. Точнее, для Икс ≠ 0

Многие результаты о круговом преобразовании Гильберта могут быть получены из соответствующих результатов для преобразования Гильберта из этого соответствия.

Еще одна более прямая связь обеспечивается преобразованием Кэли C(Икс) = ( Икся ) / ( Икс + я ) , который переносит действительную прямую на окружность и верхнюю полуплоскость на единичный диск. Он индуцирует унитарное отображение

из L2(Т) на L2(ℝ). Оператор U несет пространство Харди ЧАС2(Т) на пространство Харди ЧАС2(ℝ).[2]

Преобразование Гильберта в обработке сигналов

Теорема Бедросяна

Теорема Бедросяна утверждает, что преобразование Гильберта продукта сигнала нижних частот и сигналов верхних частот с неперекрывающимися спектрами задается произведением сигнала нижних частот и преобразования Гильберта сигнала верхних частот, или

куда жLP и жHP - сигналы нижних и верхних частот соответственно (Schreier & Scharf 2010, 14).

Сигналы с амплитудной модуляцией моделируются как продукт ограниченный диапазон форма волны "сообщения", тым(т), и синусоидальный «носитель»:

Когда тым(т) не имеет частотного содержания выше несущей частоты, затем по теореме Бедросяна:

(Бедросян 1962 )

Аналитическое представление

В контексте обработки сигналов интерпретация сопряженной функции преобразования Гильберта, обсужденная выше, дает аналитическое представление сигнала ты(т):

который является голоморфная функция в верхней полуплоскости.

Для узкополосной модели (см. Выше) аналитическое представление:

Формула Эйлера )

 

 

 

 

(Уравнение 1)

Этот комплекс гетеродин операция сдвигает все частотные составляющие тым(т) выше 0 Гц. В этом случае мнимая часть результата является преобразованием Гильберта действительной части. Это косвенный способ получения преобразований Гильберта.

Угловая (фазовая / частотная) модуляция

Форма:

называется угловая модуляция, который включает как фазовая модуляция и модуляция частоты. В мгновенная частота является Для достаточно больших ω, в сравнении с :

и:

Модуляция с одной боковой полосой (SSB)

Когда тым(т) вУравнение 1 является также аналитическое представление (формы волны сообщения), то есть:

результат односторонняя полоса модуляция:

переданный компонент которого:

Причинно-следственная связь

Функция час с час(т) = 1/ π т это беспричинный фильтр и поэтому не может быть реализован как есть, когда ты сигнал, зависящий от времени. Если ты является функцией вневременной переменной (например, пространственной), отсутствие причинности может не быть проблемой. Фильтр тоже бесконечный поддерживать, что может быть проблемой в некоторых приложениях. Другая проблема связана с тем, что происходит с нулевой частотой (DC), чего можно избежать, если s не содержит DC-компонента.

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

Дискретное преобразование Гильберта

Рисунок 1: Фильтр, частотная характеристика которого ограничена примерно до 95% частоты Найквиста.
фигура 2: Фильтр преобразования Гильберта с высокочастотной характеристикой
Рисунок 3.
Рисунок 4. Преобразование Гильберта cos (ωt) является грех (ωt). Этот рисунок показывает грех (ωt) и два приближенных преобразования Гильберта, вычисленные библиотечной функцией MATLAB, hilbert (·)
Рисунок 5. Дискретные преобразования Гильберта косинусной функции с использованием кусочной свертки

Для дискретной функции с преобразование Фурье с дискретным временем (DTFT), и дискретное преобразование Гильберта ДВПФ в регионе π <ω < π дан кем-то:

Обратный ДВПФ с использованием теорема свертки, является:

куда

который представляет собой бесконечную импульсную характеристику (БИХ). Когда свертка выполняется численно, FIR приближение заменяется на час[п], как показано в Рисунок 1. КИХ-фильтр с нечетным числом антисимметричных коэффициентов называется типом III, который по своей природе демонстрирует отклики нулевой величины на частотах 0 и Найквиста, что в данном случае приводит к форме полосового фильтра. Схема типа IV (четное число антисимметричных коэффициентов) показана на фигура 2. Поскольку амплитудная характеристика на частоте Найквиста не пропадает, он немного лучше приближает идеальный преобразователь Гильберта, чем фильтр с нечетным ответвлением. тем не мение

  • Типичный (т.е. правильно отфильтрованный и отобранный) ты[п] последовательность не имеет полезных компонентов на частоте Найквиста.
  • Для импульсной характеристики IV типа требуется12 сдвиг образца в час[п] последовательность. Это приводит к тому, что коэффициенты с нулевым значением становятся ненулевыми, как показано на фигура 2. Таким образом, конструкция типа III потенциально вдвое эффективнее, чем тип IV.
  • Групповая задержка конструкции типа III - это целое число отсчетов, что облегчает выравнивание с создать аналитический сигнал. Групповая задержка типа IV находится на полпути между двумя сэмплами.

В MATLAB функция гильберт (u, N), сворачивает последовательность u [n] с периодическое суммирование:[3]

    [4]

и возвращает один цикл (N выборки) периодического результата в мнимой части комплексной выходной последовательности. Свертка реализуется в частотной области как произведение массива с образцами я sgn (ω) распределение (чьи действительные и мнимые компоненты равны 0 или±1). Рисунок 3 сравнивает полупериод часN[п] с частью эквивалентной длины час[п]. Учитывая приближение FIR для обозначается замена для я sgn (ω) samples дает FIR-версию свертки.

Действительная часть выходной последовательности - это исходная входная последовательность, так что комплексный выход представляет собой аналитическое представление из ты[п]. Когда вход представляет собой сегмент чистого косинуса, результирующая свертка для двух разных значений N изображен в Рисунок 4 (красные и синие участки). Краевые эффекты не позволяют результату быть чистой синусоидальной функцией (зеленый график). С часN[п] не является FIR-последовательностью, теоретическая степень воздействия - это вся выходная последовательность. Но отличия от синусоидальной функции уменьшаются по мере удаления от краев. Параметр N - длина выходной последовательности. Если он превышает длину входной последовательности, входные данные изменяются путем добавления элементов с нулевым значением. В большинстве случаев это уменьшает величину различий. Но их продолжительность определяется временем нарастания и спада час[п] импульсивный ответ.

Понимание краевых эффектов важно, когда метод, называемый перекрытие-сохранение используется для выполнения свертки на длинном ты[п] последовательность. Сегменты длины N свёртываются с периодической функцией:

Когда длительность ненулевых значений является выходная последовательность включает NM + 1 образцы M − 1 выходы отбрасываются из каждого блока N, и входные блоки перекрываются на это количество, чтобы предотвратить пробелы.

Рисунок 5 является примером использования функции Гильберта (·) БИХ и приближения КИХ. В этом примере синусоидальная функция создается путем вычисления дискретного преобразования Гильберта косинусной функции, которая была обработана в четырех перекрывающихся сегментах и ​​снова собрана вместе. Как показывает результат КИХ (синий), искажения, заметные в результате КИХ (красный), не вызваны разницей между час[п] и часN[п] (зеленый и красный в Рисунок 3). Дело в том, что часN[п] конический (оконный) действительно полезен в этом контексте. Настоящая проблема в том, что у него недостаточно окон. Фактически, M = N , тогда как метод сохранения перекрытия требует M < N .

Теоретико-числовое преобразование Гильберта

Теоретико-числовое преобразование Гильберта является расширением (Как 1970 ) дискретного преобразования Гильберта к целым по модулю подходящего простого числа. Здесь следует обобщение дискретное преобразование Фурье к теоретико-числовым преобразованиям. Теоретико-числовое преобразование Гильберта можно использовать для генерации наборов ортогональных дискретных последовательностей (Как 2014 ).

Смотрите также

Рекомендации

  1. ^ Видеть:
  2. ^ Розенблюм и Ровняк 1997, п. 92
  3. ^ видеть Теорема свертки
  4. ^ Для четных значений N, эквивалентная закрытая форма:
    Видеть http://www.rle.mit.edu/dspg/documents/HilbertComplete.pdf экв. (17), (18) и немеченое уравнение ниже (18).

Источники

  • Баргманн, В. (1947). «Неприводимые унитарные представления группы Лоренца». Анна. математики. 48 (3): 568–640. Дои:10.2307/1969129. JSTOR  1969129.
  • Bracewell, R. (2000). Преобразование Фурье и его приложения (3-е изд.). Макгроу – Хилл. ISBN  0-07-116043-4.
  • Карлсон; Крилли и Ратледж (2002). Системы связи (4-е изд.). ISBN  0-07-011127-8.
  • Duoandikoetxea, J. (2000). Фурье-анализ. Американское математическое общество. ISBN  0-8218-2172-5.
  • Дурен, П. (1970). Теория -Пространства. Нью-Йорк, Нью-Йорк: Academic Press.
  • Графакос, Лукас (1994). «Элементарное доказательство суммируемости квадратов дискретного преобразования Гильберта». Американский математический ежемесячный журнал. Математическая ассоциация Америки. 101 (5): 456–458. Дои:10.2307/2974910. JSTOR  2974910.
  • Графакос, Лукас (2004). Классический и современный анализ Фурье. Pearson Education. С. 253–257. ISBN  0-13-035399-X.
  • Кинг, Фредерик В. (2009a). Преобразования Гильберта. 1. Кембридж, Великобритания: Издательство Кембриджского университета.
  • Кинг, Фредерик В. (2009b). Преобразования Гильберта. 2. Кембридж, Великобритания: Издательство Кембриджского университета. п. 453. ISBN  978-0-521-51720-1.
  • Кресс, Райнер (1989). Линейные интегральные уравнения. Нью-Йорк, штат Нью-Йорк: Springer-Verlag. п. 91. ISBN  3-540-50616-0.
  • Ланг, Серж (1985). SL (2, ℝ). Тексты для выпускников по математике. 105. Нью-Йорк, штат Нью-Йорк: Springer-Verlag. ISBN  0-387-96198-4.
  • Панди, Дж. (1996). Преобразование Гильберта распределений Шварца и приложения. Wiley-Interscience. ISBN  0-471-03373-1.
  • Пихоридес, С. (1972). «О наилучшем значении констант в теоремах Рисса, Зигмунда и Колмогорова». Studia Mathematica. 44 (2): 165–179. Дои:10.4064 / см-44-2-165-179.
  • Розенблюм, Марвин; Ровняк, Джеймс (1997). Классы Харди и теория операторов. Дувр. ISBN  0-486-69536-0.
  • Шварц, Лоран (1950). Теория распределений. Париж, Франция: Германн.
  • Schreier, P .; Шарф, Л. (2010). Статистическая обработка комплексных данных: теория несобственных и некруглых сигналов. Кембридж, Великобритания: Издательство Кембриджского университета.
  • Сугиура, Мицуо (1990). Унитарные представления и гармонический анализ: введение. Математическая библиотека Северной Голландии. 44 (2-е изд.). Эльзевир. ISBN  0444885935.
  • Зигмунд, Антони (1988) [1968]. Тригонометрический ряд (2-е изд.). Кембридж, Великобритания: Издательство Кембриджского университета. ISBN  978-0-521-35885-9.

внешняя ссылка