1096

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

Лекция

Математика и математический анализ

Математические модели цифровых рекурсивных фильтров. Методика синтеза цифровых рекурсивных фильтров. Численное исследование методики синтеза цифровых рекурсивных фильтров.

Русский

2013-01-06

1.61 MB

130 чел.

Лекция

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

Рассматриваются следующие вопросы:

- математические модели цифровых рекурсивных фильтров;

- методика синтеза цифровых рекурсивных фильтров.

2.1. Математические модели цифровых рекурсивных фильтров

Свойство рекурсивный фильтр обозначает, что выходной сигнал фильтра зависит как от значения входного сигнала, так и от предшествующих значений выходного сигнала. Рекурсивный фильтр имеет импульсную характеристику бесконечной длины и поэтому еще называется фильтром с импульсной характеристикой бесконечной длины (или БИХ-фильтром). АЧХ рекурсивного фильтра точнее аппроксимирует АЧХ аналогового прототипа, чем нерекурсивные.

Рекурсивные ФНЧ, ФВЧ, ФП, ФР выглядят следующим образом

,  (2.1)

,  (2.2)

,  (2.3)

,  (2.4)

где   - порядок фильтра, .

Связь между передаточной функцией (частотной характеристикой)  и переходной функцией (импульсной характеристикой)  цифрового рекурсивного фильтра, например рекурсивного ФНЧ, может быть представлена в виде (2.5) или с помощью  z-преобразования в виде (2.6).

, (2.5)

,   (2.6)                     

где   - шаг квантования,   - частота дискретизации.

Амплитудно-частотная характеристика (АЧХ) (2.7) и фазо-частотная характеристика (ФЧХ) (2.8) рекурсивного фильтра выглядит в виде

(2.7)

 (2.8)

Для выделения частотных составляющих звуков используются каскады фильтров, т.е. последовательное соединение фильтров, например ФНЧ 1-го и 2-го порядков.

 (2.9)

где  - передаточная функция i-го ФНЧ 1-го или 2-го порядка,

– количество фильтров в каскаде.

Если порядок заданного фильтра N четный, то =M/2 и используется M/2 фильтров 2-го порядка, иначе =(M+1)/2  и используется (M-1)/2 фильтров 2-го порядка и один фильтр 1-го порядка.

Обработку дискретного сигнала рекурсивным фильтром, например ФНЧ (рис.2.1), можно представить в виде параметрической структуры, в которой используются коэффициенты фильтра , , и его порядок M .

 

Рис. 2.1. Структура фильтрации дискретного сигнала

В качестве входных переменных блока фильтрации выступают дискретный сигнал x(n), его  длина N, функция окна w(n), частота среза , частота дискретизации , а в качестве выходных – фильтрованный сигнал xН(n).

Цифровая фильтрация рекурсивным фильтром требует  умножений.

2.2 Методика синтеза цифровых рекурсивных фильтров

Цифровые рекурсивные фильтры можно реализовывать с различными характеристиками: отсутствие пульсации в полосе пропускания, отсутствие пульсации в полосе непропускания (запирания), крутизна АЧХ при переходе от полосы пропускания в полосу непропускания, пригодность к передаче импульсов. Ни один из рекурсивных фильтров не отвечает этим характеристикам. Согласно современным исследованиям, в преобразованиях аналоговых и цифровых сигналов наиболее распространенными являются фильтры Баттерворта, Бесселя, Чебышева, Кауэра.

Алгоритм синтеза цифрового фильтра приведен на рис.2.2.

Фильтр Баттерворта. Оптимизирован для максимально плоской АЧХ в полосе пропускания. Крутизна АЧХ увеличивается с порядком фильтра. Является наиболее распространенным и требующим меньшей вычислительной сложности при реализации по сравнению с остальными фильтрами.

Рис. 2.2. Алгоритм синтеза цифровых фильтров

