Рассмотрите эту теоретическую, правильно-стохастическую матрицу перехода стохастического процесса.

Элемент $$P_{ij}$$ вероятность, что переходы процесса, чтобы утвердить j во время t + 1, учитывая, что это находится в состоянии i во время t для всего t.

Создайте Цепь Маркова, которая характеризуется матрицей P перехода.

P = [0.5 0.5 0 0; 0.5 0 0.5 0; 0 0 0 1; 0 0 1 0];
mc = dtmc(P);

mc является объектом dtmc, который представляет Цепь Маркова.

Отобразите количество состояний в Цепи Маркова.

numstates = mc.NumStates

numstates = 4

Постройте ориентированного графа Цепи Маркова.

figure;
graphplot(mc);

Заметьте, что утверждает 3, и 4 формируют абсорбирующий класс, в то время как состояния 1 и 2 являются переходными.

Рассмотрите эту матрицу перехода в который элемент $$(i,j)$$ наблюдаемое число раз, утверждают i переходов, чтобы утвердить j.

Например, $$P_{\mathrm{}32}=7$$ подразумевает что состояние 3 перехода, чтобы утвердить 2 семь раз.

P = [16 2  3  13;
     5  11 10 8;
     9  7  6  12;
     4  14 15 1];

Создайте Цепь Маркова, которая характеризуется матрицей P перехода.

mc = dtmc(P);

Отобразите нормированную матрицу перехода, сохраненную в мГц. Проверьте, что элементы в строках суммируют к 1 для всех строк.

mc.P

ans = 4×4

    0.4706    0.0588    0.0882    0.3824
    0.1471    0.3235    0.2941    0.2353
    0.2647    0.2059    0.1765    0.3529
    0.1176    0.4118    0.4412    0.0294

sum(mc.P,2)

ans = 4×1

     1
     1
     1
     1

Постройте ориентированного графа Цепи Маркова.

figure;
graphplot(mc);

Считайте деловой цикл с двумя состояниями США действительным валовым национальным продуктом (ВНП) в [3] p. 697. Во время t, действительный GNP может быть в состоянии расширения или сокращения. Предположим, что следующие операторы верны в демонстрационный период.

Создайте матрицу перехода для модели.

p11 = 0.90;
p22 = 0.75;

P = [p11 (1 - p11); (1 - p22) p22];

Создайте Цепь Маркова, которая характеризуется матрицей P перехода. Маркируйте два состояния.

mc = dtmc(P,'StateNames',["Expansion" "Contraction"])

mc = 
  dtmc with properties:

             P: [2x2 double]
    StateNames: ["Expansion"    "Contraction"]
     NumStates: 2

Постройте ориентированного графа Цепи Маркова. Укажите на вероятность перехода при помощи цветов обводки.

figure;
graphplot(mc,'ColorEdges',true);

Чтобы помочь вам исследовать функции объекта dtmc, mcmix создает Цепь Маркова из случайной матрицы перехода использование только конкретного количества состояний.

Создайте Цепь Маркова с пятью состояниями из случайной матрицы перехода.

rng(1); % For reproducibility
mc = mcmix(5)

mc = 
  dtmc with properties:

             P: [5x5 double]
    StateNames: ["1"    "2"    "3"    "4"    "5"]
     NumStates: 5

mc является объектом dtmc.

Постройте собственные значения матрицы перехода на комплексной плоскости.

figure;
eigplot(mc)

Этот спектр определяет структурные свойства Цепи Маркова, такие как периодичность и смешивание уровня.