50686

Моделирование дискретной случайной величины

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

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

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

Русский

2014-01-28

267 KB

3 чел.

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

Моделирование дискретной

случайной величины

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

Clc

clear all

format compact

format long

%--------------------------------------------------

k=input('Enter k=');    % ввод числа членов полинома

     

Enter k=30     

%--------------------------------------------------

vsv=1:k;     

%--------------------------------------------------

dz= vsv.^6.*pi^6;

p= 945./dz;      % ввод полинома

%--------------------------------------------------

cp=cumsum(p);    %сумма полинома

%--------------------------------------------------

figure(1)     

   subplot(2,1,1)    %график плотности распределения

   plot(p,'b.')     

   title('Density')    

  

  subplot(2,1,2)    

   bar(vsv+0.5,cp,1,'or')   % график ф-ции распред.

   title('Function of distribution')   

         

%--------------------------------------------------

n=input('Enter n=');     % кол-во случайных величин

     

Enter n=70     

%--------------------------------------------------

for t=n:-1:1,     

   x(t)=sum(cp<=rand)+1;   % генератор счлуч величин

end      

%--------------------------------------------------

figure(2)     

   plot(x,'m*')     

   title('discrett chance value')  % возможные значения случ вел.

   xlabel('N');     

   ylabel('value');    

%--------------------------------------------------

m=mean(x);     

sko=std(x);     

dissv=sko*sko;    

mt=sum(vsv.*p);    

dissvt=sum(vsv.*vsv.*p)-mt*mt;  % вывод теоретич. и эксп. величин

skot=sqrt(dissvt);    % матожид, дисперсии

     % квадратич. отклонения

Theoretic           Experimental

mean=1.0193        mean=1.03

disp=0.024998        disp=0.049596

sko=0.15811          sko=0.2227      

%--------------------------------------------------

disp('Theoretic           Experimental')

disp(['mean=',num2str(mt),'        mean=',num2str(m)])

disp(['disp=',num2str(dissvt),'        disp=',num2str(dissv)])

disp(['sko=',num2str(skot),'          sko=',num2str(sko)])

%--------------------------------------------------%

g=input('Enter level of reliability:  '); % пераразбиваем карманы для более точного

disp('Theoretic reliable interval')  % рассчета hi2

     

Enter level of reliability:  0.9   

%--------------------------------------------------

z=erfinv(g)*sqrt(2);     %рассчет доверительного интервала(теор)

delta=z*skot/sqrt(n);    

display(['At ',num2str(m-delta),' to ',num2str(m+delta)])

Theoretic reliable interval   

At 0.99901 to 1.061    

%--------------------------------------------------

disp('Experimental reliable interval')  

q=tinv((g+1)/2,n-1);    % рассчет доверительного интервала(практ.)

delta1=q*sko/sqrt(n);    

display(['At ',num2str(m-delta1),' to ',num2str(m+delta1)])

     

Experimental reliable interval  

At 0.98581 to 1.0742    

%--------------------------------------------------

v=hist(x,vsv);     

%--------------------------------------------------

figure(3)     

   bar(vsv,v,1,'m')    

   title('Histogram on k karmanov')  

   ylabel('N*freq');    

   xlabel('n')     

%--------------------------------------------------

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

karm=ceil(k/gr);

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

v1=sum(reshape([v zap],gr,karm));  

p1=sum(reshape([p zap],gr,karm))*n;

%--------------------------------------------------

figure(4)     

   xc=(gr+1)/2:gr:k+gr/2;   

   bar(xc,v1,1,'g')    

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

   xlabel('Intervals')    

   ylabel('N*freq')    % вывод распределения случ величины

    hold on      

       bar(xc,p1,0.8,'b')    

   hold off     

%--------------------------------------------------

hi2=sum(((v1-p1).^2)./p1);  

stsv=karm-1;     

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))])

disp(['80% ot   ',num2str(chi2inv(0.10,stsv)),'      do  ',num2str(chi2inv(0.90,stsv))])

Enter n=50

Vvedite k=10000

Theoretic           Experimental

mean=1.0193        mean=1.0173                 Enter level of reliability:  0.95

disp=0.024998        disp=0.024403             Theoretic reliable interval At 1.0142 to 1.0204

sko=0.15811          sko=0.15622                  Experimental reliable interval  At 1.0142 to 1.0204

How many karmanov in group:  2                hi2 =0.2152     Degrees of freedom=24

50% ot   19.0373      do  28.2412                  60% ot   18.0618      do  29.5533

70% ot   16.9686      do  31.1325                  80% ot   15.6587      do  33.1962


 

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

64117. Структура управления и ее влияние на эффективность работы торгового предприятия (на примере ЗАО «Связной Логистика») 1.74 MB
  Роль структуры управления в эффективном управлении организацией Анализ организационной структуры и структуры управления организацией Анализ управления спросом на товары и услуги предприятия...
64118. Технологія зборки та монтажу плати блоку біжучого рядка на світлодіодах 363.5 KB
  Яскрава картинка яка динамічно змінюється привертає увагу і допомагає донести велику кількість інформації використовуючи мінімум простору. принцип управління схемний; кількість входів 23; Кількість виходів 21; входи з клавіатури та блоку живлення...
64119. Оcновные нaпpaвления деятельноcти оpгaнов пpедвapительного pacледовaния в cовpеменный пеpиод 1.92 MB
  Пpо‏во‏димые‏ в Pо‏ccии cуде‏бно‏ пpaво‏вые‏ pе‏фо‏pмы нaпpaвле‏ны нa эффе‏ктивную paбо‏ту пpaво‏о‏xpaните‏льныx о‏pгaно‏в по‏ бо‏pьбе‏ c пpе‏cтупно‏cтью в cо‏че‏тaнии c гумaнным о‏тно‏ше‏ние‏м к че‏ло‏ве‏ку, по‏пaвше‏му в cфе‏pу уго‏ло‏вно‏- пpо‏це‏ccуaльныx о‏тно‏ше‏ний.
64120. Потребительский кредит 257.5 KB
  История возникновения потребительского кредита Принципы и методы потребительских кредитования. Все факты экономического развития России говорят о том что необходимо уделять большое внимание проблеме кредита в том числе и потребительского так как экономическое состояние страны...
64121. ФЛОРА ЛЕСОВ МАЙСКОЙ ГОРЫ ЗАВОЛЖСКОГО РАЙОНА ГОРОДА УЛЬЯНОВСКА И ЕЁ КОМПЛЕКСНЫЙ АНАЛИЗ 456.5 KB
  Богат в лесах видовой состав кустарников: бересклет бородавчатый орешник жимолость лесная крушина. По видовому составу к кустарниковым степям близки луговые степи с преобладанием костреца берегового занимающие чаще всего северные склоны.
64122. Проектирование технологического процесса работы отделения топливной аппаратуры грузовых автомобилей 1.04 MB
  Самое простое отнести форсунки в ближайший дизель-сервис и оставить там некоторое количество денег. Для начала надо разобраться в устройстве дизельной форсунки и в том что там происходит и что на что может влиять. Сечение форсунки представлено на рис.
64124. Автоматизация процессов сбора, передачи, обработки и накопления информации в системе управления персоналом предприятия 4.33 MB
  Цель исследования: обосновать пути автоматизации процессов сбора, передачи, обработки и накопления информации в системе управления персоналом предприятия (на примере ООО «Ковровый двор»).
64125. КРИМІНАЛЬНА ВІДПОВІДАЛЬНІСТЬ ФІЗИЧНИХ ОСІБ ЗА МІЖНАРОДНІ ЗЛОЧИНИ 167 KB
  Розкрити поняття міжнародного злочину та встановити, які саме діяння містять його ознаки; визначити компетенцію Міжнародного кримінального суду; охарактеризувати міжнародний та регіональний рівні співробітництва держав у боротьбі з міжнародної злочинністю.