Корректность задачи адаптивного управления диаграммой направленности антенны в неоднородной среде Текст научной статьи по специальности «Физика»

Корректность задачи адаптивного управления диаграммой

направленности антенны в неоднородной среде

В.К. Игнатьев, А.В. Никитин

Волгоградский государственный университет

Аннотация: В работе представлены результаты исследования корректности задачи адаптивного управления диаграммой направленности антенны в неоднородной среде. С помощью спектральных методов найдено решение уравнения Гельмгольца с переменными параметрами в приближении замороженной среды и впервые обоснована корректность задачи формирования заданной диаграммы направленности в неоднородной среде по измеренному распределению комплексной амплитуды волны в плоскости приемников. Обоснована возможность использования антенны с синтезированной апертурой для адаптивной коррекции диаграммы направленности излучателя, оптимальной для обеспечения максимальной пропускной способности канала связи. Ключевые слова: канал связи, диаграмма направленности, комплексная амплитуда, синтезированная апертура, метод углового спектра, адаптивное управление.

Создание как оптического, так и акустического высокоскоростного канала двухсторонней связи с мобильным объектом является в настоящий момент не решенной, но в соответствии с современными тенденциями развитии информационно-телекоммуникационных навигационных систем практически актуальной задачей [1]. Ее решение требует создания антенн с синтезированной апертурой [2], обладающих одновременно широкой полосой и высокой чувствительностью, например - фазированных решеток [3]. Адаптивное управление диаграммой направленности передающей антенны позволяет фокусировать излучение на приемнике, а адаптивное управление приемной антенной позволяет осуществлять пространственную фильтрацию помех, направление на источник которых неизвестно [4].

Управление диаграммой направленности на практике осуществляется коррекцией волнового фронта излучателя, а основной причиной, ухудшающей фокусировку как оптической [5], так и акустической [6] волны на приемнике, являются неоднородности среды распространения, как регулярные, так и случайные. В оптической связи для компенсации

крупномасштабных флуктуаций скорости ветра, тепловой дефокусировки в нерегулярной среде и др. с успехом используется метод апертурного зондирования, численному исследованию которого посвящен ряд работ [7]. Применимость в гидроакустических системах адаптивных фазированных антенных решеток [8] и параметрических антенн, которые формируют направленное излучение в пределах небольших телесных углов при практически полном отсутствии бокового поля [9] в известной мере ограничивается неидентичностью фазовых характеристик гидроакустических преобразователей, размещаемых в узлах антенных решеток. Несмотря на значительные усилия, предпринятые по совершенствованию конструкций гидроакустических антенн [10] и исследованию влияния конструктивных параметров на их эффективность [6, 11], обеспечение идентичности пьезокерамических датчиков является в настоящее время нерешаемой проблемой.

Разброс фазовых характеристик элементов антенной решетки в определенной мере эквивалентен наличию дополнительной неоднородности среды распространения и может быть компенсирован методами апертурного сканирования, если в качестве приемника использовать антенную решетку, аналогичную излучающей, и измерять распределение комплексной амплитуды волны в плоскости матричного приемника. Современные методы оценивания угла фазового сдвига между двумя квазигармоническими сигналами [12] обеспечивают необходимую для устойчивой работы системы точность измерений.

Комплексная амплитуда гармонической волны в неоднородной среде с

фазовой скоростью с(г) удовлетворяет уравнению Гельмгольца:

2

АЛ + А = 0. (1)

с (г)

Положим ю/с(г) = к [1 + «(г)]. Для параксиальных волн, распространяющихся вблизи продольной оси г, положим

