СПИСОК ЛИТЕРАТУРЫ
1. Рожкова С.В. Оптимальная непрерывно-дискретная передача сигнала по каналам с памятью при наличии бесшумной обратной связи // Известия Томского политехнического университета. - 2012. - Т. - №5. - С. 6-10.
2. Липцер Р.Ш., Ширяев А.Н. Статистика случайных процессов.
- М.: Наука, 1974. - 696 с.
3. Демин Н.С., Рожкова С.В., Рожкова О.В. Обобщенная скользящая экстраполяция стохастических процессов по совокуп-
ности непрерывных и дискретных наблюдений с памятью // Известия РАН. Теория и системы управления. - 2000. - № 4. -С. 39-51.
4. Gallager R.G. Information theory and reliable communication. -New York: Wiley, 1968 - 346 p.
Поступила 03.10.2012 г.
УДК 681.51
КОМПЬЮТЕРНАЯ МОДЕЛЬ СОПРЯЖЕНИЯ ЭЛЕМЕНТОВ В СТРОГО-ИЕРАРХИЧЕСКОЙ СТРАТИФИЦИРОВАННОЙ СИСТЕМЕ
С.Ш. Баласанян
Государственный инженерный университет Армении, г. Капан, Армения E-mail: suni-com@syunik.am
Разработана компьютерная модель сопряжения элементов в строго-иерархической стратифицированной системе, позволяющая свести имитацию сопряжения элементов смежных страт, т. е. формирование входного сигнала каждой страты на основании выходного сигнала предыдущей страты, к стандартной операции умножения матриц. Показано, что предложенная модель допускает распараллеливание процесса имитации передачи сигналов отэлементов каждой страты элементам последующей страты, что предоставляет возможность при использовании многопроцессорных компьютеров значительно ускорить процесс имитации. Применение предложенной модели сопряжения позволяет при проведении компьютерных экспериментов ввести необходимые изменения в матрицы сопряжения смежных страт в соответствии с конкурирующими вариантами структур исследуемой системы в процессе экспериментов, без изменения базовой имитационной программы.
Ключевые слова:
Страта, стратифицированный, контакт, имитация, эксперимент, сопряжение, структура, канал. Key words:
Strata, stratified, contact, simulation, experiment, link, structure, channel.
Одной из основных проблем, возникающих при программной реализации компьютерных стратифицированных моделей сложных технологических систем со многими состояниями [1], ориентированных на исследование эффективности их функционирования с учетом надежности элементов, является разработка процедуры имитации сопряжения элементов смежных страт в системе, позволяющей обеспечить удобство и простоту проведения компьютерных экспериментов и уменьшить при этом объем перепрограммирования.
Суть этой проблемы заключается в том, что при рассмотрении на компьютерной модели конкурирующих вариантов структуры исследуемой системы возникает необходимость в модификации исходной (базовой) имитационной программы, обусловленной изменением структуры связей между элементами смежных страт.
В настоящей работе предлагается компьютерная модель сопряжения элементов смежных страт в строго-иерархической стратифицированной системе (СИСС) [1], позволяющая успешно решить указанную выше проблему.
Предлагаемая модель базируется на математической модели сопряжения элементов в СИСС [1, 2],
рассматривающей взаимодействие между элементами в рамках механизма обмена сигналами [3, 4], который включает следующие составляющие:
1) процесс формирования выходного сигнала элементом, выдающим сигнал;
2) определение адреса передачи для каждой характеристики выходного сигнала;
3) прохождение сигналов по каналам связи и компоновка (формирование) входных сигналов элементов;
4) функционирование элемента, принимающего входной сигнал.
Первая и четвертая составляющие не рассматриваются в рамках модели взаимодействия, поскольку относятся к построению моделей функционирования элементов системы.
Третья составляющая механизма обмена сигналами связана с прохождением сигналов через реальные (неидеальные) каналы связи, формально рассматриваемые как самостаятельные элементы системы, функционирование которых сводится к соответствующим задержкам и искажениям сигналов. В этом случае фиктивные каналы, соединяющие элементы, передают сигналы мгновенно и без искажений, т. е. формально являются идеальными.
Адресация характеристик выходных сигналов и их компоновка во входные сигналы элементов осуществляется схемой сопряжения, формальное описание которой опирается на следующие предположения, вытекающие из наблюдений над формализованными в виде СИСС реальными сложными системами.
1. Взаимное влияние элементов внутри СИСС, а также взаимодействие между СИСС и внешней средой реализуется только посредством обмена сигналами.
2. Входной сигнал х'(/)еХ' (X'- множество входных сигналов страты 8'), поступающий к страте 8' в момент времени /, можно рассматривать как упорядоченную совокупность сигналов х'(/), /е/'={1,2,...р1}, одновременно поступающих на входы ее элементов 8', /е/':
х'(г) = (х'(г), X'«X ; (г)).
В свою очередь, входной сигнал х'(/), поступающий к элементу 8' в момент времени /, можно представить в виде упорядоченной совокупности элементарных сигналов х/(/), неИ'={1,2,...,н'}, одновременно поступающих на его вход:
х'(г) = (х£(г), х^(г),..., х^ (г)).
Аналогично выходной сигнал у'(/)еУ' (У' -множество выходных сигналов страты 8'), выдаваемый стратой 8' в момент времени /, можно рассматривать как упорядоченную совокупность сигналов уДО, /е/', одновременно выдаваемых элементами страты 8', а выходной сигнал у'(/) элемента 8' - как упорядоченную совокупность элементарных сигналов у;?(0, /еХ'={1,2,...,/'}, одновременно выдаваемых элементом 8'.
Предположения 1 и 2 основываются на том, что взаимодействие элементов в процессе функционирования реальной системы рассматриваются как результат совокупности воздействий (сигналов) каждого элемента на другие элементы. Сигнал представляет собой воздействие, представленное набором характеристик.
Для общности описания внешняя среда представлена в виде совокупности фиктивных страт 8"и 8м, причем страта 83 состоит из элементов 8-, /е/0={1,2,...р0}. Таким образом, описание взаимодействия СИСС с внешней средой включает лишь формальное описание структуры связей между стратой 8'СИСС и фиктивной стратой 83. Сигнал у°(0е У0 (У0 - множество выходных сигналов внешней среды), поступающий в СИСС из внешней среды, является выходным сигналом страты 80 и представляет собой упорядоченную совокупность выходных элементарных сигналов у,7°=у;0(0 е У° элементов 80, /е/0.
3. Элементарные сигналы передаются в СИСС независимо друг от друга по элементарным каналам связи, причем каждый элементарный канал, подключенный к выходу элемента 8', способен передавать только элементарные сигналы, принадлежащие одному из множеств У', не И'.
Предположение 3 сформулировано на основе того, что для функционирования реальной системы существенно лишь наличие на входе каждого элемента в данный момент времени соответствующего набора характеристик, описывающих сигнал. Данное предположение допускает следующую интерпретацию.
Вход каждого элемента СИСС представляет собой упорядоченную совокупность элементарных входов или входных контактов. Множество входных контактов элемента 8' обозначено К'. Каждый входной контакт кнеК' предназначен для приема только элементарных сигналов хнеХ'. Аналогично выход каждого элемента СИСС представляет собой упорядоченную совокупность выходных контактов. Каждый выходной контакт с/е С' (С' - множество выходных контактов элемента 8 ) элемента 8' предназначен только для выдачи элементарных сигналов у 'е У'.
Вход каждой страты 8' представляет собой упорядоченную совокупность К входных контактов всех элементов 8 , е / данной страты:
к' = у К'.
1е1'
Соответственно выход страты 8 представляет собой упорядоченную совокупность С выходных контактов всех элементов 8 , е / :
с' = у с'.
1е1'
Элементарные сигналы, выдаваемые данным выходным контактом страты 8'-1, передаются некоторому входному контакту страты 8 лишь в том случае, если эти контакты соединены между собой элементарным каналом связи.
4. К каждому входному контакту элемента или страты подключается не более одного элементарного канала связи; к каждому выходному контакту страты 8'-1 может быть подключено конечное число элементарных каналов при условии, что к входу одного и того же элемента страты 8 подключается лишь один из этих каналов. Следовательно, число элементарных каналов, подключенных к любому выходному контакту страты 8'-1, не может превышать числа элементов страты 8 , а число выходных контактов страты 8'-1 не может быть больше числа входных контактов страты 8 .
Данное предположение сформулированно с целью исключения неоднозначности входного сигнала за счет возможного появления на его входе в данный момент времени нескольких (несовпадающих) сигналов, поступающих из разных источников.
5. Элементы страты 8' могут принять только сигналы, выдаваемые элементами страты 8'-1. Предположение 5 вытекает из строго-иерархической конфигурации реальных систем, формализованных в виде СИСС.
Предположения 4 и 5 следует рассматривать как ограничение на структуру сети элементарных каналов связи в СИСС.
6. Сигналы передаются по элементарным каналам связи мгновенно и без искажений.
Данное предположение вытекает из сути третьей составляющей механизма обмена сигналами, допускающей рассмотрение реального канала связи как совокупности самостаятельного элемента системы, функционирование которой сводится к соответствующим задержкам и искажениям сигналов, и фиктивного (идеализированного) канала, передающего сигналы мгновенно и без искажений.
Необходимо заметить, что аналогичные предположения лежат в основе математической модели сопряжения элементов в сложной системе, разработанной Н.П. Бусленко [3], широкое практическое применение которой свидетельствует об ее адекватности исследуемым системам.
Математическая модель сопряжения любых двух элементов 8?е8? и 8/-1е8“?-1 представляет собой пару множеств К??, С??-1 и однозначный оператор Я?, реализующий отображение
сТ = Щ (К), і
• Т?-1 : Ч ,
:АТ-
Оператор сопряжения Я? данному входному контакту кк?єК? элемента 8? страты 8? ставит в соответствие единственный выходной контакт с/-1^ С?-1 элемента С?-1 страты 8?-1, связанный с ним элементарным каналом связи. Следовательно в область определения оператора Я? входят лишь те входные контакты элемента 8?, которые связаны с выходными контактами элемента С?-1 элементарными каналами связи.
Оператор Я? задается в виде бинарной матрицы Л? размера (? хп,?), столбцы и строки которой соответствуют входным и выходным контактам элементов 8? и С??-1.
Структура связей между каждым элементом 8? и стратой 8?-1 описывается оператором сопряжения Я?:
сТ = Щ(К), і є I?-1, і
■ Г?-1
■Ч ,
N,
реализующим отображение (инъекцию) К^С?-', ставящим в соответствие каждому входному контакту кп?еК? элемента 8? страты 8? единственный выходной контакт с/-1^С?-1страты 8?-1, соединенный с контактом к/ элементарным каналом связи.
В отличии от оператора Я??, реализующий функцию К?^ С??-1 оператор Я? всюду определен на множестве контактов элемента 8?. Область значений оператора Я? представляет собой совокупность всех тех выходных контактов страты 8?-1, которые соединены с входными контактами элемента 8?.
Оператор Я?
ЛМЛц?,Лй?,...,Л%\т размера
‘х «г
задается бинарной матрицей
Г л
I;-1
Для описания структуры связей между смежными стратами 8? и 8?-1 введен оператор сопряжения Я?:
с?-1 = Щ?(к?), і ЄI?, і є I?-1, и
Ч-
:Ц? ,
реализующий отображение (сюрьекцию) К?-^С?-1.
В силу предположения 5 область значений оператора Я? совпадает с множеством выходных контактов С?-1 страты 8?-1. Оператор Я? задается матри-
цей Л?=\Л1?, Л2?,..-, Л/\ размера
Я:
ІєІ ?-1
?-1
І
хЕп?
іЄІ ?
Необходимо отметить, что, несмотря на громоздкость матрицы А?, она обладает тем премуще-ством, по сравнению с более компактной таблицей сопряжения [3, 4], что позволяет применять методы теории матриц и теории графов при изучении структуры сопряжения в СИСС, а также существенно упрощает процедуру имитации сопряжения смежных страт.
Контактная модель сопряжения элементов в СИСС
Формально сопряжение любых смежных страт 8?-1и 8?, т. е. формирование входного сигнала %?({) страты 8? на основе выходного сигнала у?-1(0 страты 8?-1с использованием модели сопряжения (К?,С?,Я?), можно описать с помощью фиктивной страты 8? следующим образом.
Согласно принятым предположениям 2 и 5, входной сигнал любой страты формируется исключительно из выходных элементарных сигналов страты 8?-1. Следовательно, в силу предположения
6, каждому выходному сигналу у?-1(/)еРастраты 8?-1 можно ставить в соответствие единственный,
Рис. 1. Структура связей между элементами страт Б? 1 иБ?
Рис. 2. Сопряжение страт Б* иБ* с помощью фиктивной страты Б*
вполне определенный входной сигнал Х^ОеХ* страты 8*, и наоборот. Тогда, с формальной точки зрения можно считать, что между любыми смежными стратами 8“-1и 8* расположена некоторая фиктивная страта 8*, преобразующая выходные сигналы страты 8*-1 в соответствующие им входные сигналы страты 8* в соответствие^ матрицей А*.
Введение фиктивных страт 8*, *е/={1,2,...,к} позволяет из рассмотренных моделей сопряжения элементов и страт в СИСС перейти к эквивалентной, более удобной в теоретическом и прикладном отношениях контактной модели сопряжения, которая предполагает передачу и прием сигналов страт непосредственно через их контакты, без элементарных каналов связи.
В рамках механизма обмена сигналами каждая фиктивная страта 8* описывается множествами входных контактов К*= С*-1 и выходных контактов 8*=К*.
^ Функционирование каждой фиктивной страты 8*, как временной статической системы без памяти [5], детерминировано и описывается оператором V*:
х1*(г) - V * ( у*1 (I)),
реализующим однозначное отображение У1-1^Х*.
С целью выяснения сути оператора V* и выявления его аналитического вида рассмотрим следующий пример.
Пусть структура связей между элементами страт 8* и 8*-1 имеет конфигурацию, изображенную на рис. 1. Путем введения фиктивной страты 8* рассмотренную структуру предствавим в виде, показанном на рис. 2.
В силу предположения 3, каждый входной контакт к^еК может принимать только входной элементарный сигнал х,/(0 еХ* страты 8*, а каждый выходной контакт с/-1 е С*-1 может выдавать только выходной элементарный сигнал у/-1(0 е у*-1 страты 8*-1. Следовательно, данный входной элементарный сигнал х,/(0 представляет собой некоторый
выходной элементарный сигнал у/-1(0, если контакт к/ соединен с контактом с/-1 элементарным каналом, т. е. на пересечении соответствующего столбца и строки матрицы А* стоит единица. Учитывая это, с помощью матрицы А* для рассмотренной структуры сопряжения получим (/упущено):
х* = у*-1; х* = у*-1; х*з = у*-1; х*4 - у*-1;
х2,1 — у1,1 ; х2,2 — ^2,1 ; ^з — у2,2 ; х2,4 — У,з ;
х* = у*-1- х* - у*-1
л2,5 У3,1 > л2,6 У3,2 '
Нетрудно показать, что входной сигнал х*(0 страты 8* представляет собой результат операции умножения выходного сигнала (вектор-строки) страты 8*-1на матрицу А*
(х*, х*л, х*з, х*4, х* 1, х*2, х*,з, х*4, х*5, х*б) -
Действительно, в результате умножения вектор-строки у*-1на матрицу А* получается вектор-строка х*, каждый элемент хы* которой равен сумме произведений выходных элементарных сигналов у/-1 страты 8*-1 на соответствующие элементы столбца (¡,н) матрицы А*.
В силу предположения 4, каждый столбец А* содержит только один единичный элемент, расположенный на пересечении данного столбца (¡,н) и той строки (/,/), соответствующий контакт с/-1 которой соединен с контактом к* элементарным каналом связи. Следовательно, каждый элемент х,/ вектор-строки х* равен тому элементу у*-1 вектор-строки у*-1, который при вычислении хы* умножается на единичный элемент столбца (¡,н) (остальные элементы вектор-строки у -1 умножаются на нуль).
Таким образом, установлено, что оператор V*, ставящий в соответствие каждому выходному сигналу у*-1(0 е У*-1 страты 8*-1 входной сигнал х*(0 еХ*
ХМ = Л’М~1 1,4 У 2,2
▲
I
______I_________
Vй 8 “ ] К Л1,1 1,2 Л1,3 1,4 1 2,1 К 2,2 км 2 3 Ка к2,5 км 2,6
сі,і 1 1 1 1 1
сц~1 1,2 1 1 1
с^[ 2,1 1 1 1 1 1
СМ^ 2,2 ► -'1 1
с^[ С2,3 1
С3,1 1
С3,2 1
Рис. 3. Матрица А? сопряжения страт Э? 1 иЭ?
страты Б», реализует операцию умножения вектор-строки у»-1(0 на матрицу Л»:
х »(*) = у»-1 (г) • Л» =
= ОТЧО, уПО,..., О)) ■ а? =
( >'^1(?УГЧО, УД 0,..., ^
Г-1
Р?-іЛ
(0
(1)
При этом, в силу принятых предположений о структуре сети каналов связи в СИСС, вектор-строка у»-1(0 и матрица Л» всегда согласуются
^1х I ¡Г
к .Ш»-1
по форме (у?-1(0 имеет размер
а ма-
трица А? -
Л
Ґ
І ? * I «г).
«
Уіє!» у
Одно из привлекательных с практической точки зрения свойств фиктивных страт заключается в том, что каждый оператор V? можно представить в виде совокупности независимых друг от друга операторов V?, ієІ?, каждый из которых реализует отображение х?= V?(y?-1(t)).
Для определения вида каждого оператора V? представим матрицу А? через блоки V., ієІ? и подставим в (1):
х?(1) = у?-1(1)[ А?, А., —, А?? ] =
= [ у?-1(1) ■ А?' уГ-\1) ■ А?'-' у?-1 (I) ■ А?? ]. (2)
С другой стороны, в силу предположения 2, имеем:
Х?(*) = (Х?(*),X?(*),...,X? а)).
(3)
При сопоставлении соотношений (2) и (3) становится очевидным, что каждый оператор V» реализует операцию умножения выходного сигнала (вектор-строки) у^1({) страты ¿»^насоответствующую матрицу Л»:
х»(1) = у»\1) • л».
То обстоятельство, что каждый оператор V» допускает декомпозицию на независимые друг от друга операторы V», ^!», предоставляет возможность распараллеливать процесс имитации передачи сигналов от любой страты $“-1 к страте ¿».
Таким образом, процедура имитации сопряжения элементов любых смежных страт ¿» и $“-1 в СИСС сводится к реализации стандартной операции умножения матриц, допускающей распараллеливание процесса имитации сопряжения.
Выводы
Разработана имитационная модель сопряжения элементов в СИСС, позволяющая свести процедуру имитации сопряжения элементов ее смежных страт, т. е. формирование входных сигналов элементов »-й страты на основе выходных сигналов элементов (»-1)-й страты, к реализации стандартной операции умножения матриц, допускающей распараллеливание процесса имитации сопряжения. Применение разработанной модели дает возможность при проведении компьютерных экспериментов ввести необходимые изменения в матрицы сопряжения смежных страт в соответствии с конкурирующими вариантами структур исследуемой системы в ходе экспериментов, без изменения базовой имитационной программы, а также существенно ускорить моделирование при использовании многопроцессорных компьютеров.
Предложенная модель сопряжения использована при компьютерном моделировании процесса функционирования технологической системы измельчения руды Зангезурского медно-молибденового комбината (Армения) с учетом надежности оборудования. Результаты проведенных имитационных экспериментов свидетельствуют об адекватности разработанной модели сопряжения исследуемой технологической системе.
СПИСОК ЛИТЕРАТУРЫ
1. Баласанян С.Ш. Стратифицированная модель для оценки и анализа эффективности функционирования сложных технологических систем со многими состояниями // Известия Томского политехнического университета. - 2011. - Т. 318. - № 5.
- С. 25-30.
2. Баласанян С.Ш. Метод стратифицированной формализации сложных систем с учетом надежности // Вестник ГИУА. Серия «Моделирование, оптимизация, управление». - 2007. - Т. 1. -Вып. 10. - С. 22-32.
3. Бусленко Н.П. Моделирование сложных систем. 2-е изд. - М.: Наука, 1978. - 399 с.
4. Сирота А. Компьютерное моделирование и оценка эффективности сложных систем. - М.: Техносфера, 2006. - 280 с.
5. Месарович М., Такахара И. Общая теория систем: математические основы. - М.: Мир, 1978. - 311 с.
Поступила 19.06.2012 г.
УДК 519.872
МНОГОМОДАЛЬНОСТЬ ВЫСОКОИНТЕНСИВНОГО ПОЛУМАРКОВСКОГО ПОТОКА В УСЛОВИИ ПРЕДЕЛЬНО РЕДКИХ ИЗМЕНЕНИЙ ЕГО СОСТОЯНИЙ
А.А. Назаров1, В.З. Ямпольский2, М.А. Яценко3
Томский государственный университет 2Томский политехнический университет 3ООО«ИНКОМ»
E-mail: yvz@tpu.ru
Построена математическая модель телекоммуникационного потока в виде высокоинтенсивного полумарковского потока в условии предельно редких изменений его состояний, который обладает свойством многомодальности распределения вероятностей числа его событий, наступивших за единицу времени. Показана линейная зависимость его характеристик от времени наблюдения за потоком.
Ключевые слова:
Телекоммуникационные потоки, высокоинтенсивные полумарковские потоки, предельно редкие изменения состояний потока, многомодальность распределения.
Key words:
Telecommunication streams, high-intensity semi-Markov streams, extremely rare stream state changes, multimodal distributions.
Введение
Статистические исследования реальных телекоммуникационных потоков показывают, что некоторые из них обладают свойством многомодаль-ности гистограмм числа их событий, наступивших за единицу времени. Применение методов проверки статистических гипотез не отклоняют гипотезу о том, что теоретическое распределение генеральной совокупности является смесью нормальных распределений.
Ставится задача построения математической модели таких потоков, обладающих свойством многомодальности распределений числа событий, наступивших в этих потоках за единицу времени. В работе показано, что этим свойством обладает высокоинтенсивный полумарковский поток в условии предельно редких изменений его состояний.
Высокоинтенсивный полумарковский поток
в условии предельно редких изменений его состояний
Пусть задана диагональная матрица А(х) размерности L с элементами Д(х), 1=1,, и стохастическая матрица
Р = I + 50 (1)
размерности ,х, с элементами
р1112 = 5у2 + 5<?у2>
где I - единичная матрица; 5^ - символ Кронеке-ра; 5 - некоторый малый параметр. Также задан положительный большой параметр N.
Рассмотрим полумарковский процесс 1(7), принимающий значения 1(7)=1=1,, с условно независимыми компонентами [1], полумарковская матрица которого факторизуется и равна произведению
А(ЛХ)Р =А(ЛХ)(1+50). (2)
Обрабатывающий его процесс п(7) определяет число событий, наступивших за время 7 в полумар-ковском потоке (8М-потоке), заданным полумар-ковской матрицей (2).
Полумарковский процесс 1(7) будем называть процессом, управляющим полумарковским потоком, а значения 1(7)=1 этого процесса будем называть состояниями 8М-потока.
В состоянии I длина ^ интервала между моментами наступления событий рассматриваемого потока определяется функцией распределения