Лайонс Р. Цифровая обработка сигналов. Второе издание. Пер. с англ. (2006) (1095937), страница 75
Текст из файла (страница 75)
Затем мы используем действительную и мнимую части коэффициентов фильтра Ьвр (Ь) в двух отдельных КИХ-фильтрах с действительными коэффициентами, как показано на рисунке 9.15. На языке ЦОС фильтр, выдающий последовательность х)(п), называется 1-каналом для синфазнойсоставляющей, а фильтр, выдающий х~(л), называется Я-каналом для квадратурной составляющей. Эта схема генерации аналитического сигнала с использованием смешивания, приведенная на рисунке 9.15, обладает рядом интересных особенностей: 1. Смешивающее преобразование коэффициентов Ьгл (Ь) путем умножения на е)о'о"г~ приводит к уменьшению амплитудно-частотных характеристик двух действительный КИХ-фильтров в 2 раза.
Удвоение значений Ь~р (Ь) перед смешиванием устранит эти потери. 2. Мы можем взвесить коэффициенты Ьгр(Ь) окном перед смешиванием для уменьшения величины пульсаций в полосе пропускания и разности АЧХ двух действительных фильтров. (Нам хотелось бы сохранить коэффициенты передачи в полосе пропускания обоих фильтров как можно ближе друг к другу.) Конечно, взвешивание несколько ухудшает крутизну АЧХ в переходной полосе, поэтому может потребоваться больше коэффициентов.
3. В этом случае можно использовать фильтры как нечетной, так и четной длины. 4. Если либо Ь,,(Ь), либо Ь„„(Ь) симметричны, можно использовать сложеннув КИХ-структуру (см. раздел 13.7) для уменьшения в два раза количества умножений на один выходной отсчет. 5. Если центральные частоты исходного полосового сигнала и комплексного полосового фильтра равны четверти частоты дискретизации (~,/4), вам крупно повезло, потому что примерно половина коэффициентов каждого из действительных фильтров в этом случае равна нулю, что снижает количество умножений еще в два раза.
9.6. С аанениеметодоагене циианалитическихсигналоа 379 6. Особенно эффективен полуполосный КИХ комплексный полосовой фильтр, центральная частота которого перенесена на~,=/', /4. В этом случае половина коэффициентов Ьгох(1) принимают нулевые значения и все, кроме одного, коэффициенты Ь,.;„(а) равны нулю! 7. При аппаратурной реализации должны быть реализованы оба действительных КИХ-фильтра, изображенные на рисунке 9.15.
Если же вы реализуете генерацию аналитического сигнала просто для тренировки на каком-нибудь языке высокого уровня, например, в рамках МАТЮКАВ, используемый язык может дать возможность реализации бар(Ф) в виде одного комплексного фильтра. Запомните, что последовательность х0(п) на рисунке 9.15 не является преобразованием Гильберта входной последовательности х„(п).
Это не бьио нашей целью. Мы намеревались получить аналитический сигнал хг(п), квадратурная составляющая которого х~(п) представляет собой преобразование Гильберта синфазной последовательности х)(п). 9.6. Сравнение методов генерации аналитических сигналов Проектирование преобразователя Гильберта во временной области представляет собой, по существу, упражнение в проектировании фильтров нижних частот.
И как таковой, идеальный дискретный КИХ-преобразователь Гильберта, как и идеальный КИХ-фильтр нижних частот, не может быть реализован практически. К счастью, мы обычно можем обеспечить полосу пропускания преобразователя Гильберта, перекрывающую полосу частот исходного сигнала, который мы сдвигаем по фазе на + 90'. Использование большего количества ответвлений улучшает характеристики преобразователя (минимизируя пульсации в полосе пропускания), а выбор нечетного или четного количества ответвлений определяется тем, хотим ли мы, чтобы коэффициент передачи на частоте Найквиста был равен О, или нет, а также тем, должна ли задержка равняться целому количеству периодов дискретизации.
Эти замечания справедливы также и в случае генерации комплексного сигнала по входному действительному сигналу методом комплексной фильтрации (рисунок 9.15). Один возможный недостаток этого метода состоит в том, что действительная часть получаемого комплексного сигнала не совпадает с исходным входным вещественным сигналом (т. е.'хоп) не равна х,(п).) При использовании вычислений с плавающей точкой схема комплексной фильтрации обеспечивает разность между двумя вещественными фильтрами длиной в 100 отсчетов каждый не более 0.15 дБ, тогда как КИХ-фильтр Гильберта требует примерно 30 ответвлений для обеспечения такой же величины пульсаций в полосе пропускания.
Эти числа имеют место, когда не используется взвешивание с помощью окон. Таким образом, КИХ-фильтр Гильберта более эффективен с вычислительной точки зрения, чем схема комплексной фильтрации. Для высококачественной генерации аналитических сигналов (когда ширина полосы пропускания преобразователя Гильберта должна быть очень близка к полной полосе 380 Глава 9. Диск етноеп еоб взоввниеГильбе гн сигналов 1; /2, и пульсации в полосе пропускания должны быть очень низкими) следует рассмотреть возможность применения метода ДПФ. Схема комплексной фильтрации вносит потери по амплитуде в два раза, а методы ДПФ и КИХ-фильтрации не дают таких потерь.
Метод ДПФ обеспечивает самую точную генерацию комплексного аналитического сигнала по вещественному сигналу, и при использовании ДПФ действительная часть х,(п) совпадает с вещественным входным сигналом х,(п). Выбор используемого метода в каждом конкретном приложении требует моделирования на целевой аппаратуре с использованием ожидаемой формы представления чисел (целых или с плавающей точкой) при типовых входных сигналах, которые вы собираетесь обрабатывать. Библиография 1.
РгоаЫз, 1., апс1 Мапо!а!с1з, 1). П8(га! 5!8па! Ртосезяп8с Рт!пар!ез, А!нот!Йспз аль Арр!!саг!опе, Ргепс1се-На!1, 1)ррег 5айИе К1чег, Ь(етт )егзеу, 1996, рр. 618, 657. 2. КаЬ(пег, 1., апс1 Со!Й, В. 77се Т7сеотуапс7Аррйсабоп о/!7!8!го!5!8па7Ртосея!пя, Ргепс1се-На!1, Еп81еъ оос1 С!1!!з, Мечт )егзеу, 1975, рр. 67, 168 (есть русский перевод: Рабинер Л., Голд В. «Теория и применение цифровой обработки сигналов», Мс Мир, 1978, доступен по адресу Ыр://8ео81п.пагос).ги/агЬ1ч/сЬр/сЬрЗ.Ьст). 3.
ОррепЬепп, А., апс1 5сЬа!ег, К. Э!зстеге-Т!те 5!8па! РтосеЫпя, Ргепс(се-На!1; Еп8!есчоос! С!11!з, Ыечс )егзеу, 1зг Ес). 1989, 2пс( Ес(. 1999 (имеется русский перевод одного из предыдущих изданий: Оппенгейм А. В., Шафер Р. В. «Цифровая обработка сигналов», пер. с англ. / под ред. С. Я. Шаца, М.: Связь, 1979, доступен по адресу с)зр-Ьоо)с.пагос).гп/Ор5ЬП5РЛ1!чи).
4. ааааа)со, Т., ес. а1. «5р!1сс(п8 "сЬе 11пй Пе1ау»,1ЕЕЕ5!8па!РтосеззтйМайаз(пе ! ап. 1996. 8. Сотез, 1., апс1 Регга81(а, А. «Ап Апа1о8 5атр1ес1-1)аса П5В со 55В Сопчегсег Пяпй Кесигяче Н11Ъегс Тгапз!огщег 1ог Ассигаге 1 апс1 (4 СЬаппе1 МассЬ)п8», 1ЕЕЕ Ттапх оп Соси!гз апй 5уя-11, Апа!о8 апй Р(8(га! 5!8па! Ртосезт!п8, Чо1. 39. Ыо. 3, Маг. 2002, рр. 177-187. 8. Апзаг1, К. «ПК Н11Ьегг ТгапЫогщегз», 1ЕЕЕ Тта~и. Асоизг., 5реес)с 5!8па1Ртосезяп8, Уо!. 35, Ап8. 1987, рр. 1116-1119. 7. Магр1е 1г., 5. «Соарес)пй сЬе Г)1зсгесе-гппе 'Апа1ус1с' 518па1 ч(а ГЕТ», 1ЕЕЕ Ттапз, оп 5!8па! Ртос., Ъ'о1.
47, Ь(о. 9, 5ерп 1999, рр. 2600-2603. 8. Ке)11у, А., ек а1. «Апа1ус1с 518па1 бепегас(оп — Т1рз апс1 Тгарз», 1ЕЕЕ Ттапх оп 5!8па! Ртос., Уо1. 42, Ь(о. 11, Ыоч. 1994. эвг Глава!О. П еоб азованиечастотыдиск етизации 10.1. Прореживание Изменение частоты дискретизации встречается в двух разновидностях: как понижение и повышение. Понижение частоты дискретизации известно как децимация'. (Термин децимация не совсем подходит в этом случае, т. к.
изначально децимация понималась как уменьшение в десять раз. В настоящее время термин децимация используется для обозначения понижения частоты дискретизации в любое целое число раз.) Процесс повышения частоты дискретизации известен как интерполяция, т. е. получение оценок значений промежуточных отсчетов. Поскольку прореживание является самым простым случаем изменения частоты дискретизации, сначала рассмотрим именно его. Мы можем проредить, или понизить частоту дискретизации, последовательность дискретных отсчетов в Р раз, выбрав каждый Ь-й отсчет и отбросив все остальные отсчеты.
По отношению к исходной частоте дискретизации, 1 ы, новая частота дискретизации равна Уяеа УоЫ /Р (10-1) ф Хи«(П) 1 ° а ааа ааа а аа аа а а а а Н-Н-«~-~- «~ «««««««« ~-~. ««~-Ф (а) О -1 а а ° а а а 'а»» ~х (и) 1 °- -а а а а с (Ь) О ° а '"а Рис. 10.1. Преобразование частоты дискретизации: (а) исходная последовательность; (Ь) прореженная в три раза последовательность Спектральные эффекты прореживания вполне предсказуемы, как показано на рисунке 10.2, где ограниченный спектр непрерывного сигнала показан сплошной линией.
На рисунке 10.2 (а) показан периодический спектр Х )~(т) дискретной 1 Хотя в русскоязычной литературе понижении частоты дискретизации чаще называют прореживанием, здесь мы используем для обозначения этого приема транслитерацию английского термина дсснпзбоп, Причины тому две: в последнее время у нас получил распространение и термин децимация, а, кроме того, при использовании термина «прореживзние» следующее замечание автора было бы непонятно читателя> — (ирим.
иере»). Например, чтобы проредить последовательность хм)(л) в три раза (Р-З), мы оставляем хм,)(О) и отбрасываем х„)д(1) и х„)д(2), сохраняем х )д(3) и отбрасываем хмд(4) и х„)д(5), сохраняем х„)д(6) и так далее, как показано на рисунке 10.1. Таким образом,х„„(п) = х )д(Зп), где и - О, 1, 2 и т.д. Результат такого процесса прорежива- ниЯ идентичен РезУльтатУ исходной дискРетизации с частотой/я =~,)д/3длЯ получения х„,м(п).
звз 10.1. П о еживание последовательности хоы(п). Спектр Хаа (тл) при х„, (и) = хоь((Зп) показан на рисунке 10.2 (Ъ). Рисунок 10.2 иллюстрирует две важные особенности. Во-первых, Ха,„(т) можно было бы получить непосредственно, дискретизируя исходный непрерывный сигнал с частотой~оси, а не путем прореживания хос((п) в три раза. И, во-вторых, значение коэффициента прореживания ограничено шириной спектра В исходного сигнала. Мы должны выдерживать соотношение/'„„, > 2В, чтобы не допустить возникновения наложений после прореживания.