Аппроксимирующая функция (функция Баттерворта) представлена в виде (2.10).

, где  M – порядок фильтра (2.10)

Для ФНЧ ,  для ФВЧ  .

Ослабление АЧХ вычисляется в виде (2.11).

 (2.11)

Для значительного ослабления (затухания) АЧХ в полосе непропускания необходимо брать  M>>3.

Если , то полюса передаточной функции K(), лежащие в левой полуплоскости, вычисляются в виде (2.12).

,   (2.12)

ФНЧ Бесселя. Имеет лучшую передаточную характеристику для прямоугольных импульсов. АЧХ спадает менее круто, чем в фильтре Баттерворта. Отсутствует пульсация.

, (2.13)

где  - функция Бесселя M-го порядка, .

Для ФНЧ ,  для ФВЧ  .

ФНЧ Чебышева. В полосе пропускания либо полосе непропускания имеет пульсацию с максимальной амплитудой (в отличие от фильтра Баттерворта). Крутизна АЧХ больше, чем у фильтра Баттерворта.

Аппроксимирующая функция представлена в виде (2.14).

,  (2.14)

где   - полином Чебышева  M-го порядка

Для ФНЧ ,  для ФВЧ .

Ослабление АЧХ вычисляется в виде (2.15).

 (2.15)

Амплитуда пульсации АЧХ определяется в виде (2.16)

 (2.16)

Если , то полюса передаточной функции K(), лежащие в левой полуплоскости, вычисляются в виде (2.17).

. (2.17)

ФНЧ Кауэра (эллиптический ФНЧ). Имеет пульсацию в полосе пропускания и полосе непропускания. Это цена за большую крутизну фронта импульса фильтра.

Аппроксимирующая функция представлена в виде (2.18).

,  (2.18)

где   - эллиптическая функция

Для ФНЧ , для ФВЧ .

Ослабление АЧХ вычисляется в виде (2.19).

 (2.19)

В блоке 2 определим с помощью нормированных полюсов, например с помощью (2.12), коэффициенты передаточной функции аналогового ФНЧ 1-го порядка в виде (2.20), а коэффициенты передаточной функция аналогового ФНЧ 2-го порядка в виде (2.21).

,   (2.20)

, (2.21)

.

Цифровой ФНЧ получают (синтезируют) по аналоговому ФНЧ. В   блоке 3 коэффициенты цифрового ФНЧ определены следующим образом. Используя билинейное преобразование, произведем в (2.20) и (2.21) замену  

, ,

где  - шаг квантования,  - частота дискретизации.

Тогда передаточная функция, определенная в (2.6), для ФНЧ 1-го (2.22) и 2-го (2.23) порядка будет представлена в виде

, (2.22)

, (2.23)

Цифровой ФВЧ получают из цифрового ФНЧ. В блоке 4 коэффициенты цифрового ФВЧ определены следующим образом. Используя частотное преобразование, произведем в (2.22) и (2.23) замену  

,

где .

Тогда передаточная функция для ФВЧ 1-го (2.24) и 2-го (2.25) порядка будет представлена следующим образом

, (2.24)

, (2.25)

Цифровой ФП получают из цифрового ФНЧ. В блоке 4 коэффициенты цифрового ФП определены следующим образом. Используя частотное преобразование, произведем в (2.24) и (2.25) замену  

,

где , , , .

Тогда передаточная функция, для ФП 1-го (2.22) и 2-го (2.23) порядка будет представлена следующим образом

, (2.26)

, (2.27)

Цифровой ФР получают из цифрового ФНЧ. В блоке 4 коэффициенты цифрового ФР определены следующим образом. Используя частотное преобразование, произведем в (2.22) и (2.23) замену  

,

где , , , .

Тогда передаточная функция, для ФР 1-го (2.28) и 2-го (2.29) порядка выглядит следующим образом

, (2.28)

, (2.29)

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

,  (2.30)

