Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Лаба 3 (Цепи Маркова и СМО) / Красивый цикл

.txt
Скачиваний:
23
Добавлен:
25.12.2019
Размер:
620 б
Скачать
N=[200 1000 10000];
z=zeros(N(3),1,3);
r=rand(N(3),1,3);

for i=1:3
z(1,1,i)=1;
for n=2:N(i)
k=1;
while r(n-1,1,i)>=P_cum(z(n-1,1,i),k)
k=k+1;
end
z(n,1,i)=k;
end
%Построение графиков
figure(2)
subplot(3,1,i)
plot(z(1:N(i),1,i))
txt=['Поведение цепи Маркова в течение ',num2str(N(i)),' итераций'];
title(txt)
xlabel('t')
ylabel('States')
end


for n=2:N
k=1;
while r(n-1,1) >= P_cum(a(n-1,1),k)
k=k+1;
end
a(n,1)=k;
end
figure();
stairs(a');
Соседние файлы в папке Лаба 3 (Цепи Маркова и СМО)