50690

Моделирование потока Пуассона

Лабораторная работа

Информатика, кибернетика и программирование

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

Русский

2014-01-28

158 KB

11 чел.


Лабораторная работа № 4

Моделирование потока Пуассона

Цель работы. Практическое освоение алгоритма программной генерации стационарного потока Пуассона и методов статистической проверки разработанного генератора.

clc

clear all

format compact

format long

W1=1;W2=48; %интервал моделирования

la=1.5;

N=100;

%количество экспериментов

amo=input('Enter quantity of ехperiments: ');

t=ones(1,N);

t(1)=0;

i=2;

delta=la*(W2-W1)

l0=exp(-delta);

l=[1,(ones(1,N)*delta)./(1:N)];

for k=N+1:-1:1

   P(k)=prod(l(1:k))*l0;

   

end;

%стационарный поток

mas_of_n=rand(1,amo);

%---------массив экспериментальных значений n----

for o=1:amo

i=2;  

t(1)=W1;

%----------экспериментальное вычисление n-----------

 while t(i-1)<W2

  z=rand;

  t(i)=t(i-1)-log(z)/la; %рекуррентное вычисление t

  i=i+1;

end

%i-3 - это и есть n

mas_of_n(o)=i-3;

end;

%+++++++++++++++ сам поток ++++++++++++++++++

figure(1)

%---v t realiz puas potoka---

t=t(2:i-3);

h=ones(2,i-3);

plot(t,t,'*');%,h,'*');

title('Puasson''s stream');

xlabel('t - time');

d=rand(1,N+1);

for i=1:N+1

   d(i)=i;

 end

%------------формирование nu-------------

%-----массив экспериментальных количеств--

%---всего amo экспериментов-----

E=hist(mas_of_n,d);

E=E/amo;%---частоты----

%%%----sum(E)=1---sum(P)=1----

%-----E - экспериментальный поток----

%-----P - теоретический поток--------

%-------(частоты)--------------------

%+++++++++++++++ plotnosti ver +++++++++

figure (2)

bar(d,E,1,'b')

hold on

bar(d,P,0.1,'k')

title('Density of distribution. Blue - experimentel, Black - theoretic')

xlabel('n - amount of events')

ylabel('p - probablity')

hold off

%+++++++++++++++ Массив экспериметнов ++++++++++++

figure (3)

plot(mas_of_n,'*r');

title('Array of amounts');

xlabel('Number of experiment');

ylabel('Amount of events');

%+++++++++++++++ exp i teor hist++

figure (4)

mi=min(mas_of_n);

ma=max(mas_of_n);

hist(mas_of_n,(mi:ma));

title('Histograms')

xlabel('Amount of events(theor)');

ylabel('Amount of experiments');

hold on

bar((mi:ma),P((mi+1):(ma+1))*amo,0.6,'r');

hold off

k=N+1;

gr=input('How many karmanov in group:  ');

%new value of karmanov

karm=ceil(k/gr);%----round value-----

zap=zeros(1,karm*gr-k); %nuli

%--------merge in groups---------

%--------sum is the row vector---

E1=sum(reshape([E zap],gr,karm))*amo;%-------resize of v--

P1=sum(reshape([P zap],gr,karm))*amo;%-------resize of p--

%---попробуем просуммировать----

o=1;h=0;

while E1(o)<5

   h=h+E1(o);

   o=o+1;

end;

E1(o)=E1(o)+h;

o2=size(E1);

o1=o2(2);

h=0;

while E1(o1)<5

   h=h+E1(o1);

   o1=o1-1;

end;

E1(o1)=E1(o1)+h;

h=0;

for i=1:o-1,

   h=h+P1(i);

   P1(i)=0;

end;

P1(o)=P1(o)+h;

h=0;

for i=o2(2):-1:o1+1,

  h=h+P1(i);

  P1(i)=0;

end;

P1(o1)=P1(o1)+h;

figure(5)

xc=(gr+1)/2:gr:k+gr/2;%-----centers of bars-------

bar(xc(o:o1),E1(o:o1),1,'g')

title('Histogram on KARM karmanov (Exp-blue, Theor-green)')

xlabel('Intervals')

ylabel('N*freq')

hold on

bar(xc(o:o1),P1(o:o1),0.8,'b')

hold off

PP=ones(1,o1-o+1);

EE=PP;

for i=o:o1,

   PP(i-o+1)=P1(i);

   EE(i-o+1)=E1(i);