где   - передаточная функция реализованного фильтра,

- передаточная функция идеального фильтра, .

Приведем пример определения среднеквадратичной ошибки  для ФНЧ.

Если заменить

,

то для ФНЧ 1-го порядка передаточные функции (2.22) примут вид (2.31).

,  (2.31)

Отсюда

 (2.32)

Для фильтра 2-го порядка передаточные функции (2.23) примут вид (2.33).

, (2.33)

Отсюда

 (2.34)

На основании (2.32) и (2.34) вычисляется среднеквадратичная ошибка , которая затем сравнивается с порогом    и в случае превышения порога производится уточнение коэффициентов ФНЧ. Аналогично вычисляются , , .

2.3 Численное исследование методики синтеза цифровых рекурсивных фильтров

В работе [1] осуществлялось подавление высокочастотных составляющих дискретного сигнала. Для этого выбран цифровой ФНЧ восьмого порядка, который синтезируется из аналогового ФНЧ, составленного из четырех Г-образных ФНЧ второго порядка. При этом в табл.2.1 заданы величины емкостей , индуктивностей , сопротивлений .

Таблица 2.1

Заданные значения элементов аналоговых ФНЧ

№ фильтра

C (мкФ)

L (Гн)

R (Ом)

1

1.632

0.062

500

2

0.286

0.354

1000

3

0.096

1.059

2000

4

0.041

2.498

4000

Согласно (2.15)-(2.17), приведенным в части 1, определены параметры этих фильтров, представленные в табл.2.2.

Таблица 2.2

Расчетные значения параметров аналоговых ФНЧ

№ фильтра

(мкс)

Q

1

612.894

1.632

2.563

2

1745.375

0.573

0.9

3

2612.139

0.383

0.601

4

3081.228

0.325

0.51

Расчет фильтров произведен по алгоритму, приведенному на рис.2.2. Согласно (2.12), в блоке 1 рис.2.2 вычислены нормированные полюса аналоговых ФНЧ, их значения приведены в табл.2.3.

Таблица 2.3

Значения нормированных полюсов ФНЧ

p1

p2

p3

p4

-0.195+ 0.981j

-0.556+ 0.831j

-0.831+ 0.556j

-0.981+ 0.195j

p5

p6

p7

p8

-0.981-0.195j

-0.831-0.556j

-0.556-0.831j

-0.195-0.981j

В блоке 2, согласно (2.21), осуществлен расчет коэффициентов каскада аналоговых ФНЧ, которые приведены в табл.2.4

Таблица 2.4

Расчетные значения параметров аналоговых ФНЧ

№ фильтра

a0на = b2на

b1на

1

9869604.401

1225.789

2

9869604.401

3490.751

3

9869604.401

5224.278

4

9869604.401

6162.456

Для значений нормированных полюсов, приведенных в табл.2.3, в блоке 3, согласно (2.23), вычисляются коэффициенты каскада цифровых ФНЧ, которые приведены в табл.2.5.

Таблица 2.5

Расчетные значения коэффициентов цифровых ФНЧ

фильтра

a0нц = a2нц 

a1нц

b1нц

b2нц

1

0.005

0.010

-1.926

0.946

2

0.005

0.009

-1.835

0.854

3

0.005

0.009

-1.771

0.789

4

0.004

0.009

-1.738

0.756

Для обоснования правомочности выбора ФНЧ восьмого порядка, в блоке 5 вычисляется среднеквадратичная ошибка (2.30), значения которой приведены в табл.2.6. Данные этой таблицы позволяют сделать вывод, что с возрастанием кратности наборов фильтров среднеквадратичная ошибка резко уменьшается для ФНЧ восьмого порядка (почти вдвое) и незначительно уменьшается для ФНЧ кратности 10. На основании анализа среднеквадратичных ошибок делается вывод о целесообразности применения ФНЧ восьмого порядка.

Таблица 2.6

