Добавил:
Студент, если у тебя есть завалявшиеся работы, то не стесняйся, загрузи их на СтудентФайлс! Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

ЦИФРОВЫЕ КИХ-ФИЛЬТРЫ

.pdf
Скачиваний:
6
Добавлен:
18.04.2022
Размер:
359.99 Кб
Скачать

 

+

ωд / 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 ( )|

|Hd ( 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 ( )

с шагом ∆f, в результате чего получаем ДЧХ Hd ( k); k = 0, 1, ..., N − 1. Определяем номера k единичных, нулевых и варьируемых частот-

ных выборок.

Задаем начальные значения Hi.нач оптимизируемых частотных выборок в каждой переходной полосе, например путем линейной интерполяции АЧХ между ее граничными частотами среза и задерживания.

Шаг 4. Рассчитываем частотную характеристику Н( ) и находим значения 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