- •2. Основные понятия математического моделирования
- •2. Принципы построения математических моделей
- •3.Классификационные признаки и классификация моделей
- •1.Основные этапы математического моделирования
- •2. Понятие о вычислительном эксперименте
- •3.Оценка свойств моделей
- •2.1. Метод нелинейного преобразования, обратного функции распределения
- •2.2. Метод суперпозиции
- •2.3. Некоторые специальные методы моделирования случайных величин
- •2.3.1. Метод Неймана
- •3.2. Метод кусочной аппроксимации
- •2.4. Моделирование векторных случайных величин
- •Раздел 2. Модели и методы исследования структурных свойств сетей связи
- •2. Изоморфность графов
- •1. Подграфы и дополнения
- •2. Деревья, разрезы, циклы
- •1. Матрица циклов и ее связь с матрицей инцидентности.
- •2.Матрица разрезов и ее связь с матрицей циклов
- •Матрицы смежности и перечисления путей
- •2. Связность
- •1. Степенная последовательность вершин графа
- •2. Алгоритм синтеза графов с максимальной связностью
- •3. Алгоритмы синтеза графов с максимальной связностью при заданном числе вершин и ребер.
- •4.Однородные графы. Мера неоднородности.
- •Модели и методы оценки разведзащищенности сетей связи
- •1.Показатели разведзащищенности сетей связи.
- •2.Модель оценки разведзащищенности сети связи.
- •Матрица смежности графа а
- •Матрица смежности графа в
- •Зависимость разведзащищённости от степени неоднородности сетей
- •5. Синтез сетей связи, оптимальных по показателям структурной устойчивости и разведзащищенности.
- •3.1. Моделирование марковских случайных процессов
- •3.2. Разностные и дифференциальные стохастические уравнения
- •2. Понятие эквивалентной функции стохастической сети
- •1.Метод двухмоментной аппроксимации.
- •3.Разложение на простые дроби.
- •3. Общие правила моделирования исследуемых процессов
1.Метод двухмоментной аппроксимации.
Особенностями этого метода является:
Придание вероятностного смысла эквивалентной функции стохастической сети: Q(s) – вероятность того, что процесс передачи сообщений по системе связи будет завершен успешно до случайного момента времени старения этих сообщений.
Принятие ограничений на вид функции распределения времени старения передаваемой информацииZ(t).
В общем случае, вероятность того, что случайное время передачи tП, распределенное по закону H(t), будет меньше случайного времени старения информации tc, распределенное по закону Z(t) равна
.
(23)
При использовании преобразования Лапласа - Стилтьеса
Z(t)= 1 – exp[-t/ tc],
а параметр преобразования s = 1/tcимеет смысл интенсивности старения передаваемой информации. Отсюда, вероятность своевременной передачи сообщений

и является, по существу, характеристической функцией, позволяющей определить k-е начальные моменты случайного времени реализации стохастической сети по формуле:
.
(24)
Отсюда среднее время передачи сообщения, определяемое как начальный момент первого порядка равно:

а дисперсия времени передачи D[tП], определяемая как второй центральный момент, вычисляется по формуле:
.
Вычисление математического ожидания tПи дисперсии D[tП] позволяет приближенно определить функцию распределения времени передачи как неполную гамма – функцию
(25)
где
,
- параметры формы и масштаба, соответственно.
Таким образом, реальная функция распределения времени передачи была аппроксимирована неполной гамма - функцией, с параметрами, определяемыми из эквивалентной функции стохастической сети, что и дало название этому методу. Следует заметить, что если исследователю априорно известен закон распределения, наиболее полно характеризующий исследуемый процесс, то аналогичным образом можно выполнить аппроксимацию этим законом реальной функции распределения.
Описанный метод позволяет получать достаточно точные результаты при моделировании систем массового обслуживания и обуславливает его широкое использование специалистами для моделирования систем связи общего назначения. Однако при моделировании систем связи ВМФ, из-за несоответствия закона старения передаваемой информации реальному, применение метода двухмоментной аппроксимации может привести к существенным ошибкам. Поэтому для получения более точных результатов моделирования и при решении исследовательских задач необходимо использовать второй подход.
2.Метод обращения.
В отличие от первого, второй подход к определению вероятностно-временных характеристик систем связи не предусматривает ограничений на вид закона распределения времени старения передаваемой информации. Поэтому, эквивалентная функция стохастической сети рассматривается не как вероятность своевременной передачи, а как изображение исследуемого процесса доставки сообщений по Лапласу. Действительно, функция
, (26)
является классическим представлением интегрального преобразования Лапласа функции–оригинала f(tП) - функции плотности вероятностей времени tП. В свою очередь, определение изображений подпроцессов исследуемого процесса передачи и сам алгоритм формирования эквивалентной функции является типичным для операционного исчисления методом, широко применяемым для исследования сложных функций. Что касается параметра s, то в данном случае он представляет собой не интенсивность старения передаваемой информации, а является комплексной переменной:
s = + jy, - < y < +.
Таким образом, интеграл (26) является преобразованием Лапласа, ставящим в соответствие каждой функции (оригиналу) f(tП), tП> 0, единственную функцию Q(s) (изображение) комплексной переменной, причем для существования изображений необходима сходимость этого интеграла. Отсюда вытекает первый важный шаг анализа – проверка сходимости интеграла Лапласа. Напомним, что признаком абсолютной и равномерной сходимости (26) является существование предела

и аналитичность функции Q(s), приRe(s) = 0 Точная нижняя грань a действительных чисел 0,для которых выполняется это условие называется абсциссой абсолютной сходимости преобразования Лапласа. При этом в полуплоскости Re(s)0 изображение Q(s), определяемое из стохастической сети с использованием уравнения Мейсона, как было доказано выше, является голоморфной (однозначной аналитической) функцией. Областью определения функции Q(s) является полнаяs-плоскость, за исключением особых точек, расположенных слева от прямой =a.
После проверки сходимости интеграла Лапласа, следует второй шаг, состоящий в переходе из пространства изображений в пространство оригиналов, т.е. осуществление обратного преобразования Лапласа, основанного на теореме обращения (см. Преобразования Лапласа).
Достаточными условиями существования обратного преобразования являются:
Если Q(s) аналитична для a и имеет порядок меньше –1, то обратное преобразование Лапласа L-1[Q(s)] существует; оно непрерывно для всех t>0 и f(t) асимптотически сходится к функции exp[at] при t с абсцисой абсолютной сходимости a, а при t 0 L-1[Q(s)] 0;
Если Q(s)=[Qi(s), i=1,n], при (si, i=1,n)аналитичной относительно каждого si, равна нулю при si=0, i=1,n и Qi(s) = L[fi(t)] (ai,i=1,n), то L-1[Q(s)] существует и соответствующее преобразование Лапласа имеет абсциссу абсолютной сходимости.
Задача определения оригинала f(t), соответствующего эквивалентной функции стохастической сети Q(s) является достаточно сложной и для ее решения используются способ контурного интегрирования и методы разложения Q(s) (с целью упрощения ее вида), с последующим почленным переходом в пространство оригиналов.
Способ контурного интегрирования подразумевает полное определение функции Q(s) для всех значений комплексной переменной s=α+jy, -<y<+, а поэтому его использование зачастую весьма трудоемко и на практике применяется достаточно редко.
Значительно упростить процедуру обращения (перехода от Q(s) к f(t)) позволяют методы разложения, причем в зависимости от вида эквивалентной функции используются различные разложения. Рассмотрим их более подробно.
Эквивалентная функция Q(s) является рациональной алгебраической функцией.
Разложение Хевисайда.
В этом, самом распространенном случае, эквивалентная функция представляется в виде отношения двух многочленов f(s) иφ(s), степеней m и n соответственно, т.е.
,
deg[f(s)]=m;
deg[φ(s)]=n;
n>m
где deg[*] – степень соответствующего многочлена.
Условие n>m позволяет использовать теорему о вычетах и при
φ(s)
= a0(s
- s1)(s
- s2)…(s
– si)…(s
- sn),
si
si+1
,
i = 1,n;
представить эквивалентную функцию в виде суммы:
,
где φ’(si) –значение производной многочлена знаменателя в точке si. А значит
.
В случае кратных
корней, т.е. если
φ(s) = a0(s - s1)a1 (s - s2)a2 …(s – si)ai … (s - sn)an , si si+1 , i = 1,n;
оригинал определяется по формуле:
.
Если φ(s) имеет комплексно сопряженные корни s = αjy, то пары членов, соответствующих этим корням можно преобразовать

,
гдеR
= 2A2+B2;
= -arctg[A/B];
= arctg[B/A].
При этом,
если A, B, R,
и
- являются действительными, то и L-1[Q(s)]
= f(t) является
действительной функцией.
