Радиотехника и электроника, 2019, T. 64, № 2, стр. 144-151

Цифровые методы линеаризации фазочастотных характеристик аналоговых фильтров

Ю. А. Гребенко 1, Р. И. Поляк 1*, А. И. Стариковский 2, Г. В. Куликов 2

1 Национальный исследовательский университет “МЭИ”
111250 Москва, ул. Красноказарменная, 14, Российская Федерация

2 МИРЭА – Российский технологический университет
119454 Москва, просп. Вернадского, 78, Российская Федерация

* E-mail: PoliakRI@mpei.ru

Поступила в редакцию 12.04.2018
После доработки 01.06.2018
Принята к публикации 11.06.2018

Полный текст (PDF)

Аннотация

Показана возможность линеаризации фазочастотной характеристики (ФЧХ) аналогового фильтра, включенного перед аналого-цифровым преобразователем с помощью цифрового фильтра с конечной импульсной характеристикой (КИХ). Путем схемотехнического моделирования подтверждено, что сквозная ФЧХ становится более линейной, а сквозная амплитудно-частотная характеристика (АЧХ) равна квадрату АЧХ аналогового фильтра. Также показана возможность линеаризации ФЧХ полосовых комплексных аналоговых фильтров. Приведены результаты моделирования последовательного соединения вещественного фильтра нижних частот, а также комплексного полосового фильтра, рассчитанных с использованием низкочастотного прототипа Баттерворта 5-го порядка, с соответствующими линеаризующими КИХ-фильтрами. Найдено, что неравномерность группового времени запаздывания в полосе пропускания может быть уменьшена более чем в 100 раз.

ВВЕДЕНИЕ

Система цифровой обработки сигналов (ЦОС) обычно содержит аналоговый входной фильтр, который используется для ограничения полосы частот входного аналогового сигнала перед его оцифровкой, чтобы уменьшить эффект наложения спектров при дискретизации. При аналого-цифровом преобразовании широкополосных сигналов часто требуется использовать аналоговый фильтр нижних частот (ФНЧ) высокого порядка, имеющий заметную нелинейность фазочастотной характеристики (ФЧХ) [1].

В радиоприемных устройствах с переносом спектра на низкую промежуточную частоту использование аналогового комплексного полосового фильтра (ПФ) высокого порядка после квадратурного преобразователя частоты позволяет обеспечить существенное подавление помехи по зеркальному каналу приема и других внеполосных помех [2]. В этом случае влияние нелинейности ФЧХ аналогового комплексного полосового фильтра на спектр сигнала промежуточной частоты может оказаться недопустимо большим.

В рассмотренных случаях возникает необходимость в разработке методов линеаризации ФЧХ аналогового вещественного ФНЧ или аналогового комплексного ПФ.

Линейной ФЧХ соответствует постоянство группового времени запаздывания (ГВЗ) в интересующем нас диапазоне частот. Неравномерность зависимости ГВЗ от частоты в рабочем диапазоне частот характеризует нелинейность ФЧХ и обычно приводится в качестве контролируемого параметра в техническом задании на устройство обработки сигнала. Предлагается для уменьшения неравномерности частотной зависимости ГВЗ аналогового вещественного ФНЧ или аналогового комплексного ПФ использовать соответственно вещественный или комплексный цифровой фильтр с конечной импульсной характеристикой (КИХ).

1. ЦИФРОВАЯ КОМПЕНСАЦИЯ НЕЛИНЕЙНОСТИ ФЧХ АНАЛОГОВОГО ФНЧ

Идея предлагаемого метода базируется на возможности компенсации нелинейности ФЧХ произвольного цифрового КИХ-фильтра, путем последовательного соединения с ним КИХ-фильтра, имеющего обратный порядок следования коэффициентов (назовем его обратным КИХ-фильтром). Продемонстрируем такую возможность.

Пусть передаточная функция КИХ-фильтра с нелинейной ФЧХ имеет вид

$T(z) = {{a}_{0}} + {{a}_{1}}{{z}^{{ - 1}}} + {{a}_{2}}{{z}^{{ - 2}}} + {{a}_{3}}{{z}^{{ - 3}}},$
где ${{a}_{0}},{{a}_{{13}}},{{a}_{2}},{{a}_{3}}$ – вещественные коэффициенты.

