CC BY
493
УДК 621.391.26 Дата подачи статьи: 28.02.2014
МЕТОДИКА ОПРЕДЕЛЕНИЯ ПОКАЗАТЕЛЕЙ КАЧЕСТВА ОБНАРУЖЕНИЯ РАДИОЛОКАЦИОННЫХ СТАНЦИЙ
А.Ю. Платонов, начальник группы; О.Н. Помазуев, зам. начальника управления (Главное управление научно-исследовательской деятельности МО РФ, ул. Профсоюзная, 84/32, г. Москва, 117997, Россия, [email protected], [email protected]); А.М. Миронов, референт зам. министра обороны РФ (МО РФ, ул. Знаменка, 19, г. Москва, 119160, Россия, [email protected]); Ф.Н. Абу-Абед, к.т.н., доцент (Тверской государственный технический университет, наб. Аф. Никитина, 22, г. Тверь, 170026, Россия, [email protected]); А.Р. Рюмшин, к.т.н., с.н.с., старший научный сотрудник (НИИ «Центрпрограммсистем», просп. 50 лет Октября, 3а, г. Тверь, 1 70024, Россия,
rvdopira@yandex. ru)
В настоящее время из-за широкого использования технологий снижения заметности излучаемый радиолокационный сигнал существенно отражается и поглощается. Таким образом, значительно уменьшается мощность поступающего на вход приемного устройства РЛС отраженного сигнала, что затрудняет обнаружение самолета либо другого объекта, выполненного по подобной технологии. В этих условиях сложно оценить показатели качества обнаружения малозаметных и малоразмерных целей.
Во многих случаях амплитуда и фаза сигналов РЛС являются случайными неизвестными параметрами. В условиях достаточно быстрых по отношению к тактовой частоте сигналов флуктуаций амплитуды и фазы, когда отсутствует возможность определения их значений по результатам предварительного статистического оценивания на предыдущих интервалах, часто используют модели сигналов.
В работе предложена методика, позволяющая определять показатели качества обнаружения РЛС для широкого класса моделей сигналов, в том числе малозаметных и малоразмерных целей, в условиях стационарных гауссовских, шумовых импульсных помех, а также в беспомеховой обстановке. Методика использует полученные аналитические выражения и аппроксимации распределений рядом Грама-Шарлье в приближении Эджворта. Значение порога обнаружения для заданных вероятностей ложной тревоги определяется методом Ньютона, что позволяет повысить точность и упростить вычисления. Получены оценки показателей качества обнаружения в рассматриваемых условиях. Разработанная методика и полученные оценки могут быть использованы для определения показателей качества систем более высокого иерархического уровня, например показателей вторичной обработки сигналов, вероятности обнаружения (захвата) сигнала в зоне поиска.
Ключевые слова: показатели качества обнаружения, модели сигнала, накопление пачки импульсов, обнаружитель сигналов, вероятности правильного обнаружения и ложной тревоги.
THE METHOD OF QUALITY FACTORS ESTIMATION FOR RADAR DETECTION Platonov A Yu., Chief of Group; Pomazuev O.N. , Deputy Chief of Department (Central administrative board of research activity the Ministry of Defence the Russian Federation, Profsojznaya St. 84/32, Moscow, 117997, Russia, [email protected], [email protected]);
Mironov A M., Deputy Assistant of Defense of the Russian Federation (Ministry of Defense, Znamianka St. 19, Moscow, 119160, Russia, [email protected]);
Abu-Abed F.N, Ph.D. (Engineering), Associate Professor (Tver State Technical University, Quay Nikitin 22, Tver,170026, Russian Federation, [email protected]);
RyumshinAR., Ph.D. (Engineering), Senior Researcher (Research Institute «Centerprogramsystem», 50 let Octyabrya Av. 3a, Tver, 170024, Russia, [email protected])
Received 28.02.2014
Abstract. Nowadays stealth technologies are widely used. That is why the radiated radar-tracking signal is essentially reflected and absorbed. Thus, the power of a reflected signal on an intake radar entry is decreased. It complicates airplane detection, or other object created according to similar technology. Under these conditions it is difficult to estimate quality factors of radar detection of stealth and small targets.
In many cases the radar signals amplitude and a phase is unknown parametres. The models of signals are used often when there are fast (in relation to clock frequency) fluctuations amplitude signals and a phase when there is no possibility to estimate their values according to the results of preliminary statistical estimation on the previous intervals.
We propose a technique which allows determining the quality of radar detection for a wide class of models of signals (including stealth and small targets) in terms of stationary Gaussian noise, impulse noise and in non-jamming environment. The technique uses the obtained analytical expressions and approximation of distributions with Gram-Charlier series in Edge-worth Expansion. Detection threshold for given probabilities of false alarm is defined by Newton's method. It improves the accuracy and simplifies the calculations. The article shows the estimates for quality detection under these conditions. The developed methodology and the resulting estimates can be used to determine quality factors of higher hierarchical level systems. For example, indicators of secondary signal processing, probability of signal detection (capture) in the search area.
Keywords: detection quality factors, signal model, accumulation of a pack of impulses, signals detector, correct detection and false alarm probabilities.
Развитие технологий снижения заметности [1, 2] привело к тому, что излучаемый радиолокационный сигнал существенно отражается и поглощается. Таким образом, значительно уменьшается мощность поступающего на вход приемного устройства радиолокационной станции (РЛС) отраженного сигнала, что затрудняет обнаружение самолета либо другого объекта, выполненного по подобной технологии. В этих условиях сложно оценить показатели качества обнаружения малозаметных и малоразмерных целей.
Необходимо разработать методику М: X х У ^ ^•{Д позволяющую определять показатели качества обнаружения РЛС для широкого класса моделей сигналов, в том числе малозаметных и
малоразмерных целей. X - вектор характеристик принимаемых сигналов, определяемый целевой и
помеховой обстановками; У - вектор характеристик используемых обнаружителей в РЛС; {В,
- вероятности правильного и ложного обнаружения.
Во многих случаях амплитуда и фаза сигналов РЛС являются случайными неизвестными параметрами. В условиях достаточно быстрых по отношению к тактовой частоте сигналов флуктуаций амплитуды и фазы, когда отсутствует возможность определения их значений по результатам предварительного статистического оценивания на предыдущих интервалах, часто используют модели сигналов. Характеристики принимаемых сигналов задаются релеевским распределением, которое правильно описывает отраженный сигнал от целей с большим числом случайно расположенных независимых отражающих элементов, а также пятью моделями Сверлинга [3]:
- первая соответствует цели, которую можно представить как совокупность независимых случайных рассеивателей при флуктуации от цикла к циклу сканирования;
- вторая - модель флуктуирующей цели от импульса к импульсу, которую можно представить как совокупность независимых случайных рассеивателей;
- третья соответствует цели, которую можно свести к одному большому и прочим малым отражателям (малозаметные и малоразмерные цели на фоне отражений) при флуктуации от цикла к циклу сканирования;
- четвертая - модель флуктуирующей цели от импульса к импульсу, которую можно представить как один большой и прочие малые отражатели;
- пятая соответствует цели с одним рассеивающим элементом в отсутствие флуктуаций.
Накопление нескольких реализаций, соответствующих отражениям от одной совокупности целей, существенно улучшает характеристики РЛС. Обработка множества независимых наблюдений значительно увеличивает вероятность обнаружения цели на фоне мешающих отражений. Определим показатели качества обнаружения (вероятности правильного обнаружения D и ложной тревоги —) для РЛС, использующих когерентное обнаружение сигнала в квадратурах и некогерентные обнаружители.
Рассмотрим когерентное обнаружение пачки М принимаемых импульсов сигнала, флуктуирующего по закону Релея, в квадратурах на выходе линейного детектора в условиях стационарных шумовых гауссовских помех. Показатели качества обнаружения можно определить следующим образом: D = ехр(- S¡/(2M (1 + и2) + М У), (1)
^ = ехр(- Б20/ (2М (1 + и2)), (2)
2 2
где и , q - нормированные к мощности внутрен-
2
него шума аш значения мощности помехи и сигнала соответственно; 50 - значение порога обнаружения, нормированное к стш [4].
Во многих случаях отражения от целей могут иметь протяженный характер, тогда их обнаружение может быть сведено к модели многоканального некогерентного весового обнаружителя сигналов различной мощности. Рассматривается О-канальный весовой обнаружитель с весовыми коэффициентами Ъ,г, осуществляющий в каждом
,^-м канале (g = 1, G) некогерентное (на выходе квадратичного детектора) накопление пачки из Я импульсов различной мощности ст^, г = 1, R . Каналы этого обнаружителя идентичны (дисперсия
2
внутренних шумов аш во всех каналах одинакова), но отличаются весами Ъ,г. Обнаружение цели осуществляется в условиях помех, представляющих собой стационарный гауссовский шум посто-
2
янной мощности стП^..
Для заданных условий и флуктуирующего по закону Релея сигнала выражения для определения
показателей качества обнаружения [5] имеют следующий вид:
л=а £ £ ехр(2)
__О Я '
1 "П П (А -«г)
1=1 }=1 V* &
F = D
(3)
(4)
а * а;
,/ = 1,О и г, у = 1, Я,
где Д = ПП а
Я=1 г=1
а = В2 /2Ь
gг Pgг/ е
2
; Р^- = У(1);
9 9/9 9 9/9
Цт =стпг/ ; =<5^! стш ; 2 - порог обнаружения.
На рисунке 1 приведены кривые обнаружения D=f(qос) при фиксированном значении -=10-8, где qос - ослабление сигнала за счет дальности до объекта локации. Характеристики получены для следующих исходных данных: 1 2 1,5
=
4,5 50 10 2,5 4 3,3
2 Л
1 + ^
2 2
кривая 1 - Ь =
V " У
кривая 2 - Ъ,= 1; кривая 3 - Ъ,= О' ;
К 2 Л ' с
1 + -
кривая 4 - Ьёг = У
2
На практике часто приходится сталкиваться с ситуацией, когда помеха представляет собой случайный стационарный поток взаимно неперекрывающихся шумовых радиоимпульсов постоянной
2
мощности ст^, (g = 1, G, г = 1, R), в общем случае
различной для каждого из каналов, с гауссовским распределением амплитуды напряжения внутри импульса, заданный средними значениями длительности тП и ТП периода следования импульсов шумовой импульсной помехи. При этом обеспечивается А/п~А/с, где А/щс) - ширина спектра помехи (сигнала) по уровню половинной мощности.
Для сигнала, флуктуирующего по закону Ре-лея, и рассматриваемых условий показатели качества обнаружения имеют следующий вид:
а я
ПП К -Р*)
D--
ПП а
F = D
2 G R
Ъ Ъ Ъ
k=1 ¡=1 j=1
exp
(-zPkj)
Pk
q„ = 0 ^
п п п (p-g --pkj)
(5)
1=1 g =1 r=1
lg rikij
agr ф atJ, Vl, k = 1,2 , g, i = 1, G и r, j = 1, R,
где а2gr + (1 p2gr,
Pgr 1
Pigr =
2bgr (1 + 2)
P2 gr =
1grr 2gr
a,
P2 =■
sr 1 +u2„
(1 + 2)'
ц = тп /Тп - вероятность появления импульса помехи в элементе разрешения обнаружителя [6].
На рисунке 2 приведены кривые обнаружения D=f(qос) при фиксированном значении
^=10-4. Характеристики получены для следующих исходных 1 2 3
qgr
4 5 6 7 8 9
(сплошные линии) и
qgr =
1 2 1,5 4,5 50 10 2,5 4 3,5
(штриховые линии) при по-
стоянной нормированной средней мощности помехи та=и2ц=1. Проведенный анализ эффективности шумовых импульсных помех [6] показывает, что в условиях таких помех качество функционирования обнаружителя снижается по сравнению со стационарными шумовыми помехами, это происходит вследствие увеличения вероятности ложных тревог. Существует оптимальное значение ц0, соответствующее наихудшему качеству обнаружения полезного сигнала, то есть минимуму D для фиксированных К Кривые 1 и 2 на рисунке 2 отображают воздействие стационарных шумовых (ц=1) и шумовых импульсных (ц=ц0) помех соответственно.
Рассмотрим обнаружение других моделей сигнала в беспомеховой обстановке.
D
0.8
0.6
0.4
0.2
/
1 1 I 2
1 1
// 1
Jj /
10
15
20
25
30 35
q , дБ
Рис. 2. Кривые обнаружения для релеевской модели сигнала при фиксированной F=10'4
Fig. 2. Detection curves for Rayleigh signal model when F=10'4 is fixed
Плотность распределения вероятности случайного напряжения сигнала А>0, представленного пятой моделью Сверлинга, в случае некогерентного накопления пачки N импульсов на выходе квадратичного детектора в условиях отсутствия помех имеет следующий вид:
ч( N -I)! 2
р (U)=
2 U
v NAV аш ,
х ехр
-U -
NA
с
(6)
A
2NU ✓
V I ш у
где /А-1 (•) - модифицированная функция Бесселя N-1 порядка [7].
Точных выражений для определения показателей качества обнаружения в этом случае получить не удается, значение вероятности правильного обнаружения можно определить при использовании в качестве аппроксимации распределения ряд Грама-Шарлье в приближении Эджворта [8, 9]. Тогда вероятность правильного обнаружения может быть определена следующим образом:
D =
erfc (у/42) exp (-V 2/ 2) 2
[с (У2 -1)+
+C4V (V2 - з) - C6V (V4-10V2 +15)] , (7) Q+13
и • Q =
2a2
где С3 = -
VN (2Q +1)
С = ^О+Ь ; Сб = C3/2 ; V = ^^ ;
N (2Q +1) B N (2Q + 1); erfc ( z ) = 1- erf ( z ), 2 z 2
erf (z) = —= i e-9 dф - функция ошибок; ST - по-
vn о рог обнаружения.
0
5
-5
0
g=1 r=1
g=1 r=1
X
I
N-1
Для рассматриваемого случая вероятность ложной тревоги можно определить, используя следующее выражение:
F = 1 - Г,
Sj
N -1
(8)
где Г (•) - неполная гамма-функция [7].
Используя метод Ньютона, можно получить значение порога обнаружения для заданной
О (^, т-! ) ,
ST = ST , ■
T, m T, m-1
m = 1, 2,3,
(9)
О (^, )
где О (, т) = 0,5^ - Г7 (, N);
(^ . _ ехр () ^(10)
^) _ (N -1)! ' (10)
пР - количество ложных тревог [8].
Начальное значение порога обнаружения [7] определяется как
^ = N+ у—ЁР -1). (11)
Количество итераций определяется следующим выражением:
- ^< ^/ 104 . (12)
На рисунке 3 представлены кривые обнаружения D=f(Q) при фиксированном значении ^=10-9 для рассматриваемого случая, кривая 1 соответствует N=1, а кривая 2 - N=10.
Плотность распределения вероятности случайного напряжения сигнала, представленного четвертой моделью Сверлинга, в случае некогерентного накопления пачки N импульсов на выходе квадратичного детектора в условиях отсутствия помех имеет вид
* (U ) =
UN-' exp(-U/(1 + Q/2)) ■ (1 + Qf (n -1)! 1
-N, N;--QU
1 + Q
(13)
где 1Е1 (•) - вырожденная гипергеометрическая
функция [7].
Для случая, когда N<50, значение вероятности правильного обнаружения [7] определяется как
' |) N * +(| V ^ ■
D =1 - [ Т. + ( Q )
2!
-Y 2 +.
Y N
где у,. = Г
1+Q
St
(14)
, N + i
I + Ql 2'
Для случая, когда N>50 точных выражений для определения показателей качества обнаружения получить не удается, значение вероятности правильного обнаружения можно определить, используя в качестве аппроксимации распределения ряд Грама-Шарлье в приближении Эджворта [8, 9], выражением (7), коэффициенты в котором вычисляются следующим образом:
1 2Р3-1 1 2Р4-1
Сз =
3N (2в2-if 4N (2в2-1)2
B = yjN (2p2 -1); p = 1 + Q.
C2
С = —; в = 6 2 ;
1 = 1+—.
2
Значение порога обнаружения для заданной ^ можно получить, используя выражения (9)-(12).
На рисунке 4 представлены кривые обнаружения D=f(Q) при фиксированном значении ^=10-9 для рассматриваемого случая, кривая 1 соответствует N=1, кривая 2 - N=10, кривая 3 - N=50, а кривая 4 - N=100.
D
0.8
0.6
0.4
0.2
У! I 2
[ ' 1
г \ 4 1
11 1 Л
1 п
Q^
-10
10
15
20
25
30
Рис. 4. Кривые обнаружения для сигнала, соответствующего четвертой модели Сверлинга при F=10'9
Fig. 4. Detection curves for a signal corresponding to the forth Swerling model when F=10'9 is fixed
Плотность распределения вероятности случайного напряжения сигнала, представленного третьей моделью Сверлинга, в случае некогерентного накопления пачки N импульсов на выходе квадра-
- N
2
2
0
5
5
0
тичного детектора в условиях отсутствия помех [7] имеет вид
f \
Р (U) =
= UN-1 exp (-U)
1 + NQ |(N -1)!
i F
2, N
U
1 + -
NQ
.(15)
В этом случае для N<2 значение вероятности правильного обнаружения [7] определяется следующим образом:
D = 11 +
exP I
2
NQ ST
1+
St , -—( N - 2 ) 1 + NQ/2 NQV ;
(16)
1 + N<2/
Для случая, когда N>2, выражение для определения вероятности правильного обнаружения [7] будет иметь вид
Л = , ^ехр (\ +1 - г (^, N -1) +
(1 + NQ/2)(N - 2)!
+ |1+Л| .11 +
NQ J
S , - — (N - 2) 1 + NQ/2 NQV ;
x exp
St
Г
1 + NQ/ 2 ) у 1 + 2/NQ
S
, N-1 I.
(17)
Значение порога обнаружения для заданной ^ можно получить, используя выражения (9)-(12).
На рисунке 5 представлены кривые обнаружения D=f(Q) для рассматриваемого случая при фиксированном значении ^=10-9, кривая 1 соответствует N=1, кривая 2 - N=10, кривая 3 - N=50, а кривая 4 - N=100.
D
0.8
0.6
0.4
0.2
2
/ / -У1 / / ] / 1
/ • г
11 /
/
QдБ
-10 -5
10
15
20 25
30
Рис. 5. Кривые обнаружения для сигнала, соответствующего третьей модели Сверлинга при F=10'9
Fig. 5. Detection curves for a signal corresponding to the third Swerling model when F=10'9 is fixed
Плотность распределения вероятности случайного напряжения сигнала, представленного второй моделью Сверлинга, в случае некогерентного на-
копления пачки N импульсов на выходе квадратичного детектора в условиях отсутствия помех [7] имеет вид
Р (U) =
1
(1 + Q)N (N -1)!
U N exp
U 1+Qy
(18)
Для случая, когда N<50, значение вероятности правильного обнаружения [7] определяется следующим образом:
D = 1 - Г,
Sr
.1+Q
, N
(19)
Для случая, когда N>50 точных выражений для определения показателей качества обнаружения получить не удается, значение вероятности правильного обнаружения можно определить, используя в качестве аппроксимации распределения ряд Грама-Шарлье в приближении Эджворта [8, 9], выражением (7), коэффициенты в котором вычисляются следующим образом:
С =-
1 1 С1
—С4 = —; С6 = ; В = (1 + 0).
Ъ^Ы 4 4Ы' 6 2 ' У
Значение порога обнаружения для заданной ^ можно получить, используя выражения (9)-(12).
На рисунке 6 представлены кривые обнаружения D=f(Q) для рассматриваемого случая при фиксированном значении ^=10-9, кривая 1 соответствует N=1, кривая 2 - N=10, кривая 3 - N=50, а кривая 4 - N=100.
Плотность распределения вероятности случайного напряжения сигнала, представленного первой моделью Сверлинга, в случае некогерентного накопления пачки N импульсов на выходе квадратичного детектора в условиях отсутствия помех [7] имеет следующий вид:
2
N -2
X
0
5
0
Р NQ l,1 + NQ) exp [ 1+ NQ
хГ,
U
, N -1
1+
NQ
(20)
Для этих условий вероятность правильного обнаружения [7] можно определить как
С л
D = 1 - Г, (Sr, N -1) + ^ + J
хГ,
St
, N -1
1+
NQ
exp
St
1 + NQ
(21)
Значение порога обнаружения для заданной ^ можно получить, используя выражения (9)-(12).
На рисунке 7 представлены кривые обнаружения D=f(Q) для рассматриваемого случая при фиксированном значении ^=10-12, кривая 1 соответствует N=1, кривая 2 - N=10, кривая 3 - N=50, а кривая 4 - N=100.
D
0.8
0.6
0.4
0.2
Q^
-10 -5 0
5 10 15 20 25 30
Рис. 7. Кривые обнаружения для сигнала, соответствующего первой модели Сверлинга при F=10'9
Fig. 7. Detection curves for a signal corresponding to the first Swerling model when F=10'9 is fixed
Разработанная методика позволяет рассчитать показатели качества обнаружения РЛС для широкого класса моделей принимаемых сигналов, в том числе малозаметных и малоразмерных целей, которые могут быть использованы для определения показателей качества систем более высокого ие-
рархического уровня, например, показателей вторичной обработки сигналов, вероятности обнаружения (захвата) сигнала в зоне поиска.
Литература
1. Куприянов А.И., Сахаров А.В. Теоретические основы радиоэлектронной борьбы. М.: Вузовская книга, 2007.
2. Современная радиоэлектронная борьба. Вопросы методологии; [под ред. В.Г. Радзиевского]. М.: Радиотехника, 2006.
3. Справочник по радиолокации; [под ред. М. Сколника]. М.: Сов. радио, 1977. Т. 1.
4. Кузьмин Г.В., Иванов А.Н., Рюмшин А.Р., Ягольни-ков С.В. Методика оценки эффективности РЛС обнаружения и сопровождения траекторий целей в условиях помех // Радиотехника (Журнал в журнале). 1996. № 11.
5. Кузьмин Г.В., Рюмшин А.Р., Шевчук В.И., Ягольни-ков С.В. Рабочие характеристики многоканального некогерентного весового обнаружителя сигналов различной мощности в условиях стационарных шумовых помех // Радиотехника (Журнал в журнале). 1997. № 11.
6. Кузьмин Г.В., Рюмшин А.Р., Шевчук В.И., Ягольни-ков С.В. Оценка эффективности и обоснование помех многоканальному некогерентному весовому обнаружителю сигналов различной мощности // Радиотехника (Журнал в журнале). 1998. № 2.
7. Mеikle Н. Modern Radar Systems. 2nd ed. Artech House, 2008.
8. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1984.
9. Справочник по специальным функциям; [под ред. М. Абрамовица и И. Стиган]. М.: Наука, 1979.
References
1. Kupriyanov A.I., Sakharov A.V. Teoreticheskie osnovy radioelektronnoy borby [Theoretical basics of radioelectronic struggle]. Moscow, Vuzovskaya kniga, 2007, 356 p.
2. Radzievskiy V.G., ed. Sovremennaya radioelektronnaya borba. Voprosy metodologii [Modern radioelectronic struggle. Methodology issues]. Moscow, Radiotexnika Publ., 2006, 424 p.
3. Skolnik M. Spravochnik po radiolokatsii [Radiolocation guide]. Moscow, Sоv. radw Publ., vol. 1, 1977, 456 p.
4. fauzmin G.V., Ivanov A.N., Ryumshin A.R., Yagolni-kov S.V. A technique for efficiency assessment of radar detection and moving target tracking under noise conditions. Radiotekhnika [Radioengineering]. 1996, no. 11.
5. fauzmin G.V., Ryumshin A.R., Shevchuk V.I., Yagolni-kov S.V. Operating characteristics of multy-channel incoherent weighted detector for signals with different power under permanent noise conditions. Radiotekhnika [Radioengineering]. 1997, no. 11.
6. fauzmin G.V., Ryumshin A.R., Shevchuk V.I., Yagolni-kov S.V. An efficiency assessment and noise validation for multy-channel incoherent weighted detector for signals with different power. Radiotekhnika [Radioengineering]. 1998, no. 2.
7. Mеikle Н. Modern radar systems. 2nd edition, Artech House Publ., 2008, 563 p.
8. Кот G.A., Кот T.M. Spravochnik po matеmаtikе dlya naushnykh rabotnikov i inzhеnеrоv [Guide on mathematics for scientists and engineers]. Russ. ed.: Moscow, Nauka Publ., 1984, 832 p.
9. Abramovits M., Stigan I. eds. Spravochnik po spetsialnym funktsiyam [Guide on special functions]. Moscow, Nauka Publ., 1979, 832 p.
X
1
X
1
Вниманию авторов!
Пристатейный список литературы (не менее 5 пунктов) необходимо транслитерировать в соответствии со стандартом BGN/PCGN. Самоссылки в списке не должны составлять более 20 %. Новые редакционные требования на сайте журнала www.swsys.ru.
