
Лекции по ТПС
.pdfДля получения передаточной функции рекурсивного фильтра принято использовать Z-преобразование:
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X (z) xn z n . |
|
|
|||
|
|
|
|
|
n |
|
|
|
|
Умножая y(n) на b0 |
1 и подвергая обе части уравнения (33.1) Z- |
||||||||
|
|
|
преобразованию, получим |
|
|||||
|
|
|
M |
|
|
|
|
N |
|
|
|
|
bk y(n k ) z n |
|
|
|
ak x(n k ) z n . |
||
n k 0 |
|
n |
k 0 |
|
|||||
Поскольку в уравнении (33.1) принято |
|
|
|
|
|||||
|
|
|
ak = 0 |
при k<0 и k>N, |
|
|
|
|
|
|
|
|
|
bk = 0 |
при k<0 и k>M, |
|
|||
можно расширить границы суммирования до : |
|
||||||||
|
|
|
|
|
|
|
|
||
|
bk y(n k ) z n |
ak x(n k ) z n |
|||||||
n k |
|
n k |
|
||||||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
bk z k y(n k)z (n k ) |
ak z k x(n k)z |
(n k ) |
|||||||
k |
n |
|
k |
|
|
n |
|
||
Обозначив m = (n - k), получим |
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
bk z k |
ym z m an z k xm z m |
||||||
|
|
k |
m |
k |
|
m |
|
или в компактном виде
B(z) Y(z) = A(z) X(z),
где A(z), B(z), X(z), Y(z) - Z-преобразования соответствующих числовых последовательностей.
Отсюда следует, что Z-преобразование передаточной функции фильтра (т.е. отношения выходной реакции ко входному воздействию) имеет вид:
271

H д ( z) |
Y ( z) |
|
A( z) |
|
a a z 1 a z 2 |
a z N |
|||
|
|
0 |
1 |
2 |
|
N |
|||
|
|
|
|
|
|
||||
|
X ( z) |
|
B( z) |
|
b b z 1 |
b z 2 |
b z M |
||
|
|
|
|
|
0 |
1 |
2 |
M |
|
(33.2)
После подстановки в (33.2) z exp( j Tд ) получим передаточную функцию в явном виде, т.е. в виде зависимости коэффициента передачи от частоты:
H д ( j ) H д (z exp( j Tд )) .
Коэффициент |
передачи - периодическая функция частоты с периодом |
2 Tд . |
|
В частном случае нерекурсивного фильтра |
|
|
0 |
B(z) bn z n bn z n b0 1 |
|
n |
n 0 |
и
H д ( j ) A(z exp( j Tд )) ak exp( j nTд ) .
n
Расчет рекурсивных цифровых фильтров. Расчет (проектирование) ре-
курсивных фильтров существенно сложнее расчета нерекурсивных фильтров. Существует большое количество различных методик, основанных на методиках проектирования аналоговых фильтров и применении математического анализа, связанного с преобразованием Лапласа, Z-преобразованием, теории вычетов.
Рассмотрим наиболее простую и популярную методику, называемую «ча-
стотное преобразование».
Сущность этой методики состоит в трансформации передаточной характеристики H д (z* ) некоего ФНЧ, именуемого «ФНЧ-прототип», в передаточную характеристику H д (z) нужного фильтра (НЧ, ВЧ, полосового, режекторного), с последующей заменой z exp( j Tд ) .
Представим такую методику в виде последовательности “шагов”:
Шаг 1: Исходя из требований к качеству проектируемого фильтра (гладкость в полосе пропускания, фазовая характеристика, допустимый уровень пульсаций в полосах пропускания и запирания), выбирают тип фильтра и порядок фильтра N. При этом следует знать:
272

ФНЧ Баттерворта обеспечивает максимально плоскую характеристику в зоне пропускания, однако плохо передает фронты прямоугольных импульсов. Крутость переходной зоны (из полосы пропускания в полосу запирания) растет с увеличением N.
ФНЧ Бесселя обеспечивает линейность фазовой характеристики в зоне пропускания, что хорошо для фильтрации прямоугольных импульсов (малые выбросы на фронтах отфильтрованных импульсов), крутость переходной зоны меньшая, чем у фильтров Баттерворта.
ФНЧ Чебышева обеспечивает фиксированную пульсацию в зоне пропускания (фильтры типа I) или в зоне запирания (фильтры типа II). Плохо передает фронты прямоугольных импульсов. Крутость переходной зоны большая, чем у фильтров Баттерворта и Бесселя.
ФНЧ Кауэра имеет пульсацию как в зоне пропускания, так и в зоне запирания. Плохо передает фронты прямоугольных импульсов. Крутость переходной зоны большая, чем у фильтров Баттерворта и Бесселя.
Шаг 2. После выбора типа фильтра обращаются к одной из таблиц филь- тров-прототипов, из которых выбирают числовые значения коэффициентов ak и
bk . К таким таблицам «привязан» также числовой параметр fд* fc* - отношение частоты дискретизации к частоте среза фильтра-прототипа. В [1] предлагаются таблицы для fд*
fc* 4 , имеющие вид:
Таблица 33.1. Коэффициенты для прототипов ФНЧ Баттерворта N-го порядка
a0i |
a1i |
a2i |
b1i |
b2i |
0,5 |
0,5 |
0 |
0 |
0 |
0,2929 |
0,5858 |
0,2929 |
0 |
0,1716 |
0,5 |
0,5 |
0 |
0 |
0 |
0,3333 |
0,6667 |
0,3333 |
0 |
0,3333 |
|
|
|
|
|
Таблица 33.2. Коэффициенты для прототипов ФНЧ Бесселя N-го порядка
273
a0i |
a1i |
a2i |
b1i |
b2i |
0,5 |
0,5 |
0 |
0 |
0 |
0,335 6 |
0,671 2 |
0,3356 |
0,2564 |
0,0806 |
0,5695 |
0,5695 |
0 |
0,1390 |
0 |
0,4037 |
0,8075 |
0,4037 |
0,4222 |
0,1928 |
|
|
|
|
|
Таблица 33.3. Коэффициенты для прототипов ФНЧ Чебышева I-го типа N-го порядка с пульсацией 11% (1дБ)
a0i |
a1i |
a2i |
b1i |
b2i |
0,5 |
0,5 |
0 |
0 |
0 |
0,2595 |
0,5190 |
0,2595 |
-0,2862 |
0,3242 |
0,3110 |
0,3110 |
0 |
-0,3780 |
0 |
0,3636 |
0,7273 |
0,3636 |
-0,1496 |
0,6042 |
|
|
|
|
|
Таблица 33.4. Коэффициенты для прототипов ФНЧ Кауэра N-го порядка
с пульсацией 15% (1,5дБ)
a0i |
a1i |
a2i |
b1i |
b2i |
0,3306 |
0,3306 |
0 |
-0,3389 |
0 |
0,4898 |
0,6612 |
0,4898 |
-0,0219 |
0,6627 |
0,1851 |
0,2397 |
0,1851 |
-0,5678 |
0,1778 |
0,4347 |
0,8693 |
0,4347 |
-0,0215 |
0,7602 |
|
|
|
|
|
Рассматривая таблицы, нетрудно видеть что коэффициенты зависят также от некоего индекса i - номера каскада. Дело в том, что рекурсивные фильтры порядка N>2 принято получать путем последовательного соединения фильтров
274

2-го порядка. Поэтому если передаточная функция фильтра-прототипа 2-го порядка имеет вид:
H д (z* ) |
a* a* z* 1 |
a* z* 2 |
|
||
0 |
1 |
2 |
, |
||
b* b* z* 1 |
b* z* 2 |
||||
|
|
||||
|
0 |
1 |
2 |
|
то передаточная функция фильтра-прототипа более высокого порядка получается как произведение передаточных характеристик P частных фильтров:
P |
a* |
a* z* 1 |
a* |
z* 2 |
|
||
H д (z* ) |
0i |
1i |
2i |
|
|
. |
|
b* |
b* z* 1 |
b* |
z* 2 |
||||
i 1 |
|
||||||
0i |
1i |
2i |
|
|
|
||
Шаг 3. Исходя из требуемого соотношения |
fс |
fд |
|
(при расчете НЧ или ВЧ |
|||
фильтра) либо требуемой пары соотношений fс1 |
fд |
и |
fc2 fд , fc1 fc2 (при рас- |
чете полосового или режекторного фильтра), заменяют переменную z* на функцию переменной z , для чего используют соотношения из табл 33.5.
Таблица 33.5
ФНЧ ФНЧ |
|
|
|
|
z* 1 |
|
|
|
z 1 |
d |
, |
|
|
d |
sin( fс* |
f д* fс |
|
|
f д |
) |
|
|
|
|
|
||||||||||
|
|
|
|
|
1 dz 1 |
|
|
sin( fс* |
f д* fс |
|
|
f д |
) |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
ФНЧ ФВЧ |
|
|
|
|
z* 1 |
|
z 1 |
d |
, |
|
d cos( fс* |
fд* fс |
|
f д |
) |
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
1 dz 1 |
|
|
|
|
|
|
|
sin( fс* |
f д* fс |
|
fд |
) |
|
|
|
|
||||||||||
ФНЧ поло- |
|
|
|
|
z 2 |
|
2dk z 1 |
k 1 |
|
|
|
|
|
cos( f |
|
f |
|
f |
|
fд ) |
|||||||||||||||
совой фильтр |
z |
* 1 |
|
|
|
|
|
k 1 |
|
|
|
k |
1 |
|
, |
|
d |
с1 |
д |
с2 |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
k 1 |
|
2 |
|
2dk |
|
|
1 |
1 |
|
cos( f |
|
f |
|
f |
|
f ) |
|||||||||||||||||
|
|
|
|
|
k 1 z |
|
|
|
k 1 z |
|
|
|
|
|
|
|
с2 |
|
д |
|
|
|
|
с1 |
д |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
k ctg( f |
с2 |
f |
д |
f |
с1 |
f |
д |
)tg f * |
f * |
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
с |
|
д |
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
ФНЧ ре- |
|
|
|
|
z 2 |
|
2d |
z 1 |
|
1 k |
|
|
|
|
|
cos( f |
|
|
|
f |
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f |
|
|
fд ) |
|
||||||||||||||
жекторный фильтр |
z |
* 1 |
|
|
|
|
|
k 1 |
|
|
|
k |
1 |
|
, |
|
d |
с1 |
д |
с2 |
|||||||||||||||
|
|
1 k |
|
2 |
|
2d |
|
|
|
1 |
|
|
|
|
cos( fс2 |
f |
д |
fс1 |
fд ) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
k 1 z |
k 1 z |
1 |
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
k tg( fс2 |
fд fс1 |
fд )tg fс* fд* |
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
275
Шаг 4. Полученное аналитическое выражение передаточной характеристики H д (z) упрощают так, чтобы в числителе и знаменателе оказались полиномы переменной z 1 .
Шаг 5. Числитель и знаменатель функции H д (z) делят на такое число, чтобы выполнялось условие b0 1. Результирующие коэффициенты в числителе и знаменателе (после такой нормировки) и образуют искомые множества коэффи-
циентов ak |
и bk . |
|
|
|
|
|
Шаг |
6. В |
аналитическом |
выражении |
для H д (z) производят замену |
||
z exp( j Tд ) , |
получая, |
таким |
образом, |
частотную |
характеристику |
|
H д ( j ) H д (z exp( j Tд )) |
синтезированного фильтра. АЧХ |
фильтра получают |
как корень квадратный из суммы квадратов действительной и мнимой частей комплексной функции H д ( j ) .
В случае, когда АЧХ не устраивает пользователя, производят пересчет коэффициентов для другого типа фильтра либо увеличивают порядок фильтра.
КОНТРОЛЬНЫЕ ВОПРОСЫ:
1.Какой импульсной характеристикой обладают рекурсивные цифровые Фильтры (РЦФ)?
2.Какое преобразование используют для получения передаточной функции РЦФ?
3.В чем состоит сущность методики расчета РЦФ, называемой «частотное преобразование»?
4.Какими особенностями обладают РЦФ нижних частот Баттерворта, Чебышева, Кауэра, Бесселя?
5.Как получить рекурсивные фильтры порядка N>2 из фильтров 2-го по-
рядка?
276

