- •Лекция 13
- •Здесь t - шаг дискретизации, который связан с частотой Найквиста F формулой.
- •Значок « d » относится к дискретному фильтру, а значок « a »
- •Поэтому можно говорить только о приближенном совпадении частотных характеристик дискретного и аналогового фильтров.
- •На втором рисунке частота Найквиста в три раза больше, она
- •Поэтому условие (6) в полосе пропускания выполняется лучше, отсюда неплохое совпадения АЧХ дискретного
- •Однако, хорошее совпадение частотных характеристик дискретного и аналогового фильтров при малых частотах в
- •Красный цвет – АЧХ аналогового фильтра , синий цвет - АЧХ
- •Геометрические свойства билинейного преобразо
- •Для любой точки в левой полуплоскости комплексной
- •В дроби (16) числитель меньше знаменателя, поэтому дробь (16) меньше единицы. Отсюда получаем
- •Другими словами, точки мнимой оси комплексной
- •Далее вспоминаем, что условием устойчивости дискретного фильтра. Дискретный фильтр является устойчивым фильтром, если
- •Метод инвариантной импульсной характерист
- •Схему синтеза дискретного фильтра можно представить в следующем виде.
- •На самом деле, чтобы получить частотную характеристику Kd ( ) дискретного фильтра по
- •Здесь S(f) – спектр непрерывного сигнала, SD( f ) - спектр
- •Вспомним из прошлой лекции, как выражаются частотные характеристики аналогового и дискретного фильтров через
- •Если импульсную характеристику аналогового фильтра рассматривать как непрерывный сигнал, а импульсную характеристику дискретного
- •Частотные характеристики аналоговых фильтров, которые мы рассмотрели выше, были нормированы на единицу. Дальше
- •Формулы (25), (26) являются основными в методе
- •На следующем рисунке частота Найквиста выбрана в два раза большей 48 кГц .
- •Синтез нерекурсивных фильтров с использование
- •Поэтому передаточная функция нерекурсивного фильтра имеет следующий вид.
- •Кроме того, элементы импульсной характеристики h(n) совпадают с коэффициентами bn основного уравнения фильтра.
- •Таким образом, если мы имеем нерекурсивный фильтр порядка N, и нам известна импульсная
- •Увеличение порядка фильтра N означает увеличение электрических элементов в конструкции фильтра. Поэтому N
- •Метод окон является одним их таких методов синтеза нерекурсивных фильтров.
- •На рисунке показан график спектра.
- •Так как мы собираемся рассматривать действительные сигналы (импульсные характеристики), то спектр, как известно,
- •На рисунке показан график расширенного спектра.
- •Подставляем (35) в интеграл (32), интегрируем и получаем элементы дискретного сигнала
- •Для удобства мы положили множитель в (35) равным.
- •На следующем рисунке показаны значения усеченного сигнала.
- •Этот усеченный дискретный сигнал подставляем в ряд Фурье (32), который теперь является конечной
- •На горизонтальной оси выводится частота в единицах
- •Теперь нам надо связать дискретный сигнал с импульсной характеристикой фильтра, а спектр дискретного
- •Тогда частотная характеристика фильтра и спектр дискретного сигнала будут связаны соотношение.
- •Весовые функции, окна
- •Весовая функция прямоугольного окна имеет вид.
- •Это преобразование будем обозначать следующим символом.
- •Импульсная функция фильтра будет находится из нового ограниченного дискретного сигнала.
- •войства некоторых популярных весовых функци
- •Треугольное окно
- •Весовая функция треугольного окна
- •Спектр треугольного окна. Уровень бокового лепестка
- •Применим треугольное окно для синтеза дискретного ФНЧ 32 порядка.
- •Окно Бартлетта
- •w( n ) Bartlett
- •S( f ) Bartlett
- •Применим окно Бартлетта для синтеза дискретного ФНЧ 32
- •Окно Ханна
- •w( n ) Hann
- •S( f ) Hann
- •Применим окно Ханна для синтеза дискретного ФНЧ 32
- •Окно Хэмминга
- •w( n ) Hamming
- •S( f ) Hamming
- •Применим окно Хэмминга для синтеза дискретного ФНЧ 32
- •Окно Блэкмена
- •w( n ) Blackman
- •S( f ) Blackman
- •Применим окно Блэкмена для синтеза дискретного ФНЧ 32
- •Окно Кайзера
- •w( n ) Kaiser
- •S( f ) Kaiser
- •Применим окно Кайзера для синтеза дискретного ФНЧ 32 порядка. Результат показан на следующем
- •Окно Чебышева
- •w( n ) Chebyshev
- •S( f ) Chebyshev
- •Применим окно Чебышева для синтеза дискретного ФНЧ 32 порядка. Результат показан на следующем
На горизонтальной оси выводится частота в единицах |
|
частоты Найквиста. |
41 |
|
Теперь нам надо связать дискретный сигнал с импульсной характеристикой фильтра, а спектр дискретного сигнала с частотной характеристикой фильтра. Для этого надо сравнить формулу (31), которая связывает импульсную характеристику фильтра с частотной характеристикой, и формулу (38), которая связывает дискретный сигнал со спектром дискретного сигнала.
Свяжем дискретный сигнал и импульсную характеристику следующим соотношением.
h(n) s |
N , |
n 0,1, , N (39) |
|
n |
|
|
|
2 |
|
|
|
|
|
|
|
42
Тогда частотная характеристика фильтра и спектр дискретного сигнала будут связаны соотношение.
i N f
K( f ) 2FSD ( f )e 2 F (40)
Таким образом, переход от характеристик дискретного сигнала к характеристикам фильтра осуществляется с помощью формул (39), (40).
43
Весовые функции, окна
При переходе от дискретного сигнала к импульсной характеристики фильтра, мы бесконечную последовательность усекли до конечной последовательности с числом элементов N+1 . По сути дела, мы использовали весовую функцию (weighting function) прямоугольного окна.
Весовая функция это некоторый дискретный сигнал w(n) , в точности равный нулю за пределами заданного интервала.
w(n) 0, n [ M , M ]
(41)
w(n) 0, n [ M , M ]
44
Весовая функция прямоугольного окна имеет вид.
1, |
n [ M , M ] |
(42) |
w(n) |
n [ M , M ] |
|
0, |
|
|
В рассмотренном примере, мы полагали число M равным. |
||
|
M |
N |
|
|
|
|
|
|
||
|
2 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
Весовая функция используется по следующему правилу. |
||||||||||
Дискретный сигнал и его спектр связаны преобразованием |
||||||||||
Фурье. |
|
1 |
|
|
|
i |
n |
f |
||
SD ( f ) |
|
|
e |
|||||||
|
|
|
||||||||
|
|
sn |
|
F |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2F n |
|
|
|
(43) |
||||
F |
|
|
|
|
n |
|
|
|
||
|
|
|
|
|
|
|
|
|||
sn SD ( f ) e i |
|
f d f |
|
|||||||
F |
45 |
|||||||||
F |
|
|
|
|
|
|
|
|
|
|
Это преобразование будем обозначать следующим символом.
sn SD ( f ) |
(44) |
У весовой функции имеется свой Фурье образ.
w(n) W ( f ) |
(45) |
Перемножаем дискретный сигнал и весовую функцию, в результате получаем новый ограниченный дискретный сигнал
~ |
M n M (46) |
sn sn w(n), |
46
Далее мы находим спектр нового сигнала по формулам (43). Как известно~из теории преобразований Фурье, новый полученный спектр SD ( f ) будет равен свертке спектра исходного сигнала SD( f ) и спектра весовой функции W ( f ) .
~
SD ( f ) SD ( f ) W ( f ) (47)
Новый спектр будет, по сути дела, частотной характеристикой фильтра, которая связана с ним соотношением.
i N f
K( f ) 2FSD ( f )e 2 F (48)
47
Импульсная функция фильтра будет находится из нового ограниченного дискретного сигнала.
~ |
N , |
n 0,1, , N |
(49) |
|
h(n) s |
||||
n |
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
Весовых функций существует довольно много. Многие весовые функции дают гораздо лучший результат, чем прямоугольное окно. Это в первую очередь определяется характером спектра весовой функции. При небольшом навыке, по виду спектра окна можно сделать выводы о частотных свойствах синтезируемого фильтра.
48
войства некоторых популярных весовых функци
Прямоугольное окно
Весовая функция прямоугольного окна (42) была рассмотрена выше. На рисунке показан спектр прямоугольного окна.
Уровень первого бокового лепестка составляет – 15.0 дБ. 49
Треугольное окно
Отсчеты треугольного окна рассчитываются по следующей формуле, для нечетный n.
|
2k |
|
|
, |
|
1 k |
n 1 |
, |
|
|
|
|
|
|
2 |
||||
|
|
|
|
||||||
|
|
|
|
|
|
|
(50) |
||
w(k) n 1 |
|
|
|
|
|||||
|
2(n k 1) |
, |
n 1 |
k n |
|
||||
|
|
n 1 |
2 |
|
|||||
Для четных n это будет другая формула.
|
2k 1 |
, |
|
|
1 |
k |
n |
, |
|
|
n |
|
|
|
2 |
||||
|
|
|
|
|
|
(51) |
|||
w(k) |
|
|
|
n |
|
|
|||
|
2(n k 1) |
, |
1 k n |
||||||
|
|
n |
2 |
||||||
50
