Лайонс Р. Цифровая обработка сигналов. Второе издание. Пер. с англ. (2006) (1095937), страница 73
Текст из файла (страница 73)
Уравнение (9-11) выглядит совсем не так, как уравнение для импульсной характеристики, которое я нашел в другой книге по ЦОС. Которое из них правильно? ь И читатель будет прав, потому что в литературе популярно следующее выражение для 11(п); Ответ состоит в следующем. Во-первых, вывод уравнения (9-12) основан на предположении, что частота дискретизации/«равна единице. Кроме того, если вы стряхнете пыль с вашего старого справочника по математике, то внутри него вы найдете тригонометрическое тождество: яп2(а) = [1 — соз(2а))/2.
Если в (9-11) мы подставим 1 вместо/л и подставим 2яп2(лп/2) вместо ) 1 — сов(лп)), то мы увидим, что (9-11) и (9-12) идентичны. л п(п) Зл 2 ° ' ° ° Ф ° «в -1 ь -15 -10 -5 0 и 10 15 Рис. 9.9. Дискретная импульсная характеристика преобразования Гипьберта при 1, = 1 1, 0.5)- 0, Глава9.Диск етноеп еоб азованиеГильбе га 11(п) = ( 1/лпт«))1 — соз(2л/т п/2~,)) =~, /лп 11 — соз(лп)) .
(9-11) для п э~ О и ~Ь(п) = О, для п=О) . Альтернативная форма: п(п) = 2яп2(лп/2)/(лп) . (9-12) 9.4. П оекти ованиедиск етногоп еоб азователяГильбе та 371 Посмотрев снова на рисунок 9.9, мы можем усилить обоснованность нашего вывода Ь(л). Заметим, что при п > О значения Ь(л) отличны от нуля только для нечетных значений п. Кроме того, зти ненулевые значения уменьшаются как ряд чисел 1/1, 1/3, 1/5, 1/7 и т.
д. Что это вам напоминает? Правильно, ряд Фурье для периодического прямоугольного сигнала! Это не случайно, т. к. Ь(л) представляет собой обратное преобразование Фурье прямоугольной функции Н(вз), показанной на рисунке 9.2. Более того, Ь(п) антисимметрична, и это согласуется с тем, что Н(гз) — чисто мнимая. (Бели бы нам требовалось сделать Ь(п) симметричной, проинвертировав ее значения при п < О, новая последовательность была бы пропорциональна ряду Фурье действительного периодического прямоугольного колебания.) Теперь, когда мы имеем выражение для импульсной характеристики ПГ Ь(п), используем его для построения дискретного преобразователя Гильберта.
9.4. Проектирование дискретного преобразователя Гильберта Дискретное преобразование Гильберта может быть реализовано как во временной, так и в частотной области. Сначала рассмотрим преобразователь Гиль- берта во временной области. 9.4.1. Преобразование Гильберта во временной области: реализация в виде КИХ-фильтра Взгляните снова на рисунок 9.4. Имея Ь(п), мы хотим узнать, как получить дискретный сигнал х;(и). Вспомнив, что в (9-1) используется умножение в частотной области, мы можем сказать, что х;(л) образуется в результате свертки х„(л) с Ь(Ь). Аналитически зто выглядит так: х;(л) =,,~ Ь(Ь)х,(п — Ь).
(9-13) к -ао Это значит, что мы можем реализовать преобразователь Гильберта в виде дискретного нерекурсивного фильтра с конечной импульсной характеристикой (КИХ); как показано на рисунке 9.10. Рис.9.10. Реализация преобразователя Гильберта с К ответвлениями в виде КИХ-фильтра Зтг Глава 9. иск егное и еоб азование Гульбе га Проектирование традиционного КИХ-преобразователя Гильберта сводится к вычислению значений Ь(Ь), благодаря чему становится возможной реализация схемы, приведенной на рисунке 9.4. А, может быть, просто взять коэффициенты Ь(л) согласно (9-11) или рисунку 9.9, и использовать их в качестве Ь(Ь) на рисунке 9.10? Это почти правильный подход.
К сожалению, последовательность Ь(п) на рисунке 9.9 имеет бесконечную длину. И все усилия при проектировании сосредоточиваются на выяснении того, как следует усекать Ь(л). Для начала мы должны решить, должна ли усеченная Ь(л) иметь четную или нечетную длину. Мы принимаем это решение, вспомнив, что реализации КИХ- фильтров, имеющих антисимметричные коэффициенты и нечетную или четную длину, называются системами П1 или 1Ъ' типа соответственно [1 - 3]. Частотные характеристики антисимметричных фильтров этих типов удовлетворяют непреодолимым ограничениям: Длина Ь(я) иечетпв (тпп 111) Длина Ь(я) четпа (тпп 1У) ) Н(0)) - 0 .) Н(0)) - 0 ) Н(ш,/2)) = 0 ! Н(ш, /2)) без ограничений Эта маленькая табличка говорит нам, что амплитудно-частотная характеристика преобразователей Гильберта с нечетным числом ответвлений равна нулю как на частоте 0 Гц, так и на частоте, равной половине частоты дискретизации.
Преобразователи Гильберта с четным количеством ответвлений всегда имеют АЧХ, равную 0 на нулевой частоте. Рассмотрим несколько примеров. На рисунке 9.11 показана частотная характеристика КИХ-преобразователя Гильберта длиной в 15 ответвлений (Тип П1, нечетная длина), коэффициенты которого обозначены как Ь1(Ь). Эти графики могут многому научить нас.
64 ° О 2 4 6 « 6 10 12 14 -16 -122 124 -124 з 1 Ф -2 Ьи -124 Рис. 9. 11. Частотная характеристика Н,(ш) преобразователя Гильберта с 15 ответ- влениями, имеющего импульсную характеристику л 1(к) 9.4. П оекти оеаниедиек етногоп еоб азоеателя Гильбе та 373 1. Например, КИХ-реализация нечетной длины действительно имеет нулевую АЧХ на частотах О Гц и +/,/2 Гц. Это значит, что КИХ-реализации нечетной длины (П1-го типа) ведут себя как полосовые фильтры.' 2. Нт(щ) имеет пульсации в полосе пропускания.
Нам следовало этого ожидать, т. к. мы не можем использовать последовательность Ь ~(Ь) бесконечной длины. Здесь, точно так же, как и при проектировании стандартных КИХ-фильтров нижних частот, усечение импульсной характеристики во временной области приводит к появлению пульсаций в частотной области. (Когда мы резко обрываем функцию в одной области, Природа мстит нам, вводя в действие явление Гиббса, результатом которого являются пульсации в другой области.) Как вы, видимо, уже догадались, мы можем уменьшить пульсации (Н~(м) (, умножив усеченную последовательность Ьт(Ь) на окно. Однако умножение коэффициентов на окно приведет к некоторому сужению полосы пропускания Нт(ш), поэтому может потребоваться увеличение количества коэффициентов.
Попробуйте убедиться в том, что взвешивание усеченной последовательности Ь ~(Ь) окном дает положительный эффект. 3. Очень трудно рассчитать ПГ для низкочастотных сигналов. Мы можем несколько расширить полосу пропускания и уменьшить ширину переходной полосы Н~(о>), но это потребует существенного увеличения длины фильтра. 4. Фазо-частотная характеристика Н~(ш) линейна, как и должно быть при симметричной импульсной характеристике. Наклон ФЧХ (который в нашем случае постоянен) пропорционален задержке, которую испытывает сигнал, проходя через КИХ-фильтр. Ниже мы поговорим об этом подробнее.
Разрыв ФЧХ на частоте О Гц соответствует и радианам, как и должно быть в соответствии с рисунком 9.2. Отлично, это как раз то, чего мы добивались в первую очередь! В нашем неустанном стремлении к точным результатам при генерации аналитического сигнала х,(и) мы должны компенсировать линейный фазовый сдвиг, присутствующий в Йт(тл) — а именно, постоянную задержку во времени, равную групповой задержке фильтра. Мы делаем это, задерживая во времени исходный сигнал х,(и) на время, равное групповой задержке КИХ-преобразователя Гильберта с импульсной характеристикой Ьт(Ь).
Напомним, что групповая задержка 6 КИХ-фильтра длиной в К ответвлений, измеряемая в отсчетах, равна 6 = (К вЂ” 1)/2 отсчетов. Так что блок-схема процесса генерации комплексного сигнала х (и) с помощью КИХ-фильтра приобретает вид, показанный на рисунке 9.12 (а). Мы задерживаем х„(и) на 6 = (7-1) / 2 = 3 отсчета, получая при этом задержанную последовательность х'т(и). Эта задержанная последовательность теперь корректно выровнена по времени с х;(и). Если вы реализуете КИХ-преобразователь Гильберта с нечетным количеством ответвлений в аппаратуре, простейший способ получения х'„(и) состоит в том, чтобы взять копию задержанной исходной последовательности х„(и) со среднего ответвления лийии задержки КИХ-преобразователя Гильберта, как 374 Глава 9.
иск етноеп еоб азованиеГильбе та показано на р)(сунке 9.12 (Ь). Если вы моделируете схему, приведенную на рисунке 9.12 (а), программно, то последовательность х'„(и) можно получить, введя С = 3 нулевых отсчета в начало исходной последовательности х,(п). «',(и) х,(п) = «',(и) +/х,(п) — Ф х,(п) х,(п) (ь) ь(о) Рис.
9.12. Генерация последовательности х (и) с помощью КИХ-фильтра Гильбер- та, имеющего 7 ответвлений: (а) схема обработки; (Ь) структура фильтра Мы можем, например, реализовать КИХ-преобразователь Гильберта, используя структуру!'т/ типа с четным количеством ответвлений. На рисунке 9.13 изображен этот случай, причем коэффициенты обозначены как /)2(к).
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
