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


 

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

8624. Философия Гераклита и Демокрита 33 KB
  Философия Гераклита и Демокрита Афоризмы Гераклита: Эту-вот Речь (Логос) сущую люди не понимают и прежде, чем выслушать её, и выслушав однажды. Ибо, хотя все люди сталкиваются напрямую с этой-вот Речью (Логосом), они подобны незнающим её, даро...
8625. Сократ и софисты 113 KB
  Сократ и софисты Вопросы: 1. Какие определения благочестия даёт Евтифрон. Как опровергает эти определения Сократ? 3. Какова цель Сократовских бесед, если ответа в конце каждого диалога не даётся? 3. Можно ли вообще дать определение благочестию? Е...
8626. Философия Аристотеля. Аристотель кается перед Платоном 92.5 KB
  Философия Аристотеля Вопросы: 1. Что Аристотель называет началом, которым «должен владеть всякий, кто постигает какую-либо вещь»? При помощи какого метода строит свою философию Аристотель? 2. В чём состоит критика Аристотелем Платоновской теории иде...
8627. Философия Средних веков (богословие) и эпохи Возрождения 130 KB
  Философия Средних веков (богословие) и эпохи Возрождения Вопросы: 1. В чём отличие богословия от философского мышления? 2. Какие ступени познания выделяет Плотин? В чём отличие учения Плотина от богословия? 3. Какие два монотеизма противостоят друг ...
8628. Философия Р.Декарта и Ф.Бэкона 249 KB
  Философия Р.Декарта и Ф.Бэкона Р.Декарт РАЗЫСКАНИЕ ИСТИНЫ ПОСРЕДСТВОМ ЕСТЕСТВЕННОГО СВЕТА, который сам по себе, не прибегая к содействию религии или философии, определяет мнения, кои должен иметь добропорядочный человек относительно всех предметов, ...
8629. Философия И.Канта 102.5 KB
  Философия И.Канта Вопросы: 1. В чём отличие эмпирических знаний от априорных? 2. Какие признаки априорного знания выделяет И.Кант? Какова задача философии, по Канту? Чем обусловлена постановка задачи? 3. Какую классификацию суждений предлагает Кант?...
8630. Русская религиозная философия 19-20 веков 68 KB
  Русская религиозная философия 19-20 вв. Вопросы: 1. Какой должна стать русская философия? С каким предметом она должна иметь дело? 2. Как И.Ильин аргументирует необходимость отказа от намеренного выдумывания философских систем? 3. Каков закон иссл...
8631. Сознание как философская категория 56.5 KB
  Сознание как философская категория Вопросы: 1. Какие правила мышления выделяет И.Кант? Как понимается мышление в европейской традиции? 2. Что такое сознание? Как формулирует своё предназначение человек сознающий? 3. Каковы свойства гармонии, носител...
8632. Проблема познания в философии 105 KB
  В каком смысле акты научного познания - свободные явления. Можно ли рассматривать знание законов как актуализацию готовых смыслов и сущностей. Можно ли рассматривать знание как превращение...