Слайд 1Лекция 12
Эллиптический фильтр
Эллиптический фильтр в некотором смысле
объединяет в себе свойства фильтров Чебышева первого и второго рода,
поскольку АЧХ эллиптического фильтра имеет пульсации заданной величины, как в полосе пропускания, так и в полосе задерживания.
На трех рисунках показаны АЧХ фильтра Чебышева первого рода, второго рода и АЧХ эллиптического фильтра.
Слайд 2На первом рисунке у фильтра Чебышева первого рода уровень пульсации
в полосе пропускания взят равным 0.5 дБ (Rp = 0.5).
Слайд 3На втором рисунке, у фильтра Чебышева второго рода, уровень пульсации
в полосе задерживания взят равным 20 дБ (Rs = 0.5).
Слайд 4На третьем рисунке, у эллиптического фильтра, уровень пульсации в полосе
пропускания взят равным 0.5 дБ, в полосе задерживания взят равным
20 дБ .
Слайд 5 Эллиптический фильтр в советской и российской литературе называют
также фильтром Золотарева – Кауэра. В английской литературе называют фильтром
Кауэра (Cauer filter).
Эллиптический фильтр имеет максимально возможную крутизну ската АЧХ, то есть минимально возможную переходную область Δω между полосами пропускания и задерживания.
Можно доказать, что с точки зрения минимальной ширины переходной полосы эллиптический фильтр является оптимальным, т.е. для заданных порядка фильтра и уровня пульсаций, не существует других фильтров с более быстрым переходом от полосы пропускания к полосе задерживания.
Слайд 6 Эллиптический фильтр имеет как полюсы, так и нули.
Число полюсов Np определяется порядком фильтра. Число нулей Nz равно
максимальному четному числу, не превосходящему порядок фильтра
На рисунке показаны полюса и нули эллиптического фильтра.
Слайд 7В рассматриваемом примере число полюсов равно 5, число нулей равно
Слайд 8 АЧХ эллиптического фильтра, в общем случае описывается следующей
аналитической формулой.
(1)
Здесь - ω0 частота среза,
n – порядок фильтра, ε и L – параметры, определяющие величину пульсаций в полосах пропускания и задерживания. Функция Rn (x , L) называется рациональной функцией Чебышева. Эта функция выражается через эллиптические функции Якоби. Отсюда и название фильтра. Эти функции достаточно сложны, и на наше счастье MATLAB выполняет за нас все трудоемкие расчеты.
Для анализа эллиптических фильтров нам достаточно знать несколько основных параметров. Это n – порядок фильтра, ω0 - частота среза, Rp- уровень пульсации АЧХ в полосе пропускания и Rs - уровень пульсации АЧХ в полосе задерживания.
Слайд 9Для примера, на этом рисунке показана АЧХ эллиптического фильтра 4-го
порядка. Уровни пульсации имеют следующие значения.
Слайд 10 Уровни пульсации позволяют найти следующие характерные значения АЧХ.
Наименьшее значение пульсаций Ap в полосе пропускания находится по формуле.
(2)
Здесь значение A0 = 1. Наибольшее значение As пульсаций в полосе задерживания находится по формуле
(3)
Слайд 11Используя формулы (2), (3) для выбранных значений уровней пульсации получаем
Действительно такие значения показаны на рисунке. Теперь построим АЧХ
в логарифмическом масштабе, где АЧХ будет измеряться в децибелах. От относительных единиц к децибелам АЧХ преобразуется по следующей формуле
(4)
На следующем рисунке показана АЧХ в децибелах.
Слайд 12Здесь выбран интервал от –30 дБ до 5 дБ. Из
этого рисунка хороши видно, что уровень пульсации в полосе задерживания
равен 20 дБ.
Слайд 13 Уровень пульсации в полосе пропускания на этом рисунке
виден плохо. Поэтому построим АЧХ в интервале от –1 дБ
до 1 дБ. На следующем рисунке показана АЧХ в этом интервале.
Из этого рисунка хороши видно, что уровень пульсации в полосе пропускания равен 0.5 дБ.
Слайд 15Отсюда получаем общее правило для эллиптического фильтра - порядок фильтра
равен суммарному числу минимумов и максимумов АЧХ в полосе пропускания.
Слайд 16Преобразование фильтров
Мы рассмотрели несколько известных фильтров –
фильтр Баттерворта, фильтр Чебышева первого рода, фильтр Чебышева второго рода,
эллиптический фильтр.
Оказывается, что с помощью нескольких рассмотренных фильтров можно получить огромное количество других фильтров, с иными свойствами. За основу берется передаточная функция H(s) исходного фильтра. Затем производится некоторое математическое преобразование над комплексным аргументом s .
(5)
Слайд 17 Здесь L(s) – некоторое преобразование. В результате появляется
другая передаточная функция.
(6)
Эта новая передаточная функция описывает
новый фильтр со свойствами отличными от исходного фильтра.
Поэтому одним из способов получения новых фильтров, является нахождение нужного преобразования L.
Слайд 18Изменение частоты среза низкочастотного фильтра ФНЧ
Частотную характеристику
ФНЧ можно сжимать и растягивать вдоль оси частот простым масштабированием
частотной оси. Для этой цели подойдет простое линейное преобразование.
(7)
Здесь ω0 - задаваемая частота среза. На рисунке показаны АЧХ двух эллиптических фильтров со следующими частотами среза ω0 = 1 рад/с и ω0 = 3 рад/с.
Слайд 19Из рисунка видно, что вторая АЧХ получена из первой АЧХ
путем растяжения в три раза.
Слайд 20Первый фильтр имеет следующие коэффициенты основного уравнения.
b =
0.30 0
0.88 0 0.59
a = 1.00 1.15 2.22 1.76 1.20 0.59
Второй фильтр имеет уже другие коэффициенты
b1 = 0.92 0 23.79 0 145.39
a1 = 1.00 3.45 20.01 47.71 97.22 145.39
Так как меняются коэффициенты bn , an фильтра, то меняются так же нули и полюса фильтра. На следующих двух рисунках показаны в комплексной s - плоскости нули и полюсы первого и второго фильтров.
Слайд 23Преобразование ФНЧ в фильтр высокой частоты ФВЧ
Из
фильтра низкой частоты можно получить фильтр высокой частоты с помощью
следующего преобразования.
(8)
Здесь ω0 - задаваемая частота среза.
На следующих двух рисунках показаны АЧХ двух фильтров, связанных преобразованием (8).
На первом рисунке показана АЧХ исходного эллиптического ФНЧ. На втором рисунке показана АЧХ фильтра высоких частот ФВЧ полученного из первого фильтра с помощью преобразования (8). Частота среза взята равной ω0 = 40 рад/с.
Слайд 26 Из рисунков видно, что действительно простое преобразование (8)
позволило получить из фильтра низких частот фильтр высоких частот.
Посмотрим, как изменяются нули и полюсы передаточной функции при этом преобразовании.
На следующих двух рисунках показаны в комплексной s - плоскости нули и полюсы ФНЧ а также нули и полюсы ФВЧ, полученного из фильтра низких частот.
Слайд 29Первое, на что следует обратить внимание, это то, что число
нулей и полюсов осталось прежним, 4 нуля и 5 полюсов.
Второе,
численные значения нулей и полюсов изменились.
Третье, взаимное расположение нулей и полюсов в комплексной плоскости так же изменилось.
Четвертое, полюса обоих фильтров лежат в левой части комплексной плоскости. Напомним, что для аналоговых фильтров условие устойчивости сводится к тому, чтобы полюсы передаточной функции имели отрицательную мнимую часть.
(9)
Таким образом, преобразование (8) из устойчивого фильтра низких частот, позволяет получить устойчивый фильтр ФВЧ.
Слайд 30Преобразование ФНЧ в полосовой фильтр
Если использовать более
сложное преобразование, чем преобразования (7) и (8), то можно получить
более сложный фильтр. Например, чтобы получить передаточную функцию полосового фильтра из передаточной функции ФНЧ, нужно использовать следующее преобразование.
(10)
Слайд 31 Полосовой фильтр характеризуется нижней границей ω1 полосы пропускания,
и верхней границей полосы ω2 пропускания. Параметры преобразования (10) связаны
с нижней и верхней границами полосы пропускания полосового фильтра следующими соотношениями.
(11)
На следующих двух рисунках показаны АЧХ двух фильтров, связанных преобразованием (10).
Частоты нижней и верхней границ полосы пропускания полосового фильтра взяты равными следующим величинам ω1 = 10 рад/с и ω2 = 30 рад/с.
Слайд 32На первом рисунке показана АЧХ исходного эллиптического ФНЧ.
Слайд 33На втором рисунке показана АЧХ полосового фильтра.
Слайд 34 Из рисунков видно, что преобразование (10) позволило получить
из ФНЧ полосовой фильтр.
Посмотрим, как изменяются нули и
полюсы передаточной функции при этом преобразовании.
На следующих двух рисунках показаны в комплексной s - плоскости нули и полюсы ФНЧ, а также нули и полюсы полосового фильтра, полученного из фильтра низких частот.
Слайд 36Первое, на что следует обратить внимание, это то, что изменилось
число нулей и полюсов. Если внимательно рассмотреть рисунок, то можно
увидеть 10 полюсов и 9 нулей.
Второе, численные значения нулей и полюсов изменились.
Слайд 37Преобразование ФНЧ в режекторный фильтр
Для преобразования фильтра
низких частот в режекторный фильтр надо совершить преобразование до некоторой
степени обратное к преобразованию (10), когда из ФНЧ получается полосовой фильтр. Дело в том, что полосовой и режекторный фильтры обратные друг другу. Там где у полосового фильтра полоса пропускания, у режекторного фильтра полоса задерживания, и наоборот.
Поэтому, чтобы получить передаточную функцию режекторного фильтра из передаточной функции ФНЧ, нужно использовать следующее преобразование.
Слайд 38(12)
Параметры преобразования (12) связаны с нижней ω1
и верхней ω2 границами полосы задерживания полосового фильтра следующими соотношениями.
(13)
На следующих двух рисунках показаны АЧХ двух фильтров, связанных преобразованием (12).
Частоты нижней и верхней границ полосы задерживания полосового фильтра взяты равными следующим величинам ω1 =10 рад/с и ω2 = 30 рад/с .
Слайд 39На первом рисунке показана АЧХ исходного эллиптического ФНЧ.
Слайд 40На втором рисунке показана АЧХ режекторного фильтра.
Слайд 41 Из рисунков видно, что преобразование (12) позволило получить
из ФНЧ режекторный фильтр.
Посмотрим, как изменяются нули и
полюсы передаточной функции при этом преобразовании.
На следующих двух рисунках показаны в комплексной s - плоскости нули и полюсы ФНЧ, а также нули и полюсы режекторного фильтра, полученного из фильтра низких частот.
Слайд 43Первое, на что следует обратить внимание, это то, что изменилось
число нулей и полюсов. Если внимательно рассмотреть рисунок, то можно
увидеть 10 полюсов и 10 нулей.
Второе, численные значения нулей и полюсов изменились.
Слайд 44Проектирование дискретных фильтров
Под проектированием (или синтезом) цифровых
фильтров понимается выбор таких наборов коэффициентов {an} ,{bn} основного уравнения
фильтра, при которых характеристики получающегося фильтра удовлетворяют заданным требованиям.
Существует большое количество различных методов синтеза дискретных фильтров. Мы рассмотрим некоторые наиболее популярные методы.
Слайд 45Синтез рекурсивного фильтра по аналоговому прототипу
Напомним, что
если хотя бы один коэффициент an отличен от нуля, то
фильтр называется рекурсивным фильтром. Рекурсивный фильтр – это фильтр с обратными связями.
Рекурсивные фильтры могут быть как КИХ – фильтрами, так и БИХ – фильтрами. КИХ – фильтры это фильтры с конечной импульсной характеристикой, БИХ – фильтры это фильтры с бесконечной импульсной характеристикой.
При проектировании дискретного фильтра с заданными характеристиками одним из методов является метод, когда за основу берется похожий аналоговый фильтр.
Аналоговый фильтр в этом случае принято называть фильтром-прототипом (на английском языке Referenced Filter).
Слайд 46 При проектировании дискретного фильтра по аналоговому прототипу необходимо
совершить переход из s - области в z - область.
Другими словами, необходимо преобразовать функцию передачи Ha (s) аналогового фильтра в функцию передачи Hd (z) дискретного фильтра.
Вспомним некоторые основные соотношения теории аналоговых и дискретных фильтров.
Сначала рассмотрим аналоговые фильтры.
Слайд 47Основным уравнением аналогового фильтра является следующее дифференциальное уравнение.
(14)
Аналоговый
фильтр описывается также функцией времени h(t) - импульсной характеристикой. Импульсная
характеристика определяется из уравнения, которое связывает входящий x(t) и выходящий y(t) сигналы.
(15)
Слайд 48Уравнение (15) означает, что выходящий сигнал является сверткой импульсной характеристики
и входящего сигнала.
Передаточная функция H(s) аналогового фильтра является изображением Лапласа
импульсной характеристики h(t).
(16)
Здесь s - комплексное число
Слайд 49Если заданы коэффициенты bn , an основного уравнения, то передаточная
функция определятся следующей формулой.
(17)
Разложив числитель и знаменатель передаточной
функции (17) на элементарные множители, мы получаем передаточную функцию в следующем виде.
(18)
Здесь k = bm /an - коэффициент усиления, zi – нули передаточной функции, pi - полюсы передаточной функции.
Слайд 50Комплексный коэффициент передачи фильтра является преобразованием Фурье импульсной характеристики.
(19)
АЧХ и ФЧХ аналогового фильтра определяются формулами.
(20)
Анализ формул (16)
и (19) показывает, что комплексный коэффициент передачи фильтра и передаточная функция H(s) связаны простым соотношением.
(21)
Слайд 51Теперь рассмотрим дискретные фильтры.
Основное разностное уравнение линейного дискретного фильтра имеет
вид.
(22)
Дискретный фильтр описывается также импульсной характеристикой h(n). Импульсная характеристика
определяется из уравнения, которое связывает входящий дискретный сигнал x(n) и выходящий дискретный сигнал y(n).
(23)
Слайд 52Уравнение (23) означает, что выходящий сигнал является сверткой импульсной характеристики
и входящего сигнала.
(24)
Передаточная функция H(z) дискретного фильтра является
Z – образом импульсной характеристики h(n).
(25)
Слайд 53 Если заданы коэффициенты bn , an основного уравнения,
то передаточная функция определятся следующей формулой
(26)
Разложив
числитель и знаменатель передаточной функции (26) на элементарные множители, мы по аналогии с аналоговым фильтром получаем передаточную функцию в следующем виде.
Слайд 54
(27)
Здесь k = b0 - коэффициент усиления, zi –
нули передаточной функции, pi - полюсы передаточной функции.
Частотная характеристика фильтра
является преобразованием Фурье импульсной характеристики
Поэтому можно представить в виде ряда.
(28)
Слайд 55Здесь Δt - шаг дискретизации, который связан с частотой Найквиста
F формулой.
(29)
АЧХ и ФЧХ дискретного фильтра определяются формулами.
(30)
Анализ
формул (25) и (28) показывает, что частотная характеристика фильтра и передаточная функция H(z) связаны простым соотношением.
(31)
Слайд 56 Сравнение формул для аналогового фильтра и дискретного фильтра
показывает, что они во многом близки. Но имеются и различия,
например, частотная характеристика дискретного фильтра, является периодической функцией с периодом 4πF .
В дальнейшем, чтобы отличать характеристики дискретного и аналогового фильтра, будем обозначать аналоговый фильтр значком « a » , а дискретный фильтр значком « d ». Так как частотная характеристика дискретного фильтра обладает свойством периодичности, то это будет обозначаться следующим образом
(32)
Слайд 57 При проектировании дискретного фильтра по аналоговому прототипу необходимо
совершить переход из s - области в z
- область. Другими словами, необходимо преобразовать функцию передачи Ha(s) аналогового фильтра в функцию передачи Hd(s) дискретного фильтра.
Получающийся дискретный фильтр не может быть полностью идентичен аналоговому прототипу по своим характеристикам. Можно говорить только об определенном соответствии характеристик аналогового и дискретного фильтров.
Слайд 58Метод билинейного - преобразования
Данный метод позволяет синтезировать
рекурсивный дискретный фильтр по частотной характеристике аналогового прототипа.
В
данном методе предлагается следующее билинейное преобразование (bilinear transformation).
(33)
Далее будем предполагать, что передаточные функции дискретного и аналогового фильтров связаны соотношением.
Слайд 59где s определяется билинейным преобразованием (33). Отсюда получаем искомую формулу
перехода от аналогового фильтра-прототипа к дискретному фильтру.
(35)
Частотная характеристика дискретного
фильтра выражается через передаточную функцию следующим образом.
(36)
(34)
Слайд 60Объединяя формулы (35) и (36) получаем.
Используя формулу Эйлера, получаем.
(38)
(37)
Подставляя выражение (38) в формулу (37) записываем частотную характеристику дискретного
фильтра в следующем виде.
Слайд 61
(39)
Наконец используем связь между комплексным коэффициентом передачи и передаточной
функцией аналогового фильтра.
(40)
Окончательно, объединяя формулы (39) и (40) находим
связь между частотными характеристиками дискретного и аналогового фильтров.
(41)
Слайд 62На двух рисунках показаны АЧХ аналогового фильтра Чебышева первого рода
с частотой среза ω0 = 1 рад/с , и АЧХ
дискретного фильтра рассчитанная по формуле (41). Частота Найквиста взята равной F = 0.5 Гц .
Слайд 64Увеличим частоту Найквиста в два раза, т.е. положим F =
1 Гц . На следующих рисунках показаны результаты расчета. Видно
что период АЧХ дискретного фильтра возрос в два раза.
Слайд 66Посмотрим на аргумент комплексного коэффициента передачи аналогового фильтра в формуле
(41), который обозначим Ω .
(42)
При малых значения аргумента тангенса,
можно использовать приближенную формулу.
(43)
Поэтому при низких частотах величина Ω будет равна.
(44)
Слайд 67Подставляя (44) в (41) получаем.
(45)
Поэтому в области
низких частот частотные характеристики аналогового и дискретного фильтров почти совпадают.
Далее, по мере ускорения роста функции тангенса, частотная характеристика дискретного фильтра все сильнее сжимается по горизонтали (по сравнению с аналоговым прототипом).
На следующих рисунках, показаны АЧХ аналогового и дискретного фильтров, вблизи полосы пропускания . Красный цвет – АЧХ аналогового фильтра , синий цвет - АЧХ дискретного фильтра. На первом рисунке частота Найквиста равна F = 0.5 Гц, на втором равна F = 0.3 Гц .
Слайд 70На двух следующих рисунках показаны АЧХ аналогового фильтра Чебышева второго
рода с частотой среза ω0 = 1 рад/с , и
АЧХ дискретного фильтра рассчитанная по формуле (41). Частота Найквиста взята равной F = 0.5 Гц .