УДК 621.396
СПОСОБ РАЗРЕШЕНИЯ ПО ДАЛЬНОСТИ ДВУХ СОСЕДНИХ СЖАТЫХ ЧМ СИГНАЛОВ НА ВЫХОДЕ ВЕСОВОЙ ОБРАБОТКИ
Иванников Анатолий Петрович
доктор технических наук, профессор; ведущий научный сотрудник ОАО «Федеральный научно-производственный центр «Нижегородский научно-исследовательский институт радиотехники».
Адрес: г. Н. Новгород, ул. Шапошникова, 5. Иванников Дмитрий Анатольевич кандидат технических наук, доцент; технический директор ОАО «Специальное конструкторское бюро радиоизмерительной аппаратуры».
E-mail: ccs5000@mail.ru. Адрес: г. Н. Новгород, пр. Гагарина, 174.
Аннотация: Известно, что уменьшение боковых лепестков автокорреляционных функций (АКФ) сигналов с ЛЧМ обеспечивают с помощью весовой обработки на базе функций Хемминга, Наталла, Чебышева, Тейлора и других, а также синтеза сигналов с нелинейной частотной модуляцией (НЧМ). Это, естественно, приводит к расширению главных лепестков АКФ ЛЧМ и НЧМ сигналов и перекрытию соседних сигналов во времени. В связи со значительным увеличением амплитудного динамического диапазона приемного тракта и потенциала современных РЛС происходит подавление слабых сигналов сильными, что ухудшает разрешающую способность РЛК по дальности. Поэтому устранение этого недостатка является актуальным. В работе проводится анализ известных способов уменьшения боковых лепестков АКФ импульсных сигналов с частотной модуляцией (ЧМ), на основание которых предлагается алгоритм разрешения по дальности двух соседних сжатых сигналов после весовой обработки. Оценка эффективности рассматриваемого способа проводится с помощью математического моделирования в среде MathCAD.
Ключевые слова: сигналы с линейной частотной модуляцией (ЛЧМ), сигналы с нелинейной частотной модуляцией (НЧМ), АКФ импульсных сигналов с ЛЧМ и НЧМ, разрешающая способность РЛК по дальности, функция неопределенности.
В настоящее время импульсные сигналы с частотной модуляцией (ЧМ) широко используются в РЛС, так как имеется возможность формировать и обрабатывать их быстродействующими цифровыми системами с использованием программных средств в реальном масштабе времени.
Одними из наиболее важных критериев их эффективности являются разрешающая способность по дальности и скорости. В общем случае для одновременного разрешения по каждому из этих параметров используется функция «разрешения» (неопределенности), предложенная Вудвордом Ф.М. [11].Кроме того, большой интерес представляет эффективность системы «передающий модуль фазированной решетки (ФАР) + согласованный фильтр (СФ)»в зависимости от вида излучаемого сигнала, режима работы передающего модуля и характеристик СФ.
На первом этапе проведем исследование АКФ и спектральных плотностей ЛЧМ и НЧМ сигнала после весовой обработки.
Автокорреляционная функция комплексной огибающей ЛЧМ сигнала с прямоугольным импульсом и нечетной симметрией закона модуляции частоты ю(0 относительно несущей частоты о)0 на выходе СФ в дискретном виде имеет вид [10]:
N—1
«лчм (т) = — ^ cos [Т Шл(пД0], (1)
п=0
где мл(пД£) =Д^-пД£;ти- длительность импуль-
ти
са; Дю - девиация частоты; At длительность выборки по времени.
Соотношение (1) является четной функцией типа sin(x)/x, которая приведена на рис. 1.
Особенностью АКФ излучаемых ЛЧМ сигналов является то, что после приема и об-
Рис.1. Нормированная АКФ ЛЧМ сигнала
работки их в согласованном фильтре амплитуда первого бокового лепестка АКФ ниже уровня главного на 13,2 дБ (0,217), а последующие боковые лепестки уменьшаются от лепестка к лепестку в соотношении 0,128; 0,091; 0,071; 0,058; 0,049; 0,042;0,037; 0,032; 0,028; 0,025 и т.д. Распределение мощности по времени между главным и боковыми лепестками составляет 1:1. То же имеет место у сигналов с ФКМ.
Известно, что уменьшение уровня боковых лепестков ЛЧМ сигналов осуществляется за счет применения весовой обработки. В качестве примера выбирается функция весовой обработки дискретном виде:
cos I-1 при nAt< —
w(nAt)=< 4 ти у 2 при m>1, (2)
^0 при пМ > ^
тогда автокорреляционная функция комплексной огибающей ЛЧМ сигнала (1) с учетом весовой обработки (ВО) определяется так:
Длчм+во(т) Zn=dw(nAi)cos [т ^л(пД0], (3)
В качестве исходных данных для расчета АКФ и спектральной плотности ЛЧМ и НЧМ сигналов выбраны следующие параметры: Дм= 2п 10б Гц - девиация частоты; ти =200 мкс; Дтг= 1 мкс - длительность сигнала на выходе СФ; Д t = 10_1мкс.
Нормированную спектральную плотность ЛЧМ сигнала в дискретном виде, используя соотношения (1) и (3), можно представить в виде:
SnJf) J? cos[т /чм(пДГ))] cos(f^t)dx.(4)
Известно [4], что cos(x)cos(y)= -[cos(x-
-y) + cos(x + y)], тогда при ю > 0 соотношение для нормированной спектральной плотности импульсного ЛЧМ сигнала имеет вид:
SяЧм(/) =Пл Ж sine [/ - /чм(пДС)]}, (5)
а
R2(x)
- 1х
10 5 - 7.5x10 6 -5*10 6 -2.5x10 6 0 2.5x10 6 5x10 6 7.5x10 6 1x1 т Рис.2. Нормированная АКФ ЛЧМ сигнала с весовой обработкой при m = 2
1 / \ 1
.....и ■л.
-2x10е —1.5 к -МО6 -З;-;^5 о 5ч:»5 1>ю(
Рис. 3. Спектральная плотность ЛЧМ сигнала с весовой обработкой (сплошная линия) при т = 2
где sinc(x)= —-—; пл - нормирующий коэффициент.
Нормированную спектральную плотность ЛЧМ сигнала после весовой обработки, используя соотношение (3) и (5), можно представить в виде:
£ячм+во(/)=Пл Щ w(nдt)sinc [/ - /чмСпАГ)]}(6)
На базе соотношений (1), (3), (5) и (6) созданы математические модели в среде MathCAD, с учетом приведенных исходных данных проведена оценка влияния весовой обработки на АКФ и спектральную плотность ^лчм+во(/) ЛЧМ сигнала. В качестве примера на рис. 2 и 3 показаны соответственно АКФ Ялчм+во (т) и спектральная плотность ^лчм+во(/) ЛЧМ сигнала после весовой обработки при т = 2.
Как и следовало ожидать, при уменьшении
боковых лепестков АКФ сигналов ЛЧМ за счет весовой обработки на 20 дБ увеличилась длительность главного лепестка почти в 2 раза, а форма спектральной плотности значительно отличается от прямоугольной, причем тем больше, чем больше подавление боковых лепестков.
Аналогично рассмотрим особенности характеристик НЧМ сигналов, синтезированных во временной области [10]:
МНчм(0 = О^лОО+а - а)шн0:) ПриН < (7) где а - весовой коэффициент, причем 0< а < 1; мл(1) -линейная часть закона изменения частоты; нелинейная часть закона изменения частоты.
Выражения автокорреляционных функций, а также спектральной плотности5Нчм (/) НЧМ
Рис. 5. Спектральная плотность НЧМ сигнала tg(^) при ц=1,35 и а =0,45
сигнала можно представить аналогично выражениям (1), (5) и (6) в виде [10]:
«Нчм (т) = cos [т Шнл(пД0], (8)
^нчм+Во(т) = -^Zn=o w(nДí)cos [т ^нл(пД£)],(9) 5н(/)=Пн Ж w(^t)sinc [/ - /нл(пДС)]}. (10)
На базе соотношений (8), (9) и (10) разработаны математические модели в среде MathCAD^ учетом приведенных исходных данных проведена оценка влияния весовой обработки на АКФ и спектральную плотность ^нчм (/) НЧМ сигнала.
В качестве примера на рис. 4 и 5 приведены соответственно АКФ йнчм (т) и спектральные плотности 5нчм(/) комплексных огибающих НЧМ сигнала, сформированного на базе нелинейной функции типа (ц).
Из приведенных графиков следует, что при уменьшении уровня боковых лепестков АКФ сигналов с ЧМ главный лепесток расширяется во времени и происходит искажение формы огибающей спектральной плотности от прямоугольной. Очевидно, что применение ЛЧМ и НЧМ сигналов после весовой обработки в качестве зондирующих сигналов может привести к снижению КПД передающего модуля ФАР (работающего в большинстве случаях в ключевом режиме) и дополнительным потерям в СФ ввиду отличия формы огибающей отраженного ЧМ сигнала от прямоугольной.
На втором этапе рассмотрим возможность повышения разрешающей способности ЧМ сигналов по дальности с учетом выше приведенных результатов. При этом предполагается:
■ огибающие АКФ сигналов с ЧМ на выходе фильтров с весовой обработкой, определяемые соотношениями (6) и (10),полностью известны, что позволяет создать эталонные АКФ ЧМ сигналов;
■ два сжатых соседних ЧМ сигналов, сдвинутых относительно друг от друга на время, равное их длительностит1=1/Д£, ортогональны во времени.
Рассмотрим предложенный способ на примере ЛЧМ сигналов. Для различия двух сигналов щф и и2(^т) предлагается использование дискретных трансверсальных фильтров с импульсными характеристиками и которые согласованы с огибающими эталонной АКФ ЛЧМ сигналов
Мк) =^=<,"(^0^ [т шл(&Д0], (11) ^(М^Х^о^ - я)Дф05 [т 0>л[(£ - П)Д£]]. (12)
Оптимальная схема различения двух сигналов с использованием СФ приведена на рис.6.
Обобщенное уравнение для дискретной фильтрации с конечной импульсной характеристикой (КИХ) можно представить в виде:
у^Х^ЩМп-^, (13)
где х(п) и у(п) - соответственно сигналы на входе и вход и выходе фильтра, - импуль-
сная характеристика дискретного фильтра, ^ коэффициенты фильтра, «=0,1,2, 3,4, ....
Дискретный фильтр 1
hi(t) - «i(k) ut(n)
2(П)
h2(fc — n) = u2(k — n) u2(k
Дискретный фильтр 2
Рис.6. Оптимальная схема для различения двух сигналов
Для двух данных последовательностей х(«) и у(п) длины N оценка их взаимной корреляции равна
^Я^МЮуа + п) при п = 0,1, 2,3,4,5 ^
Кху(")= 1™1 (14)
У (. ^Ж?"1*» —п)у(к) при п = 0,-1, —2 4 7
Кху(0)=^Хк="оП_1^(^)2 при п = 0,1,2... (15)
На рис.7 показана блок-схема дискретного трансверсального фильтра, реализующего рассматриваемый метод.
Практическая реализация рассматриваемого метода заключается в обеспечении свертки АКФ сигнала, отраженного от цели, с импульсной характеристикой дискретного трансвер-сального фильтра во временных координатах ^ф(я), весовые коэффициенты и длительность которого определяются параметрами эталонной АКФ:
и1(п) = ^х^х^о^е^-поо, (16)
и2(п — !Т1) = ^-¡ЙЙ Й=0т«(т — ^К^Сг — iT1), (17)
где тп(т)= w(nДí) cos[т млпД^|;
тк_п (т) = w[(fc — п)Д£] cos[т ыл(к — п)Д£];
тп(т — 1тх) =w(nДt) cos[(т — 1т1)^лпД£]; тк_п(т — — я)Д£] соб[(т — 1г1) ыл(& — и)Д(:];
- нормирующий коэффициент.
Итак, практическая реализация предложенного метода возможна с помощью трансвер-сальных дискретных фильтров с параметрами:
- весовые коэффициенты соответствуют значениям огибающих АКФ ЛЧМ сигналов, отраженных от целей, с учетом весовой обработки;
- количество дискретных отсчетов по дальности N в оперативной памяти выбирается в пределах длительности излучаемого импульса;
- обеспечиваются свертки АКФ сигналов, отраженных от цели, с импульсной характеристикой дискретного трансверсального фильтра во временных координатах ^ф(я);
- ортогональные во времени ЛЧМ сигналы с центральных выходов дискретных фильтров поступают на входы двух ортогональных каналов;
- информация на входе дискретного фильтра синхронизируется с разверткой по дальности РЛС.
Учитывая, что разрешаемые ЛЧМ сигналы имеют одинаковые эталонные АКФ, можно реализовать рассматриваемый метод с помощью скользящего окна, состоящего из транс-версального дискретного фильтра, информация с которого поступает на два ортогональных канала: четные на первый канал, а нечетные -на второй канал (рис.8).
Известно, что свертка во временных координатах, эквивалентная умножению в частотной области, обеспечивает увеличение от-
Дискретный фильтр 1 Управ 11 I I I I I I I I ■! ,1111.......—н, Коэффициенты
т
Pl
Управ
Pi
Регистры
Управ
П ÍI I I I I I I I I I I I I I I I I гп
И:' :,ixi и1 —>
Рис.7. Структурная схема дискретного фильтра одного канала
Дискретный фильтр Управ
Рис. 8. Структурная схема дискретного фильтра со скользящим окном
ношения сигнал/шум. Алгоритм вычисления нормированной спектральной плотности ЛЧМ сигналов на выходе трансверсального дискретного фильтра, используя соотношение (5), (16) и (17), можно представить в виде:
S(f)=Чs Ж Ж w(nДt)w[ВД sinc {|[/ -
/чм (пДО]} йпс [/ - /чм (кДС)]}, (18) где нормирующий коэффициент.
На базе соотношений (16), (17) и (18) созданы математические модели в среде MathCAD, с помощью которых возможно провести анализ параметров ЧМ сигналов и оценку эффективности предлагаемого способа.
В качестве примера выбраны следующие исходные данные: Дм= 2п 10б Гц - девиация частоты; ти =200 мкс; А тг = 1 мкс - длительность сигнала на выходе СФ; А t = 10_1мкс; т =1 и 2мкс - относительное смещение двух сигналов во времени.
С учетом приведенных выше исходных данных на рис. 9 и рис.10 показаны АКФ комплексной огибающей ЛЧМ сигнала с ВО при т = 2 и его спектральная плотность.
С учетом приведенных выше исходных данных на рис.11-20 показано разрешение двух ЛЧМ сигналов, сдвинутых относительно друг от друга на время, равное их длительности т1=1/Дf (т=2) и имеющих разные амплитуды.
На рис.11-15приведены соответственно сигналы 1 и 2 на входе дискретного фильтра при т=2 и т=1мкс, их смесь на входе дискретного фильтра и сигналы 1 и 2 на выходе дискретного фильтра в каналах 1 и 2.
На рис. 16-20 приведены соответственно сигналы 1 и 2 на входе дискретного фильтра при т=2 и т=2 мкс, их смесь на входе дискретного фильтра и сигналы 1 и 2 на выходе дискретного фильтра в каналах 1 и 2.
В качестве примера выбраны следующие исходные данные:
Дм= 2п 10б Гц - девиация частоты; ти =200 мкс; Атг = 1 мкс - длительность сигнала на выходе СФ; Аt = 10_1мкс; т = 1 и 2 мкс - относительное смещение двух сигналов во времени.
С учетом приведенных выше исходных данных на рис.9 и рис.10 показаны АКФ комплексной огибающей ЛЧМ сигнала с ВО при т=2 и его спектральная плотность.
С учетом приведенных выше исходных данных на рис. 11-20 показано разрешение двух ЛЧМ сигналов, сдвинутых относительно друг от друга на время, равное их длительности т1=1/Д/" (т=2) и имеющих разные амплитуды.
На рис. 11-15 приведены соответственно сигналы 1 и 2 на входе дискретного фильтра при т=2 и т=1 мкс, их смесь на входе дискретного фильтра и сигналы 1 и 2 на выходе дискретного фильтра в каналах 1 и 2.
На рис.16-20 приведены соответственно сигналы 1 и 2 на входе дискретного фильтра при т=2 и т=2 мкс, их смесь на входе дискретного фильтра и сигналы 1 и 2 на выходе дискретного фильтра в каналах 1 и 2.
При наличии смеси сигнала и шума с гауссовым распределением вероятность полной ошибки обнаружения двух объектов определяется как
Рис.9. АКФ комплексной огибающей ЛЧМ с ВО при т=2
5:; г .рЭ
/
¡1.
-5 101
5 1С1
|ч(0|
Рис.10. Спектральная плотность ЛЧМ сигнала с ВО при т=2
- 1x10 : -7.5x10 6 -5x13 6 -2.5x13 6
2.5x10
-6
5x13
-6
7.5x10
Рис.11. Сигнал 1 на входе дискретного фильтра при т =2и т=1мкс
кз О)
0.6
0.4
ч
1x13
1 8
| 0 .6'
я 0.4"
/ 0.2^
- 1хЮ
Рис.12. Сигнал 2 на входе дискретного фильтра при т=2 и т=0
=Щис1 )И/лт(ис1 ) +Ж(ис2 )И^т(ис2 ),(19) где Щис1),Щис2 ) - вероятность правильного
обнаружения 1 и 2 объектов, И^т(ис1 ),
И/лт(ис1 ) - вероятность ложного обнаружения 1 и 2 объектов.
Для априорно равновероятных сигналов вероятность полной ошибки обнаружения двух объектов при наличие шума равна [12]:
и/0ш =1 - Ф(7(Е/^)(1 — г12), (20) где Е/Л^0 - отношение сигнал/шум, г12 - коэффициент взаимной корреляции между сигналами, Ф(х) - интеграл вероятности.
R4(T)
/0.8
/ 0.6
/ 0.4
/ 0.2
ImmI
Рис. 13. Смесь сигналов 1 и 2 на входе дискретного фильтра
1
0.£
0.- /О
г
-1x10 " — 7.5 10 5 - 5 10 5 -2.5<10 5
2_sao
,-в
S<10
7.SilO й 1x10 "
Рис.14. Сигнал 1 на выходе дискретного фильтра при т=2 и т=1 мкс в канале 1
\
Ю-8" \
1 0.6" \
( 0.4" \
г 0.2" V
R3M
2.5x10
- S
IricoI
Рис.15. Сигнал 2 на выходе дискретного фильтра при т=2 и т=0 в канале 2
0.&
0.4
/ 0.2
-2.5x10
2.5x10
1x10
Рис. 16. Сигнал 1 на входе дискретного фильтра при т=2 и т=2 мкс
Например, при 2Е/М0=20 и 2Е/М0=8, г12=0 вероятность Wош соответственно равна 3 * 10"2 и 10_3, что приемлемо на практике.
Выводы
1. При выборе уровня боковых лепестков АКФ излучаемых ЧМ сигналов необходимо учитывать влияние весовой обработки на КПД
ЕЗ Ст)
\
\
\
/ 0.4" \
V
Рис.17. Сигнал 2 на входе дискретного фильтра при т=2
114(0
1
ш 0.6"
/ 0.4"
¿г 0.2
2.5x10
Рис.18. Смесь сигналов на входе дискретного фильтра при т=2 и т=2 мкс
|К2(0 ист
0.4"
—7^*4
- IX 13
- 7.5x10
-5x13 " -2.5x15
2.5x13
5x10
7.5x10
Рис.19. Сигнал 1 на входе дискретного фильтра при т=2 и т=2 мкс в канале 1
из (О
2.5x10
- 6
1х 13'
\
р.%- \
I 0.6" \
/ 0.4" \
ё 0.2" \
Рис.20. Сигнал 2 на входе дискретного фильтра при т=2 и т=2 мкс в канале 2
транзисторных усилительных модулей ФАР и потери в СФ.
2.Предложенный способ оптимизации разрешающей способности РЛС по дальности является эффективным, что подтверждено его макетированием.
Литература
1.Ширман Я.Д. Основы теории обнаружения радиолокационных сигналов и измерение их пара-метров.-М: Сов. Радио,1963. - 416с.
2. Варакин Л.Е. Теория сложных сигналов.-М.Сов. Радио.1970.-372
3. Кук Ч., Бернфельд М. Радиолокационные сигналы. - М. Сов. Радио.1974, 568с.
4. Современная радиолокация. Анализ, расчет и проектирование систем. Перевод с английского под редакцией Кобзарева Ю.Б, изд-во «Советское радио», 1969. 702с.
5. Френкс Л. Теория сигналов.- М: Сов. Ра-дио,1974.-344 с.
Поступила 05 июля 2016 г.
6.Слока В.К. Вопросы обработки радиолокационных сигналов-М: Сов. Радио. 1970, -256 с.
7. Вакман В.Е. Сложные сигналы и принцип неопределенности, Изд. Сов. Радио. 1965.-303с.
8.Оконечников В.С., Кочемасов В.Н. Сжатие частотно-модулированных сигналов с наибольшим произведением девиации частоты на длительность импульса, Зарубежная электроника, №1, 1987,с. 82-93.
9. Денисенко А.Н. Сигналы. Теоретическая радио - техника. Справочное пособие.- М: Горячая линия-Телеком, 2005.-704с
10. Иванников А.П., Иванников Д.А. Метод синтеза импульсных сигналов с нелинейной частотной модуляцией \\ Радиотехнические и телекоммуникационные системы, 2015, №2.с.11-17.
11.Вудворд Ф.М. Теория вероятностей и теория информации с приложениями в радиолокации. Пер. с англ., под редакцией Г.С. Горелика. Изд. «Советское радио», 1955.-126с.
12.Тихонов В.И., Харисов В.Н. Статистический анализ и синтез радиотехнических устройств и систем. Учебное пособие для вузов.- 3-е изд.- М: Горячая линия.-Телеком, 2015.-608 с.
English
Range resolution method of two adjacent compressed FM signals in the weighting output
Anatoliy Petrovich Ivannikov - Doctor of Engineering, Professor, Leading Researcher Research and Development Center JSC Nizhny Novgorod Radio Engineering Research Institute.
Address: N. Novgorod, Shaposhnikov Str., 5.
Dmitry Anatolyevich Ivannikov - Candidate of Technical Sciences, Associate Professor Technical Director of JSC Special Design Office of Radio- Measuring Devices.
E-mail: ccs5000@mail.ru.
Address: N. Novgorod, Gagarin Ave., 174.
Abstract: It is known that reduction of side lobes of LFM signal autocorrelation functions (ACF) is provided by weighting based on functions of Hamming, Nuttal, Chebyshev, Taylor and others, and also nonlinear frequency modulation (NFM) signal synthesis. This certainly involves the main lobes expansion of LFM and NFM signal autocorrelation functions (AKF) and overlapping of the adjacent signals in time. Suppression of weak signals by strong ones occurs due to significant increase in amplitude dynamic range of receive path and due to today's radar capability that worsens range resolution capacity of radar facility. Therefore, elimination of this shortcoming is vital. The paper work gives analysis of known ways of reduction of side lobes of side lobes of FM signal autocorrelation functions (ACF). It is shown that application of weighting provides side lobes reduction of FM signal autocorrelation functions with simultaneous main lobe expansion in time. Besides, if using probing LFM and NFM signals after weighting it leads to a non-rectangular distortion of spectrum envelope that can cause decrease in efficiency factor of transmitting phased antenna array operating mainly in key mode. To avoid the specified shortcomings the range resolution algorithm of two adjacent compressed FM signals in the weighting output is proposed and the algorithm efficiency is confirmed by MathCAD mathematical modeling.
Key words: linear frequency modulated signals, nonlinear frequency modulated signals, AKF of LFM and NFM pulsed signals, radar facility range resolution, uncertainty function.
References
1. Shirman Ya.D. Fundamentals of radar signal detection and their parameters measurement. - M: Sov. Radio, 1963. - 416 p.
2. Varakin L.E. Aggregate signal theory. - M. Sov. Radio.1970.-372 p.
3. Cook Ch., Bernfeld M. Radar signals. - M. Sov. Radio.1974, 568 p.
4. Modern Radar. Analysis, Evaluation, and System Design. Translated from English. Edited by Kobzarev Yu.B., publishing house Sovetskoye radio, 1969. 702 p.
5. Franks L. Signal Theory. - M: Sov. Radio, 1974. - 344 p.
6. Sloka V. K. Aspects of radar signal processing -M: Sov. Radio. 1970, - 256 p.
7. Vakman V. E. Aggregate signals and uncertainty principle, Publ. Sov. Radio. 1965. - 303 p.
8. Okonechnikov V. S., Kochemasov V. N. Compression of frequency modulated signals with utmost frequency deviation for pulse width, Zarubezhnaya elektronika, No. 1, 1987, p. 82-93.
9. Denisenko A.N. Signals. Radio technology in theory. Reference book. - M: Goryachaya liniya -Telecom, 2005. - 704 p.
10. Ivannikov A.P., Ivannikov D. A. Synthesis method of pulse nonlinear frequency modulated signals \\ Radio-tekhnicheskiye i telekommunikatsionnye sistemy, 2015, No. 2. p. 11-17.
11. Woodward F.M. Probability and Information Theory, with Applications to Radar. Transl. from Engl., Ed. by G.S. Gorelik. Publ.h. Sovetskoye radio, 1955. - 126 p.
12. Tikhonov V. I., Harisov V. N. Statistical analysis and synthesis of radio engineering devices and systems. Manual for higher education institutions. - 3rd iss. - M: Goryachaya liniya. - Telecom, 2015. - 608 p.