А(х, у, г) = Ат(х, у, г)ехр(г'кг), (2)

где функция Ат(х, у, г) меняется в направлении г, то есть в продольном, медленно в сравнении со скоростью изменения в сечении (х, у), то есть в

поперечном направлении. Тогда, подставляя выражение (2) в уравнение (1) и

2 2

пренебрегая слагаемым д Ат/дг , получаем параболическое уравнение:

2 А А А

—Г + -дАт + 2^к—гт + к 2«(2 + п )Ат = 0, Ат (х, у,0) = Ао (х, у), (3) дх ду дг

где А0(х, у) - комплексная амплитуда волны в сечении г = 0. Для слабо

неоднородной среды при |п(г)| << 1 это уравнение можно решать методом

последовательных приближений, положив Ат(г) = А1(г) + А2(г) + ..., где |А2(г)|

~ |п(г)|-^^1(г)|. Уравнение (3) распадается на два уравнения с

соответствующими начальными условиями:

+ + 21кдА = 0, А1 (х, у,0)= А00 (х, у), (4)

дх 2 ду 2 дг

^ + + 2^ = -2к2пА1, А2 (х, у,0)= 0. (5)

дх ду дг

Для решения уравнений (4) и (5) методом углового спектра положим

1 о о

А1,2 (x, У, г) = Ца1,2 ((, ку, г)ехр(/кхх + 1куу Умку , (6)

—о—о

1,2 ((х, ку, г) = Ц А1,2 (x, У, г)еХР(— — 1куу )dxdУ,

—о—о

о о

(, ку) = Л 4)(x, у )ехР(— 1кхх — &уу .

(7)

а0

Выполняя над уравнениями (4) и (5) преобразования Фурье вида (7), получим обыкновенные дифференциальные уравнения

о—о

-(к2 + К )а2 = — 2к 2Ъ (2), а2 (°)= 0,

й2

где обозначено

то то

Ъ ((, ку, г) = Л п(х, у, г ) (х, у, г )ехр(- ¡кхх — ¡куу \dxdy

—то—то

Решение однородного уравнения (8) имеет вид

/

а1 (х, ку , 2)= а° (х, ку )еХР поэтому с учетом формулы (6) получаем

(8) (9)

(10)

' X + К ^

— I—-^2

V

2к

А

1 то то

1((у,2) = 4П2 Па°((,ку)ехР

' к2 + К ^

I—-^2

2к

ехр(г'кхх + йуУУмку . (11)

Решение неоднородного уравнения (9) имеет вид

а,

2

((х, ку , 2 )= ¡к | еХР

' М + ку2

— V у

2к

(2 — 2') к (х, ку, ,

(12)

поэтому с учетом формулы (6) получаем

А (х, у, 2) =

~ / -О -О

¡к

2 тото ' к2 + к2

4п

Ш

0 —то—то

ехр

-¡

х у

2к

(2 — Ъ1 (кх, ку, 2')ехр(кхх + ¡куу)йМку^'.

(13)

Выберем оси х и у в плоскости решетки излучателей, а ось 2 нормальной к ней. Для того чтобы задать начальное распределение комплексной амплитуды в виде среза гауссового пучка:

А0 (х, у ) = А ехр

( 2 , 2 Л

х + у

2 Л

V у

ехр

22 ¡кх +у л

2 Р

V у

на решетке излучателей следует задать следующее распределение:

то—то

Ао (х, у ) = Лт ехр

( 2 , 2 Л ' х + у л

2 В2Т у

ехр

22 1кх +у л

V

2 Я,

ехр (¡к (- рх + ау)). (14)

т у

Здесь Вт - ширина гауссового пучка на транспаранте передатчика, Ят -радиус кривизны волнового фронта на передатчике, Ат - комплексная амплитуда излучения в центре передатчика, а и в - полярный и нормальный углы волнового вектора, которые можно считать малыми.

Тогда, используя формулы (7) и (11), получим выражение для комплексной амплитуды гауссового пучка на расстоянии г от передатчика:

А (х, у, 2) = Ао ехр

Г (гВЯт2 - ¡г2В(Ят + г)- ¡Я]г)(х2 + у2)

к-

2 ( (Ят + г )2 + Ят 2 г2

(15)

где гВ = кВ^ - дифракционная длина пучка,

Аг

Ят (гВ (Ят + г)- ¡Ятг) ( ( + г )2 + Ят 2 г2) '

(16)

Уравнение (15) описывает гауссов пучок с шириной

В (г) = Вту/(1 + г/Ят )2 +(г/гд )2

и радиусом кривизны

Я(г )

г

В

Координата перетяжки

гА

(1 + г/Ят )2 + (г/ гр )2 г/ гВ + (1 + г/Я )гВ/ЯГ

Ят

1 + Я/ г

(17)

В

достигает максимума, равного гВ/2, при Ят = - гВ. Минимальная ширина пучка достигается в точке перетяжки и равна

Во = В(го ) = Втго1 гВ , Отклонение точки перетяжки гауссового пучка от центра приемника может быть обусловлено как погрешностью в определении расстояния от передатчика до приемника, так и неоднородностью среды распространения.

Пусть А (х, у, 2\) - комплексная амплитуда пучка на приемнике, сформированная при комплексной амплитуде А0(х, у) на передатчике, где г! -координата приемника. Найдем коррекцию ЛА0(х, у) комплексной амплитуды на передатчике, создающей на приемнике оптимальное распределение комплексной амплитуды А\(х, у, 2\), при котором точка перетяжки находится в центе приемника, то есть 2\ = 20. Обозначив отклонение комплексной амплитуды на приемнике от оптимальной ЛА (х, у, 2\) = А (х, у, 2\) - А\(х, у, 2\), получим:

то то

Ла((х, ку, 2\) = || ЛА(х, у, 2\ )ехр(- ¡кхх — ¡куу')dxdy .

—то—то

Тогда из формул (6), (7) и (!!) можно получить:

\ то то ^

ЛA(x,y,^)=— — ЦЛао((,ку)ехР -

—то V

Г к 2 + к2 л

х у

— I-— 2,

2к

ехр(кхх + ¡куу )к^ку

у

.кх2+к; ^

I- 2\

2к

Ла0 ((х , ку ) = —Ла((х , ку , 21 )еХР

V У

а затем найти коррекцию комплексной амплитуды на передатчике:

\ то то

ЛА0 ( у )=— 4^2 | |Ла((х , ку , 21 )еХР

с

К + ку2

21 + кхх + куу

2к 1 х у"

dkxdky =

к

2п2\

тото

| |ЛА(х!, у\, 2\ )ехР

к ( — х!)2 +(у — л) 22

2

(!8)

dx!dy!.

— со—с»

Из формулы (!8) следует возможность адаптивной коррекции диаграммы направленности передающей антенны по измеренному в плоскости приемника распределению комплексной амплитуды излучения. Поскольку ядро интегрального преобразования в правой части формулы (18) непрерывно и ограничено по модулю, задача адаптивного управления является корректной.

со— со

При численном моделировании процесса управления диаграммой направленности антенны в неоднородной среде использовалась модель показателя преломления в виде суммы детерминированной и случайной составляющих:

I-1

п(х, у, I) = По + пхх + пуу + п21 + ^ п (х, у, г).

1=о

В качестве детерминированных составляющих выступают линейные функции координат, а в качестве случайной - "облачка" случайного размера, расположенные в случайном месте:

П ((, у, г )= По,, ехр

(х - х0,)2 +(у - уо,)2 +(г - ^):

2 Л

О?

1 = 0,...,I -1.

Здесь п0[, х01, у01, 20и о , и I - случайные величины.

Рассмотрим область пространства, в начале координат которого расположена антенная решетка передатчика размером 2Ьхх2Ьу. Проведем дискретизацию этой области с шагами Дх = Ьх/Кх, Ду = Ьу/Ку и Дг. Область пространства ограничим размерами 2Вхх2ВухОг, где Бх = 2МхДх, Бу = 2МуДу, = М2Д2, а Мх, Му и Мг являются степенями двойки. Дискретный вариант показателя преломления примет вид:

п[х , ту , т1 ]= П0 + Пхтх Дх + Путу Ду + П1т1ДI + Х П, [тх, ту, т1]

1=0

П [мх, ту, тг ] = П0, ехр

' (тхДх - х0, )2 + (тхДх - ^0, )2 + (тхДх - 10, )2 Л

О 2

тх =-Мх,...,Мх -1, ту =-Му,...,Му -1, тг =-Мг,...,Мг -1,

1 = 0,., I -1.

В качестве модели закона работы антенной решетки передатчика используем выражение (14) с заданной шириной гауссового пучка на передатчике и радиуса кривизны волнового фронта на передатчике В пределах антенной решетки дискретный вариант выражения (14) примет вид:

А0 [х , ту ] = А еХР

(тх + 1/2)2 Л2 +(ту + \2)2 Л

2 Л 2 Л

2Ц2

х

У

х ехР

к (тх + \2)2Л2х +(ту +1/2)2Л

2

у

2ЯТ

ехР(к (— (тх + 1/2)Л хв + (ту + 1/2) у а))

тх =—Кх,..., Кх — 1, ту =— Ку,..., Ку — 1.

а в остальной области:

А0 [тх, ту ] = 0, тх =—Мх ,...,—Кх + 1 тх = КхМ —1,

ту =—Му:

Ку +1, ту = Ку,..., Му — 1.

Спектр а0(кх, ку) можно рассчитать с помощью выражения (7), используя дискретное преобразование Фурье:

а0 [/х , 1у ] = а0

Мх —1 Му-1

= I I А0 [тх, ту ]ехР

' п1х п/у Л

V мх л :му л у у

— гп

V V

с

( + 1/2)х , (ту + У2).

\\

= (— 1))+/у ехР

/

. п

— г-

2

М,

1х + 1у

М

у УУ

\\

х

М х М

I I А0 [тх — Мх, ту — Му ]ехР

х

2 М х —1 2 М у —1

V Vх у J у / г

т

тх=0 ту=0

Л I

А + ту1у

— гп

V Мг Му V V х у уу

/х = 0,...,2М х — 1, /у = 0,...,2Му — 1. Это выражение подготовлено к использованию двумерного алгоритма быстрого преобразования Фурье (БПФ).

На следующем этапе в соответствии с выражением (11) с помощью обратного БПФ рассчитаем распределение А1(х, у, 2):

тх =—Мх ту =—М у

А[тх, ту, тг ] = А ((тх + 1/2)Д х, (ту + 1/2)Д у, (тг + 1/2)Д г)

4ММуДхДу х =0 /у =0

х ехр

Мх-12 Му-1

X X [х, 1у ]еХР

- . ( +1/2)Д2п2

2к

• + •

Мх2Д2х Му2Д2у

х

( " /х + /у Л (

. п 1 — ехр 1п

2 му Му

V _ у у _ у V

/хтх + ;уту

ММ

т = -М ,...,М -1, т =-М ,...,М -1, т = 0,...,М -1.

х х? ? х ? у у' 5 у 5 I ^ ^ I

Видно, что расчет сводится к применению М2 раз алгоритма БПФ к

последовательности

а

0 [, /у]

ехр

-г

( + 1/2)Д I п

2к

2 ( / 2

1х + 1у

2

. м2Д2х м2д2 ..

V х х у у у у

ехр

п

г—

2 ч.

V V

1х + 1у

ММ

у

/х = 0,...,2Мх -1, /у = 0,.,2Му -1.

для каждого значения т2.

На следующем шаге найдем дискретный образ функции Ь1(кх, ку, I), определенной выражением (10):

Ь1[/х ,/у, т ] = Ь1

п/ п/у

VМхДх МуДу

, (т +12) =(-1)+/у ДхДу ехр

п

- г-

2

/Х . + . /у

\\

ММ

х

х

2Мх —12М у-1

X X П[тх , ту , mI кк , ту , mI ]ехр

тх =0 ту =0

- гп

V Vх у у у \\

^ + ¿т.

Мх Му V V х у уу

/х = 0,.,2Мх-1, /у = 0,...,2Му-1, mI =-Мг.,Мг-1.

Затем в соответствии с выражением (12) вычислим спектр а2(кх, ку, I) как дискретную свертку:

а

т2

[/х, /у, mI ]= 1кД I X ехр

т=0

. п2Д,

2к

/2

/

2

+ ■

V МХ Му2Д2у у

(mI - т) Ь1[/х , /у , т]

и по формуле (13) с помощью обратного БПФ найдем функцию А2(х, у, I):

2

/

/

1

Л [mx, my, mz ] = А2 ((тх + 1/2)А х, (ту + 1/2)А у , ^ + 1/2)А г) =

1

4MxMyА xА у ^

Мх-12 М у -1

X X [х, 1у, mz ]еХР

/

/

. п

г—

2 ч

V V

\\

X . +. у

М~ Шу уу

ехр

г (» 1 ^ • 1тх 1уту

V V у у уу

тх = -Шх,...,Шх -1, ту = -Шу,...,Му -1, mz = 0,...,Мг -1.

На последнем этапе вычислим комплексную амплитуду волны в соответствии с формулой (2):

А[

тх, ту, mzi

= А

[х, ту, ^ ] + А2 [тх, ту, ^ ]ехр(гк( + 12)Аz ) .

Для апробации разработанного алгоритма проведено численное моделирование диаграмм направленности оптического излучателя с помощью программы, реализующей описанный алгоритм. При численном моделировании использовались аналитические выражения (15) и (16), полученные для случая а = 0 и в = 0. Их дискретный вариант имеет вид:

ATzDRT (Л + mz Аz )- 1К-ТтгАz 2х

х ехр

А1а [тх, ту, т,] 2 ( А )2 0 2 2 А2

zD (т + тгАz) + Rт т]А2

(zDRт2 (т2 А2х + т^ А2у)-^ (( + тг А г)-( тг А г ) А2х + т^ А2уУ

£

2( (( + тг А г )2 + (т2 т22 А22)

тх =-Мх,...,Мх -1, ту =-Му,...,Му -1, т2 =-Мг,...,Мг -1.

Для того чтобы описанный метод давал верные результаты, необходимо, чтобы спектры функций А0(х, у), А1(х, у, z) и А2(х, у, z) были ограничены величинами кх тах, ку тах и кг тах, а интервалы дискретизации

выбирались в соответствии с теоремой Котельникова [13]:

А <

п

к

А у < к

п

А<

п

у тах

к

На рис. 1 показаны полученные с помощью описанной программы зависимости амплитуды поля от координат А^у^) и А^х^у полученные при следующих параметрах транспаранта и области пространства, в которой производился расчет: Ьх = Ьу = 0,05 м, Dx = Dy = 0,16 м, Ах = Ау = 0,0025 м,

0

А2 = 0,4 м, Кх = Ку = 20, Мх = Му = 64, Мг = 1024, Dт = 0,0125, (т = -350, X = 1 мкм, а = в = 0. Рассчитанная по формуле (17) координата перетяжки при таких параметрах z0 = 310,5 м. Графики слева на рис. 1 соответствуют z1 = 0, справа - Z2 = 310,8 м.

0.8 -

0.4 -

0.8 -

0.4 -

-0.16

-0.08

1—1 Г -0.16 -0.08 0

0.08

0.16 у, м

2 -

0

I 1 I ^

0.08 0.16 х, м -0.16 -0.08

2 -

-0.16

-0.08

0 0.08 0.16 х, м

0.08

0.16 у, м

Рис. 1. - Аналитические (линия) и численные (точки) зависимости амплитуды поля от координат А^ху^) и А1(х1,у.^) при х1 = 0, у1 = 0, z1 = 0 и

Z2 = 310,8 м

1

0

1

0

0

На рис. 2 показаны соответствующие диаграммы, на которых черным цветом показаны максимальные, а светлым - минимальные значения. На рис. 3 показаны аналитическая и полученная численно зависимости А^х^у^).

Рис. 2. - Диаграммы А^ху,^) и А^ху,^) при ^ = 0 и ^ = 310,8 м

Al(xl,Уl,I)l

Рис. 3. - Аналитическая (тонкая линия) и численная (толстая линия) зависимости амплитуды поля А^х^у^) при х1 = 0 и у1 = 0

Видно, что полученные численно и аналитически результаты хорошо совпадают, а координата перетяжки, определенная по графику А1(х1,у1,£) (рис. 3), с точностью до Д^ составляет 311,2 м. На рис. 2 хорошо видно улучшение фокусировки в области I0.

На рис. 4 показаны диаграммы А^у^) (^ = 310,8 м) для углов а = 1,2х10^ и в = 1,4х10-4. Расчет смещений фокуса для этих углов дает 0,037 м

и 0,044 м соответственно. На рисунках видно, что получаются именно такие смещения.

Рис. 4. - Диаграммы Л^ху^) при а = 1,2 хЮ-4, в = 0 (слева) и а = 0

в = 1,4х10-4 (справа) в точке перетяжки (при 22 = 310,8 м)

Результаты численного моделирования подтверждают, что в первом приближении задача управления диаграммой направленности передатчика с синтезированной апертурой по измеренному в плоскости приемника распределению комплексной амплитуды волны является корректной. Это позволяет создать систему адаптивной коррекции диаграммы направленности излучателя, оптимальной для обеспечения максимальной пропускной способности канала связи.

Исследование выполнено при финансовой поддержке РФФИ и Администрации Волгоградской области в рамках научного проекта 15-4702297 «р_а».

Литература

1. Asadov H., Ismaylov K. Information Method for Synthesis of Optimal Data Subsystems Designated for Positioning, Location and Remote Sensing Systems // Positioning. 2011. № 2. pp. 61 - 64.

2. Stergiopoulos S. Implementation of adaptive and synthetic aperture beamformers in sonar systems // Proc. IEEE. 1998. Vol. 86. № 2. p. 358.

3. Autrey S.W. Passive Synthetic Arrays // J. Acoust. Soc. Am. August 1988. Vol. 84, pp. 592 - 598.

4. Lasky M., Doolittle R. D., Simmons B.D., Lemon S.G. Recent progress in towed hydrophone array research // IEEE J. Oceanic Eng. 2004. Vol. 29. № 2. pp. 374 - 387.

5. Osche G.R. Optical detection theory for laser applications. Wiley-Interscience. Berlin. 2002. 418 p.

6. Кириченко И.А. Повышение эффективности направленных свойств адаптивных гидроакустических систем с параметрическими антеннами // Инженерный вестник Дона. 2016. № 4. URL: ivdon.ru/ru/magazine/archive/n4y2016/3845/.

7. Егоров К. Д., Чесноков С. С. Апертурное зондирование в задачах фокусировки мощных световых пучков при пульсациях скорости ветра на трассе // Квантовая электроника. 1987. т. 14. № 6. с. 1269 - 1273.

8. Sullivan E.J. Passive Acoustic Synthetic Processing // IEEE Ocean. Eng. Soc. Newsletter. 2003. Vol. 38. № 1. pp. 21 - 24.

9. Воронин В.А., Тарасов С.М., Тимошенко В.А. Гидроакустические параметрические системы. Ростовское областное книгоиздательство. 2004. 400 с.

10. Митько В.Н., Эсси-Эзинг А.С. Проектирование и моделирование антенны гидролокатора бокового обзора // Инженерный вестник Дона. 2016. № 2. URL: ivdon.ru/ru/magazine/archive/n1y2016/3534.

11. Кириченко И.А., Котов В.Ю., Лонкин П.В., Слуцкий Д.С. Оценка влияния конструктивных характеристик на направленность гидроакустических антенн // Инженерный вестник Дона. 2015. № 4. URL: ivdon.ru/ru/magazine/archive/n4y2015/3483.

12. Игнатьев В.К., Станкевич Д.А. Оптимизация метода измерения угла фазового сдвига между двумя квазигармоническими сигналами // Инженерный вестник Дона. 2015. № 4. URL: ivdon.ru/ru/magazine/archive/n4p2y2015/3490/.

13. Ifeachor E.C., Jervis B.W. Digital Signal Processing. A Practical Approach. Prentice Hall. Second Edition. 2002. 960 p.

References

1. Asadov H., Ismaylov K. Positioning. 2011. № 2. pр. 61 - 64.

2. Stergiopoulos S. Proc. IEEE. 1998. Vol. 86. № 2. p. 358.

3. Autrey S.W. J. Acoust. Soc. Am. August 1988. Vol. 84, pp. 592 - 598.

4. Lasky M., Doolittle R. D., Simmons B.D., Lemon S.G. IEEE J. Oceanic Eng. 2004. Vol. 29. № 2. pp. 374 - 387.

5. Osche G.R. Optical detection theory for laser applications. Wiley-Interscience. Berlin. 2002. 418 p.

6. Kirichenko I.A. Inzenernyj vestnik Dona (Rus). 2014. № 4. URL: ivdon.ru/ru/magazine/archive/n4y2016/3845/.

7. Egorov K. D., Chesnokov S.S. Soviet Journal of Quantum Electronics. 1987. Vol. 17. № 6. pp. 808 - 810. doi: 10.1070/QE1987v017n06ABEH009365.

8. Sullivan E.J. IEEE Ocean. Eng. Soc. Newsletter, Winter 2003, pp. 21 - 24.

9. Voronin V.A., Tarasov S.P., Timoshenko V.I. Gidroakusticheskie parametricheskie sistemy. [Parametric sonar systems]. Rostov-na-Donu: Rostizdat. 2004. 400 р.

10. Mit'ko V.N., Essi-Esing A.S. Inzenernyj vestnik Dona (Rus). 2016. № 2. URL: ivdon.ru/ru/magazine/archive/n1y2016/3534/.

11. Kirichenko I.A., Kotov V. Yu, Lonkin P.V., Slutskii D.S. Inzenernyj vestnik Dona (Rus). 2015. № 4. URL: ivdon.ru/ru/magazine/archive/n4y2015/3483/.

12. Ignatjev V.K., Stankevich D.A. Inzenernyj vestnik Dona (Rus). 2015. № 4. URL: ivdon.ru/ru/magazine/archive/n4p2y2015/3490/.

13. Ifeachor E.C., Jervis B.W. Digital Signal Processing. A Practical Approach. Prentice Hall. Second Edition. 2002. 960 p.