Тогда передаточная функция обратного КИХ-фильтра будет иметь вид

$\begin{gathered} {{T}_{{{\text{о б р }}}}}(z) = {{a}_{3}} + {{a}_{2}}{{z}^{{ - 1}}} + {{a}_{1}}{{z}^{{ - 2}}} + {{a}_{0}}{{z}^{{ - 3}}} = \\ = {{z}^{{ - 3}}}\left( {{{a}_{0}} + {{a}_{{1\quad}}}{{z}^{1}} + {{a}_{2}}{{z}^{2}} + {{a}_{3}}{{z}^{3}}} \right). \\ \end{gathered} $

Запишем частотную характеристику КИХ-фильтра с нелинейной ФЧХ

$\begin{gathered} T\left( {\exp \left( {j{\Phi }} \right)} \right) = {{a}_{0}} + {{a}_{{1\quad}}}\exp \left( { - j{\Phi }} \right) + \\ + \,\,{{a}_{2}}\exp \left( { - 2j{\Phi }} \right) + {{a}_{3}}\exp \left( { - 3j{\Phi }} \right) = \\ = {{a}_{0}} + {{a}_{1}}\left( {{cos\Phi } - j{sin\Phi }} \right) + \\ \end{gathered} $
$ + \,\,{{a}_{2}}\left( {{\text{cos}}2{\Phi } - j{\text{sin}}2{\Phi }} \right) + {{a}_{3}}\left( {{\text{cos}}3{\Phi } - j{\text{sin}}3{\Phi }} \right),$
где ${\Phi }$ – цифровая частота (в радианах).

Соответственно, частотная характеристика обратного КИХ-фильтра имеет вид

$\begin{gathered} {{T}_{{{\text{о б р }}}}}\left( {{\text{exp}}\left( {j{\Phi }} \right)} \right) = {\text{exp}}\left( { - 3j{\Phi }} \right)({{a}_{0}} + {{a}_{{1\quad}}}{\text{exp}}\left( {j{\Phi }} \right) + \\ + \,\,{{a}_{2}}{\text{exp}}\left( {2j{\Phi }} \right) + {{a}_{3}}{\text{exp}}\left( {3j{\Phi }} \right)) = \\ = {{a}_{0}} + {{a}_{1}}\left( {{cos\Phi } + j{sin\Phi }} \right) + \\ + \,\,{{a}_{2}}\left( {{\text{cos}}2{\Phi } + j{\text{sin}}2{\Phi }} \right) + {{a}_{3}}\left( {{\text{cos}}3{\Phi } + j{\text{sin}}3{\Phi }} \right). \\ \end{gathered} $

Найдем выражение для АЧХ КИХ-фильтров

$\begin{gathered} \left| {T\left( {{\text{exp}}\left( {j{\Phi }} \right)} \right)} \right| = \sqrt {{{{\left( {{{a}_{0}} + {{a}_{1}}{cos\Phi } + {{a}_{2}}{\text{cos}}2{\Phi } + {{a}_{3}}{\text{cos}}3{\Phi }} \right)}}^{{2\quad}}} + {{{\left( {{{a}_{0}} + {{a}_{1}}{sin\Phi } + {{a}_{2}}{\text{sin}}2{\Phi } + {{a}_{3}}{\text{sin}}3{\Phi }} \right)}}^{2}},} \\ \left| {{{T}_{{{\text{о б р }}}}}\left( {{\text{exp}}\left( {j{\Phi }} \right)} \right)} \right| = \sqrt {{{{\left( {{{a}_{0}} + {{a}_{1}}{cos\Phi } + {{a}_{2}}{\text{cos}}2{\Phi } + {{a}_{3}}{\text{cos}}3{\Phi }} \right)}}^{{2\quad}}} + {{{\left( {{{a}_{0}} + {{a}_{1}}{sin\Phi } + {{a}_{2}}{\text{sin}}2{\Phi } + {{a}_{3}}{\text{sin}}3{\Phi }} \right)}}^{2}}} . \\ \end{gathered} $

Можно сделать вывод, что АЧХ рассматриваемых КИХ-фильтров одинаковые. Найдем выражение для ФЧХ КИХ-фильтров:

