CC BY
159
УДК 519.8
ВЕРОЯТНОСТНАЯ МОДЕЛЬ БОЕВЫХ ДЕЙСТВИЙ ПРИ УПРЕЖДАЮЩЕМ УДАРЕ ОДНОЙ ИЗ СТОРОН
В.Ю. Чуев, И.В. Дубограй
МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]; [email protected]
Разработаны вероятностные модели двухсторонних боевых действий для различных начальных численностей противоборствующих группировок на основе теории непрерывных марковских процессов. Исследовано влияние упреждающего удара одной из противоборствующих сторон на исход боя и его основные показатели. Получены расчетные формулы для вычисления основных показателей боя. Приведены результаты расчетов, показавшие, что это влияние существенно в бою близких по силам группировок. Нанесение одной из противоборствующих сторон упреждающего удара может на 30 % уменьшить ее потери и на 30% увеличить потери противника. Показано, что при большом (в 3 раза и более) начальном превосходстве одной из сторон это влияние незначительно. Отмечено увеличение влияния упреждающего удара на исход боя и на его основные показатели при возрастании начальных численностей противоборствующих группировок.
Ключевые слова: вероятностная модель двухсторонних боевых действий, эффективная скорострельность, непрерывный марковский процесс, параметр соотношения сил.
STOCHASTIC MODEL OF BOTH-SIDED BATTLE ACTIONS DURING THE PRE-EMPTIVE ATTACK BY ONE OF WARRING PARTIES
V.Yu. Chuev, I.V. Dubogray
Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]; [email protected]
The stochastic models of both-sided battle actions for various initial numerosity of opposing groupings have developed on the basis of continuous markov processes. The influence ofpre-emptive attack of one of the opposing force on the battle outcome and its main indicators have investigated. Calculation formulas for counting main indicators of battle were obtained. The results of calculations showed that this influence is essentially in battle of groupings similar in forces are given. Preventive strike by one of the opposing forces can reduce its losses by 30 % and increase the enemy's casualties up to 30 %. It has been shown that this effect is negligible at large (3-fold or more) initial superiority by one of the parties. An increase in the influence of the pre-emptive attack on the battle outcome and its main indicators is noted with an increase the initial numerosity of opposing groupings.
Keywords: stochastic model of both-sided battle actions, effective rapidity of fire, continuous markov process, balance of forces parameter.
Введение. Важнейшие характеристики систем вооружения и военной техники — показатели их боевой активности, позволяющие определить степень их приспособленности к решению конкретных боевых задач. При проведении такой оценки необходимо учитывать противодействие противника. Наиболее полно это можно выполнить с помо-
щью моделей двухсторонних боевых действий [1-3]. В качестве основы оценки целесообразно использовать вероятностные модели боя, так как с учетом его стохастического характера указанные модели позволяют отобразить бой с большей степенью точности и полноты, чем модели динамики средних, а также достовернее исследовать его основные показатели.
Один из возможных методов моделирования боя — применение теории марковских процессов [4, 5]. Процесс, протекающий в некоторой системе, называется марковским, если для каждого момента времени вероятность состояния системы в будущем зависит только от состояния системы в настоящий момент и не зависит от того, каким образом система пришла в это состояние [6].
Описание протекания боя. Рассмотрим бой двух группировок. В начале боя первая группировка X имеет т однотипных боевых единиц, а вторая (противоборствующая) группировка У — п однотипных боевых единиц, не обязательно однородных с боевыми единицами группировки X.
Бой происходит следующим образом. Противоборствующие стороны открывают огонь по противнику одновременно и ведут бой до тех пор, пока одна из сторон не будет полностью уничтожена. При этом стороны имеют полную и незапаздывающую информацию о состоянии боевых единиц противника (поражены они или нет) и ведут огонь только по непораженным боевым единицам [7]. Последовательность выстрелов, осуществляемых каждой участвующей в бою единицей, представим в виде пуассоновского потока событий [8]. При моделировании боя перейдем к потоку успешных выстрелов, который также будем полагать пуассоновским. Выстрел назовем "успешным", если он поражает боевую единицу противника [9, 10].
Совокупность противоборствующих сторон назовем системой. Обозначим ее состояния (г; ]), что означает, что в данный момент времени сохранились г единиц группировки X и ] единиц группировки У. Пусть ^ (¿) вероятность того, что в момент времени £ система находится в состоянии (г; ]). Отметим, что состояние (0; 0) не является состоянием рассматриваемой системы, так как вероятность одновременного поражения двух и более единиц является бесконечно малой величиной. Примем, что моменту начала боя соответствует ¿0 = 0.
В этом случае процесс боевых действий описывается системой дифференциальных уравнений
F0 (t) = № (t), i = 1; m;
F0j (t) = juFij (t), j = 1; n;
Fj (t) = - (iv + ju) Fj (t) +
+ivFi,j+i (t) + juFi+ij (t) , i = 1; m - 1; j = 1; n - 1; (1)
FLj (t) = - (mv + ju) Fmj (t) + mvFm,j+i (t), j = 1; n - 1;
Fn (t) = - (iv + nu) (t) + nuFj+i,„ (t), i = 1; m - 1;
FL (t) = - (mv + nu) Fmn (t) с начальными условиями
to = 0; Fmn (to) = 1; Fj (t0) = 0 при i + j<m + n, (2)
где v = pxAx, u = pyAy — эффективные скорострельности боевых единиц группировок X и Y; ,py — вероятности поражения боевой единицы противника одним выстрелом единицы группировок X и Y; Ax, Ay — практические скорострельности боевых единиц группировок X и Y.
Отметим, что устойчивыми состояниями рассматриваемой системы являются только состояния (1; 0), (2; 0) ,..., (m; 0), (0; 1) , (0; 2) ,
...,(0; n).
Исследование влияния упреждающего удара на основные показатели боя. Теперь рассмотрим ситуацию, когда определенные обстоятельства (например, хорошая маскировка) боевых единиц одной из сторон (например, стороны X) позволяет ей в течение некоторого, достаточно небольшого времени tc вести огонь по боевым единицам другой стороны, не испытывая противодействия противника. Тогда на промежутке времени t £ [t0, tc] протекание боя можно описать системой уравнений
F^o (t) = mvFmi (t);
F'm3 (t) = mv (Fm,j+i (t) - Fmj (t)), j = 1; n - 1; FLn(t) = mvFm n (t);
Fj (t) = 0, i = 0; m - 1, j = 0; n с начальными условиями (2).
Если принять величины v и u в течение всего боя постоянными v = pxAx = const; u = py Ay = const), то получим следующие вероят-
ности состояний системы в момент времени ¿с:
п
Рт0 с) = 1 — ^^ Рт] (£с) = Со;
3=1
Рш (¿с) = ^Г = Сз, , = 1; п — 1; (4)
(п —;)!
р (/) = = С ;
тп \ьс/ ^ °п 1
Р? (¿с) = 0, г = 0; т — 1, з = 0; п.
Дальнейшее протекание боя описывается системой уравнений (1) с начальными условиями
Рт3 Ос) С3, ^ 0; п; (5)
рз (¿с) =0, г = 0; т — 1, = 0; п. (5)
Протекание боя при упреждающем ударе группировки У можно описать аналогично уравнениям и условиям (3)-(5).
При постоянных эффективных скорострельностях и и V боевых единиц сторон авторами настоящей работы получены формулы для вычисления вероятностей текущих (¿)) и окончательных (то)) состояний для боев "1 : п" (одна боевая единица группировки X против произвольного числа однотипных боевых единиц группировки У), "п:1" и "3:3".
Для боя "3:3" при упреждающем ударе группировки X вероятности окончательных состояний вычисляются как
с^и2
10 = (V + и) ^ + и) ^ + и) +
/ 9 2 3 2 6 \
+с2 V2 (V + и) — 2у + м + — V + 2и — 3^Т2Й) +
3 /_5___2___2 3
+ Сз^ + и 3^ + и) 3 (V + 2и) + 4^ + и) +
3 6 6
+
F20 (то) =
4 (v + 3u) 3v + 2u 2v + 3u 2civu
(2v + u) (3v + u)
2 2 3 6
+2c2v---1-----1--+
v + u 2v + u 3v + u 3v + 2u,
6C3V2 / 17 2 3 6
+---г + ^-г - ^-г + ^-Г- +
u V 6 (v + u) 3 (2v + u) 2 (3v + u) 3v + 2u 2v + 3u/ '
2
„ / ч 3civ 9c2v2 / 1 1 \
F30 = Co + -— + - -— - +
3v + u u \3v + u 3v + 2uy
27c3v3 / 1 2 1
2u2 \3v + u 3v + 2u 3v + 3u
F01 (то) = ^^10 (то);
^ 4c2u3 / 1 1 1 .
Fo2 (то) = ^ -vr^ + vr^u + + (6)
v2
6c3u2 / 17 2 3 2 6
+ v V 6 (v + u) + 3(v + 2u) -2 (v + 3u)+3v + 2u+ + 3u ^ 27сзм3 / 1 2 1
03 = 2v2 Vv + 3u - 2v + 3u + 3v + 3u^ ' (7)
(то) =0, г = 1, 2, 3; j = 1, 2, 3,
ci = (3v2¿c) e-3víc; c2 = 3vtce-3vtc; c3 = e-3vtc; c0 = 1 - ci - c2 - c3.
Вероятности окончательных состояний для боя "3:3" при упреждающем ударе группировки Y находятся по формулам, аналогичным формулам (6).
Для боя "1 : n" получаем
1) при упреждающем ударе стороны X:
n
V-^ ciuvt
Fok M = % (ku + v)... (iu + v), k = 1; n - 1;
nCnU
Fon (то) =
-"n
;
Co = 1 - Cfe; Cfe = ( k)! e vtc, k = 1; n - 1; Cn = e
nu + v (8)
n i
, , ^—л cív v , ,
Fio (то) = Co + > 7—■—г-/-■—г = Co + -F01 (то);
(u + v)... (nu + v) u
Í=1 4 x 4 x
Fik (то) =0, k = 1; n,
n / . \n—k
(vtc) p-VÍc l, 1. „ _ „-Vtc.
k=i (n - k)!
2) при упреждающем ударе группировки Y :
^ kdiuvn-k
Fok (то) = тг—:—ч-7-:—г, k = 1; n - 1;
(ku + v)... (nu + v)
Fon = do + n^f-; (9)
nu + v
Fio (то)= div"
(u + v)... (nu + v)' где do = 1 - e-nutc; di = e-nuíc.
Формулы для вычисления вероятностей окончательных состояний для боя "п : 1" при упреждающем ударе одной из сторон аналогичны формулам (7) и (8).
Перейдем к исследованию основных показателей боя. В первую очередь, к ним относятся вероятности победы Р0х и Р0у сторон X и У, а также математические ожидания относительного числа Мх и Му сохранившихся боевых единиц сторон к концу боя. В общем случае перечисленные показатели вычисляются по формулам
Pox = Е Fi0 (то) ;
i=1 n
poy = E Foj(то);
£
j=1 m (10)
Mx = - V iFi0 (то); m
i= 1
1
m = -E j'F0j (то)
, 03 (ю)
п
3 = 1
Отметим, что, если группировка X в начале боя имеет одну боевую единицу т = 1, то Мх = Р0х. Аналогично, Му = Р0у при п =1.
п и ™
Введем параметр соотношения сил к = —*/ — [9] и приведен-
т V V
ное время Ь = л/июЬ [11]. В этих обозначениях формулы (4) примут следующий вид:
п
Рт0 (£с) = 1 _ Рт-з (£с) = с0;
3=1
- 1 (ПЪС \П-3 _Мс
Ртз (£С) = (п _ .)Л е К = с-, 3 = 1; п - 1; (11)
_ _п1с
Ртп (£с) е к сп;
(¿С) = 0, г = 0; т _ 1, 3 = 0; п, где ¿С = ^йюЬс.
Расчеты основных показателей боя при упреждающем ударе одной из сторон. На основе формул (6)-(10) проведены расчеты основных показателей боя при различных значениях параметра соотношения сил к и приведенного времени ¿с, в течение которого группировка X не испытывает противодействия противника. Значения математических ожиданий относительного числа сохранившихся боевых единиц Мх и Му сторон к концу боя "1 : 1" (дуэльного боя)
даны в табл. 1. Для этого боя вероятности победы сторон X и У равны: Р0х = Мх и Р0у = Му. Для боя "1 : 5" (табл.2) основные показатели боя составляют Р0х = Мх, Р0у = 1 — Р0х, для боя "5 : 1" (см. табл. 1) — Р0у = Му, Р0х = 1 — Р0у. Кроме значений Мх и Му для боев "2 : 2" и "3 : 3" (табл. 3) указаны вероятности победы группировки X — Р0х. Отметим, что в этих случаях вероятность победы группировки У равна Р0у = 1 — Р0х.
Таблица 1
Значения математического ожидания относительного числа сохранившихся боевых единиц сторон к концу боя «1:1» и «5:1»
к tc tc
0 0,1 0,2 0,3 0,4 0,5 0 0,1 0,2 0,3 0,4 0,5
Бой "1 : 1" Бой "5 : 1"
0,4 0,862 0,138 0,892 0,108 0,916 0,084 0,935 0,065 0,949 0,051 0,960 0,040 0,811 0,068 0,853 0,053 0,885 0,041 0,911 0,032 0,930 0,025 0,946 0,019
0,8 0,610 0,390 0,656 0,344 0,696 0,304 0,732 0,268 0,763 0,237 0,791 0,209 0,446 0,429 0,511 0,379 0,568 0,334 0,619 0,295 0,664 0,260 0,703 0,230
1,0 0,500 0,500 0,548 0,452 0,591 0,409 0,630 0,370 0,665 0,335 0,697 0,303 0,328 0,571 0,392 0,517 0,450 0,468 0,502 0,423 0,550 0,383 0,593 0,346
1,2 0,410 0,590 0,457 0,543 0,500 0,500 0,541 0,459 0,578 0,422 0,611 0,389 0,248 0,672 0,308 0,619 0,363 0,569 0,414 0,524 0,461 0,482 0,504 0,443
1,4 0,338 0,662 0,383 0,617 0,426 0,574 0,465 0,535 0,502 0,498 0,537 0,463 0,192 0,744 0,248 0,693 0,299 0,645 0,348 0,601 0,393 0,559 0,435 0,521
3,0 0,100 0,900 0,130 0,870 0,158 0,842 0,186 0,814 0,212 0,788 0,238 0,762 0,047 0,936 0,078 0,905 0,109 0,876 0,138 0,847 0,166 0,819 0,193 0,792
Таблица 2
Значения математического ожидания относительного числа сохранившихся боевых единиц сторон к концу боя «1:5»
к t c
0 0,1 0,2 0,3 0,4 0,5
0,4 0,909 0,067 0,942 0,037 0,965 0,019 0,981 0,009 0,990 0,004 0,995 0,002
0,8 0,693 0,232 0,746 0,176 0,795 0,131 0,838 0,096 0,876 0,069 0,906 0,049
1,0 0,571 0,328 0,627 0,266 0,681 0,212 0,732 0,168 0,778 0,131 0,819 0,101
1,2 0,456 0,423 0,511 0,358 0,566 0,299 0,619 0,248 0,669 0,204 0,717 0,166
1,4 0,355 0,511 0,405 0,445 0,457 0,385 0,509 0,330 0,560 0,281 0,610 0,238
3,0 0,030 0,871 0,040 0,832 0,052 0,794 0,066 0,755 0,082 0,718 0,101 0,680
Левый столбец в таблицах соответствует одновременному открытию огня противоборствующими сторонами, остальные столбцы — ситуации, когда единицы группировки X проводят по одному-два выстрела до открытия группировкой У ответного огня.
Таблица 3
Значения математического ожидания относительного числа сохранившихся боевых единиц сторон к концу боя и вероятности победы стороны X для боя
«2:2» и «3:3»
К tc tc
0 0,1 0,2 0,3 0,4 0,5 0 0,1 0,2 0,3 0,4 0,5
Бой "2 : 2" Бой "3 : 3 "
0,4 0,871 0,045 0,943 0,909 0,029 0,963 0,937 0,018 0,975 0,957 0,012 0,984 0,971 0,007 0,990 0,981 0,005 0,993 0,881 0,018 0,974 0,921 0,010 0,985 0,950 0,005 0,992 0,969 0,003 0,996 0,981 0,001 0,998 0,989 0,001 0,999
0,8 0,559 0,281 0,656 0,624 0,228 0,714 0,682 0,185 0,763 0,732 0,150 0,804 0,776 0,121 0,839 0,813 0,098 0,868 0,542 0,228 0,688 0,619 0,172 0,755 0,687 0,130 0,810 0,747 0,097 0,854 0,797 0,072 0,888 0,839 0,054 0,915
1,0 0,417 0,417 0,500 0,481 0,355 0,563 0,542 0,302 0,620 0,598 0,256 0,671 0,648 0,217 0,715 0,693 0,184 0,755 0,378 0,378 0,500 0,454 0,307 0,578 0,526 0,248 0,647 0,593 0,200 0,708 0,654 0,160 0,759 0,708 0,128 0,803
1,2 0,305 0,534 0,371 0,364 0,469 0,432 0,422 0,411 0,490 0,477 0,361 0,544 0,529 0,316 0,592 0,577 0,277 0,637 0,253 0,512 0,344 0,320 0,436 0,419 0,387 0,370 0,490 0,452 0,312 0,556 0,515 0,263 0,616 0,573 0,220 0,669
1,4 0,222 0,627 0,273 0,274 0,564 0,329 0,327 0,507 0,383 0,378 0,454 0,435 0,428 0,407 0,483 0,475 0,364 0,529 0,166 0,620 0,230 0,221 0,546 0,295 0,279 0,478 0,360 0,337 0,417 0,423 0.395 0,363 0,483 0,452 0,315 0,539
3,0 0,025 0,911 0,031 0,039 0,875 0,048 0,056 0,840 0,066 0,074 0,806 0,086 0,094 0,773 0,107 0,115 0,741 0,130 0,008 0,920 0,011 0,014 0,883 0,020 0,024 0,846 0,032 0,036 0,809 0,047 0,050 0,773 0,064 0,066 0,738 0,083
Анализы результатов и выводы. Результаты расчетов позволяют сделать следующие выводы.
1. Нанесение упреждающего удара одной из противоборствующих сторон существенно влияет на исход боя достаточно близких по силам группировок. В ряде случаев это приводит к уменьшению ожидаемых потерь стороны, наносящей упреждающий удар, и увеличению ожидаемых потерь стороны, наносящей ответный удар (относительно начальных численностей группировок) более чем на 30 %, а также к изменению более чем на 30 % вероятностей побед противоборствующих сторон.
2. С возрастанием начальных численностей противоборствующих группировок увеличивается влияние упреждающего удара одной из сторон на исход боя и его основные показатели, причем увеличение этого влияния достаточно заметно.
3. При существенном начальном превосходстве стороны, наносящей упреждающий удар, т.е. при превосходстве более чем в 3 раза (к < 1/3) для боя "n : 1" и более чем в 2,5 раза (к < 0,4) для всех остальных боев, влияние упреждающего удара на исход боя и его основные показатели незначительно.
4. При значительном начальном превосходстве стороны, наносящей ответный удар, т.е. при превосходстве более чем в 4 раза (к > 4) для боя "n : 1" и более чем в 3 раза (к > 3) для всех остальных боев, влияние упреждающего удара на ожидаемые потери стороны, наносящей упреждающий удар, а также на вероятности побед сторон незначительно. При этом упреждающий удар может оказать достаточно весомое влияние на потери стороны, наносящей ответный удар: более 15% для боя "n : 1" и более 10% относительно начальной численности группировки для всех остальных боев.
ЛИТЕРАТУРА
1. Ткаченко ДН.Математические модели боевых действий. М.: Сов. радио, 1969. 240 с.
2. Ильин В.А. Моделирование боевых действий сил флота // Программные продукты и системы. 2006. № 1. С. 23-27.
3. Winston W.L. Operations Research: applications and algorithms. Duxbury Press, 1998. P. 128.
4. Алексеев О.Г., Анисимов В.Г., Анисимов Е.Г. Марковские модели боя. М.: Министерство обороны СССР, 1985. 85 с.
5. Чуев В.Ю. Вероятностная модель боя многочисленных группировок // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. Спец. выпуск "Математическое моделирование". 2011. С. 223-232.
6. Вентцель Е.С. Теория вероятностей. М.: Высш. шк., 1999. 576 с.
7. Jaswal N.K. Military Operations Research. Quantitative Desigion Making. Kluwer Academic Publishers, 1997. P. 388.
8. Вентцель Е.С. Исследование операций. М.: УРСС, 2006. 432 с.
9. ЧуевЮ.В. Исследование операций в военном деле. М.: Воениздат, 1970. 270 с.
10. Shanahan L. Dynamics of model battles. N.Y.: Phisics Department, State University of New York. 2003. P. 1-43.
11. Дубограй И.В., Дьякова Л.Н., Чуев В.Ю. Учет упреждающего удара при моделировании двухсторонних боевых действий // Инженерный журнал: наука и инновации. 2013. Вып. 7 [Электронный ресурс] URL: http://enggournal.ru/catalog/mathmodel/hidden/842.html (дата обращения: 21.12.2013).
REFERENCES
[1] Tkachenko P.N. Matematicheskie modeli boevykh deystviy [Mathematical models of battle actions]. Moscow, Sov. Radio Publ., 1969. 240 p.
[2] Il'in V.A. Modeling of naval combat operations. Prog. prod. i sist. [Software and Systems], 2006, no. 1, pp. 23-27 (in Russ.).
[3] Winston W.L. Operations Research: applications and algorithms. Duxbury Press, 1998, p. 128.
[4] Alekseev O.G., Anisimov V.G., Anisimov E.G. Markovskie modeli boya [Markovian battle model]. Moscow, Ministerstvo oborony SSSR Publ., 1985. 85 p.
[5] Chuev V.Yu. Probability model of a double-sided battle of numerous groups. Vestn. Mosk. Gos. Tekh. Univ. im.N.E. Baumana, Estestv. Nauki., Spetsvyp. "Matematicheskoe modelirovanie" [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci., Spec. Iss. "Mathematic simulation"], 2011, pp. 223-232 (in Russ.).
[6] Venttsel' E.S. Teoriya veroyatnostey [Probability theory]. Moscow, Vysshaya Shkola Publ., 1999. 576 p.
[7] Jaswal N.K. Military Operations Research. Quantitative Desigion Making. Kluwer Academic Publishers, 1997, p. 388.
[8] Venttsel' E.S. Issledovanie operatsiy [Operations research]. Moscow, URSS Publ., 2006. 432 p.
[9] Chuev Yu.V. Issledovanie operatsiy v voennom dele [Operations research in military affairs]. Moscow, Voenizdat Publ., 1970. 270 p.
[10] Shanahan L. Dynamics of model battles. New York, Phisics Department, State University of New York, 2003, pp. 1-43.
[11] Dubogray I.V., Dyakova L.N., Chuev V.Yu. Preemptive attack consideration when duel combat operations simulating. Jelektr. nauchno-tehn. izd. "Inzhenernyj zhurnal: nauka i innovacii", MGTU im. N.E. Baumana [El. Sc.-Tech. Publ. "Eng. J.: Science and Innovation" of Bauman MSTU], 2013, no. 7 (19). Available at: http://engjournal.ru/eng/catalog/mathmodel/hidden/842.html (accessed 21.12.2013) (in Russ.).
Статья поступила в редакцию 03.03.2014
Чуев Василий Юрьевич — канд. техн. наук, доцент кафедры "Вычислительная математика и математическая физика" МГТУ им. Н.Э. Баумана. Автор более 20 научных работ в области прикладной математики.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Chuev V. Yu. — Cand. Sci. (Eng.), assoc. professor of "Computational Mathematics and Mathematical Physics" department of the Bauman Moscow State Technical University. Author of more than 20 publications in the field of applied mathematics. Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
Дубограй Ирина Валерьевна — доцент кафедры "Вычислительная математика и математическая физика" МГТУ им. Н.Э. Баумана. Автор нескольких научных работ в области прикладной математики.
МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.
Dubogray I.V. — assoc. professor of "Computational Mathematics and Mathematical Physics" department of the Bauman Moscow State Technical University. Author of some publications in the field of applied mathematics.
Bauman Moscow State Technical University, 2-ya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.
