УДК 621.372
DOI: 10.17586/0021-3454-2017-60-7-641-646
АНАЛИЗ КОМПЛЕКСНЫХ ФИЛЬТРОВ НА ОСНОВЕ ПЕРЕХОДНЫХ ХАРАКТЕРИСТИК
С. И. Зиатдинов, Ю. В. Соколова
Санкт-Петербургский государственный университет аэрокосмического приборостроения,
190000, Санкт-Петербург, Россия E-mail: [email protected]
Разработана методика исследования комплексных фильтров с использованием переходных характеристик. Путем представления комплексного входного сигнала в виде суммы сдвинутых относительно друг друга ступенчатых воздействий получена новая форма комплексного интеграла наложения, позволяющая при заданном комплексном входном сигнале и известной комплексной переходной характеристике фильтра найти его выходной сигнал. Показано, что использование комплексных фильтров позволяет легко менять их частоту настройки без изменения формы амплитудно- и фазочастотных характеристик. Введены понятия комплексной дельта-функции и комплексной переходной характеристики, используемые при исследовании динамических свойств фильтров и анализе передачи через них разнообразных сигналов. В отличие от известного метода анализа линейных систем на основе импульсной характеристики предложенная методика позволяет синтезировать комплексные дискретные фильтры, динамические свойства которых совпадают с динамическими свойствами непрерывных фильтров—аналогов. Выдвинутые теоретические положения подтверждены конкретными примерами.
Ключевые слова: комплексный сигнал, комплексный фильтр, комплексная импульсная характеристика, комплексная переходная характеристика, интеграл наложения
Введение. При обработке информации наиболее часто используется процедура фильтрации сигналов с помощью фильтров нижних и верхних частот, полосовых и режекторных фильтров. В ряде случаев применяются адаптивные и комплексные фильтры, параметры которых могут изменяться в зависимости от спектральных характеристик обрабатываемых сигналов. В комплексных фильтрах управлять средней частотой передаточной функции фильтра значительно проще, чем в действительных фильтрах.
В настоящее время теория исследования прохождения сигнала через разнообразные фильтры достаточно хорошо разработана [1—11].
Исчерпывающей характеристикой фильтра является его частотно-передаточная функция. Для анализа фильтров используются два основных метода — частотный и временной. В случае частотного метода исследования линейной системы с частотной передаточной функцией W(Ош) комплексный выходной сигнал 2вых (t) при известной спектральной плотности z^ (jro) комплексного входного сигнала z^ (t) определяется на основании обратного преобразования Фурье [1, 2]:
1 !Х>
(t) = — J WО)Zвх (»Лш.
7
^вых х ' 0
При использовании временного метода анализа комплексных фильтров с импульсной
*
характеристикой И (V) применяется интеграл наложения (интеграл Дюамеля) [2, 3]:
ю *
7вых (0 = I 7вх (* _ т)И ((1)
0
где комплексные входной гвх (*) и выходной гвых (*) сигналы, а также комплексная импульс*
ная характеристика И (*) определяются как
^вх (*) = х« (*)+Жх (*); 2вых (*) -*вых (*) + 7Увых (О, И (г) = Их(*) + И (*). (2)
В соотношениях (2) хвх (V) и увх (V) — квадратурные (комплексно-сопряженные) составляющие входного сигнала; хвых (V) и Увых ( *) — квадратурные составляющие выходного сигнала; Их ( *) и Иу( *) — квадратурные составляющие импульсной характеристики комплексного фильтра.
В работах [4—6] показано, что при переходе от непрерывных фильтров к дискретным существующая методика расчетов коэффициентов разностных уравнений на основе отсчетов импульсных характеристик не обеспечивает правильной реализации переходных процессов. При этом методика расчетов коэффициентов с использованием отсчетов переходных характеристик позволяет получить динамические свойства дискретных фильтров, совпадающие с динамическими свойствами непрерывных фильтров как в переходном, так и в установившемся режиме.
В этих условиях возникает актуальная задача разработки методики анализа фильтров на основе их переходных характеристик. Решение данной задачи составляет основное содержание настоящей статьи.
Интеграл наложения с использованием комплексных переходных характеристик.
Для решения поставленной задачи составляющие хвх (V) и увх (V) комплексного входного
сигнала 2вх (V) представим в виде суммы последовательности ступенчатых функций, начинающихся в моменты времени ^, ^ , ..., *п (на рисунке показан сигнал хвх (V) ):
хвх (0 =[ хвх (*" Ч) - хвх (*" *0)] +[ хвх (* " Н) - хвх (* " *1)] + ..
+[ хвх (* - Ц ) - хвх (* - ь-1)] +... = АхзХ (VI) + Лхвх +... + Лхвх () +...;
Увх(*) = [Увх(* - Ч) - Увх(* - *0)] + [Увх(* - *2) -Увх(* - *1)] +... +[Увх (* - Ц) - Увх (* - Ц-1)] + ... = АУвх (Ь ) + АУвх (¿2) + ... + АУвх (*,■) +..., где Ахвх (* г ) = хвх (* - *г) - хвх (* - *г -1), АУвх (* г) = Увх (* - *г) - Увх (* - *г-1) — элементарные ступенчатые функции, возникающие в момент времени * = .
При этом комплексный входной сигнал определяется выражением
да да
2вх (*) = Е Ахвх (*,,) + ] ЕАУвх (*,). ,=1 ,=1
В результате комплексный выходной сигнал фильтра будет определяться соотношением
7вых ()
го го
ЕЛХзх (и) + ] ЕАУвх (и)
г =1 г =1
^ (V _ ),
(3)
где g (V) — комплексная переходная характеристика, являющаяся реакцией фильтра на комплексное ступенчатое входное воздействие Л7вх (V) = Лх^ (V) + у'Аувх (V).
Переходная характеристика комплексного фильтра находится из соотношения
g*(г) = | И*( х)^х.
С учетом того, что И (V) = Их (V) + ]Иу (V), комплексная переходная характеристика g (V)
*
записывается через квадратурные составляющие g (V) = gx (V) + jgy (V), где
г г
gx (0 = 1 Их (т)ат gy (0 = 1 Иу (тМт.
00
В результате выражение (3) приобретает следующий вид:
го
го
7вых (О = Е ЛХвХ (^ ^ (V _ ^ ) _ Е Лувх (V,- ^у (V _ ^ ) +
г=1 г=1
+j
го
го
(4)
Е Лхвх (^ ^у (V _ V,-) + Е Лувх (и )gх (V _ ^)
_г=1 г=1
Пусть отсчеты квадратурных составляющих входного сигнала хвх (V,) и увх (V,) при
формировании ступенчатых функций берутся через одинаковые отрезки времени, т. е. = гТ,
где Т — период следования отсчетов.
Помножив и поделив выражение (4) на Т , получим
М) = Т
ЕЛхвх (Ц )gx (V _ 4 )/Т _ Е Лувх (4 ^у (V _ ^)/Т
г=1
г=1
+jT
го го
Е ЛхВх (^ )gy (V _ ц)/ Т + Е ЛуВх (^ ^х (V _ ц)/ Т
. г=1 г=1
В предельном переходе при Т ^ 0 можно записать:
7вых (V) = I dxвв^gx (V _ &Т _ I ^^у (V _ т)^т +
0 ат о ат
+j
I (V _ т)а т+1 (V _ т)а т
0 ат о ат
(5)
Полученное выражение (5) для выходного сигнала комплексного фильтра является комплексным интегралом наложения, связывающим производные квадратурных составляющих входного сигнала и переходной характеристики фильтра.
Известно [2, 7], что импульсная характеристика линейной системы является производ-
, * . . dg (V) ^
ной от переходной характеристики И (V) =-. Тогда выражение для выходного сигнала
фильтра принимает вид
' аг (т)
7вых (') = I 7вх (' _ ат.
0 ат
Таким образом, два равноценных преобразования (5) и (6) позволяют найти выходной
сигнал комплексного фильтра при известных комплексных входном сигнале zвх (t) и пере*
ходной характеристике g (t).
Для проверки правильности полученных результатов рассмотрим конкретные примеры. Пример 1. Пусть входной сигнал фильтра является комплексным единичным ступенчатым воздействием z^ (t) = 1(t) (cos ф + j sin ф) , где угол ф может принимать любые значения и, в частности, при ф=0 входной сигнал становится вещественным.
Для данного сигнала ^вх (t) = g [.] = g[] (cos ф + j sin ф) — комплексная дельта-функция,
dt
после подстановки которой в выражение (5) получим
г , *
2вых (*) = I 5 № (* - т)Л = % (*).
0
В результате при комплексном единичном ступенчатом входном воздействии выходной сигнал комплексного фильтра представляет его комплексную переходную характеристику.
Пример 2. Рассмотрим комплексный полосовой фильтр первого порядка с частотной передаточной функцией вида [9]
Ж (=-1-,
1 + ](Ю - Ю0)/Юср
где юср — частота среза фильтра, Ю0 — средняя частота (частота настройки фильтра).
Пусть на входе фильтра действует комплексное гармоническое колебание 2вх (*) = [^(ю*) + ] ^п(ю*)] с частотой ю.
Для рассматриваемого фильтра комплексные импульсная и переходная характеристики записываются следующим образом [2, 10]:
* —ю *
И (*) = юсре ср [cos(ю0t) + ] sin(ю0t)];
* Юср / ( -ю Г )
% (*) = 2 2 (\е ср [ю0^ЧЮ00 - юср + юср [-
ю0 + ю^и
|е"ЮсрГ [Юср ^(ю,/) + Ю0 ^(ю0*)] + Ю0 ^. (7)
игнала определяется вы
= ю[cos(юt) - j'sm(юt)]. dt
После подстановки производной входного сигнала в (5) с учетом выражения (7) получим комплексный выходной сигнал рассматриваемого фильтра:
* юс
При этом производная входного сигнала определяется выражением
^вх (*)
W(t) = J ro[cos(rox) - jsm(roT)] 2 ср 2 х
ю0 + Юср
х е
L
|е ЮсрТ) [ю0^п (ю0 (* - т))"юср cos (ю0 (* - Т))] + юср | - j |е"юср -т) ^юср ^ (ю0(* - Т)) + ю0 cos ((ю0 (* - Т) )] + ю0 ^
Интеграл в данной формуле является табличным [11]. В окончательном виде для установившегося режима можно записать следующие соотношения для комплексного выходного сигнала полосового фильтра:
x(t) = ■
1
(
г cos
1+
Гю -
ю - ю.
- гюЛ
юt - arctg-
ю - ю,
'0
ю
0
ю
ср У
y(t) = ■
(
г sin
1+
С ю - ю0 ^
ср У
-ю
юt - arctg-
ю - ю,
0
ю
ю
сР У
ср У
Данные выражения строго определяют выходной сигнал комплексного полосового фильтра первого порядка при гармоническом входном сигнале.
Полученные в работе соотношения (5) и (6) для комплексного интеграла наложения дополняют известную методику анализа линейных систем.
Заключение. В результате проведенных исследований получила дополнительное развитие теория линейных систем. Представленная форма интеграла наложения позволяет при заданной производной комплексного входного сигнала и известной комплексной переходной характеристике фильтра либо при заданном комплексном входном сигнале и известной производной комплексной переходной характеристики найти выходной сигнал комплексного фильтра. Разработанная методика может быть использована при синтезе по заданной переходной характеристике комплексных дискретных линейных систем, динамические свойства которых совпадают с динамическими свойствами непрерывных линейных систем—аналогов.
список литературы
1. Воробьев С. Н. Цифровая обработка сигналов. СПб: Изд. дом „Академия", 2013. 318 с.
2. Гоноровский И. С. Радиотехнические цепи и сигналы. М.: Радио и связь, 1986. 512 с.
3. Оппенгейм А., ШафферР. Цифровая обработка сигналов. М.: Техносфера, 2006. 855 с.
4. Зиатдинов С. И. Синтез нерекурсивных дискретных фильтров во временной области // Информационно-управляющие системы. 2016. № 5. С. 98—101.
5. Зиатдинов С. И. Синтез рекурсивных дискретных фильтров во временной области // Изв. вузов. Радиоэлектроника. 2016. № 3. С. 3—6.
6. Зиатдинов С. И. Анализ линейных систем на основе переходных характеристик // Информационно-управляющие системы. 2016. № 2. С. 104—106.
7. Зиатдинов С. И. Импульсная характеристика комплексного полосового фильтра Баттерворта // Изв. вузов. Приборостроение. 2015. Т. 58, № 8. С. 167—172.
8. Зиатдинов С. И. Синтез комплексного фильтра с заданной передаточной функцией // Изв. вузов. Приборостроение. 2016. Т. 59, № 7. С. 253—259.
9. Микропроцессорные системы автоматического управления / Под общ. ред. В. А. Бесекерского. Л.: Машиностроение, 1988. 355 с.
10. Бесекерский В. А. Цифровые автоматические системы. М.: Наука, 1976. 576 с.
11. Брычков Ю. А., Маричев О. И., Прудников А. П. Таблицы неопределенных интегралов. М.: Наука, 1986. 192 с.
Сергей Ильич Зиатдинов Юлия Витальевна Соколова
Рекомендована кафедрой информационно-сетевых технологий
Сведения об авторах д-р техн. наук, профессор; СПбГУАП, кафедра информационно-сетевых технологий; E-mail: [email protected]
аспирант; СПбГУАП, кафедра информационно-сетевых технологий; E-mail: [email protected]
Поступила в редакцию 20.02.17 г.
Ссылка для цитирования: Зиатдинов С. И., Соколова Ю. В. Анализ комплексных фильтров на основе переходных характеристик // Изв. вузов. Приборостроение. 2017. Т. 60, № 7. С. 641—646.
ANALYSIS OF COMPLEX FILTERS WITH THE USE OF TRANSITION CHARACTERISTICS
S. I. Ziatdinov, Yu. V. Sokolova
St. Petersburg State University of Aerospace Instrumentation, 190000, St. Petersburg, Russia E-mail: [email protected]
A method for studying complex filter with the use of the transition characteristics is developed. Presentation of input signal as a sum of step functions shifted relative to each other leads to a new form of complex superposition integral allowing to determine the filter output signal for a given complex-valued input and a known complex transition characteristic. Application of complex filters is shown to enable variation of their frequency settings without a change in shape of the amplitude and phase-frequency characteristics. The concepts of complex delta function and complex transition characteristics are introduced; the concepts are employed in the study of dynamic properties of the filters and transmission of various signals. In contrast to the known method of analysis of linear systems based on the pulse characteristic, the proposed method may be used to synthesize complex discrete filter with dynamic properties coinciding with those of a continuous filter-analog. The developed theoretical principles are confirmed by presented examples.
Keywords: complex signal, complex filter, complex pulse characteristic, complex transition characteristic, superposition integral
Data on authors
Sergey I. Ziatdinov — Dr. Sci., Professor; St. Petersburg State University of Aerospace In-
strumentation, Department of Information Network Technologies; E-mail: [email protected]
Yulia V. Sokolova — Post-Graduate Student; St. Petersburg State University of Aerospace
Instrumentation, Department of Information Network Technologies; E-mail: [email protected]
For citation: Ziatdinov S. I., Sokolova Yu. V. Analysis of complex filters with the use of transition characteristics. Journal of Instrument Engineering. 2017. Vol. 60, N 7. P. 641—646 (in Russian).
DOI: 10.17586/0021-3454-2017-60-7-641-646