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


 

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

38844. Методические рекомендации по выполнению дипломного проекта 1.19 MB
  Дипломная работа выполняется в форме дипломного проекта. Каждому студентудипломнику назначается руководитель проекта как правило из числа преподавателей вуза. Выпускная квалификационная работа выполняется в форме проекта в соответствии с утвержденным Советом вуза Положением о выпускной квалификационной работе.
38845. Экспертная система комплексного диагностирования линейной части магистрального газопровода ЭС «Диагностика ЛЧ МГ» 3.02 MB
  В рамках проделанной работы проведено детальное обследование предметной области, составлен перечень функциональных задач, исследованы аналоги автоматизированных систем, спроектирована инфологическая модель предметной области, реализована ЭС «Диагностика ЛЧ МГ»
38846. МЕТОДИЧЕСКИЕ УКАЗАНИЯ К ВЫПОЛНЕНИЮ ДИПЛОМНОЙ РАБОТЫ 186 KB
  Плеханова Кафедра документоведения и документационного обеспечения управления МЕТОДИЧЕСКИЕ УКАЗАНИЯ К ВЫПОЛНЕНИЮ ДИПЛОМНОЙ РАБОТЫ Для студентов специальности 032001. Петраченко Методические указания к выполнению дипломной работы для студентов спец. Методические указания к выполнению дипломной работы разработаны преподавателями кафедры документоведения и документационного обеспечения управления Саратовского государственного социальноэкономического университета на основании требований государственных образовательных стандартов высшего...
38847. ДИПЛОМНА РОБОТА СПЕЦІАЛІСТА, ЇЇ МЕТА І ЗАВДАННЯ 253 KB
  Вона є кваліфікаційним документом на підставі якого Державна екзаменаційна комісія визначає рівень теоретичної підготовки випускника його готовність до самостійної роботи за фахом і приймає рішення про присвоєння кваліфікації. Працюючи над ДР студент має засвоїти навики правильної постановки проблеми та обґрунтування її актуальності формулювання мети і завдань дослідження побудови логічного плану оптимальної структури роботи з літературними джерелами і статистичною інформацією аналізу та оцінки різних аспектів діяльності організації...
38848. ДИПЛОМНОЕ ПРОЕКТИРОВАНИЕ 231.5 KB
  Оно проводится в целях выполнения квалификационной работы проекта соответствующей государственным требованиям к уровню подготовки инженера систематизации и закрепления знаний по образовательнопрофессиональной программе совершенствования умений их применения для решения задач в области мостостроения. Следует иметь в виду что основную ответственность за правильность принятых в проекте технических решений и всех данных несёт студент автор дипломного проекта. Темы дипломных проектов как правило соответствуют одному из направлений:...
38849. Курсовые и дипломные работы по теории обучения иностранным языкам: практические рекомендации для студентов вузов 592 KB
  Романовская Курсовые и дипломные работы по теории обучения иностранным языкам: практические рекомендации для студентов вузов Учебнометодическое пособие по теории обучения иностранным языкам Ульяновск 2008 ББК 81. Курсовые и дипломные работы по теории и методике обучения иностранным языкам: практические рекомендации для студентов вузов : учеб. Даются рекомендации по поиску и изучению литературы написанию теоретической части планированию и проведению эмпирического исследования подготовке и проведению защиты дипломной работы. Для студентов...
38850. Эстетическое воспитание подростков на занятиях по художественному текстилю на примере серии работ «Флора» в технике холодный батик 1.24 MB
  ИСТОРИКОКУЛЬТУРОЛОГИЧЕСКИЙ АСПЕКТ СТАНОВЛЕНИЯ И РАЗВИТИЯ ХУДОЖЕСТВЕННОЙ РОСПИСИ ТКАНИ И ЕЁ ТРАНСФОРМАЦИЯ В СОВРЕМЕННОМ ТВОРЧЕСТВЕ. История возникновения и развития художественной росписи ткани в странах Азии и Европы10 Исторические аспекты развития художественной росписи ткани России20 1.3 Роль художественной росписи ткани в современном декоративно прикладном искусстве в контексте моды как культурного феномена и ее связи с современным искусством росписи...
38851. Обучающая подсистема для лабораторного исследования характеристик замкнутых САУ в среде интернет 4.57 MB
  В данном дипломном проекте рассматривается обучающая подсистема для лабораторного исследования характеристик замкнутых САУ в среде интернет. В экономической части дается техникоэкономическое обоснование разработки Обучающей подсистемы для лабораторного исследования характеристик замкнутых САУ в среде интернет с помощью частотных критериев устойчивости проводится расчет ее сметной стоимости и стоимости эксплуатации Содержание.1 Описание предметной области по характеристикам замкнутых САУ.
38852. Фотокаталітичне очищення стічних вод від біогенних елементів 5.96 MB
  Представлено та описано технологічну схему очищення стічних вод від біогенних елементів фотокаталітичною технологією за допомогою солей Ti(IV) та Zr(IV). Виконано розрахунок і вибір основного та допоміжного обладнання у відповідності з задоною продуктивністю установки. Надано його характеристику.