В цифровую обработку

• Название: В цифровую обработку
• Дата: 01.03.2020
• Размер: 0.7 Mb.
• Имя файла: 010828.docx
• Тип: Документы; #110381
• страница: 13 из 23

Расчет начинается с выбора величины N.

Пусть N = 8. Отсюда интервал между выборками W₁ = = 0,125.

Формула для ФЧХ (3.3): j(w) = -wT . Отсюда

j (W) = -7pW, или для частот выборки j (kW₁) = -7pW₁,

Отсчеты АЧХ - по требуемой АЧХ на графике Рис. 3.8, а.

Следовательно, комплексные частотные отсчеты:

Н(jkW₁) = {1e ^j0; 1e ^-j0,875p ; 1e ^-j1,75p ; 0; 0; 0; -1e ^-j5,25p ; -1e ^-j6,125p }

Отсюда расчет импульсной характеристики по формуле обр. ДПФ

h (nT) = H (jkW₁) e ^{j (2p/N) kn} =

={0,065; -0,165; 0,025; 0,53; 0,53; 0,025; -0,165; 0,065}

что соответствует схеме фильтра на Рис. 3.8, б

Расчетная формула АЧХ такого типа фильтра - (3.5).

Поэтому Н(W) = 1,06 cos pW + 0,05 cos 3pW - 0,33 cos 5pW + 0,13 cos 7pW

Результаты расчета реализованной АЧХ приведены на графике Рис. 3.8, а (штриховая линия).

В окрестности точек разрыва требуемой АЧХ (в данном примере это частоты 0,25 и 0,75) отклонение от нормы реализованных характеристик получается значительным вследствие влияния эффекта Гиббса. Ослабить влияние эффекта Гиббса удается введением весовой функции (метод взвешивания) к импульсной характеристике.

Новая импульсная характеристика формируется по правилу:

h' (nT) = W (nt) * h (nT)

Где W (nT) - весовая функция или "сглаживающее окно".

Находят применение различные типы окон, например "окно" Хэмминга:

W(nT) = 0,54 + 0,46 cos [2p ], (3.9)

где n = 0, 1, 2, ... (N - 1)

Для рассматриваемого примера

W (nT) = {0,08; 0,244; 0,64; 0,96; 0,96; 0,64; 0,244; 0,08}

h' (nT) = {0,005; -0,04; 0,016; 0,51; 0,51; 0,016; -0,04; 0,005}

Отсюда новые коэффициенты фильтра и новая передаточная функция

H'(Z) = 0,005 - 0,04Z^-1 + 0,016Z^-2 + 0,51Z^-3 + 0,51Z^-4 + 0,016Z^-5 - 0,04Z^-6 +

+ 0,005Z^-7

График АЧХ с учетом сглаживающего окна приведен на Рис. 3.9. Расчетная функция получена из формулы для Н'(Z) после подстановки

Z = e^jwT = e^j2pW.

Сравнивая реализованные АЧХ на Рис. 3.8, а и Рис. 3.9, можно убедиться в улучшении качества аппроксимации требуемой АЧХ при введении "окна".

С ростом N положительный эффект от применения "сглаживающего окна" возрастает.



В рассмотренном примере нормы на отклонение реализованной АЧХ от требуемой не заданы. Если эти нормы не выполняются, то.... (строчка ксерокопии не влезла)
3.6. Метод частотной выборки

Коэффициенты не рекурсивного ЦФ (Рис. 3.2, а) соответствуют отсчетам импульсной характеристики. Схему не рекурсивного ЦФ можно преобразовать таким образом, чтобы коэффициенты фильтра соответствовали отсчетам другой системной характеристики - передаточной функции. Новая схема ЦФ является основой конструирования фильтров по методу частотной выборки.

3.6.1 Схема фильтра.

Схема фильтра формируется по результатам эквивалентных преобразований передаточной функции не рекурсивного ЦФ

H(Z) = a_n Z^-n

где в соответствии с формулой обратного ДПФ

a_n = h (nT) = H (jkw₁) e^j(2p/N)kn

следовательно

Н(Z) = H (jkw₁) e^j(2p/N)kn Z^-n = (e^j(2p/N)kn Z^-1)ⁿ

Применяя здесь формулу суммы N первых членов геометрической прогрессии



получаем

H(Z) = = P(Z) (3.10)

где

P(Z) = 1 - dZ^-N, F_k(Z) = 1 / (1 - b_kZ^-1), d = e^j2pk, b_k = e ^j2pk/N (3.11)



Схема фильтра, соответствующего (3.10), приведена на Рис. 3.10, а. Схемы звеньев фильтра, соответствующих (3.11), приведены на Рис. 3.10, б.

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

Нули и полюсы H(Z) (3.10), т.е. корни уравнений

1- e^j2pk Z^-N = 0, 1 - e ^j2pk/N Z^-1 = 0

Расположены на единичной окружности плоскости Z в точках

Z_k = e ^j2pk/N

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

Z_k = e ^-aT/N e ^j2pk/N, где aТ < 10^-5

что соответствует коэффициентам фильтра

d = e^-aT e ^j2pk, b_k = e^-aT e ^j2pk/N (3.12)

Небольшая поправка коэффициентов фильтра (3.12) практически не отразится на характеристиках фильтра.

3.6.2 Частотная характеристика фильтра

Частотная характеристика фильтра по методу частотной выборки получается подстановкой

Z = e^jwT,