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


 

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

43621. Управление мотивацией персонала на примере ГУП “Таткоммунпромкомплект” 943.46 KB
  Анализ деятельности предприятия ГУП â€œТаткоммунпромкомплектâ€. Вместе с тем работающие предприятия имеют свои мотивационные модели практическое использование которых позволяет им держаться на плаву и даже увеличивать объем производства. Применяются ли они на наших предприятиях кроме конечно чисто русской модели кнута и пряника К сожалению нет. Таким образом учитывая вышеизложенное нами выше можно сделать следующий вывод: весь персонал предприятия условно можно подразделить на определенные...
43622. Разработка коробки скоростей токарно-револьверного станка 979.5 KB
  К достоинствам схемы относятся высокий КПД, удобство монтажа и демонтажа. К недостаткам – неравномерное распределение нагрузки вдоль зуба и неодинаковое давление на подшипники (установленные на разных концах вала) от окружного усилия и веса колёс.
43623. Информатика. Информация. Устройство ЭВМ 1.08 MB
  Информатика – современная научная база информационной сферы деятельности людей. Не существует однозначного понятия информатики – как науки. Рассмотрим следующие: Информатика – это наука, изучающая все аспекты получения, хранения, преобразования, передачи и использования информации.
43624. Опыт осуществляемый первичной профилактикой наркомании в гимназии РГСУ 60.45 KB
  Программа наркомании несовершеннолетних в Российской Федерации. Социальная профилактики наркомании среди подростков. Методы работы социального педагога по профилактике наркомании среди несовершеннолетних Глава 2. Опыт осуществляемый первичной профилактикой наркомании в гимназии РГСУ.
43625. Формирование диалогической речи детей дошкольного возраста с задержкой психического развития по картинке 65.27 KB
  Первое условие потребность детей высказываться; второе о чём нужно сказать т. и обобщенные в монографии“Обучение детей с задержкой психического развития†1981 констатируют у детей с ЗПР неустойчивость внимания недостаточность развития фонематического слуха зрительного и тактильного восприятия оптикопространственного синтеза моторной и сенсорной стороны речи долговременной и кратковременной памяти зрительномоторной координации автоматизации движений и действий. Общение необходимо для психического развития детей...
43626. Влияние применения ближнего света фар в дневное время суток на безопасность дорожного движения 1.99 MB
  В разделе «Введение» изучается значение станций технического обслуживания, необходимость их планирования, приводится описание станций технического обслуживания, видов СТО и производимых на них работах и услугах. Раздел «Технологический расчет»: представляет собой расчет производственной программы, годового объема работ, количества постов, численности рабочих и площадей помещений производственного корпуса.
43627. Анализ и оптимизация САР частоты вращения вала двигателя постоянного тока 2.54 MB
  Даля кафедра автоматизации и компьютерноинтегрированных технологий Курсовой проект по дисциплине Теория автоматического управления на тему: анализ и оптимизация САР частоты вращения вала двигателя постоянного тока В №13.2 Описание принципа действия САР 1.3 Определение передаточных функций элементов САР 1.4 Построение модели структурной схемы САР 1.
43628. УПРАВЛЕНИЕ ПЕРСОНАЛОМ Т.В. Зайцева, А.Т. Зуб 2.12 MB
  В учебнике рассматриваются важнейшие вопросы теории и практики управления персоналом организации: история становления данной науки; принципы, методы и подходы к построению системы управления персоналом в организации; основные процессы. Проанализированы различные аспекты основных направлений деятельности в рамках управления персоналом: планирование потребности в персонале и найм, оценка результативности труда, мотивация и оплата труда персонала...
43629. Организация проведения ТО и ТР автомобилей УАЗ -31512 в условиях ОАО Ростелеком г. Ухта 1.04 MB
  Транспорт важнейший элемент инфраструктуры, под который понимают отрасли народного хозяйства, создающие общие его функционирования. Транспорт оказывает активное влияние на процесс расширенного воспроизводства, величину запасов, сырья, топлива и промышленной продукции, производственную мощность складов