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


 

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

27464. Действие нормативно-правовых актов во времени, в пространстве и по кругу лиц 27.5 KB
  Действие нормативноправового акта порождение тех юридических последствий которые в нем предусмотрены. Действие т. во времени Действие нормативноправовых актов во времени определяется двумя моментами: моментом вступления нормативноправового акта в силу и моментом утраты им юридической силы. Существуют такие понятия как футуроспективное действие закона обращен в будущее рассчитан на facta futura ретроактивное распространение на отношения возникшие до вступления данного акта в силу и ультраактивное действие акт утрачивает силу а...
27465. Действие права 28.5 KB
  Действие права это действие составляющих его разнообразных норм которые в официальной форме выражены в различных нормативноправовых и иных правоустановительных актах. Эти характеристики пределов действия норм права официально закрепляются в соответствующих нормативноправовых актах. Действие во времени определяется периодом времени от момента вступления нормативноправового акта и содержащихся в нем норм в силу до момента ее утраты. Время вступления в законную силу разных нормативноправовых актов определяется различными способами: 1...
27466. Демократический политический режим 24.5 KB
  Государственнополитический режим это элемент формы государства характеризующий совокупность приемов методов способов и средств осуществления государственной власти. Основными признаками демократического режима являются: господство закона; разделение властей; наличие у граждан реальных политических и социальных прав и свобод и их юридическая защищенность и гарантированность; выборность и сменяемость центральных и местных органов государственной власти их подотчетность перед избирателями гласность; наличие свободно формируемых...
27467. Законность, понятие и характеристика 25.5 KB
  Понятие законности неотделимо от понятия права. В результате реализации этого требования органами государства и его должностными лицами возникает особый режим система общественных отношений режим законности. Законность как принцип одно из руководящих системообразующих начал построения и деятельности правового государства выражающееся в системе принципов и требований отражающих идею законности. Режим законности такое состояние общественных отношений когда во взаимоотношениях между государственной властью и населением господствует...
27468. Законотворчество в РФ 31.5 KB
  104 Конституции РФ принадлежит Президенту РФ Совету Федерации членам Совета Федерации депутатам Государственной Думы Правительству РФ законодательным представительным органам субъектов Федерации. 50 1 голос; федеральные же конституционные законы считаются принятыми если за них проголосовало не менее 2 3 голосов от общего числа депутатов Государственной Думы; б одобрение закона Советом Федерации в соответствии с ч. 105 Конституции РФ федеральный закон считается одобренным Советом Федерации если за него проголосовало более...
27469. К каким источникам права относятся, либо с ними связаны: договор купли-продажи, договор найма жилого помещения, союзный договор между суверенными республиками 25 KB
  К каким источникам права относятся либо с ними связаны: договор куплипродажи договор найма жилого помещения союзный договор между суверенными республиками. Виды источников: правовой обычай правовой прецедент правовая доктрина договоры нормативного содержания нормативный правовой акт. Договор куплипродажи и договор найма жилого помещения типовые договоры. Типовой договор стандартная типовая отпечатанная типографским способом форма договора со стандартными условиями и формулировками.
27470. К какому виду норм можно отнести ст.7 ГК РФ 40 KB
  Гражданское законодательство и нормы международного права 1. Общепризнанные принципы и нормы международного права и международные договоры РФ являются в соответствии с Конституцией РФ составной частью правовой системы РФ. Здесь присутствуют нормы: 1. исходные нормыначала 2.
27471. К какому виду норм относятся: правила дорожного движения; правила безопасной эксплуатации домашних бытовых приборов; правила, регламентирующие порядок проведения митингов и демонстраций 30.5 KB
  1 Правила дорожного движения общие императивные запрещающие нормы. Общие нормы это предписания охватывающие своим действием как правило все правовые институты той или иной отрасли. Императивные нормы категорические строго обязательные не допускающие отступлений и иной трактовки предписания. Императивными являются большинство норм права относящихся к различным его отраслям а исходные юридические нормы будут таковыми всегда.
27472. К какому виду относятся правоотношения собственности (обще-регулятивные, относительные, абсолютные, правильного ответа нет) 26 KB
  К какому виду относятся правоотношения собственности общерегулятивные относительные абсолютные правильного ответа нет Правовые отношения могут классифицироваться по различным основаниям: l в зависимости от предмета правового регулирования отраслевого признака они подразделяются на: конституционные административные уголовные гражданские и т.; ПС: материальные 3 в зависимости от функциональной роли на регулятивные возникают на основе норм права или договора и охранительныесвязаны с государственным принуждением и реализацией...