Лабораторная работа 2 Моделирование непрерывной случайной величины методом обратных функций Цель рабо
Работа добавлена на сайт samzan.net:
Лабораторная работа № 2
Моделирование непрерывной случайной величины методом обратных функций
Цель работы. Практическое освоение алгоритма программной генерации непрерывной случайной величины и методов статистической проверки разработанного генератора.
clear all;
clc
N=200;
%-------------генератор---------------------------
z=rand(N,1);
x=-0.24+sqrt(0.48^2+12*z)/2;
%---------построение графиков -----------------------------------
a=0:0.001:1.5;
figure(1);
subplot(3,1,1);
plot(a,f(a),'b-');
title(['f(x) - plotnost raspred']);
xlabel('x');
ylabel('p(x)');
subplot(3,1,2);
plot(a,Fx(a),'b-');
title(['F(x) - f-ciya raspred']);
xlabel('x');
ylabel('P(X<x)');
Fe=((1:N)-0.5)/N;
x1=sort(x);
Ft=Fx(x1);
subplot(3,1,3);
plot(x1,Ft,'b-');
title(['F(x) - f-ciya raspred teor i praktich']);
xlabel('x');
ylabel('P(X<x)');
hold on
plot(x1,Fe,'r-');
hold off
%-----sлучайн величина-------------------
figure(2)
plot(x,'r*')
xlabel('n')
ylabel('sluch vel')
title('sluch vel')
%-------гистограмма----------------------
r=input('na skolko karmanov delit interval ')
figure(3);
hist(x,r);
ylabel('N * freq')
xlabel('values of sv')
title('Histogram sluch velich')
%----мат ожидание------------------------
y=2/3*a.^2+0.16*a;
disp(['теор мат ожид. =']);
mat_t=trapz(a,y)
disp(['практич. мат ожид. =']);
mat_pr=mean(x)
hist(x,r);
%---- dispersiya---------------------------
y=(a-mat_t).^2.*(2/3*a+0.16);
dis1=trapz(a,y)
dis2=cov(x)
%-----------доверительный интерв---------
s=std(x);
st=sqrt(0.1502);
aa=input('vvedi urov dover: a=');
za=erfinv(aa)*sqrt(2);
delta=(za/sqrt(N))*st;
disp(['Доверительный интервал от ',num2str(mat_pr-delta),' до ',num2str(mat_pr+delta)])
if N<=121,
zt=tinv((aa+1)/2,N-1);
delta2=(zt/sqrt(N))*s;
else
delta2=(za/sqrt(N))*s;
end
disp(['Уточнённый доверительный интервал от ',num2str(mat_pr-delta2),' до ',num2str(mat_pr+delta2)])
%--------хи квадрат---------------------
[nu,xn]=hist(x,r);
a=0:1.5/r:1.5;
F=Fx(a);
for i=2:r
p(i-1)=F(i)-F(i-1);
end
p(r)=1-sum(p);
nt=N*p;
x1=(nu-nt).^2;
x2=x1./nt;
xi2=sum(x2);
disp([int2str(r-1),' степеней свободы' ])
disp(['xi2 = ', num2str(xi2)])
%----------омега квадрат----------------
x_sort=sort(x);
F=Fx(x_sort);
for i=1:r
w2(i)=(F(i)-(i-0.5)/r).^2+1/(12*r);
end
w=sum(w2);
disp(['w2 = ', num2str(w)])
%---------------------------------------
kk=0:1.5/r:1.5;
Ftr=Fx(kk);
Pt=(Ftr(2:r+1)-Ftr(1:r))*N;
Ck=kk(1:r)+1/r;
Pe=hist(x,Ck);
figure(5)
bar(Ck,Pt,1,'b')
title('Rasch hist (theor(blue) and exp(red))')
ylabel('N * freq')
xlabel('values of sv')
hold on
bar(Ck,Pe,0.8,'r')
hold off
Fx
function yy=Fx(e);
yy=e.^2/3+0.16*e;
f
function xx=l2_1(q);
xx=2*q/3+0.16;
==========================================================================
Результаты:
Теор. матожид. =0,93 теор. дисперсия =0,1502
Практ. матожид. =0,9261 практ. дисперсия =0,1526
Интервал доверия а=0,95
Доверительный интервал от 0,88103 до 0,97119
Уточненный доверительный интервал от 0,88067 до 0,97155
5 степеней свободы
хи квадрат 5,9237
омега квадрат 0,4316
2. Адаптивная система VPN в распределенных компьютерных сетях
3. Національні економічні інтереси та загрози економічній безпеці України
4. Вашингтоне построен огромный официальный Музей Памяти Холокоста
5. под статистикой понимают отрасль практической деятельности которая имеет своей целью сбор обработку ана
6. і. СУ працівників розглядає як основу організації і орієнтоване на довгострокову перспективу
7. Мы это терпеливо выслушиваем попутно говоря свой пол и имя.html
8. Политическую коррупцию неизбежно порождаемую организованной преступностью можно смело отнести к той боле
9. вещество ' одна из центральных философских категорий обозначающая объективную реальность существующую
10. Лісове і садовопаркове господарство Фахове спрямування ~ Лісове господарство Садовопаркове госпо
11. История острова Пасхи
12. Анды
13. Тема- Душ человеческих добрые лекарибиблиотекари Выполнила- Новоселова В
14. варіант відповіді- I
15. ЛЕКЦИЯ ХРЯЩЕВЫЕ ТКАНИ 2002г
16. М. 1960 Теория языкового знака ныне явно или неявно принятая в большинстве сочинений по общей лингвист
17. 2 Дата проведения- Место проведения- стадион БГАФК
18. Контрольная работа по дисциплине- Теория Государства и Права Вариант 2 Выполнила- студентка 1 кур
19. Лабораторная работа 12
20. варіант ’ в журналі Обчислити- визначник матриці ; транспоновану матрицю ; обернену матр