УДК 621.306
Ю.Н. Макаров, А.А. Строцев ВЫБОР ВАРИАНТА СИСТЕМЫ ФУНКЦИОНАЛЬНОГО МОНИТОРИНГА С УЧЁТОМ ВЛИЯНИЯ НА НЕЁ СТОХАСТИЧЕСКИХ И НЕОПРЕДЕЛЁННЫХ ФАКТОРОВ
Массовое развитие различных направлений мониторинга технического состояния обусловлено переходом к стратегии эксплуатации сложных технических систем по их фактическому состоянию. Одним из таких направлений является функциональный мониторинг [1], [2]. На основе теоретико-игрового подхода предложена методика выбора структуры систем функционального мониторинга на этапе проектирования сложных технических систем с учётом влияния на них в процессе эксплуатации стохастических и неопределённых факторов.
Функциональный мониторинг; теоретико-игровая оптимизация; матричные игры с ограничениями.
Yu.N. Makarov, A.A. Strocev CHOICE OF THE VARIANT OF FUNCTIONAL MONITORING SYSTEM TAKING INTO ACCOUNT INFLUENCE ON IT OF STOCHASTIC AND UNCERTAIN FACTORS
Mass development of various directions of a technical condition monitoring is caused by transition to operation strategy of difficult technical systems on their actual state. One of such directions is functional monitoring [1], [2]. On the basis of the game-theoretical approach the technique of a structure choice of functional monitoring systems at a design stage of difficult technical systems taking into account influence on them while in service stochastic and uncertain factors is offered.
Functional monitoring; game-theoretical optimization; matrix games with contingencies.
Внедрение стратегии эксплуатации сложных технических систем (СТС) по их фактическому состоянию на основе функционального мониторинга в настоящее время затрудняется рядом ограничений системного характера. Прежде всего, таким ограничением является отсутствие соответствующих средств оценивания технического состояния в составе СТС, спроектированных до принятия новой стратегии эксплуатации, что не позволяет охватить весь необходимый спектр измерений [3].
Действительно, во время проектирования невозможно учесть все факторы, так как заранее неизвестно в каких условиях будет функционировать сама СТС и система функционального мониторинга её технического состояния. Ряд факторов могут быть учтены как детерминированные, однако большинство из них являются неопределёнными и стохастическими. Стохастический или неопределённый их характер может определяться не только с технической надёжностью и другими показателями качества, но и характеристиками, связанными с действиями конкурентов, т.е. учитывающих конкурентную среду.
С другой стороны, так как следует учитывать неопределённые факторы, то модель рассматриваемой проблемной ситуации является моделю принятия решений в условиях конфликта (игры с природой или конкурентом). Ограничимся рассмотрением двух сторон конфликтной ситуации и наличия конечного числа стратегий у них. Формализация рассмотренной вербальной постановки задачи может быть проведена на основе применения классической модели смешанного расширения матричных игр с ограничениями в виде равенств и неравенств. Примеры построения таких моделей для задач поиска и устранения неисправностей приведены, например, в [4], [5].
Среди класса СТС рассмотрим такие, для которых характерно малое (ограниченное) число повторений условий принятия решений (что может быть связано с их малым числом, разными условиями функционирования отдельных групп СТС и т.д.). Примером таких СТС является уникальные агрегаты механического оборудования стартового комплекса - компонента наземной инфраструктуры космического комплекса.
Следовательно, при построении методики теоретико-игрового синтеза систем функционального мониторинга СТС следует учесть малое число реализации игровой ситуации, что, например, возможно на основе построения моделей смешанного расширения матричных игр неклассического типа, теоретические основы которого рассмотрены в [6], [7].
Таким образом, разработка методики теоретико-игрового синтеза систем функциональной диагностики СТС на основе построения и применения моделей смешанного расширения матричных игр неклассического типа с ограничениями является актуальной задачей.
Рассмотрим процесс проектирования СТС, на котором необходимо учесть влияние стохастических и неопределённых факторов при работе системы её функционального мониторинга.
Пусть А - матрица игры, элементы которой , i = 1, n, j = 1, m характеризуют эффективность процесса функционального мониторинга.
При этом строки матрицы характеризуют стратегии лица, принимающего решение (ЛПР), в части структур системы функционального мониторинга СТС, рассматриваемых на этапе её проектирования. В частности, на этом этапе необходимо провести обоснование объёма и вида информации для адекватного оценивания технического состояния подсистем СТС. Подходы к решению этой проблемы для оценивания технического состояния уникального механического оборудования стартовых комплексов в процессе длительной эксплуатации в детерминированной постановке раскрыты, например, в [3]. Отметим, что стратегиями ЛПР на этапе проектирования могут быть связаны ограничения, учитывающие известные требуемые значения вероятностей применения ряда структур построения СТС.
Столбцы матрицы характеризуют условия функционирования СТС, определяемые действиями случайных и неопределённых факторов. Под случайным понимается фактор, который имеет вероятностное описание, а под неопределённым -фактор, не имеющий вероятностного описания. При этом ряд случайных факторов могут иметь области неопределённости вероятностного описания.
В отличие от моделей, построенных в [4], [5], в качестве базовой выберем модель смешанного расширения матричной игры неклассического типа
Г а (HLp) =< Ч’ >, для которой, в соответствии с [6]:
♦ функция выигрыша первого игрока, которую второй игрок стремиться минимизировать:
n m m
Hи?(X, Y ,p)=pZZ jj^+(1 _ p)Z Vjmax au, (1)
i=1 j = 1 j=1 i
где p - параметр критерия Ходжа-Лемана (критерий HL) (p e [0,1] , при предельных значениях весового коэффициента p модель предпочтения ЛПР соответствует или критерию Байеса-Лапласа (критерий BL), или минимаксному (ММ) (максимин-ному) критерию для каждого из игроков; выбор промежуточного значения коэффициента p обеспечивает построение модели предпочтения ЛПР, учитывающей ориентацию как на средний результат, так и на получение наилучшего гарантированного результата, что обеспечивает учёт числа реализаций игровой ситуации;
♦ множества е Я" и Ят являются множествами смешанных
стратегий X = (^ %2...%„ )Г , ^ = (Л\ Л2 ■Лт )Т в исходной матричной игре соответственно для первого и второго игроков, т.е. Ме и М пред-
ставляют собой симплексы:
М. = j X = (,1 ■ ■ . )т: , > о, і = 1,n, Z, = 1
і=1
(2)
М, = \ Y = ( j ■ ■ jm )T : jj > ° j = 1, m Z = 1 Г .
j=1
(3)
Пусть заданы вероятности функционирования СТС в ряде условий (т.е. задано вероятностное описание случайных факторов), области неопределенности вероятностного описания некоторых случайных факторов и требуемые значения вероятностей применения ряда структур построения СТС. Без ограничения общности будем полагать:
и. > j, j = m’+1,m”, j =rf* j = m” + 1,m, m'<m”, (4)
^=£? , I = "' + 1," . (5)
Задача заключается в определении для неопределенных факторов модели та-
^ **1^ Н ^ Г Я" *
стратегий: 77. = п., ] = 1, т , т < т, = д. ,
ких элементов смешанных
i = 1, nn, n f < n , которые обеспечивали выполнение условия
f n m m
pZ Z Vja^, + (i - p)Z Vjmax а
max min
. Ь Ъ V і=1 і=1
і=1,п /=1,п ' 4
j=1
г
= min max
j і=fe V і=1 j=1
PZ Z +(1 - P)Z ,jmax a
Л
j =1,m і =1, n
при ограничениях (4), (5) и
j=1
(6)
(7)
Z jj = 1, jj > о, j = 1, m”.
j=1
(8)
Построим двойственные задачи линейного программирования для решения (4)-(8). В соответствии, например, с подходом, рассмотренным в [5], сформируем двойственные задачи, имеющие одно и то же значение оптимизируемых функций. Для этого введем следующие обозначения:
nm
m
і=1
m
A =
A11 A12 A13
A21 A22 A2
зУ
X= Xі X
'J,7 = (7r YT YTJ,
где X = (, ■ ■ ■ . )т, X = (^ ■ ■ . )T, Y = (j ■ ■ ■ j )T, Y = (j+1 ■ ■ ■ jOT = Y = (jm+1 ■ ■ ■ )T, dim A11 = n' x m', dim A12 = n' x (m” - m'),
dim A13 = n' x (m - m”), dim A21 = (n - n') x m', dim A22 = (n - n') x (m” - m'), dim A3 = (n - n ”) x (m - m”),
при этом, учитывая (4), (5), (7), (8), будем полагать X = Xгр, Y > Yгр , Y = Y
гр
=гр
ETX + EIX = 1, EmY + Em-mY = 1 -E^Y , X* = (^■ ■ )
гр N T
y гр=(vm+1 .. . vz )T, y гр=v .. . vm )T.
Представим (6) в соответствии с введёнными обозначениями:
max min
X Y, Y
P\ XT X lA\ YT Y Y
+
+ (1 -P)
=гр
ET,A™s:xY + E\ ,A,maxY + ET ,A,maxY
m diag m -m diag m-m diag
= min max
Y, Y X
P\ XT X lA\ YT Y Y
+
+ (1 -P)
ET,A™s:xY + ET„ ,A,maxY + ET ,A,maxY
m diag m -m diag m-m diag
Л
У
=гр
(9)
где A,maxг - диагональная матрица с элементами А^Хг.. = max а , j = 1,m'
Adiag - диагональная матрица с элементами Adiag jj = max ajf, j = m ' +1, m ",
i j
=max i = max i
Adiag - диагональная матрица с элементами Adiag jj = max аг. m,+j,
j = m" +1, m, Eg - вектор с единичными элементами, dim Е = g. Рассмотрим левую часть равенства (9). В ней
min
Y, Y
+ (1 -P)
ET,AmiY + ET„ ,A,maxY + ET ,A,maxT
m diag m -m diag m-m diag
У
представляет целевую функцию задачи линейного программирования, в которой
коэффициенты зависят от элементов вектора X. Запишем эту задачу и двойственную к ней.
T
T
T
T
Вспомогательная прямая задача 1: найти
min
Y Y
P
~ т =адГ
X An + X A21
+(і - p) ETm ' Amxі
+
P Xta12+X- A221+(і-p)En"-n.AdI'::xі
Y+
+
P
~т =грт Л X TA13 + X A23
V
+(і -PE-m Admaxі
Y
в условиях ограничений
(10)
ET'7 + E\ ,Y = 1 -ET Y , Y > Yгр, Y > 0 ,, Y > 0.
m
T m -m
T m-m
= гр
(11)
где 0 - вектор с нулевыми элементами, dim 0 = p. Вспомогательная двойственная задача 1: найти
max Ці - ET "Yгр V,
_ y Л п-п Г 0
s.,s,X 4 '
+ Yгр S +
+
P
~т =грт Л
X TA13 + X A23
+ (1 -P)ETm m Л211
Y
(12)
при ограничениях
Em So < PAX + A^* )+ (1 - P) AdagEm ',
Em"-m' So +S < P(A’2-Y + AT2X” )+ (1 - P) ArjEm-m
s > 0,
(13)
(14)
(15)
Задача (12) - (15) должна решаться при дополнительных ограничениях, учитывающих (2) и (5).
Окончательно, с учетом выделения искомого вектора переменных X , для первого игрока получим: найти
тах / (.0, *, X) = тах { (1 - Ктт_т, 7)у0 + У+ {/ЗУ ^' ЛТ 1X +
^ГХ ^X I I )
+
P
Ґ=грт Л
-23
XA
+ (1 -P)El - m "AZl1
Y
при ограничениях
Em So - PATX < PAX + (1 - P) A2;En',
Em"-m So + S - PA1T2X < PAT2Xа + (1 - P)А^,,,-m
(16)
(17)
(18)
T
У
К
У
ЕТХ = 1 - Ет ,Хгр, X > 0 ,, ^ > О, ,.
т
''п-п
(19)
Аналогично получим задачу линейного программирования для определения У, У, двойственную к (16) - (19). Для этого теперь рассмотрим правую часть равенства (9). В ней
тах
лг
т \ ! т
= гРт | .[ ~т —т =гРт
Р\ Xі X \Л\ У1 У У
+
+(1 -Р)
ЕтЛ тх'У + Ет. Лтх’У + Ет „Л,тах'У
т dгag т -т т-т
представляет целевую функцию задачи линейного программирования, в которой коэффициенты зависят от элементов векторов У, У . Запишем эту задачу и
двойственную к ней.
Вспомогательная прямая задача 2: найти
тах
X
рхт|лпу + лу+ЛУ \+Рх +Л2у+л23у 1 +
+ (1 -Р)
в условиях ограничений
77 т у тах іУ 77т дшах іу, Т?т д тах ’у'
Ет'Лdiag У Ет"-т'Лdiag У Ет-т"Лdiag У
(20)
ЕтХ = 1 - Ет ,ХГ, X > О!.
п ™ ™ *
т ^П— П
Вспомогательная двойственная задача 2:
найти
Шп |(1 - Е„_„.Xг')д + РX
$,У ,У I
грт ( у — =гр
А,У + Л22У + Л23У \ +
+ (1 -Р)
при ограничениях
=гр
Т^т у шах іу т^т д шах гу^Т^т д шах іу Ет'Лdiag У Ет" -т'Лdiag У Ет-т"Лdiag У
Е„.д >Р\ЛиУ + Л12У + Л13У
-гр
(21)
(22)
(23)
Задача (22), (23) должна решаться при дополнительных ограничениях, учитывающих (3) и (4).
Окончательно, с учетом выделения искомых векторов переменных У , У , для второго игрока получим: найти
тМя,У,У) = тп{(\-ЕІпХ^)д+(рх^А2,Ні-Р)ЕТт.АГ;)ї +
ГП1П
+(РХгрТА22 + (1-^"-„АГ) У +(РХгрТА23 + (1-Р)Е1т,АЩ£') Угр } > (24)
при ограничениях
т
Г
Еп,д-РА?-РАг? >РАУ^ , (25)
ЕТт ? + ЕТт„-т. 7 = 1 - КТт_т У ^, (26)
7 > 7гр, ~ > От , 7 > От .-т - . (27)
Таким образом, методика теоретико-игрового синтеза систем функционального мониторинга с учётом выбора модели смешанного расширения матричной игры неклассического типа, содержательного характера матрицы игры (большие или меньшие значения её элементов предпочтительнее для первого игрока (ЛИР)) может быть представлена в виде следующей последовательности действий:
1. Выбор (построение) обобщённого показателя эффективности системы функционального мониторинга.
2. Формирование конечных множеств стратегий конкурирующих сторон.
3. Определение матрицы игры.
4. Выбор (или построение) модели смешанного расширения матричной игры неклассического типа как модели принятия решений в условиях неопределённо-сти, соответствующей п.3, с функцией игры, зависящей от числа разыгрываний игровой ситуации, т.е. формирование выражений вида (6)-(8).
5. Определение вероятностей функционирования СТС в ряде условий (т.е. задание вероятностного описания случайных факторов), области неопределённости вероятностного описания некоторых случайных факторов и требуемых значений вероятностей применения ряда структур построения системы функционального мониторинга, т.е. формирование ограничений вида (4), (5).
6. Иоследовательное формирование вспомогательных прямых и двойственных задач линейного программирования для определения оптимальных смешанных стратегий игроков.
7. Учёт дополнительных ограничений, учитывающих (2)-(5), и формирование окончательного вида двойственных задач линейного программирования для определения оптимальных смешанных стратегий сторон в игре с ограничениями-равенствами и ограничениями-неравенствами и функцией игры, зависящей от числа разыгрываний игровой ситуации.
8. Решение полученных в п.7 двойственных задач для рассматриваемых условий по числу реализаций игровой ситуации, анализ полученного решения (полученных решений).
9. Формирование на основе реализации п.8 рекомендаций по структурам СТС (структурам системы функционального мониторинга).
Рассмотрим пример, реализации п.8 методики для полученных выше двойственных задач (16)—(19) и (24)-(27). Иусть задана матрица, элементы которой представляют собой значения обобщённого показателя эффективности системы функционального мониторинга в случае /-го варианта структуры СТС I = 1,4 и]-го неблагоприятного фактора (условий функционирования, действия конкурентной среды и т.д.)
40 35 50 84 65 30"
37 65 28 72 48 64
46 85 73 59 88 53
4 6 65 48 23 54 49,
и ограничения для которых п' = т' = 2, т" = 4, Xгр = (0,15 0,05)Т,
7грТ = (0,15 0,05)Т, ¥грТ = (0,05 0,10)Т.
Построение обобщённого показателя эффективности системы функционального мониторинга и определение элементов матрицы игры может быть осуществлено на основе подходов, рассмотренных в [9].
В результате решения задач (16)—(19) и (24) - (27) при ( = 1,0 получим:
/1>0. ^ = 48,65 > ^0 = 46,23 > = 5,48 >
- . Цч гр /а^1 Цс ут
52 =30,08, X =(0,73 0,07) , д=45,82, 7 =(0,47 0,18) ,
¥ =(0,15 0,05)Т . Для полученных смешанных стратегий игроков математи-
ческое ожидание значения эффективности системы функционального мониторинга
Ґ т\ ( т\
_ . ______рг)1- ~ . ____лт
Т
=гр-
X*1 X
J
А
=гр
= 48,64,
а вероятность получения для второго игрока результатов худших, чем м
3=1,0,
Р = 0 33 1 3=1,0 0,33.
При 3 = 0,5 имеем: //(^Л15= 58,37,
ї0 = 55,45, ї = 7,73 , ї2 = 25,15, X =(0,80 0,00)г, q = 23,93
~* =(0,65 0,00)Т, ¥ * =(0,15 0,05)Т, М^ 5 = 49,19, ^=0 5 = 0,29.
Равенство целевых функций задач линейного программирования (16) и (24) являются признаком наличия седловой точки, а найденные смешанные стратегии определяют ситуацию равновесия. Кроме того, поскольку в модели игры
Г а (Н1}р) второй игрок определяет число разыгрываний, то при уменьшении значения параметра ( с 1,0 до 0,5 (т.е. при уменьшении числа разыгрываний) увеличивается его средний проигрыш (от М ^=1 0 = 48,64 до М5 = 49,19), но уменьшается вероятность в единичном разыгрывании получить результат хуже, чем его математическое ожидание (с Р^=; 0 = 0,33 до 5 = 0,29).
При ( = 0,5 в случае ¥гр = (0 0)Т , т.е. при отсутствии статистической
информации об интервалах неопределённости стохастического описания неопределённых факторов, в результате решения соответствующих задач получим:
М%°Т = 46,70.
Таким образом, учёт дополнительных ограничений позволил в условиях примера повысить математическое ожидание обобщённого показателя эффективности системы функционального мониторинга на 5,3 %, а полученное решение матричной игры с ограничениями полностью соответствует всем необходимым и достаточным условиям ситуации равновесия. Предложенная модель может быть применена для оптимизации структуры системы функционального мониторинга на этапе проектирования в условиях действия во время эксплуатации сочетания слу-
чайных и неопределённых факторов. Область применения методики может быть распространена для решения задач обоснования стратегий развития в условиях действия случайных и неопределённых факторов, в частности, для обоснования космических программ и проектов [], [].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Прохорович В.Е. Прогнозирование состояния сложных технических комплексов. - СПб.: Наука, 1999. - 158 с.
2. Петров Г.Д. Методологические аспекты обеспечения долговечности механического оборудования стартовых комплексов на основе функционального мониторинга: монография. - СПб.: ВКА им. А.Ф. Можайского, 2005. - 201 с.
3. Макаров Ю.Н., Лухвич А.А., Шипша В.Г. и др. Актуальные проблемы неразрушающего контроля качества космической техники: Монография. - СПб.: Альтеор, 2008. - 333 с.
4. Строцев А.А., Синицын С.В., Жадько А.А. Методика теоретико-игровой оптимизации алгоритма контроля на основе модели смешанного расширения матричной игры с ограничениями // Известия ЮФУ. Технические науки. - 2008. - № 11 (88). - С. 66-70.
5. Строцев А.А., Оганесян А.Л., Григорян М.А. Теоретико-игровая оптимизация алгоритмов контроля сложной системы на основе классических моделей матричных игр с ограничениями-неравенствами // Известия ЮФУ. Технические науки. - 2009. - № 2 (91). - С. 244-248.
6. Строцев А.А. Построение смешанного расширения матричной игры "неклассического" типа // Известия АН. Теория и системы управления. - 1998. - № 3.- С. 119-124.
7. Макаров Ю.Н., Строцев А.А. Методология исследования сложных организационно -технических систем, функционирующих в конкурентной среде при ограниченных ресурсах: Монография. - Ростов-на-Дону: Изд-во ЮФУ, 2010. -132 с.
8. Оуэн Г. Теория игр: Изд. 3-е. - М.: Изд-во ЛКИ, 2000. -216 с.
9. Асланов М.А., Кузнецов В.В., Макаров Ю.Н., Мальчевский А.А., Шатраков А.Ю. Системный анализ и принятие решений в деятельности учреждений реального сектора экономики, связи и транспорта. - М.: ЗАО «Издательство «Экономика», 2009. - 366 с.
10. Давыдов В.А., Макаров Ю.Н., Мальченко А.Н., Пайсон Д.Б. Новые концептуальные методологические подходы к проблемам формирования оптимального технического и технологического базиса программно-целевого планирования в создании и развитии ракетно-космической техники. - М.: ЗАО «ЭНЦИТЕХ», 2006.
11. Макаров Ю.Н. и др. Перспективы развития совмещённых наукоёмких технологий. Исследования вопросов совершенствования технико-экономического обоснования космических программ и проектов. - М.: ЗАО «ЭНЦИТЕХ», 2009.
Статью рекомендовал к опубликованию д.т.н., профессор А.Е. Панич.
Макаров Юрий Николаевич
Федеральное космическое агентство Российской Федерации.
E-mail: arm415@ roscosmos.ru.
107996, ГСП-6, г. Москва, ул. Щепкина, 42.
Тел.: 84956319322.
Начальник сводного управления организации космической деятельности; к. т. н.
Строцев Андрей Анатольевич
Ростовский военный институт РВ.
E-mail: [email protected].
344068, г. Ростов-на-Дону, ул. Нариманова 76, кв. 37.
Тел.: 88632926242.
Начальник кафедры систем автоматической подготовки и пуска ракет, д.т.н.; доцент.
Makarov Yury Nikolaevich
Federal Space Agency, Russian Federation.
E-mail: arm415@ roscosmos.ru.
42, Schepkina Street, Moscow, GSP-6, 107996, Russia.
Phone: +74956319322.
Head of the Consolidated Space Activity Organization Department; Cand. of Eng. Sc.
Strotsev Andrey Anatol’evich
Rostov Military Institute.
E-mail: [email protected].
76, Narimanov Street, Rostov-on-Don, 344068, Russia.
Phone: +78632926242
The Chief of Faculty of Systems of Automatic Preparation and Start-up of Rockets, Dr. of Eng. Sc.; Associate Professor.
УДК 625.7/.8-047.58
Д.А. Скоробогатченко
МЕТОДИКА МОДЕЛИРОВАНИЯ ЭКСПЛУАТАЦИОННОГО СОСТОЯНИЯ АВТОМОБИЛЬНЫХ ДОРОГ НА ОСНОВЕ НЕЧЕТКИХ НЕЙРОННЫХ
СЕТЕЙ
Целью работы является создание методики моделирования изменения эксплуатационного состояния автомобильных дорог с учетом информации качественного характера. Для достижения поставленной цели решается задача синтеза аппарата нечетких множеств. В результате авторами реализована объективная методика моделирования эксплуатационного состояния автомобильных дорог на основе нечетких нейронных сетей.
Автомобильные дороги; моделирование эксплуатационного состояния; нечеткие нейронные сети.
D.A. Skorobogatchenko
METHOD OF MODELING THE OPERATIONAL STATUS OF ROADS BASED ON FUZZY NEURAL NETWORKS
The aim is the creation methods of modeling changes in operational condition of roads, taking into account qualitative information. To achieve this goal we solve the problem of synthesis of the apparatus offuzzy sets. As a result, the authors implemented an objective method of modeling the operational status of roads based on fuzzy neural networks.
Roads; modeling the operational status; fuzzy neural networks.
Эксплуатационное состояние автомобильных дорог (ЭС АД) определяет развитие экономики страны. Однако в связи со снижением объемов финансирования, фактическое ЭС АД не соответствует нормативному. Выходом из сложившейся ситуации является создание систем управления, основу которых составляет моделирование ЭС АД [1]. В связи с тем, что специфика дорожной отрасли предполагает наличие большого объема информации, представленной вербально, для моделирования ЭС АД предлагается использовать теорию нечетких множеств [2].
Процесс изменения ЭС АД во времени образует следующие подмножества (рис. 1):
совокупность переменных, характеризующие дорожные работы:
xt е X, i = 1, q
совокупность переменных, характеризующих воздействие на ЭС АД природноклиматических и транспортных факторов:
V е V, l = 1, r совокупность переменных, описывающих начальное ЭС АД:
h е H, к = 1, m совокупность переменных, описывающих конечное ЭС АД: