ЦИФРОВЫЕ КИХ-ФИЛЬТРЫ
.pdf
|
+ |
ωд / 2 |
ωд / 2 |
– ωд / 2 |
– ωд / 2 |
|
– |
N – четное |
N – нечетное |
Рис. 6.4. Расположение точек дискретизации ЧХ на комплексной плоскости
Шаг дискретизации |
по частоте ∆ω выбирается из |
условия |
||
∆ω ≤ ∆ωпер / (L + 1), где L − целые числа, L = 0, 1, 2, ...; ∆ωпер − |
переход- |
|||
ная полоса фильтра. |
|
|
|
|
В результате получается дискретизированная частотная характе- |
||||
ристика фильтра (ДЧХ) Hd |
( j ×ωk ) = Hd ( j ×ω) |
|
(рис. 6.5). Так как за- |
|
|
||||
|
|
|
ω =ωk |
|
|
|
|
|
данная частотная характеристика соответствует физически нереализуемому фильтру с нулевым запаздыванием, то для ЦФ со ступенчатообразными АЧХ дискретизированная частотная характеристика отождествляется далее с их дискретизированной АЧХ.
Дискретизация частотной характеристики на рис. 6.5 выполнена с шагом ∆ω = ∆ωпер / 2(L = 1).
|Hd ( jω)|
|Hd ( jωk)|
H1 H1
0 |
2 |
/ 2 |
N – 1 |
|
|
(N – 1) |
|
||
∆ω |
∆ωпер ωд / 2 |
|
ωд |
ω |
Рис. 6.5. Дискретизированная ЧХ цифрового фильтра нижних частот
ДЧХ имеет значения, равные в полосе пропускания 1 (Hd ( j × ωk) = 1), в полосе задерживания − нулю (Hd ( j × ωk) = 0) и в переходной полосе – некоторым промежуточным варьируемым (оптимизируемым) значениям
118
Hd ( j × ωk) = H1 = var, от которых зависит качество аппроксимации заданной частотной характеристики.
ДЧХ Hd ( j × ωk) можно поставить в соответствие некоторую импульсную характеристику hp(m), определяемую с помощью обратного дискретного преобразования Фурье (ОДПФ), получаемого путем дискретизации по частоте общего выражения для импульсной характеристики hd (m). Она соответствует заданной (непрерывной) частотной характеристике Hd (j × ω):
|
|
|
Tд |
|
ωд |
H ( j ×ω) ×e j×ω×m×Tдdω . |
|
||||
|
|
h (m) = |
× |
∫ |
|
||||||
|
|
|
|
||||||||
|
|
d |
2 ×π |
|
d |
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N -1 |
|
Выполняя |
замены: |
ω → ωk ; dω → ω = ωд / N ; ∫→ ∑, получим |
|||||||||
|
|
|
|
|
|
|
|
|
|
k =0 |
|
импульсную характеристику hp(m): |
|
|
|
|
|||||||
|
1 |
N -1 |
|
|
|
|
|
1 |
N -1 |
|
|
hp (m) = |
∑ Hd ( jωk )e j×ωk ×m×Tд = |
∑ Hd ( jωk )e j×ωk |
×(m+i×N )×Tд , |
||||||||
N |
N |
||||||||||
|
k =0 |
|
|
|
|
|
k =0 |
|
|||
|
|
|
|
|
|
|
|
|
где i = 0, ±1, ±2, ±...
Индекс p означает, что эта импульсная характеристика является периодической с периодом Np = N, т. к. дискретизации в частотной области соответствует периодизация во временной области (рис. 6.6).
Рис. 6.6. Импульсная характеристика, соответствующая ДЧХ
В качестве импульсной характеристики синтезируемого методом частотной выборки НРЦФ выбирается один период импульсной характеристики hp(m), сдвинутый вправо на (N − 1) / 2 отсчетов (для обеспечения физической реализуемости) и усеченный прямоугольной весовой функцией (для получения КИХ-фильтра) (рис. 6.7):
h(m) = hp (m − N −1), m = 0,1, ..., N −1. 2
Рис. 6.7. Импульсная характеристика НРЦФ, синтезированного методом частотной выборки
119
По импульсной характеристике h(m) находится частотная характеристика фильтра H( j × ω), аппроксимирующая заданную:
N -1 |
|
|
|
|
|
|
|
|
N -1 |
|
|
|
|
|
N -1 |
|
|
N -1 |
|
|
|||
H ( jω) = ∑ h(m) ×e- j×ω×m×Tд |
= |
|
1 |
∑ Hd ( jωk )∑ e j×ωk m |
Tд e- j×ω×m×Tд = |
||||||||||||||||||
|
2 |
||||||||||||||||||||||
N |
|||||||||||||||||||||||
m=0 |
|
|
|
|
|
|
|
m=0 |
|
|
|
|
|
m=0 |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
1 |
N -1 |
|
- j×ω × |
N -1 |
×T N -1 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
= |
∑ Hd ( jωk ) ×e |
|
|
k |
|
д |
∑ e- j×(ω-ωk )×m×Tд = |
||||||||||||||||
|
|
2 |
|
|
|||||||||||||||||||
|
|
|
|
||||||||||||||||||||
|
N m=0 |
|
|
|
|
|
|
|
|
m=0 |
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
N -1 |
|
|
|
|
|
|
sin |
(ω − ωk ) |
× N ×Tд |
||||||||
|
|
|
N -1 |
|
|
|
|
|
|
|
|
||||||||||||
= e- j×ω×( |
)×Tд × |
1 |
∑ Hd ( jωk ) × |
|
|
|
2 |
|
|
|
|
|
|||||||||||
2 |
|
|
|
|
|||||||||||||||||||
|
|
|
|
(ω -ωk ) |
|
|
|
||||||||||||||||
|
|
|
|
|
N m=0 |
|
|
|
|
|
|
|
sin |
×T |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
д |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(при выводе использовано выражение для суммы конечного числа членов геометрической прогрессии).
В этом выражении множитель e- j×ω×( N2−1)×Tд определяет ФЧХ филь- тра: ϕ (ω) = -ω × N2−1 ×Tд , которая строго линейна вследствие симметрии
импульсной характеристики.
АЧХ фильтра на частотах ω = ωk: H(ωk) = Hd (ωk) точно совпадает с частотными выборками ДЧХ, а на частотах ω ≠ ωk H(ω) ≠ Hd (ω) − отличается от заданной на величину погрешности аппроксимации.
Качество аппроксимации в данном методе зависит от числа выборок частотной характеристики в переходной полосе L и их зна-
чений Hi.опт (i = 1, 2, ..., L), делающих аппроксимируемую функцию более гладкой.
Различным значениям L соответствуют следующие примерные значения максимального уровня боковых лепестков:
L = 0: δ2мах ≈ −20 |
дБ; |
|
L = 1: δ2мах ≈ −40 |
дБ; |
|
L = 2: δ2мах ≈ −50−60 |
дБ; |
|
L = 3: δ2мах ≈ −80−100 |
дБ. |
Реально методом частотной выборки можно синтезировать НРЦФ с минимальным затуханием в полосе задерживания до (90−120) дБ.
Таким образом, оптимизация фильтра заключается в выборе L − числа выборок в переходной полосе, и поиске их оптимальных значений Hi.опт, минимизирующих погрешности аппроксимации. Очевидно, что с увеличением числа варьируемых выборок существенно усложняется процедура оптимизации. Она достаточно эффективно реализуется на ЭВМ методом линейного программирования.
6.2.2. Методика синтеза НРЦФ методом частотной выборки
Шаг 1. По значению заданного затухания в полосе задерживания Аз выбирается число варьируемых отсчетов L частотной характеристики в переходной полосе. Например, при Аз ≤ 40 дБ L = 1.
120
Чем сложнее АЧХ фильтра, тем меньше затухание при данном значении L.
Шаг 2. Для принятого значения L и заданной переходной полосы
Dfпер = |
|
fз - fс |
|
находим шаг дискретизации частотной характеристики по |
|||||||
|
|
||||||||||
частоте Df = |
|
|
fпер |
и число точек дискретизации N = |
fд |
= (L +1)× |
fд |
. |
|||
|
L +1 |
f |
f |
|
|||||||
|
|
|
|
|
|
|
D |
|
D |
пер |
|
Приравниваем N к ближайшему целому числу, обычно нечетному. Шаг 3. Дискретизируем заданную частотную характеристику Hd ( jω)
с шагом ∆f, в результате чего получаем ДЧХ Hd ( jωk); k = 0, 1, ..., N − 1. Определяем номера k единичных, нулевых и варьируемых частот-
ных выборок.
Задаем начальные значения Hi.нач оптимизируемых частотных выборок в каждой переходной полосе, например путем линейной интерполяции АЧХ между ее граничными частотами среза и задерживания.
Шаг 4. Рассчитываем частотную характеристику Н( jω) и находим значения Hi.опт, при которых частотная характеристика удовлетворяет заданным требованиям.
Например, для ФНЧ
при L = 1, N = 33 значение H1опт = 0,3904, δ2max = −40 |
дБ; |
|
|
||||||||||
при L = 2, N = 65 H1опт = 0,588, H2опт = 0,1065, δ2max < −60 |
|
дБ. |
|||||||||||
Шаг 5. Рассчитываем импульсную характеристику НРЦФ с уче- |
|||||||||||||
том симметрии частотной характеристики: |
|
|
|
|
|||||||||
|
Hd (0) |
|
1 |
KB |
|
|
×cos (n |
- N -1)×ωk |
×Tд |
|
|
||
h(n) = |
+ |
∑2 × |
Hd ( jωk ) |
, |
(6.22) |
||||||||
|
|
||||||||||||
|
N |
|
|
N k =1 |
|
|
|
2 |
|
|
|
||
n = 0, 1, 2, …, N − 1; |
KB = (N − 1) / 2 |
при нечетном N, и KB = (N / 2) − 1 – |
|||||||||||
при четном. |
|
|
|
|
|
|
|
|
|
|
|
Шаг 6. Выбираем способ реализации НРЦФ (ДВС или ДПФ).
6.3.Контрольные вопросы
1.Методы проектирования КИХ- и БИХ-фильтров (классификация).
2.Принципы расчета КИХ-фильтров методом оконных функций.
3.Свойства прямоугольной весовой функции: АЧХ, импульсная характеристика.
4.Характеристики идеальных ЦФ.
5.Принципы расчета КИХ-фильтров методом частотной выборки.
121