РАЗДЕЛ 5. ПОМЕХОУСТОЙЧИВОСТЬ СИСТЕМ ПЕРЕДАЧИ
ДИСКРЕТНЫХ И НЕПРЕРЫВНЫХ СООБЩЕНИЙ
ЛЕКЦИЯ 34. ПОМЕХОУСТОЙЧИВОСТЬ РАЗЛИЧНЫХ ВИДОВ МОДУЛЯЦИИ ПРИ ПЕРЕДАЧЕ НЕПРЕРЫВНЫХ СООБЩЕНИЙ
С целью сравнения различных видов модуляции представляется естественным в качестве критерия выбрать величину вида:
,
где - средняя мощность сообщения,
- средние мощности сигнала и по-
мехи (в полосе сигнала) на входе приёмника. В случае белого гауссовского шума на входе приемника имеем:
.
Величина называется выигрыш модуляции. Однако использование
не всегда оказывается удобным. Например, рассматриваются две системы передачи
информации (два метода модуляции) с полосами и
соответственно. При
одинаковых ,
,
и
, т.е. при полной эквивалентности систем, имеем:
;
так что , если
.
В этом смысле более удобно использование обобщенного выигрыша модуляции вида:
(34.1)
где – максимальная частота полосы частот, занимаемой сигналом.
277

Рассмотрим физический смысл величины . При непосредственной передаче сообщения (без модуляции) имеем:
.
При использовании модуляции:
.
Приравняем в обоих случаях:
,
откуда:
.
Таким образом, величина показывает, во сколько раз можно снизить среднюю мощность сигнала на входе приемника при использовании данного метода модуляции в сравнении с непосредственной передачей сообщения при од-
ном и том же .
Очевидно, что если исходить из одинаковых значений и
, то ве-
личина показывает, во сколько раз возрастает величина
при использовании данного метода модуляции в сравнении со случаем непосредственной
передачи сообщения. Наконец, отношение обобщенных выигрышей для двух различных методов модуляции показывает, во сколько раз можно снизить (повысить) среднюю мощность полезного сигнала при переходе от одного
метода модуляции к другому при сохранении одной и той же величины . Таким образом, величина
характеризует энергетический выигрыш (проигрыш) одного метода модуляции по отношению к другому.
278

Представим в (34.1) в форме:
ПОМЕХОУСТОЙЧИВОСТЬ АНАЛОГОВЫХ МЕТОДОВ МОДУЛЯЦИИ
Амплитудная модуляция
;
.
Если |
то |
. |
Балансная модуляция
Однополосная модуляция
где - преобразование по Гильберту функции
Таким образом, во всех рассмотренных случаях модуляция не дает никакого выигрыша, либо (в случае АМ) приводит к энергетическим потерям вследствие бесполезной траты мощности на передачу несущей частоты. Следует также заметить, что в рассмотренных случаях полученные выражения являются
точными и справедливы при любых значениях .
Частотная модуляция
279

; ,
где - индекс ЧМ.
Это выражение справедливо в области , причем
может принимать значения
.
Однако увеличение связано с расширением полосы частот
, поскольку для ЧМ:
.
При этом, с увеличением , увеличивается и
, так что, начиная с некоторого значения
, нарушается условие, при котором справедливо полу-
ченное выражение . При этом происходит резкое уменьшение величины
(пороговый эффект).
Введем понятие пороговой мощности сигнала . Если
, то наступает пороговый эффект. Очевидно,
.
Наличие порогового эффекта при ЧМ иллюстрируется рис 34.1.
Рис.34.1
280