$\begin{gathered} \arg T\left( {{\text{exp}}\left( {j{\Phi }} \right)} \right) = \\ = - {\text{arctg}}\quad\left( {\frac{{{{a}_{1}}{sin\Phi } + {{a}_{2}}{\text{sin}}2{\Phi } + {{a}_{3}}{\text{sin}}3{\Phi }}}{{{{a}_{0}} + {{a}_{1}}{cos\Phi } + {{a}_{2}}{\text{cos}}2{\Phi } + {{a}_{3}}{\text{cos}}3{\Phi }}}} \right), \\ \end{gathered} $
$\begin{gathered} \arg {{T}_{{{\text{о б р }}}}}\left( {{\text{exp}}\left( {j{\Phi }} \right)} \right) = \\ = - 3{\Phi } + {arctg\;}\left( {\frac{{{{a}_{1}}{sin\Phi } + {{a}_{2}}{\text{sin}}2{\Phi } + {{a}_{3}}{\text{sin}}3{\Phi }}}{{{{a}_{0}} + {{a}_{1}}{cos\Phi } + {{a}_{2}}{\text{cos}}2{\Phi } + {{a}_{3}}{\text{cos}}3{\Phi }}}} \right). \\ \end{gathered} $

При последовательном соединении КИХ-фильтров их АЧХ перемножаются, а ФЧХ – суммируются. Таким образом, АЧХ последовательного соединения рассматриваемых КИХ-фильтров равна квадрату АЧХ исходного КИХ-фильтра, а ФЧХ становится линейной.

Используя метод инвариантной импульсной характеристики [3], аналоговому ФНЧ можно поставить в соответствие цифровой КИХ-фильтр, коэффициенты которого являются отсчетами усеченной дискретизованной импульсной характеристики аналогового ФНЧ.

При достаточной частоте дискретизации и порядке такого КИХ-фильтра, его АЧХ и ФЧХ будут близки к АЧХ и ФЧХ аналогового ФНЧ в области частот от 0 до fд/2. Зная коэффициенты КИХ-фильтра, можно составить и реализовать передаточную функцию обратного КИХ-фильтра.

Предлагаемая методика состоит из следующих этапов:

1) определяем импульсную характеристику аналогового ФНЧ;

2) проводим дискретизацию импульсной характеристики с определенной частотой;

3) ограничиваем число отсчетов импульсной характеристики, обеспечивая допустимый уровень искажений АЧХ;

4) используя полученные отсчеты, составляем передаточную функцию обратного КИХ-фильтра.

Поясним предлагаемую методику на примере.

Рассмотрим ФНЧ Баттерворта 5-го порядка (рис. 1), АЧХ которого приведена на рис. 2 . Видно, что полоса пропускания ФНЧ равна 158.95 кГц. Непрерывная импульсная характеристика ФНЧ предствлена на рис. 3.

Рис. 1.

Модель аналогового ФНЧ Баттерворта 5-го порядка.

Рис. 2 .

Амплитудно-частотная характеристика A(f) аналогового ФНЧ.

Рис. 3.

Отсчеты импульсной характеристики (Fд = 1 МГц).

Проведем дискретизацию импульсной характеристики ФНЧ. Для уменьшения порядка линеаризующего КИХ-фильтра желательно уменьшать частоту дискретизации. Результаты проведенных исследований показали, что эта частота должна быть как минимум в шесть раз больше значения полосы пропускания аналогового ФНЧ. При этом с увеличением частоты дискретизации более чем в шесть раз относительно полосы пропускания эффективность линеаризации ФЧХ практически не изменяется.

Выберем частоту дискретизации Fд = 1 МГц. Отсчеты дискретизованной импульсной характеристики представлены на рис. 3. Количество используемых отсчетов импульсной характеристики определяет эффективность линеаризации ФЧХ. Ограничим импульсную характеристику 18 отсчетами. В этом случае импульсная характеристика будет представлена дискретной последовательностью $h(n) = \left\{ {{{h}_{1}},{{h}_{2}},{{h}_{3}}, \ldots ,{{h}_{{18}}}} \right\}\quad.$