Среднеквадратичные ошибки оценивания ФНЧ

M=2

M=4

M=8

M=10

470.775

350.240

273.821

272.125

Результат применения ФНЧ восьмого порядка к дискретному сигналу (рис.2.3), содержащему высокочастотный звук «у», приведен на рис.2.4.

Рис. 2.3. Исходный сигнал, содержащий звук «у»

Рис. 2.4. Отфильтрованный сигнал, содержащий звук «у»


x (n),
дб

N

xн (n), дб

N


 

А также другие работы, которые могут Вас заинтересовать

2437. Расчет горно-транспортного комплекса в соответствии с эксплуатационной производительностью экскаватора ЭКГ-8,3 715.56 KB
  Развитие открытых горных разработок в ближайшие 10-15 лет определяется не только все возрастающими потребностями общества в минеральном сырье и конкурентоспособностью с подземным способом разработки, но также все увеличивающимися экологическими ограничениями.
2438. Експлуатація і обслуговування машин 902 KB
  Мета і завдання курсу Експлуатація і обслуговування машин. Планування технічного обслуговування та ремонту будівельних машин. Види і періодичність ТО тракторів і машин. Вплив степеня очищення повітря. Вплив стану експлуатаційних матеріалів на надійність роботи двигуна. Зберігання і консервація машин.
2439. Исследование характеристик тиристора и управляемого выпрямителя 1.67 MB
  Целью работы является: исследование вольтамперной характеристики и определение параметров тиристора, получение семейства статических характеристик тиристора, исследование работы регулируемого однополупериодного выпрямителя.
2440. Строение и характеристики подъемных кранов 105.05 KB
  Назначение, общее устройство неходовой рамы, способы крепления ее к раме автомобиля. Назначение и общее устройство выносных опор. Проверка уровня масла в редукторе механизма подъема груза. Назначение, устройство и работа опорно-поворотного устройства. Виды грузозахватных приспособлений и требования к ним. Грузы, запрещенные для подъема краном. Проверить уровень масла в картере редуктора механизма вращения.
2441. Оценка качества и его потребительские свойства 153.5 KB
  Качество как категория управления - интегральное понятие, которое характеризует эффективность всех сторон жизнедеятельности предприятия от процесса разработки стратегии, маркетинга и далее включая все этапы жизненного цикла продукта или услуги. Стоимость – это категория исключительного товарного (рыночного) хозяйства, представляющая собой общественное отношение товаропроизводителей, которое возникает при обмене товаров.
2442. Теория менеджмента. Менеджмент как учебная дисциплина 296.5 KB
  Экономические методы управления. Социально-психологические методы управления. Понятие и содержание коммуникации. Менеджмент как особый тип управления организацией в рыночной системе хозяйствования. Механизмы рыночной экономики. Концепция жизненного цикла организации. Актуальность управления персоналом. Государственная кадровая политика Республики Беларусь.
2443. Загальна медична психологія 23.08 KB
  Особливості психіки пацієнтів на етапах, підготовки, виконання хірургічних втручань та в післяопераційний період. Особливості психіки хворих, що страждають різними захворюваннями (серцево-судинні, інфекційні, гінекологічні, і т.д.). Особливості психічних процесів хворих з дефектами органів та систем (сліпота, глухота і т.д.).
2444. Податкова система і податкова політика 23.26 KB
  Суть, значення та роль податкової системи в системі державних доходів. Податкова політика та її характеристика. Склад податкової системи України. Порівняльна характеристика податкових систем у розвинутих державах світу
2445. Маркетинг продаж 519.07 KB
  Значение маркетинга в повышении эффективности продаж. Особенности продажи различных товаров и услуг. Управление персоналом службы продаж. Обучение и мотивация. Термин продажа обычно истолковывается очень широко и охватывает весь спектр действий, направленных на то, чтобы помочь клиенту сделать выбор и совершить покупку.