end;    

   

hi2=sum(((EE-PP).^2)./PP);

stsv=o1-o;

disp(['hi2 =',num2str(hi2),'         Degrees of freedom=',int2str(stsv)])

disp(['50% ot   ',num2str(chi2inv(0.25,stsv)),'      do  ',num2str(chi2inv(0.75,stsv))])

disp(['60% ot   ',num2str(chi2inv(0.2,stsv)),'      do  ',num2str(chi2inv(0.8,stsv))])

disp(['70% ot   ',num2str(chi2inv(0.15,stsv)),'      do  ',num2str(chi2inv(0.85,stsv))])

========================================

результат: хи квадрат= 2.021;  delta =70.5


 

А также другие работы, которые могут Вас заинтересовать

49230. Механические часы, показывающие текущее время 160.41 KB
  Целью работы также является проведение сравнительного анализа языков программирования. Сравнить языки потребуется как в общем, так и применительно поставленной задаче. Одной из подзадач является создание описание алгоритма программы и составление математической модели.
49231. Разработка программы ускоренных эквивалентных испытаний для подтверждения назначенного ресурса газотурбинных установок 70.86 KB
  Обоснование применения выбранного метода испытаний ГТУ. Расчет коэффициентов ускорения испытаний и времени испытаний. Программа испытаний на малоцикловую усталость.
49232. Гидравлический расчёт трубопровода 226.48 KB
  АЛЕКСЕЕВА НАЧНООБРАЗОВАТЕЛЬНЫЙ ИНСТИТУТ ЯДЕРНОЙ ЭНЕРГЕТИКИ И ТЕХНИЧЕСКОЙ ФИЗИКИ Кафедра ЯРиЭУ Гидравлический расчёт трубопровода наименование проекта ПОЯСНИТЕЛЬНАЯ ЗАПИСКА к курсовой работе КОНСУЛЬТАНТЫ: РУКОВОДИТЕЛЬ Содержание курсовой работы Теоретическая часть Заданная схема трубопровода 6 Исходные значения Выбор и...
49233. Игра “Зенитная установка” 191.25 KB
  Цель выполнения курсовой работы формирование представления о взаимосвязи различных стадий единого процесса проектирования сложных систем. В результате выполнения работы необходимо получить практические навыки реализации основных стадий разработки: системного анализа проектирования программирования и тестирования. В результате выполнения работы были реализованы основные стадии разработки программного средства имеющего признаки сложной системы : системный анализ ;...
49234. 20 разрядный аналого-цифровой преобразователь, изготовленный по технологии КМОП 0,9 пм 873.5 KB
  Традиционные конструкции аналого-цифровых преобразователей АЦП использовали параллельную архитектуру и биполярные технологии для получения 8битного разрешения для преобразования с частотой 20 выборок с. Для более чем 8битной конструкции АЦП параллельной архитектуры требуют больших запрещенных районов и большей потребляемой мощности. Несмотря на то что биполярные технологии используются для самых быстрых преобразований для создания многоступенчатого АЦП были так же использованы БиКМОП технологии которые обеспечивают...
49235. Застосування комп’ютерних технологій в управлінні кадрами підприємства для організації ефективної роботи служби управління персоналом 437.33 KB
  Отже питання продуманого управління персоналом повинні займати одне з ключових місць у загальній стратегії розвитку підприємства. Важлива роль в оптимізації управління співробітниками підприємства належить автоматизованим системам управління персоналом так званим Humn Resource systems. Розвиток ринкових відносин у країнах світового економічного співтовариства свідчить що успіх роботи підприємств визначається розробкою всіх компонентів їх систем управління одним з яких є управління персоналом. Одним із найбільш складних але й...
49236. Обход Шахматной Доски Конем 206.97 KB
  Для решения этой задачи на компьютере необходимо разработать правила, в соответствии с которыми компьютер будет выбирать ход. В принципе, очередной ход можно выбирать случайным образом
49238. Минимальное покрытие 58.36 KB
  Требуется рассчитать минимальное покрытие для заданной таблицы, который можно редактировать с рабочей формы Визуализация В программе реализован алгоритм вершинного покрытия. Пользователь может изменят размеры таблицы. На форме две кнопки: определить минимальное покрытие - соответственно после нажатия определит минимальное покрытие, результат будет выведен в соседней таблице. Кнопка "перемешать" заполняет исходную таблицу случайными значениями.