Импульсная характеристика обратного КИХ-фильтра имеет обратный порядок следования отсчетов $\tilde {h}\quad(n) = \left\{ {{{h}_{{18}}},{{h}_{{17}}},{{h}_{{16}}}, \ldots ,{{h}_{1}}} \right\}.$ Структурная схема обратного КИХ-фильтра представлена на рис. 4. АЧХ такого фильтра практически совпадает с АЧХ аналогового ФНЧ.

Рис. 4.

Линеаризующий КИХ-фильтр.

На рис. 5 представлена зависимость ГВЗ от частоты в полосе пропускания для аналогового ФНЧ и обратного КИХ-фильтра. Среднее значение ГВЗ для аналогового фильтра и КИХ-фильтра составляет соответственно 4.193 и 13.758 мкс, а разница между минимальным и максимальным значением ГВЗ в полосе пропускания ∆ГВЗ = 1.911 и 1.910 мкс.

Рис. 5.

Зависимость ГВЗ в полосе пропускания для аналогового ФНЧ (2) и линеаризующего КИХ-фильтра (1) от частоты.

Используя моделирование в среде Micro-Cap, найдем частотные характеристики последовательного соединения аналогового ФНЧ и линеаризующего КИХ-фильтра (рис. 6). Полоса пропускания при последовательном соединении аналогового ФНЧ и линеаризующего КИХ-фильтра равна 145 кГц (см. вертикальную прямую на рис. 6).

Рис. 6.

Амплитудно-частотная характеристика A(f) последовательного соединения аналогового ФНЧ и линеаризующего КИХ-фильтра.

На рис. 7 представлена зависимость ГВЗ от частоты для последовательного соединения аналогового ФНЧ и линеаризующего КИХ-фильтра в полосе пропускания аналогового ФНЧ. Как видно, среднее ГВЗ последовательного соединения аналогового ФНЧ и линеаризующего КИХ-фильтра в полосе пропускания составляет 17.9 мкс, а неравномерность ГВЗ ∆ГВЗ = 0.144 мкс. Таким образом, неравномерность ГВЗ этого соединения уменьшилась по сравнению с неравномерностью ГВЗ аналогового ФНЧ в 13.3 раза, но если увеличить количество отсчетов импульсной характеристики с 18 до 35, то неравномерность ГВЗ этого соединения уменьшится по сравнению с неравномерностью ГВЗ аналогового ФНЧ в 140 раз.

Рис. 7.

Зависимость ГВЗ от частоты в полосе пропускания для последовательного соединения аналогового ФНЧ и линеаризующего КИХ-фильтра.

2. МЕТОДИКА РАСЧЕТА ЦИФРОВОГО КОМПЕНСАТОРА НЕЛИНЕЙНОСТИ ФЧХ АНАЛОГОВОГО КОМПЛЕКСНОГО ПОЛОСОВОГО ФИЛЬТРА

Принципы построения и методики расчета аналоговых комплексных полосовых фильтров хорошо известны. В работе [4], например, показано, что при использовании метода смещения частотных характеристик форма частотных характеристик при изменении частоты смещения не искажается, поэтому можно использовать результаты, полученные при линеаризации ФЧХ ФНЧ. Метод смещения частотных характеристик используется и при расчете цифровых комплексных полосовых КИХ-фильтров. Наиболее удобной формой реализации этого метода является использование комплексных задержек [5]. В такой ситуации результаты линеаризации ФЧХ ФНЧ могут полностью использоваться для построения комплексных полосовых фильтров с линейными ФЧХ.

Предлагаются следующие этапы линеаризации ФЧХ аналогового комплексного полосового фильтра:

1) разрабатываем схему аналогового ФНЧ на идентичных звеньях первого порядка [6] с полосой пропускания, равной половине полосы пропускания комплексного полосового фильтра;

2) в соответствии с изложенной выше методикой находим передаточную функцию и структурную схему КИХ-фильтра, компенсирующего нелинейность ФЧХ ФНЧ;

3) используя метод смещения по частоте [4], составляем схему комплексного фильтра, заменяя вещественные идентичные звенья первого порядка на комплексные звенья, сохраняя в структурной схеме значения коэффициентов обратной связи, рассчитанные для ФНЧ;

