50686
Моделирование дискретной случайной величины
Лабораторная работа
Информатика, кибернетика и программирование
Цель работы. Практическое освоение алгоритма программной генерации дискретной случайной величины и методов статистической проверки разработанного генератора.
Русский
2014-01-28
267 KB
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
А также другие работы, которые могут Вас заинтересовать | |||
6794. | УПРАВЛЯЕМЫЙ ВЫПРЯМИТЕЛЬ | 139 KB | |
УПРАВЛЯЕМЫЙ ВЫПРЯМИТЕЛЬ Цель работы - ознакомление с принципом работы и экспериментальное исследование характеристик однофазного управляемого выпрямителя на тиристорах при различных видах нагрузки. 2.1 Основные положения Для регулирования выпря... | |||
6795. | ПОЛУПРОВОДНИКОВЫЕ СТАБИЛИЗАТОРЫ | 126 KB | |
ПОЛУПРОВОДНИКОВЫЕ СТАБИЛИЗАТОРЫ Цель работы - исследование рабочего участка вольт-амперной характеристики полупроводникового стабилитрона, изучение основных характеристик параметрического и компенсационного стабилизаторов напряжения и определен... | |||
6796. | Усилители постоянного тока | 121.5 KB | |
Усилители постоянного тока Цель работы - изучение принципа работы простейшей схемы дифференциального усилителя и экспериментальное определение его параметров ознакомление с работой операционного усилителя (ОУ) и экспериментальное определение е... | |||
6797. | Усилительные каскады на транзисторах | 135.5 KB | |
Усилительные каскады на транзисторах Цель работы - исследование одиночных усилительных каскадов RC-типа на биполярном и полевом транзисторах исследование характеристик и определение параметров изучаемых усилителей. 5.1 Основные положения Усилители ... | |||
6798. | Аналоговые схемы на операционных усилителях | 146.5 KB | |
Аналоговые схемы на операционных усилителях Цель работы - изучение некоторых схем включения операционного усилителя для обработки аналоговых сигналов определение характеристик и параметров инвертирующего и неинвертирующего усилителей, сумматор... | |||
6799. | Импульсные схемы на операционных усилителях | 123.5 KB | |
Импульсные схемы на операционных усилителях Цель работы - изучение принципа работы компаратора и триггера Шмитта на операционном усилителе исследование и определение параметров схем мультивибратора, одновибратора и генератора треугольных импульсов... | |||
6800. | Генераторы импульсов на транзисторах и интегральных микросхемах | 98 KB | |
Генераторы импульсов на транзисторах и интегральных микросхемах Цель работы - ознакомление с принципом работы мультивибратора на биполярных транзисторах и определение его параметров изучение принципа работы мультивибратора, одновибратора и генерато... | |||
6801. | Электронная таблица Microsoft Office | 40 KB | |
Электронная таблица MicrosoftOffice Задание: Построить таблицу значений координат двух графиков F1 и F2. По полученной таблице создать диаграмму с тремя графиками: F1, F2 и их суммы. Найти среднее значение на заданном диапазоне для каждо... | |||
6802. | Определение параметров электрических колебаний | 208 KB | |
Цель работы: В лабораторной работе ставится целью определение параметров электрических колебаний по их осциллограммам. Содержание отчета. Измерение параметров синусоидального напряжения с помощью осциллографа. Схема исследуемой электр... | |||