Vol. 19, № 02, 2016
Civil Avition High TECHNOLOGIES
УДК 629.7.07
КОМБИНИРОВАННАЯ ИМИТАЦИОННО-АНАЛИТИЧЕСКАЯ МОДЕЛЬ ПРОПУСКНОЙ СПОСОБНОБНОСТИ АЭРОДРОМНОЙ ВОЗДУШНОЙ зоны
В.Л. КУЗНЕЦОВ, А.И. ГОРЧАКОВ
В работе развивается подход к задаче оценки пропускной способности воздушной аэродромной зоны, основанный на идее «квантования» воздушного пространства, перехода в дискретное время при наложении ограничений на заполнение квантовых состояний воздушными судами [4-6]. Показано, что такой формализм позволяет существенно снизить трудоемкость задачи имитационного моделирования движения воздушных судов в аэродромной зоне, учесть ограничения, накладываемые на ее пропускную способность пересечениями маршрутов прилета и вылета ВС. В статье приводятся результаты проведенного имитационного эксперимента, иллюстрирующие зависимость пропускной способности зоны от количества пересечений маршрутов воздушных судов.
Ключевые слова: пропускная способность аэродрома, квантование воздушного пространства, имитационный эксперимент.
ВВЕДЕНИЕ
В настоящее время вопросам моделирования систем организации воздушного движения (ОрВД) и их компонент - процессов планирования и управления воздушным движением (УВД), уделяется большое внимание. Одним из важнейших параметров, ограничивающих возможности УВД, является пропускная способность аэродрома. Следует отметить, что утвержденные в настоящий момент методики [1] и международные рекомендации [2] по этому вопросу носят сугубо инженерный характер и базируются на излишне упрощенных моделях, не позволяющих глубоко анализировать вопросы оптимизации пропускной способности аэродрома. В частности, пропускную способность аэродрома отождествляют с пропускной способностью взлетно-посадочных полос, не принимая во внимание ограничения, возникающие в воздушном пространстве аэродромной зоны, связанные с пересечением траекторий движения взлетающих и заходящих на посадку воздушных судов (ВС).
В силу сложности аналитического решения задач УВД в настоящее время в научных исследованиях используют возможности и методы имитационного моделирования. Сейчас в мире выделяют уже более десятка различных комплексов имитационного моделирования, различающихся отдельными, порой очень существенными аспектами решаемых задач [3]. Как правило, аналитика в таких комплексах используется для описания динамики ВС с целью максимального приближения имитируемой картины воздушной обстановки к реально складывающейся. По существу, применение такой аналитики приводит как к росту адекватности модели (желаемый результат), так и к резкому увеличению потребных ресурсов вычислительной техники и, соответственно, времени, затрачиваемою на расчет (плата за результат). Компромисс между требованиями адекватности и простоты модели решается в пользу адекватности. Однако в конкретной задаче УВД о пропускной способности аэродромной воздушной зоны, предполагающей максимально плотную «упаковку» ВС в воздушном пространстве, этот компромисс можно сместить в сторону требования простоты модели.
В представляемой работе развивается подход к задаче о пропускной способности аэродрома, основанный на предложенном ранее способе формализации - абстрактной структуризации воздушного пространства [4-6], в результате которой аналитический этап моделирования позволяет существенно снизить уровень сложности модели, подвергаемой имитационному эксперименту.
Civil Avition High TECHNOLOGIES
Vol. 19, № 02, 2016
1. ПОСТАНОВКА ЗАДАЧИ И ОСНОВНЫЕ ПРИБЛИЖЕНИЯ МОДЕЛИ
Рассмотрим аэродром, работающий по пропускной способности на пределе своих возможностей. Будем полагать, что структура воздушного пространства аэродромной зоны задана. Схематично она представлена на рис. 1. Это значит, что известны все геометрические параметры маршрутов прилета и вылета ВС - их SID (Standard Instrument Departure) и STAR (Standard Terminal Arrival Route), а также точки на карте и высотные эшелоны их пересечений.
Заданными считаются также статистические данные о частоте прибытия и отлета ВС по воздушным трассам, подходящим к рассматриваемой аэродромной зоне, а следовательно, и все соответствующие вероятности, описывающие входные и выходные потоки ВС.
Целью имитационного эксперимента является оценка максимальных интенсивностей встречных потоков ВС, бесконфликтно проходящих воздушное пространство аэродромной зоны. Эти интенсивности и определяют пропускную способность аэродрома.
В качестве основных приближений модели будем полагать выполнение следующих условий.
1. Все ВС технически исправны, а экипажи строго выдерживают предписанный план полета, двигаясь по маршрутам друг за другом (векторение не допускается).
2. Все ВС принадлежат одному классу, поэтому продольные интервалы эшелонирования в каждой данной точке пространства одинаковы для всех ВС.
3. Все ВС обладают равным приоритетом, обслуживаются в порядке прибытия и, в случае необходимости, задерживаются в зоне ожидания на границе района аэродрома.
Рис. 1. Схематичное представление структуры воздушного пространства в зоне аэродрома
2. МЕТОД КВАНТОВАНИЯ АЭРОДРОМНОГО ВОЗДУШНОГО ПРОСТРАНСТВА
2.1. Основные положения метода
Основная идея развиваемого метода заключается в том, что при расчете пропускной способности воздушной зоны аэродрома плотность потоков взлетающих и заходящих на посадку ВС должна быть максимально высокой при условии соблюдения норм эшелонирования. При этом необходимо учитывать, что ВС, находящимся на определенном эшелоне и удаленности от ВПП, предписывается некоторая рекомендуемая скорость движения, а на SID/STAR должно ре-ализовываться стационарное движение.
Мысленно очертим вокруг всех ВС, движущихся в воздушном аэродромном пространстве, области в виде прямоугольных параллелепипедов, ориентированных по направлению скорости ВС. Договоримся сразу, что линейные размеры этих областей должны быть не меньше
Vol. 19, № 02, 2016
Civil Avition High TECHNOLOGIES
установленных норм эшелонирования. Конфигурация мгновенного расположения ВС в воздушном пространстве должна быть такова, что построенные параллелепипеды не накладываются друг на друга. Это ограничение связано с тем, что факт наложения параллелепипедов однозначно приводит к нарушению норм эшелонирования, ВС оказываются ближе друг к другу, чем это предписано Федеральными авиационными правилами [7].
Поскольку все ВС, согласно допущению модели, двигаются строго вдоль своих SID/STAR, то естественно сделать вывод о том, что и выделенные прямоугольные параллелепипеды двигаются вдоль этих линий. В предельном случае они двигаются непосредственно друг за другом.
Однако такая динамическая картина представления потоков ВС может быть существенно упрощена, если перейти в дискретное время с периодом квантования, равным времени, за которое некоторый выделенный параллелепипед займет положение предыдущего. Это время т определяется из условия:
т = max
{Твзл посад '{/V }} '
Здесь тезл ,тпосад - времена занятости ВПП при взлете и посадке, Lj - значения минимумов продольного эшелонирования на i-м участке траектории (1<i<P, P - число квантовых состояний от ВПП до точки выхода на трассу), а Vi - соответствующая плановая скорость ВС. Теперь линейный размер выделенного объема (параллелепипеда) в направлении SID/STAR на i-м участке определяется по формуле L j = Vj ■ т > Lj.
Переход в дискретное время позволяет «остановить» движение параллелепипедов, сводя изменение конфигурации ВС в воздушном пространстве к скачкообразным переходам из одного «кванта пространства» в другой, следующий по ходу движения.
Следующий шаг в унификации модели связан с введением т.н. путевых координат - расстояний, измеренных вдоль соответствующих SID/STAR от ВПП до точек нахождения ВС.
Если теперь высоту положения ВС связать с занимаемым им вертикальным эшелоном, то топологически сложную картину распределения ВС в пространстве аэродрома можно заменить простой прямоугольной матрицей (техн.), где в горизонтальной строке последовательно зафиксированы положения квантованных элементов воздушного пространства, находящихся на разном удалении от ВПП, а столбцы описывают вертикальное эшелонирование. Добавив теперь некоторые ограничения на эволюцию заполнения матрицы взлетающими и заходящими на посадку ВС, получаем простую и наглядную картину воздушной обстановки в зоне аэродрома (рис. 2). Каждая «строка» соответствует своему эшелону по высоте, а «столбец» - удалению от ВПП. Кружками помечены заполненные в данный момент времени квантовые состояния, т.е. состояния, в которых находятся ВС. Другие состояния свободны. Стрелками указаны переходы ВС в новые состояния на следующем временном такте.
По мере приближения к ВПП заходящие на посадку ВС будут занимать все более низкие эшелоны по высоте. При этом в каждом квантовом состоянии по-прежнему может находиться
0. >■ ; ■ Б О
•—
•
* г *—
•
—• _
—•
•
Заходящие на посадку ВС
Рис. 2. Иллюстрация картины воздушной обстановки в зоне аэропорта в формализме квантования воздушного пространства
Civil Avition High TECHNOLOGIES
Vol. 19, № 02, 2016
не более одного ВС. Принимая во внимание синхронность перемещений всех ВС в квантованном воздушном пространстве, приходим к заключению о том, что в каждом «столбце» нашего пространства может находиться не более одного заходящего на посадку (или, аналогично, взлетающего) ВС. Такой столбец ассоциируется с вектором у+ = (0,0,1,0,...,0)Г в K-мерном пространстве (K - полное число эшелонов по высоте в воздушной зоне аэродрома), компоненты которого в основном нули (соответствуют отсутствию ВС на данном эшелоне), и лишь одна компонента вектора у+ может равняться единице (состояние занято). Вектор у описывает заходящие на посадку, а вектор у+ - взлетающие ВС.
При рассмотрении всего квантованного пространства удобно перейти к описанию с помощью векторных полей Ц (t',р'). Здесь t' - текущий момент времени, t' <= Z, р - порядковый номер столбца (удаленность от ВИН, 1<р<P ). Индекс n фиксирует рассматриваемую траекторию SID или STAR. Здесь и далее мы перешли к безразмерному времени - t'=, а расстояние р до ВИЛ измеряем числом квантовых состояний, как это отмечено на рис. 3. Далее, для простоты обозначений, штрихи у t' и р' будем опускать.
Пересечения возможных траекторий полета взлетающих и заходящих на посадку ВС и ограничения на реализации таких полетов, обусловленные наличием минимумов эшелонирования, в предложенном методе реализуются с помощью т.н. стационарного матричного поля критических перекрытий («минного поля» в работах [4-6]). Это поле - R(p,p), выделяет пары столбцов состояний с координатами (р,р2) Для взлетающих и заходящих на посадку ВС соответственно, в которых могут возникать конфликты между ВС. Для таких пар R(p,p) - единичная матрица, для остальных - нулевая.
2.2. Моделирование ветвлении и слиянии маршрутов прилета и вылета ВС
Описанное выше и изображенное на рис. 2 графическое представление квантованного воздушного пространства справедливо для простого случая, когда в зоне аэродрома реализуется только один маршрут вылета (с соответствующим SID) и один маршрут прилета (с соответствующим STAR). Понятно, что такое упрощение имеет право на существование лишь как иллюстративный пример, разъясняющий первые основные шаги метода. Основные проблемные ситуации в воздушном пространстве возникают при пересечении траекторий взлетающих и заходящих на посадку ВС, двигающихся по разным маршрутам. Для получения адекватной в отношении описания эффектов ветвления и слияния маршрутов ВС модели воздушного пространства, при сохранении преимуществ описанной выше структуры, используем следующий прием.
Объединим квантованные элементы пространства, соответствующие как всем разрешенным маршрутам прилета и STAR, так и маршрутам вылета и SID, сформировав новое пространство, подобно тому как это изображено на рис. 3. Здесь вертикальные пунктирные линии - границы слияния (разделения) маршрутов. Левее этих линий объединенные маршруты тождественно совпадают. Темные полосы - иллюстрация для областей поля критических перекрытий. ВС, попавшие в эти полосы на одинаковом высотном эшелоне, конфликтуют. В таком расширенном пространстве векторы для каждого сечения по «дальности» будут представлять собой блочный столбец, блоками которого будут выступать векторы у",i = 1,m, где m - общее количество различных разрешенных маршрутов. Однако, как и в простейшем случае с одним маршрутом прилета и одним маршрутом вылета, модуль вектора у* = (у±,...,цт) может принимать только два значения - нуль или единицу, поскольку ограничения, связанные с пропускной способностью ВИЛ, сохраняются. Дру-
Vol. 19, № 02, 2016
Civil Avition High TECHNOLOGIES
гими словами, иа расстоянии рк = к (в выбранных единицах дальности) на всей совокупности разрешенных маршрутов вылета (прилета) может находиться не более одного ВС.
• V БПП
я *
SB о я
о " ^
н sa
С.- щ W
3 О4*
О и
у.
!
А,».
к ВПП
и ° *
О" °
s ü i
ю til
С) и
ПОСАДКА
Р\
р2 = 1 ръ = 3 рл= 4 ръ = 5 р6 = 6 р7 = 7
l"
т
ц 1
■ - маршрут -3
!-маршрут -2
'-маршрут -1
'■-маршрут -3
"-маршрут -2
>-маршрут -1
Рис. 3. Схематическое изображение расширенной квантованной воздушной зоны аэродрома с тремя воздушными трассами
Отметим еще одну особенность расширенной квантованной модели воздушной зоны. В предложенном подходе на первоначальном этапе теряется информация о точках слияния SID и STAR, соответствующих разным трассам. Так, ВС, удаленные на расстояния р и р2 от ВПП, могут находиться физически на одной линии STAR, хотя на схеме они принадлежат разным маршрутам. На эту особенность можно было бы и не обращать внимания, если бы не необходимость учета пересечений SID и STAR при расчете пропускной способности воздушного пространства. Точки объединения STAR (ветвления SID) обозначим вертикальными пунктирными линиями (рис. 3), полагая, что левее этих линий соответствующие траектории движения ВС совпадают. Ниже мы покажем, как вопрос учета слияния (ветвления) решается за счет расширения матричного поля критических перекрытий квантовых состояний ВС в аэродромной воздушной зоне.
2.3. Матричное поле критических перекрытии квантовых состоянии в аэродромной воздушной зоне
В расширенной квантовой модели наличие разных маршрутов ВС можно учесть, введя блочную матрицу R с блоками Rnmр,р2) = RR(p1,p2) для маршрутов n и m , чьи SID и STAR
на интервале указанных расстояний совпадают. Так, в примере, проиллюстрированном на рис. 3, STAR первого и второго маршрутов сливаются при р = р4, а соответствующие SID ветвятся при р = р3. Поэтому матричное поле критических перекрытий для рассматриваемого случая в точке (р2,р4) будет представимо матрицей R(p2,p4) вида:
Civil Avition High TECHNOLOGIES
Vol. 19, № 02, 2016
RP2P4) =
fR R
r RR 0
0 0 0
v у
Здесь R - матрица перекрытий пары маршрутов, а 0 - нулевая матрицы.
STAR третьего маршрута пересекается с SID в точке (р2,р3). Для нее R(p2,p3) имеет вид:
R(P2' Р4) =
^0 0 0 ^ 0 0 0
v 0 0 R ,
На рис. 3 такие критические перекрытия соответствуют серым столбцам.
Существенно, что матричное поле критических перекрытий зависит лишь от геометрии линий SID и STAR в обычном координатном пространстве и не зависит от времени.
3. УРАВНЕНИЯ СОВМЕСТНОЙ ДИНАМИКИ ПОТОКОВ ВЗЛЕТАЮЩИХ И ЗАХОДЯЩИХ НА ПОСАДКУ ВС
В расширенной модели воздушного пространства со временем ВС не только смещаются, приближаясь или удаляясь от ВПП, но и меняют свой высотный эшелон, оставаясь, однако, в полосе своего маршрута. Для описания этого процесса, рассматриваемого как результат управления, введем матричные поля U± (t,р), такие, что:
у ±(г+1,р±1) = и ±( ,р)-¥±{* (2)
Для описания движения конкретного ВС необходимо учитывать все предшествующие моменту t управления. Это удобно реализовать, если перейти к новым переменным:
£ = (-р/т) - для взлетающих и п = [?-(Ь-р)/т] - для заходящих на посадку ВС. Здесь Ь - полная длина квантованной зоны, а т - период дискретизации времени (1). Тогда (2), например, для взлетающих ВС можно переписать в виде:
¥+ (£,Р +1) = и + (£,р) (£, р), (3)
где матрица и + (£,р) при фиксированном значении £ задает все управления в воздушной зоне аэродрома для ВС, взлетевшего с ВПП в момент времени £ .
Требование бесконфликтности пролета ВС во встречных потоках можно записать в виде
(+ (£,р), Ир! ,р2 (п,рг)) = 0, (4)
где параметры £ и п связаны простым соотношением £ + р1 = п + Ь -р2, вытекающим из условия одновременности входа конфликтующих ВС в зону критического перекрытия. Задание пар
Том 19, № 02, 2016_Научный Вестник МГТУ ГА
Vol. 19, № 02, 2016 Civil Avition High TECHNOLOGIES
(Е,,Ц) является необходимым при корреляции моментов взлета/посадки и номеров трасс прилета/вылета.
При невыполнении условия (4) одно из ВС конфликтующей пары должно быть исключено из рассмотрения, а соответствующее управление, например, U + - обнуляется.
4. ИМИТАЦИОННЫЙ ЭКСПЕРИМЕНТ
Описанные в первой части работы элементы абстрагирования, идеализации и формализации позволили существенно снизить уровень сложности этапа имитационного моделирования.
Для его реализации на основе анализа навигационных карт конкретного аэродрома надо построить расширенное квантованное пространство воздушной зоны и задаться статистикой прилета и вылета ВС по различным трассам. Будем полагать далее все это известным.
Расширим далее количество ограничений, накладываемых на движение ВС. Будем полагать следующее.
1. Заходящие на посадку (взлетающие) ВС не могут увеличивать (уменьшать) свой высотный эшелон.
2. При смене конкретным ВС эшелона по высоте оказываются «занятыми» два находящиеся друг над другом состояния: hi и hi_1 - для заходящего на посадку и hi, hi+1 - для взлетающего, где hi - эшелон, соответствующий началу маневра.
3. Случайными во встречных потоках ВС являются: трассы прилета и вылета, моменты смены высотного эшелона и эшелоны входа и выхода ВС в зону критических перекрытий.
4. Будем полагать для простоты, что аэродром имеет две независимые ВИН для взлета и посадки.
5. Будем полагать что к аэродромному пространству подходят четыре воздушных трассы.
6. Будем полагать ветровое поле фиксированным, не меняющимся во времени. Другими словами, направления взлета и посадки на ВИН не меняются и используется лишь одна структура взаимного расположения SID и STAR.
Последнее допущение связано с тем, что в зависимости от направления ветра может меняться и направление взлета/посадки на ВПП, а следовательно, и геометрия всей картины пересечения маршрутов ВС (матричное поле критических перекрытий). Для корректного учета влияния этого эффекта на статистические результаты имитационного эксперимента необходимо использование розы ветров аэродрома.
Имитационный эксперимент заключатся в следующем. Случайным образом формируются два потока ВС, входящих в зону критических перекрытий. Для каждого ВС из этих потоков случайным образом выбираются (из допустимых значений): высотные эшелоны влета и вылета из зоны, трассы прибытия и убытия ВС (в модели их 4). Строятся случайные «траектории» подъема и спуска ВС. По матрице критических перекрытийТ1(рх,р2) с учетом выбранных ранее трасс прибытия и убытия восстанавливаются пары р,р2) - «координаты» точек конфликтов. С учетом соотношения £ + р1 =п + L _р2, связывающего моменты влета ВС в анализируемую зону, устанавливается факт наличия или отсутствия конфликта между ними. При наличии конфликта одно из конфликтующих ВС (либо взлетающее, либо садящееся) случайным образом удаляется из потоков. Траектория оставшегося ВС проверяется далее на наличие критических перекрытий. Если в конечном итоге это ВС достигает противоположной границы зоны перекрытий, то в счетчик позитивных исходов записывается единица, в противном случае - нуль.
Civil Avition High TECHNOLOGIES
Vol. 19, № 02, 2016
Эксперимент проводится N раз. Если показания счетчика - N1, то пропускная способность аэродрома на взлет/посадку или 3 - максимальная интенсивность встречных потоков ВС,
бесконфликтно проходящих воздушное пространство аэродромной зоны, определяется по простой формуле
где Т - период квантования (1), а Т - единица измерения времени. Описанный имитационный эксперимент позволяет искусственно менять в модели q - количество точек критических перекрытий, анализируя зависимость пропускной способности аэродромной воздушной зоны от сложности топологии трасс - 3 = 3(д). На рис. 4 представлены результаты описанного эксперимента.
Число критических перекрытий квантованных состояний в аэродромной воздуш-нойзоне менялось в эксперименте с помощью вариации отдельных элементов матрицы И. На рисунке приведены усредненные (по наборам исключаемых точек перекрытия) результаты. Следует отметить, что приведенные расчеты носят лишь иллюстративный характер, поскольку не отражают топологию реальных воздушных трасс конкретного аэродрома.
Рис. 4. Зависимость коэффициента N^1N от числа критических перекрытий в аэродромной воздушной зоне
ЗАКЛЮЧЕНИЕ
В работе рассматривается задача об оценке пропускной способности воздушной аэродромной зоны, т.е. об оценке максимальной интенсивности потоков взлетающих и заходящих на посадку ВС, бесконфликтно проходящих воздушную зону. Решение задачи разбивается на два этапа. На первом этапе с использованием методов абстрагирования, идеализации и формализации строится «квантованная» модель воздушного пространства аэродромной зоны [4-6]. В представляемой работе проведено расширение этой модели, учитывающее наличие нескольких трасс прилета и вылета ВС для анализируемой зоны, что приводит к соответствующему усложнению топологической картины траекторий движения ВС.
Этот первый этап существенно упростил дальнейшее решение задачи - переход к имитационному моделированию, который позволил дать оценку для пропускной способность аэродромной воздушной зоны с учетом пересечений маршрутов прилета и вылета ВС. В статье приводятся результаты проведенного имитационного эксперимента, иллюстрирующие зависимость пропускной способности зоны от количества критических перекрытий маршрутов ВС.
Vol. 19, № 02, 2016
Civil Avition High TECHNOLOGIES
СПИСОК ЛИТЕРАТУРЫ
1. Приказ Министерства транспорта РФ от 24.02.2011 № 62. Об утверждении методики расчета технической возможности аэропортов и порядка применения методики расчета технической возможности аэропортов. - М.: РГ, 2011.
2. IATA Airport Development Reference Manual. 9th Edition. 2004.
3. Дектярев O.B., Зубкова И.Ф. Ускоренное моделирование процессов организации потоков и управления ВД. Аналитический обзор. - М.: НИЦ ЕосНИИАС, 2008.
4. Кузнецов В.Л., Чепурина А.А. Об одной аналитической модели пропускной способности аэропорта // Научный вестник МЕТУ ЕА. - 2011. - № 169. - С. 94-98.
5. Кузнецов В.Л., Филонов П.В. Моделирование пропускной способности аэродрома с учетом конфликтных ситуаций в воздушном пространстве в формализме ферми-систем // Известия РАН. Теория и системы управления. - 2012. - № 3. - С. 136-144.
6. Кузнецов В.Л., Филонов П.В., Чепурина А.А. Метод построения системы квантовых состояний в модели расчета пропускной способности аэропорта // Научный вестник МЕТУ ЕА. - 2012. - № 184. - С. 23-28.
COMBINED SIMULATION-ANALYTICAL MODEL OF AERODROME AIRSPACE CAPACITY
Kuznetsov V.L., Gorchakov A.I.
An approach to the problem of estimating an aerodrome airspace capacity based on the idea of airspace «quantization» and transition to discrete time while imposing restrictions on the filling of quantum states by aircrafts has been developed in this paper. It is shown that such a formalism can significantly reduce the complexity of aircrafts movement simulation in the airfield zone, to take into account the restrictions imposed on the aerodrome airspace capacity by intersections of arrival and departure routes. The paper presents the results of the simulation experiment, illustrating the dependence of the aerodrome airspace capacity upon the number of aircraft routes intersections, [email protected].
Key words: aerodrome airspace capacity, airspace quantization, simulation experiment.
REFERENCES
1. On approval of the methodology for calculating the technical capacity of the airport and the application of methods for calculating technical ability airports. Order of Transport Ministry of the Russian Federation of 24.02.2011 № 62. M.: RG. 2011. (In Russian).
2. IATA Airport Development Reference Manual. 9th Edition, 2004.
3. Dektjarev O.V., Zubkova I.F. Faster modeling of flow management and air traffic control. Analytical overview. - M.: NIC-GosNIIAS, 2008. (In Russian).
4. Kuznetsov V.L., Chepurina A.A. About one analitical model of the airport throughput. Scientific Herald of MSTUCA. 2011. No. 169. Pp. 94-98. (In Russian).
5. Kuznetsov V.L., Filonov P.V. Simulating Aerodrome Handling Capacity in Terms of Fermi Systems Given Airspace Conflicts. Journal of Computer and System Sciences International. 2012. Vol. 51. No. 3. Pp. 468-476.
6. Kuznetsov V.L., Filonov P.V., Chepurina A.A. The quantium states system constructing method in aiport throughput computation model. Scientific Herald of MSTUCA. 2012. No. 184. Pp. 23-28. (In Russian).
СВЕДЕНИЯ ОБ АВТОРАХ
Кузнецов Валерий Леонидович, 1949 г.р., окончил МГУ им. М.В. Ломоносова (1972). Доктор технических наук, заведующий кафедрой прикладной математики МЕТУ ЕА, автор более
Civil Avition High TECHNOLOGIES
Vol. 19, № 02, 2016
100 научных публикаций. Область научных интересов - методы математического моделирования в задачах распространения излучения в пространственно неоднородных случайных и периодических средах, УВД, безопасность полетов, электронный адрес: [email protected].
Горчаков Антон Игоревич, 1996 г.р., окончил МГТУ ГА (2015). Область научных интересов - УВД, имитационное моделирование.