0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 10
1

Стационарное распределение марковской цепи считается неправильно, почему-то

05.06.2016, 00:50. Показов 869. Ответов 3
Метки нет (Все метки)

Моделирую марковский процесс. Конкретно Матрицу переходных вероятностей. Но вот неудача, при расчете стационарного распределения(пи) получается ерунда!

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
N=10^5; %количество шагов(тактов)
MPV=zeros(4,4);
s1=0;
s2=0;
s3=0;
s4=0;
for j=1:4
    for i=2:N
        X=rand;
        if X<Ms(j,2)            s1=s1+1; % Ms-матрица кумулятивных вероятностей
        elseif    X<=Ms(j,3)   s2=s2+1;
        elseif    X<=Ms(j,4)   s3=s3+1;
        else    s4=s4+1;
        end
    end
    MPV(j,:)=[s1 s2 s3 s4];
end
PII=MPV/N;

вот результат Этого:
0.6194 0.1223 0.0307 0.2276
0.6255 0.8052 0.1507 0.4185
0.7785 1.0299 0.1578 1.0338
0.7806 1.4434 0.5710 1.2050

Но Пи не может быть больше 1!!!!!!
Кто-нибудь может понять и помочь?))))
Спасибо)
0

Помощь в написании контрольных, курсовых и дипломных работ здесь.

Programming
Эксперт
94731 / 64177 / 26122
Регистрация: 12.04.2006
Сообщений: 116,782
05.06.2016, 00:50
Ответы с готовыми решениями:

Считается неправильно, почему?
Подскажите пожалуйста! я перемножаю числа которые находятся в текст боксах, для этого написал...

Почему неправильно считается результат?
Помогите,пожалуйста:создал два вычисляемых поля Ставка и Сумма в таблице &quot;Платежи&quot;,однако Access их...

Почему неправильно считается формула?
Вопрос в заголовке.Программа написана вроде правильно, ошибок не вижу. P.S. Так как переменные...

Почему-то неправильно считается кол-во элементов массива
В последнем цикле считывается текст из файла и вместе с этим подсчитывается кол-во символов. Текст...

3
Модератор
1274 / 1189 / 357
Регистрация: 13.09.2015
Сообщений: 4,137
05.06.2016, 05:28 2
XYRMA, посмотрите внимательней, не стоит ли перед выведенной матрицей строкой выше что-нибудь вроде "1е-3*"?
0
0 / 0 / 0
Регистрация: 01.06.2015
Сообщений: 10
05.06.2016, 12:33  [ТС] 3
Нет.Ничего не стоит.

Matlab M
1
2
3
4
5
>> Spec1
    0.6194    0.1223    0.0307    0.2276
    0.6255    0.8052    0.1507    0.4185
    0.7785    1.0299    0.1578    1.0338
    0.7806    1.4434    0.5710    1.2050
Создание матрицы MS
Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
n=11;
rand('seed',132*n);
 
N=4;%количество состояний
%1-строка, 2-столбец
for i=1:N
    M=rand(1,3);
    M=sort(M);
    Mas=[0, M, 1];
    for j=1:N
        S(i,j)=Mas(j+1)-Mas(j);
    end
    Ms(i,:)=Mas;
end
0
Модератор
1274 / 1189 / 357
Регистрация: 13.09.2015
Сообщений: 4,137
06.06.2016, 12:01 4
XYRMA, скажите мне, зачем вам цикл
Matlab M
1
2
3
for j=1:N
        S(i,j)=Mas(j+1)-Mas(j);
end
если значения S нигде не используются?
Почему не получается у вас, я не стал сильно разбираться, написал свой код, где всё как надо получается:
Matlab M
1
2
3
4
5
6
7
8
9
10
n=11;
rand('seed',132*n);
N=4;
M=sort(rand(N-1,N));
K=10^5;
MS=repmat(M,1,1,K);
X=repmat(rand(1,N,K),N-1,1,1);
Y=[X<MS;ones(1,N,K)];
MPV=sum(cumsum(Y)==1,3);
PII=(MPV/K)'
0
IT_Exp
Эксперт
87844 / 49110 / 22898
Регистрация: 17.06.2006
Сообщений: 92,604
06.06.2016, 12:01

Почему неправильно считается простое математическое выражение?
Всем привет. Вопрос будет очень банальным: я правильно записал выражение? double R = (2 *...

Стационарное распределение температуры в брусе
Найти стационарное распределение температуры в брусе прямоугольного сечения (0&lt;=x&lt;=a;0&lt;=y&lt;=b). К...

Почему в однокритериальной оптимизации считается минимум для сопротивления, но не считается для тока?
В самом конце работы идёт однокритериальная оптимизация по мощности P2 для R2 и I. Надо сделать с...

Неправильно считается факториал
Добрый день,написал рекрсивную функцию. function fact(n:integer):integer; begin if n &lt; 1...

Неправильно считается функция if
MathCAD неправильно считает функцию if

Неправильно считается факториал
Добрый день,написал рекрсивную функцию. function fact(n:integer):integer; begin if n &lt; 1...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.