Вероятностный автомат является примером дискретно-стохастической функциональной модели исследуемой системы. Сущность дискретизации времени при этом подходе остается аналогичной рассмотренным ранее конечным автоматам, но появляется возможность учета влияния случайных факторов на развитие моделируемых процессов.
В общем виде вероятностный автомат (англ. probabilistic automat) можно определить как дискретный потактный преобразователь информации с памятью, функционирование которого в каждом такте зависит только от состояния памяти в нем и может быть описано статистически [10].
Применение схем вероятностных автоматов имеет важное значение для моделирования дискретных систем, проявляющих статистически закономерное случайное поведение, для выяснения алгоритмических возможностей таких систем и обоснования границ целесообразности их использования.
Введем математическое определение вероятностного автомата, используя понятия, принятые ранее для определения конечного автомата. Рассмотрим множество G, элементами которого являются всевозможные пары (xi, zk), где xi и zk - элементы входного подмножества Х и подмножества состояний Z соответственно. Если существуют две такие функции j и y, что с их помощью осуществляются отображения G®Z и G®Y, то говорят, что F=<Z, X, Y, j, y, z0 > определяет автомат детерминированного типа.
Введем в рассмотрение более общую математическую схему. Пусть Ф - множество всевозможных пар вида (zk, yi), где yi - элемент выходного подмножества Y. Потребуем, чтобы любой элемент множества G индуцировал на множестве Ф некоторый закон распределения следующего вида:
Элементы из Ф:
(z1, y1)
…
(zk, yj)
…
(zK, yJ);
Элемент (xi,zs) из G:
b11
…
bkj
…
bKJ.
При этом где bkj - вероятности перехода автомата в состояние zk и появления на выходе сигнала yj, если автомат был в состоянии zs и на его вход в этот момент времени поступил сигнал xi. Число таких распределений, представленных в виде таблиц, равно числу элементов множества G. Обозначим множество этих таблиц через В. Тогда четверка элементов P=<Z, X, Y, В> называется вероятностным автоматом.
Пусть элементы множества G индуцируют некоторые законы распределения на подмножествах Y и Z, что можно представить соответственно в следующем виде
Элементы из Y:
y1
…
yi
…
yJ
;
Элемент (xi,zs) из G:
q1
…
qj
…
qJ
;
Элементы из Z:
z1
…
zk
…
zK
;
Элемент (xi,zs) из G:
s1
…
sk
…
sK
.
При этом и , где sk и qj- вероятности перехода автомата в состояние zk и появления выходного сигнала ykпри условии, что автомат находился в состоянии zs и на его вход поступил входной сигнал xi.
Если для всех k и j имеет место соотношение , то такой автомат называется вероятностным автоматом Мили. Это требование означает выполнение условия независимости распределений для нового состояния вероятностного автомата и его выходного сигнала.
Пусть теперь определение выходного сигнала вероятностного автомата зависит лишь от того состояния, в котором находится автомат в данном такте работы. Другими словами, пусть каждый элемент выходного подмножества Y индуцирует распределение вероятностей выходов, имеющее следующий вид:
Элементы из Y:
y1
…
yk
…
yK;
Элемент zk из Z:
q1
…
qi
…
qI .
Здесь , где qi - вероятность появления выходного сигнала yi при условии, что автомат находился в состоянии zk.
Если для всех k и i имеет место соотношение qksi=bki, то такой автомат называется вероятностным автоматом Мура. Понятие вероятностных автоматов Мили и Мура введено по аналогии с детерминированным конечным автоматом, задаваемым F=<Z, X, Y, j, y, z0 >.
Частным случаем вероятностного автомата, задаваемого как P=<Z, X, Y, В>, являются автоматы, у которых либо переход в новое состояние, либо выходной сигнал определяются детерминированно. Если выходной сигнал автомата определяется детерминированно, то такой автомат называется Y-детерминированным вероятностным автоматом. Аналогично Z-детерминированным вероятностным автоматом называется автомат, у которого выбор нового состояния является детерминированным.
Пример 5.4. Рассмотрим Y-детерминированный вероятностный автомат, который задан таблицей переходов (табл. 5.6) и таблицей выходов (табл.5.7).
В этих таблицах pij - вероятность перехода автомата из состояния zi в состояние zj. При этом, как и ранее, .
Первую из этих таблиц можно представить в виде квадратной матрицы размерностью , которую будем называть матрицей переходных вероятностей или просто матрицей переходов вероятностного автомата. В общем случае такая матрица переходов имеет вид
.
Для описания Y-детерминированного автомата необходимо задать начальное распределение вероятностей следующего вида:
Элементы из Z:
z1
…
zk
…
zK .
Элементы из D:
d1
…
dk
…
dK .
Здесь dk - вероятность того, что в начале работы автомат находится в состоянии zk. При этом .
Будем считать, что до начала работы (до нулевого такта времени) автомат всегда находится в состоянии z0 и в нулевой такт времени меняет состояние в соответствии с распределением D. Дальнейшая смена состояний автомата определяется матрицей переходовPp. Информацию о начальном состоянии вероятностного автомата удобно внести в матрицуPp, увеличив ее размерность до (K+1)(K+1). При этом первая строка такой матрицы, сопоставляемая состоянию z0, будет иметь вид (0, d1, d2 ,..., dk), а первый столбец будет нулевым.
Описанный Y-детерминированный автомат можно задать в виде ориентированного графа, вершины которого сопоставляются состояниям автомата, а дуги — возможным переходам из одного состояния в другое. Дуги имеют веса, соответствующие вероятностям перехода pij, а около вершин графа пишутся значения выходных сигналов, индуцируемых этими состояниями.
Очевидно, что с точки зрения математического аппарата задание Y-детерминированного автомата эквивалентно заданию некоторой дискретной марковской цепи с конечным множеством состояний. Поэтому аппарат марковских цепей является основным для аналитических расчетов.
Пример 5.5. Пусть задан Y-детерминированный вероятностный автомат
На рис. 5.3 показан граф переходов этого автомата. Требуется оценить суммарные финальные вероятности пребывания этого автомата в состояниях z2 и z3.
При использовании аналитического подхода можно записать известные соотношения из теории марковских цепей и получить систему уравнений для определения финальных вероятностей. При этом начальное состояние z0 можно не учитывать, так как начальное распределение не оказывает влияния на значения финальных вероятностей.
Тогда имеем
, где ck — финальная вероятность пребывания автомата в состоянии zk. Получаем систему уравнений
Добавим к этим уравнениям условие нормировки с1+ с2+ с3+ +с4=1.Тогда, решая систему уравнений, получим с1=5/23, с2=8/23, с3=5/23, с4=5/23.
Таким образом, с2+ с3=13/23=0,5652. Другими словами, при бесконечной работе заданного в этом примере Y-детерминированного автомата на его выходе формируется двоичная последовательность с вероятностью появления единицы, равной 0,5652.
Подобные автоматы могут использоваться как генераторы марковских последовательностей, которые необходимы при моделировании процессов функционирования систем или воздействий внешней среды.
Для оценки различных характеристик исследуемых систем, представляемых в виде конечных и вероятностных автоматов, кроме рассмотренного случая аналитических моделей можно применять и имитационные модели, реализуемые, например, методом статистического моделирования.