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


 

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

67026. Формула професійного успіху (рольова гра «Зустріч у шкільному центрі зайнятості») 236.5 KB
  Ви мабуть замислювались над питанням, що робить людину успішною у наш стрімкий час. Одним із показників успішності є вибір професії. Я вважаю, що саме вибір професії – тяжкий та відповідальний крок у житті кожної людини – знайти своє місце у майбутній долі.
67027. Літературний вечір «Украдене щастя» (інсценізація однойменної драми Івана Франка) 64 KB
  Мета: розкрити глибину філософського змісту драми Івана Франка «Украдене щастя»; розвивати вміння аналізувати художній твір, творчі здібності учнів, навички акторської майстерності; виховувати любов до українського слова, традицій рідного краю.
67028. ТРАНСФОРМАЦИЯ КОРПОРАТИВНОЙ КУЛЬТУРЫ ДЛЯ ДОСТИЖЕНИЯ СТРАТЕГИЧЕСКИХ ЦЕЛЕЙ 645.5 KB
  Корпоративная культура может играть решающую роль в мобилизации всех ресурсов организации на достижение целей, но также может стать тормозом ее дальнейшего развития, поэтому особую актуальность в настоящее время приобрел вопрос возможности воздействия на процессы формирования, поддержания или изменения корпоративной культуры...
67029. Шкільний проект духовно – морального виховання 256.5 KB
  Вічним джерелом цієї проблеми є нестримне прагнення кожної людини визначитись з основами своєї діяльностіпоказати те чому завдячує людинароблячи той чи інший крок у своєму житті. Невід’ємною ланкою духовності людини є її культура. Поняття культура нерозривно пов’язане з усіма видами діяльності людини: політичнимекономічним правовим релігійним етичним художнім тощо.
67030. Путь к гармонии 100 KB
  Оборудование: мини-проекты толковый словарь правила Шаги к успеху и гармонии эмоциональный термометр индикаторы настроения ступенчатая шкала Этапы пути к гармонии. А чтобы достигнуть всего этого человек прежде всего должен находиться в гармонии с самим собой.
67031. Не руйнуй гармонії земної 107.5 KB
  Глобальною проблемою стали кислотні дощі, які гублять рослинність, знищують життя в прісних водах, це і негативні наслідки науково-технічної революції, неконтрольоване зростання населення Землі, кількість якого вже перевищила критичну межу, дедалі більше забруднення атмосфери, гідросфери.
67032. Що таке гендерна рівність? Всі ми – різні, але всі ми – однакові 4.19 MB
  Мета: ознайомити дітей з поняттями “гендер”, “гендерна рівність” та основами гендерно-правового законодавства України. Формувати в учнів інтерес до навчання з даної теми. Розвивати мислення, пам’ять, увагу та активізувати словниковий запас. Виховувати людяність, гуманність та толерантне ставлення один до одного, бажання приносити користь у суспільстві.
67033. Я люблю географію 100.5 KB
  Організаційний момент Сьогодні між командами змагання Нехай образ не буде серед вас Бо переможцем або переможеним Сьогодні може бути хтось із вас. За кожну правильну відповідь команда отримує 1 бал а за музичний конкурс 2 бали. 1 раунд Географічний лікбез командам прпонується однакове завдання: необхідно на слух правильно написати географічну...
67034. Географічна гра «Щасливий випадок» 69 KB
  Мета гри – дає можливість перевірити рівень засвоєних знань,виявити вміння застосовувати набуті знання, розвивати інтерес і пізнавальну діяльність учнів,сприяє розвитку навичок спілкування між підлітками. В грі приймають участь 2 команди по 5-6 учнів в кожній.