УДК 621.396:681.323 DOI: 10.17586/0021-3454-2022-65-1-19-27
СИНТЕЗ ДИСКРЕТНЫХ ПОЛОСОВЫХ И РЕЖЕКТОРНЫХ ФИЛЬТРОВ С ИСПОЛЬЗОВАНИЕМ ИНВАРИАНТНЫХ ИМПУЛЬСНЫХ И ПЕРЕХОДНЫХ ХАРАКТЕРИСТИК
С. И. Зиатдинов
Санкт-Петербургский государственный университет аэрокосмического приборостроения
Санкт-Петербург, Россия ziat.53@mail.ru
Аннотация. Рассмотрен вопрос построения дискретных полосовых и режекторных фильтров с использованием импульсных и переходных характеристик непрерывных фильтров-аналогов нижних и верхних частот. Разработана методика синтеза перестраиваемых полосовых и режекторных фильтров с неизменной шириной амплитудно-частотных характеристик независимо от частоты настройки фильтров. Найдены весовые коэффициенты разностных уравнений, определяющих работу дискретных полосовых и режекторных фильтров, обладающих неизменной шириной амплитудно-частотных характеристик в условиях переменной частоты настройки фильтров. В основу синтеза положено представление математической частотной передаточной функции фильтра отдельно в области положительных и отрицательных частот, для каждой из которых применяется смещенное на частоту настройки выделенной области билинейное z-преобразование. Получены передаточные функции в z-плоскости дискретных фильтров отдельно для областей положительных и отрицательных частот. Путем суммирования и перемножения передаточных функций получены итоговые передаточные функции в z-плоскости для дискретных полосовых и режекторных фильтров. Переходом от z-параметра к параметру jœ получены выражения для частотных передаточных функций синтезированных полосовых и режекторных фильтров. Произведен расчет амплитудно-частотных характеристик синтезированных фильтров и показано, что их ширина остается постоянной при изменении частоты настройки фильтров. Разработанная методика полезна при построении адаптивных систем и устройств обработки сигналов, таких как системы обнаружения и фильтрации сигналов, допле-ровские измерители скорости, системы селекции движущихся целей.
Ключевые слова: дискретные фильтры, частотные характеристики, импульсные и переходные характеристики, z-преобразование, весовые коэффициенты
Ссылка для цитирования: Зиатдинов С. И. Синтез дискретных полосовых и режекторных фильтров с использованием инвариантных импульсных и переходных характеристик // Изв. вузов. Приборостроение. 2022. Т. 65, № 1. С. 19—27. DOI: 10.17586/0021-3454-2022-65-1-19-27.
SYNTHESIS OF DISCRETE BANDPASS AND NOTCH FILTERS USING INVARIANT IMPULSE AND TRANSIENT CHARACTERISTICS
S. I. Ziatdinov
St. Petersburg State University of Aerospace Instrumentation, St. Petersburg, Russia
ziat.53@mail.ru
Abstract. The problem of constructing discrete bandpass and notch filters using impulse and transient characteristics of continuous analog filters of low and high frequencies is considered. A technique for the synthesis of tunable bandpass and notch filters with a constant width of the amplitude-frequency characteristics, regardless of the filter tuning frequency, is developed. The weight coefficients of the difference equations that determine the operation of discrete bandpass and notch filters with a constant width of the amplitude-frequency characteristics under conditions of a variable filter tuning frequency are found. The synthesis is based on the presentation of the mathematical frequency transfer function of the filter separately in the region of positive and negative frequencies, for each of which a bilinear z-transform is applied shifted to the tuning frequency of the selected region. The transfer functions in the z-plane of discrete filters are obtained separately for the regions of positive and negative frequencies. By summing and multiplying the transfer functions, the final transfer functions in the z-plane for discrete band pass and notch filters are obtained. By passing from the z-parameter to the jw parameter, expressions are obtained for the frequency transfer functions of the synthesized bandpass and notch filters. The amplitude-frequency characteristics of the synthesized filters are evaluated, and their width is shown to remain constant when the filter tuning frequency is changed. The developed technique is reported to be of use in building adaptive systems and signal processing devices, such as signal detection and filtering systems, Doppler velocity meters, and moving target selection systems.
© Зиатдинов С. И., 2022 JOURNAL OF INSTRUMENT ENGINEERING. 2022. VOL. 65, N 1
Keywords: discrete filters, frequency characteristics, impulse and transient characteristics, z-transform, weight coefficients
For citation: Ziatdinov S. I. Synthesis of discrete bandpass and notch filters using invariant impulse and transient characteristics. Journal of Instrument Engineering. 2022. Vol. 65, N 1. P. 19—27 (in Russian). DOI: 10.17586/00213454-2022-65-1-19-27.
Введение. При разработке различных радиоэлектронных систем и устройств обнаружения и оценки параметров сигналов, систем селекции движущихся целей, систем автоматического сопровождения объектов по таким параметрам, как дальность, угловые координаты и скорость, широко используются разнообразные фильтры нижних и верхних частот, полосовые и режекторные фильтры [1—4].
В общем виде фильтры исчерпывающе описываются их частотными передаточными функциями. В то же время в зависимости от тактических задач фильтры могут быть заданы однозначно взаимосвязанными импульсными и переходными характеристиками, дифференциальными, интегральными и разностными уравнениями. В зависимости от того, какими параметрами определен конкретный фильтр, существуют разнообразные методы синтеза, которые достаточно широко и подробно изложены в отечественной и иностранной литературе.
При построении дискретных фильтров задача синтеза сводится к нахождению весовых коэффициентов разностного уравнения, определяющего алгоритм работы. Так, в [5—7] представлена методика вычисления весовых коэффициентов нерекурсивных и рекурсивных дискретных фильтров нижних и верхних частот по отсчетам их импульсной характеристики. Методика синтеза дискретного фильтра по отсчетам переходной характеристики изложена в [8, 9]. В [10, 11] рассмотрен синтез дискретного фильтра на основе использования дифференциального и интегрального уравнений, определяющих работу фильтра.
Синтез дискретных фильтров в большинстве случаев осуществляется по заданной частотной передаточной функции непрерывного фильтра-аналога с использованием билинейного ¿-преобразования [12—15]. Метод билинейного ¿-преобразования наиболее эффективен при разработке фильтров нижних и верхних частот. Однако при синтезе полосовых и режек-торных фильтров возникают значительные трудности.
Постановка задачи. Полосовые и режекторные фильтры, как правило, являются узкополосными и характеризуются частотой настройки. Применение данного типа фильтров, например, в адаптивных системах, где частота настройки изменяется, может привести к нежелательным явлениям, связанным с изменением полосы пропускания или зоны режекции фильтров. В аналоговом варианте полосовые и режекторные фильтры выполняются обычно на основе реактивных элементов в виде катушек индуктивности и конденсаторов. Анализ показывает, что изменение в два раза резонансной частоты (частоты настройки) последовательного или параллельного колебательного контура приводит к практически двойному изменению полосы пропускания резонансной системы. Представляется актуальной разработка методики синтеза полосовых и режекторных перестраиваемых дискретных фильтров на базе заданных непрерывных фильтров-аналогов с неизменной шириной полосы пропускания или зоны режекции. Решение данной задачи составляет основное содержание настоящей статьи.
Синтез дискретных полосовых фильтров с использованием инвариантных импульсных характеристик. При известной импульсной характеристике непрерывного фильтра-аналога весовые коэффициенты дискретного фильтра определяются на основе следующих соотношений [5]:
А = Щ;
Г , ^
А = т ъг +Хвд.} ,
I ■=1 )
,=1, 2, ..., п; (1)
и
ВД +,--]■ =^п+г, ^ ,2, (2)
]=1
где Т — период следования отсчетов входного сигнала; й,=й(/Т) — отсчеты импульсной характеристики непрерывного фильтра-аналога; А,, В] — весовые коэффициенты нерекурсивной и рекурсивной частей передаточной функции дискретного фильтра-аналога; п — порядок фильтра.
По системе п линейных уравнений (2) определяются коэффициенты В,, коэффициенты А, рассчитываются последовательно по формуле (1). В результате передаточная функция дискретного фильтра-аналога в ¿-плоскости записывается следующим образом:
п
X Аг-
Ж (г ) =-^=0-, (3)
1 + Х ВгГ" ,=1
где г" = е]Т.
Известно [12,15], что частотная передаточная функция действительных фильтров включает в себя две ветви, расположенные симметрично в области положительных и отрицательных частот. Для полосового непрерывного фильтра с постоянной шириной полосы пропускания частотную передаточную функцию запишем как
Ж(]Ю) = ат[](ю + Ю0)Г + ат-1[■ (ю + Юр)]"""1 + ...а1](ю + Юр) + ар т <п (4)
Ьп[](ю + Юр)]п + Ьп-1[](Ю + Юр)]п-1 + ..Ъх](ю + Юр) + Ър ' ' где Юр — частота настройки фильтра; а,, Ъ, — коэффициенты, определяющие вид амплитудно-частотной и фазочастотной характеристик фильтра.
Для синтеза дискретного полосового фильтра частотную передаточную функцию (4) представим в виде
Жп(/Ю) = Ж1( ]ю) + Ж2(]Ю),
где
Ж(]Ю) = ат [](ю " юр)]т + ат-1[](ю - Юр)]'"1 + . а1](ю -Юр) + ар (5)
Ъп [](ю - Юр)]п + Ъп-1[](ю - Юр)]п-1 +.. А](ю - Юр) + ър Ж (/Ю) = ат[](Ю+Юр)]т + ат-1[](ю+юр)]ТО^1 + ...^Л^^ + ар 2 Ъп [ ] (Ю+Юр )]п + Ъп-1 [] (Ю+Юр )]п-1 +.. ЬЪ1] (Ю+Юр ) + Ър
Соотношение (5) определяет частотную передаточную функцию полосового фильтра в области положительных частот, а соотношение (6) — в области отрицательных частот. Частотные передаточные функции (5) и (6) соответствуют двум комплексным фильтрам с частотами настройки Юр и -Юр соответственно. Передаточные функции данных фильтров в ¿-плоскости определяются следующими выражениями:
— для фильтра с частотной передаточной функцией (5) и частотой настройки Юр
т
X Агр г-
Ж1( г ) = -; (7)
1+ Х В,гр г -,=1
— для фильтра с частотной передаточной функцией (6) и частотой настройки -ш
о
m
I A1z-1z-'
Щ>(z) = -. (8)
1 + ¿3*-*
I=1
В выражениях (7) и (8) z0 = е']ШоТ; весовые коэффициенты Ai, Bi определяются соотношениями (1) и (2).
В результате передаточная функция синтезируемого дискретного полосового фильтра в z-плоскости принимает вид
Wп(z) = Щ^) + Щ^г). (9)
Пример. Синтез дискретного полосового фильтра с постоянной полосой пропускания. В качестве непрерывного фильтра-аналога возьмем фильтр нижних частот Баттерворта второго порядка с частотной передаточной функцией [6, 7]
2
юсп
W О)=—2—г2-- (10)
(7'ю) + И2 • шсрш + шср
где юср — частота среза фильтра.
Импульсная характеристика данного фильтра имеет вид
• ш^р^ВШ^ . (11)
\2
Система уравнений (2) для нахождения весовых коэффициентов Bi записывается следующим образом:
Blh2 + B2 / = -/3; ]
B1h3 + В2 /2 = - /4,
(12)
где /1, /2, /4 — отсчеты импульсной характеристики (11) в моменты времени t, равные Т, 2Т, 3Т, 4Т соответственно.
Решение системы уравнений (12) имеет вид
в = /2 - /2/4 В = /з + В2// (13)
в2 = ,2 , , , В1 =—;-. (13)
/2 - /3 / п2
Выражения для весовых коэффициентов Аг-, определяемых из уравнения (1), с учетом (13) записываются следующим образом:
Ао = Т/о, А = Т(/ + В1/о), А2 = Т/ + В1/1 + В2/о). (14)
После подстановки (13) и (14) в (9) результирующая передаточная функция синтезируемого полосового фильтра второго порядка в z-плоскости будет определяться выражением
4
) = -^=0-, (15)
1+1
1=1
где
Со = 2Ао, С1 = 2(А1 + АоВ1)сов(шоТ), С2 = АХВХ + 2(А2 + АоА2)сов(2шоТ), С3 = 2(А2Я + А^созНТ), С4 = 2А2В2, А = 2В соэ(шоТ), А2 = 2В2 сов(2шоТ) + В?, А3 = 2В1В2 сов(шоТ), А4 = В|.
соотношении
z-1 =е-JaT,
z"2 =е-j2юТ,
z -3 = е - j 3»Т и
После подстановки в выражение (15)
г-4 = е-34юТ получим частотную передаточную функцию синтезированного дискретного полосового фильтра ( Жп( /)) в непрерывном варианте, модуль которой для трех
частот настройки/0= ю0/2л, равных 100, 200 и 300 Гц (кривые 1—3 соответственно), в области положительных частот показан на рис. 1. При расчетах были приняты следующие значения: частота среза /ср= юср/2л=8 Гц, период повторения
Г=10 с. Как видно из графиков, полоса пропускания и форма амплитудно-частотной характеристики фильтра не зависят от частоты настройки.
200 Рис. 1
300
Я Гц
Синтез дискретных полосовых фильтров с использованием инвариантных переходных характеристик. В случае когда непрерывный фильтр-аналог задан не импульсноИ, а переходной характеристикой, необходимо воспользоваться методикой, изложенной в [8].
Известно [15], что импульсная характеристика фильтра является производной от его переходной характеристики
ц,) = М) = lim gfiWfi-0
dt т
С учетом дискретности синтезируемого фильтра данное соотношение запишем в виде конечной разности
h(t) = Т [g (t) - g (t - T )]= T Ag (t),
(16)
где Ag (г) = g (г) - g (г - Т).
Тогда с учетом соотношения (16) выражения (1) и (2) для весовых коэффициентов дискретного фильтра примут следующий вид:
А = М0];
г
А = Ag[/ ] + X В3 Ag[г - 3 ], г=1, 2, ..., п;
3=1
п
-X В3Дg[n + г - 3] =М1 + г], г=1, 2, ...,
3=1
где Л?[г] = - g[i-1].
Пример. Синтез дискретного полосового фильтра. В качестве непрерывного фильтра-аналога возьмем фильтр нижних частот Баттерворта второго порядка, частотная передаточная
функция которого определяется выражением (Ю). Данный фильтр имеет переходную характеристику [8]
1 -Юср^Г • Юсрг , Юср^ gЦ) = 1 -е ср' sln—-¡=- + cos—рг-
I у/2 л/2
Тогда система уравнений (12) примет вид
В^[2] + В2 Ag[1] = -Ag[3]; 1 В^[3] + В2 Ag[2] = -А&[2],| а решение данной системы уравнений — следующий вид:
в = Ag 2[3] - Ag[2]Ag[1] в = Ag[3] + B2Ag[1] (17)
2 Ag 2[2] - Ag[3]Ag[1]' 1 Ag[2] "
Выражения для весовых коэффициентов А,, определяемых из соотношения (1), с учетом (17) записываются следующим образом:
Ар = А1 = Д^[1] + Blg[0], А2 = Д^[2] + В^[1] + В2^[р]. (18)
Синтезированный дискретный полосовой фильтр будет иметь передаточную функцию в г-плоскости, определяемую соотношением (15), с весовыми коэффициентами
Ср = 2А), С1 = 2(А1 + А)В1)соэ(юрТ), С2 = А1В1 + 2(А2 + АрА2)соб(2ЮрТ), С3 = 2(А2В + А1В2)соз(юрТ), С4 = 2А2В2,
Д = 2В1 СОБ(ЮрТ), Б2 = 2В2СОБ(2ЮрТ) + В?, £>3 = 2В1В2соз(юрТ), £>4 = в|.
Частотная передаточная функция данного фильтра совпадает с частотной передаточной функцией ранее рассмотренного полосового фильтра (см. рис. 1).
Метод синтеза с использованием импульсной характеристики непрерывного фильтра-аналога неприемлем для синтеза режекторных фильтров, так как импульсная характеристика режекторного фильтра содержит дельта-функцию. Поэтому для синтеза дискретных режек-торных фильтров используются инвариантные переходные характеристики [8].
Синтез дискретных режекторных фильтров с использованием инвариантных переходных характеристик. Для синтеза дискретных режекторных фильтров частотную передаточную функцию (4) запишем в виде
Жр( ]ю) = Ж1( ]ю) -Ж2С/ю),
где, как и ранее,
Ж1(]Ю) = ат [](Ю - Юр)]т + ат-1[](Ю -Ю+ . . .а1](Ю -Юр) + ар (19)
Ъп [](Ю - Ю^Г + Ъп-1[](Ю - Юр)]п-1 + . ..Ь1](ю- Юр) + Ьр Ж (]Ю) = ат [](ю+юр )]т + ат- 1 [](ю+юр )]ТО ^ + ...а1](ю+юр ) + ар (2р)
2 Ъп []■ (Ю+Юр)]п + Ъп-1^./ (Ю+Юр)]п -1 + ..Ъ1] (Ю+Юр) + Ър
Соотношение (19) определяет частотную передаточную функцию режекторного фильтра в области положительных частот, а соотношение (2р) — в области отрицательных частот.
Тогда на основании вышеизложенной методики результирующая передаточная функция синтезируемого дискретного режекторного фильтра в г-плоскости записывается следующим образом:
Жр(г) = Ж1(г)-Ж2(г). (21)
Пример. Синтез дискретного режекторного фильтра с постоянной шириной зоны режекции. В качестве непрерывного фильтра-аналога возьмем фильтр верхних частот второго порядка с частотной передаточной функцией
Ж (ую) =
(ую)2
(ую)2 + уЧ/2 • юсрю + ю;
ср
2 ср
Данный фильтр имеет переходную характеристику
* ()=ю2р 72
^ . юср^ юср^ ^
Sin—+ cos
V
I- I V ' л .-
у!2 л/2
у
Весовые коэффициенты Ви А^ синтезируемого дискретного режекторного фильтра определяются из соотношений (17) и (18).
С учетом выражения (21) результирующая передаточная функция рассматриваемого фильтра в г-плоскости будет определяться выражением (15), где весовые коэффициенты имеют вид
С0 = А2, с1 = 2А)А1 соэ(ю0Т), С2 = А2 + 2АоА2 соэ(2ю0Т), С3 = 2 А1А2 СОБ(ЮОТ), С4 = 2 А|, Д = 2В1 соэ(ю0Т), Б2 = 2В2 соб(2Ю0Т) + В?, £>3 = 2В1В2соэ(ю0Т), = В|.
После подстановки в выражение (15) соотношений г = е
-1 _ „-зюТ
г-2 =е-у2юТ,
г 3 = е 3 3юТ и г 4 = е 3 4юТ получим частотную передаточную функцию синтезированного дискретного режекторного фильтра (Жр(/)) в непрерывном варианте, модуль которой для
трех частот настройки/0= ю0/2п, равных 100, 200 и 300 Гц (кривые 1—3 соответственно), в области положительных частот показан на рис. 2. При расчетах, как и ранее, принято: /ср= =юср/2п=8 Гц, Г=10-3с. Анализ графиков показывает, что ширина зоны режекции и форма амплитудно-частотной характеристики фильтра не зависят от частоты настройки.
Жр(/
100
200 Рис. 2
300
/, Гц
Обсуждение результатов. Представленные на рис. 1 и 2 амплитудно-частотные характеристики действительных дискретных полосовых и режекторных фильтров имеют постоянную ширину при различных значениях частоты настройки фильтров. Данный результат является следствием того, что в основу синтеза положено представление действительных фильтров в виде двух самостоятельных узкополосных комплексных фильтров, частотные передаточные функции которых расположены симметрично относительно нуля частот как в области положительных, так и в области отрицательных частот с последующим применением для каждой частотной передаточной функции билинейного г-преобразования, смещенного на частоту настройки. Сложение частотных передаточных функций комплексных фильтров для синтеза полосовых фильтров или их перемножение для синтеза режекторных фильтров приводит к увеличению в два раза порядка результирующих полосовых и режекторных фильтров.
1
Заключение. Представленная методика разбиения частотной передаточной функции действительного непрерывного фильтра-аналога на две составляющие позволяет синтезировать перестраиваемые по частоте действительные полосовые и режекторные фильтры с постоянной шириной амплитудно-частотной характеристики независимо от частоты настройки фильтров.
Предложенная методика синтеза действительных дискретных полосовых и режекторных фильтров, рассмотренные конкретные примеры построения подобных фильтров будут полезны при создании адаптивных и когерентных систем обработки информации, таких как устройства селекции движущихся целей, доплеровские измерители скорости движения разнообразных объектов, обнаружители, устройства оценки параметров, местоположения объектов и т. д.
СПИСОК ЛИТЕРАТУРЫ
1. Бакулев П. А., Стенин В. М. Методы и устройства селекции движущихся целей. М.: Радио и связь, 1986. 286 с.
2. Бакулев П. А. Радиолокационные системы. М.: Радиотехника, 2004. 319 с.
3. Попов Д. И. Адаптивное подавление пассивных помех II Цифровая обработка сигналов. 2014. № 4. С. 32—37.
4. Котоусов А. С., Морозов А. К. Оптимальная фильтрация и компенсация помех. М.: Горячая линия-Телеком, 2008. 166 с.
5. Воробьев С. Н. Цифровая обработка сигналов. СПб: Изд. дом „Академия", 2013. 318 с.
6. Зиатдинов С. И. Синтез рекурсивных дискретных фильтров во временной области II Изв. вузов. Радиоэлектроника. 2016. № 3. С. 3—6.
7. Зиатдинов С. И. Синтез нерекурсивных дискретных фильтров во временной области II Информационно-управляющие системы. 2016. № 5. С. 98—101.
8. Зиатдинов С. И. Анализ линейных систем на основе переходных характеристик II Информационно-управляющие системы. 2016. № 2. С. 104—106.
9. ХеммингР. В. Цифровые фильтры. М.: Недра, 1987. 217 с.
10. Зиатдинов С. И. Синтез дискретных фильтров методами инвариантных дифференциальных и интегральных уравнений II Изв. вузов. Приборостроение. 2019. Т. 62, № 5. С. 424—432.
11. Микропроцессорные системы автоматического управления I Под общ. ред. В. А. Бесекерского. Л.: Машиностроение, 1988. 355 с.
12. Баскаков С. И. Радиотехнические цепи и сигналы. М.: URSS, 2016. 915 с.
13. Оппенгейм А., Шафер Р. Цифровая обработка сигналов. М.: Техносфера, 2006. 855 с.
14. КуприяновМ. С., Матюшкин Б. Д. Цифровая обработка сигналов. СПб: Политехника, 2000. 592 с.
15. Гоноровский И. С. Радиотехнические цепи и сигналы. М.: Радио и связь, 1969. 512 с.
Сведения об авторе
Сергей Ильич Зиатдинов — д-р техн. наук, профессор; Санкт-Петербургский государственный
университет аэрокосмического приборостроения, кафедра информационно-сетевых технологий; профессор; E-mail: ziat.53@mail.ru
Поступила в редакцию 07.06.2021; одобрена после рецензирования 02.07.2021; принята к публикации 02.12.2021.
REFERENCES
1. Bakulev P.A., Stenin V.M. Metody i ustroystva selektsii dvizhushchikhsya tseley (Methods and Devices for Selection of Moving Targets), Moscow, 1986, 286 p. (in Russ.)
2. Bakulev P.A. Radiolokatsionnyye sistemy (Radar Systems), Moscow, 2004, 319 p. (in Russ.)
3. Popov D.I. Digital signal processing, 2014, no. 4, pp. 32-37. (in Russ.)
4. Kotousov A.S., Morozov A.K. Optimal'naya fil'tratsiya i kompensatsiya pomekh (Optimal Filtering and Noise Compensation), Moscow, 2008, 166 p. (in Russ.)
5. Vorobiev S.N. Tsifrovaya obrabotka signalov (Digital Signal Processing), St. Petersburg, 2013, 318 p. (in Russ.)
6. Ziatdinov S.I. Journal of the Russian Universities. Radioelectronics (Izvestiya Vysshikh Uchebnykh Zavedenii Rossii. Radioelektronika), 2016, no. 3, pp. 3-6. (in Russ.)
7. Ziatdinov S.I. Information and Control Systems, 2016, no. 5, pp. 98-101. (in Russ.)
8. Ziatdinov S.I. Information and Control Systems, 2016, no. 2. pp. 104-106. (in Russ.)
9. Hamming R.W. Digital Filters, Englewood Cliffs, N J, 1977.
10. Ziatdinov S.I. Journal of Instrument Engineering, 2019, no. 5(62), pp. 424-432. (in Russ.)
11. Besekersky V.A., ed., Mikroprotsessornyye sistemy avtomaticheskogo upravleniya (Microprocessor-Based Automatic Control Systems), Leningrad, 1988, 355 p. (in Russ.)
12. Baskakov S.I. Radiotekhnicheskiye tsepi i signaly (Radio Engineering Circuits and Signals), Moscow, 2016, 915 p. (in Russ.)
13. Oppenheim A.V., Schafer R.W. Discrete-Time Signal Processing, Georgia Institute of Technology, Pearson, 2010.
14. Kupriyanov M.S., Matyushkin B.D. Tsifrovaya obrabotka signalov (Digital Signal Processing), St. Petersburg, 2000, 592 p. (in Russ.)
15. Gonorovsky I.S. Radiotekhnicheskiye tsepi i signaly (Radio Engineering Circuits and Signals), Moscow, 1969, 512 p. (in Russ.)
Data on author
Sergey I. Ziatdinov — Dr. Sci., Professor; St. Petersburg University of Aerospace Instrumentation, Department
of Information and Network Technologies; Professor; E-mail: ziat.53@mail.ru
Received 07.06.2021; approved after reviewing 02.07.2021; accepted for publication 02.12.2021.