4) применяя метод комплексной задержки, преобразуем структурную схему КИХ-фильтра, линеаризующего ФЧХ ФНЧ, в структурную схему комплексного цифрового фильтра с центральной частотой, равной центральной частоте аналогового комплексного полосового фильтра.

Проиллюстрируем предложенную методику на конкретном примере.

Возьмем в качестве исходного вариант ФНЧ с НЧ-прототипом Баттерворта 5-го порядка, с полосой пропускания ∆f = 158.95 кГц (см. рис. 1). Структурная схема обратного КИХ-фильтра, использующего 18 отсчетов импульсной характеристики и линеаризующего ФЧХ ФНЧ, представлена на рис. 4.

Структурная схема комплексного полосового аналогового фильтра, полученная из структурной схемы ФНЧ (см. рис. 1) с использованием комплексных звеньев, показана на рис. 8. Принципиальная схема комплексного базового звена с центральной частотой, равной 318 кГц, реализованная в модели в виде макроса, приведена на рис. 9, АЧХ комплексного полосового фильтра, полученная путем моделирования, – на рис. 10.

Рис. 8.

Модель аналогового комплексного полосового фильтра (КЗ – комплексное звено).

Рис. 9.

Модель базового комплексного звена.

Рис. 10.

Амплитудно-частотная характеристика A(f) комплексного полосового фильтра.

Из рис. 10 видно, что центральная частота комплексного ПФ равна ${{f}_{0}} = 318\,\,\quad{\text{к Г ц ,}}$ а полоса пропускания $\quad\Delta f = 318\,\,\quad{\text{к Г ц ,}}$ что соответствует выбранным величинам. Неравномерность ГВЗ в полосе пропускания составляет ∆ГВЗ = 1.911 мкс.

Используя метод комплексной задержки, преобразуем структурную схему вещественного обратного КИХ-фильтра, линеаризующего ФЧХ ФНЧ, в структурную схему комплексного цифрового полосового фильтра с центральной частотой ${{f}_{0}} = 318\quad\,\,{\text{к Г ц }}$ (рис. 11). Структурная схема комплексной задержки, реализованная в виде макроса и используемая в модели, показана на рис. 12. Моделируем последовательное соединение комплексного аналогового полосового фильтра и комплексного линеаризующего цифрового полосового КИХ-фильтра. Неравномерность ГВЗ в полосе пропускания для такого соединения равна ∆ГВЗ = 0.138 мкс (рис. 13). Степень уменьшения неравномерности ГВЗ в полосе пропускания такая же, как при линеаризации ФЧХ аналогового ФНЧ, и равна 13.3.

Рис. 11.

Функциональная модель цифрового комплексного полосового КИХ-фильтра.

Рис. 12.

Структурная схема комплексной задержки.

Рис. 13.

Зависимость ГВЗ в полосе пропускания для последовательного соединения комплексного полосового фильтра и комплексного полосового линеаризующего КИХ-фильтра от частоты.

ЗАКЛЮЧЕНИЕ

Продемонстрирована возможность использования для ограничения спектра сигнала до аналого-цифрового преобразования аналоговых ФНЧ или аналоговых комплексных полосовых фильтров высокого порядка, имеющих нелинейные ФЧХ. Нелинейность ФЧХ аналоговых фильтров компенсируется за счет введения после АЦП цифровых вещественных или комплексных линеаризующих КИХ-фильтров. Неравномерность ГВЗ в полосе пропускания может быть уменьшена более чем в 100 раз.

Список литературы

  1. Айфичер Э., Джервис Б., Цифровая обработка сигналов: практический подход. М.: ИД “Вильямс”, 2008.

  2. Галкин В.А. Основы программно-конфигурируемого радио. М.: Горячая линия–Телеком, 2015.

  3. Г. Лэм. Аналоговые и цифровые фильтры. Расчет и реализация. М.: Мир, 1982.

  4. Гребенко Ю.А. // Вестн. МЭИ. 2003. № 1. С. 80.

  5. Гребенко Ю.А., Акар Мьо // Вестн. МЭИ. 2009. № 1. С. 70.

  6. Гребенко Ю.А. Однородные устройства обработки сигналов. М.: ИД МЭИ, 2009.

Дополнительные материалы отсутствуют.