Дискретно-стохастические модели или P- cхемы.

Р-схемы принято называть вероятностными автоматами, с помощью которых можно определить дискретные преобразования информации с памятью. Функция в каждом такте зависит только от состояния памяти, поэтому такие автоматы могут быть описаны статистическими методами. Для вероятностных автоматов используется метод, основанный на F-автоматах. Пусть задан P-автомат в виде матрицы.

clip_image002

clip_image004

Рассмотрим граф для этого Р - автомата.

clip_image006

В данном примере рассмотрен Y - детерминированный Р - автомат. Такое представление эквивалентно заданию дискретной марковской цепи с конечным множеством состояний. При исследовании вычислительных систем часто приходится сталкиваться с дискретными случайными процессами, которые определены на конечном множестве состояний. Для анализа F - автоматов широко используются дискретные марковские цепи, которые характеризуются:

1. Множеством состояний S Î {S1, S2, S3, ... , Sn}.

2. Матрицей вероятностных переходов или переходных вероятностей.

clip_image008

3. Вектор начальных распределений.

clip_image010-этот вектор определяет вероятности рi0, соответствующие моменту времени t0 при нахождении процесса в состоянии Si.

Марковская цепь изображается в виде графа, вершины при этом соответствуют состояниям, а дуги - переходам.

clip_image012

Марковская цепь порождает множество реализаций вероятностных процессов f(t), которые можно представить в виде некоторой последовательности (множества) {f1, f2, ...}. Каждый из таких процессов происходит в соответствующие моменты времени. Начальное состояние f0 определяется вектором начальных вероятностей p, а следующие состояния f1, f2, ... определяются в соответствии с матрицей переходов Р. Таким образом, состояние f1 будет определяться как Sj, которое определяется из i-й строки матрицы Р. Процесс f(t) переходит в состояние f1 равное Sj с вероятностью pi j. Затем процесс переходит в состояние f2 = Sk c вероятностью pik. В результате N таких шагов процесс попадет в состояние SN с вероятностью piN.

Марковские цепи классифицируются в зависимости от возможности перехода из одного состояния в другое. Различают поглощающие и эргодические марковские цепи. Поглощающая марковская цепь содержит состояние, достигнув которого, процесс прекращается и это состояние в марковских цепях принято обозначать как S0. Вероятность перехода p00 = 1. Матрица вероятностей имеет следующий вид:

clip_image014

В этом смысле марковская поглощающая цепь обладает таким свойством: из какого бы состояния не начинался процесс, при N®¥ c вероятностью равной 1 этот процесс окажется в поглощающем состоянии S0. Основная характеристика случайного процесса в марковских цепях - число пребываний процесса в состояниях S1, ..., Sk до поглощающего состояния. Нахождение в промежуточных состояниях обычно минимизируют. Число пребываний в каждом из состояний Si в марковских цепях является случайной величиной, которая может быть представлена в виде средних значений (математических ожиданий), дисперсий, функций распределения. Поглощающие марковские цепи широко используются в качестве временных моделей программ и вычислительных процессов. При моделировании программ состояния цепи отождествляются с блоками программ (модулями), а матрица вероятностей переходов определяет порядок переходов между блоками, т.е. состояниями. В результате моделирования можно вычислить число обращений к блокам программы и время выполнения программы, которые оцениваются средними значениями, дисперсиями, а при необходимости - функциями распределения.

Эргодическая марковская цепь также представляет собой множество состояний, которые связаны матрицей переходных вероятностей, но отличается тем, что из какого бы состояния процесс не исходил, после некоторого числа шагов он может оказаться в любом состоянии. Другими словами, он может из любого состояния перейти в некоторое другое за n шагов. По этой причине состояния называют эргодическими, т.е. возвратными. Процесс, порождаемый эргодической цепью, никогда не завершается, а последовательно переходит из одного состояния в другое. При этом переход из одного состояния в другое может осуществляться с разной частотой, что определяется переходными вероятностями, поэтому для эргодической цепи - основная характеристика - вероятность пребывания в некотором состоянии Sj, где j=1,2, ..., k, - это есть относительная частота попадания процесса в состояние Sj и одновременно доля времени, которое процесс проводит в каждом из таких состояний эргодической цепи. Эргодические цепи также широко используются для моделирования ВС, в частности, при исследовании вопросов надежности, а также в качестве базовых моделей взаимодействия устройств с задачами, поступающими на обработку. Марковский процесс с дискретными состояниями S1, ..., Sk, переходы между которыми разрешаются в любой момент времени, принято называть непрерывными марковскими цепями. Однородную непрерывную марковскую цепь, поведение которой в любой момент времени подчиняется одному и тому же закону, принято задавать матрицей интенсивностей

переходов Q = [qi j]. Интенсивности переходов в общем случае определяются следующими формулами:

clip_image016

где pii (Dt) - вероятность перехода процесса из состояния Si за время Dt в состояние Si. Это означает, что если процесс находился в состоянии Si, то интенсивность перехода в состояние Sj в течение промежутка времени Dt отличается от интенсивности перехода в состояние Si на величину -qii Dt. Аналогичным образом, для второго случая интенсивность перехода будет определяться как qi j Dt. В качестве примера рассмотрим граф марковской цепи, состоящий из трех состояний.

clip_image018

clip_image020

clip_image022clip_image024

Основная характеристика непрерывной марковской цепи - это стационарное распределение вероятностей состояний: a = (a1, ..., ak), где a1, ..., ak - вероятности пребывания процесса в состояниях S1, ... , Sk. Распределение задается вероятностным решением системы уравнений:

clip_image026

К системе, состоящей из (k-1) уравнений, добавляется нормирующее уравнение: a1 + a2 + ... + ak = 1 (3), которое, в свою очередь определяется значениями вероятностей a = {a1, a2, ..., ak} . Уравнения (1), (2), (3) принято называть уравнениями равновесия с помощью которых определяются условия равновесия входных и выходных параметров. В частности, для нашего примера таковыми уравнениями равновесия будет являться следующая система:

clip_image028

Cостояние

Интенсивность входящего потока

Интенсивность выходного потока

S1

ma2

(l/2 + l) a1

S2

la1 + ma3

(m + l)a2 

S3

l/2a1 + la2

ma3

Этот подход широко используется при исследовании СМО. В соответствии с марковскими свойствами все предыдущие процессы сказываются на поведении в будущем только через текущее состояние, т.е. определяя дальнейший ход событий не требуется знать, как долго процесс находился в текущем состоянии. Отсюда следует вывод, что распределение оставшегося времени пребывания процесса в некотором состоянии Si, должно зависеть только от самого состояния, а не от времени пребывания в нем. Этим свойством обладает только одно распределение - экспоненциальное, поэтому уравнения равновесия можно применять для СМО с экспоненциальным распределением, однако, если допустить произвольное распределение времени пребывания процесса в том или ином состоянии, то этот процесс становится полумарковским и ведет себя только относительно моментов изменения состояний, т.е. как обычная дискретная марковская цепь и принято говорить, что в эти моменты времени имеют место вложенные марковские цепи. Характеристики полумарковских цепей и процессов определяются значительно сложнее, чем марковских, поскольку вероятности состояний связаны с параметрами процесса системами дифференциальных уравнений с частными производными.

Предлагаю ознакомиться с аналогичными статьями: