Лаба 1 (Моделирование непрерывных и дискретных случайных величин) / Lab1_1

clc;
clear all;
close all;


k=0:1:29;
p=2.^(-k-1);
p1=zeros(1,30);

for a=1:1:30
p1(1,a)=2^(-k(1,a)-1);
if (sum(p1)>0.999) && (sum(p1)<0.9995)
disp('Выберем лишь несколько элементов');
disp(a-1);
end
end

clear p1 a k p

a(1,1)=0.1;
a(2,1)=0.05;
a(3,1)=0.01;

N=[50; 200; 1000;];

KSID=zeros(3,3);
KSIH=zeros(3,3);
for j =[1 2 3]
for i =[1 2 3]
KSID(i,j)=chi2inv((1-(a(j,1)/2)),(N(i,1)-1)); %p=1-a; for a/2 = p=1-a/2 for 1-a/2 p=1-(1-a/2)=a/2;
KSIH(i,j)=chi2inv((a(j,1)/2),(N(i,1)-1));
end
end


r=rand(1,50);
i1=50;
r50=SV(r,i1); %Наша ДСВ

r=rand(1,200);
i1=200;
r200=SV(r,i1); %Наша ДСВ

r=rand(1,1000);
i1=1000;
r1000=SV(r,i1); %Наша ДСВ

clear r i1

M(1,1)=mean(r50);
M(2,1)=mean(r200);
M(3,1)=mean(r1000);

S(1,1)=std(r50);
S(2,1)=std(r200);
S(3,1)=std(r1000);

D(1,1)=var(r50);
D(2,1)=var(r200);
D(3,1)=var(r1000);

%Коэффициенты для альфа 0.5
St(1,1)=0.683;
St(2,1)=0.667;
St(3,1)=0.675;
%Коэффициенты для альфа 0.1
St(1,2)=1.711;
St(2,2)=1.671;
St(3,2)=1.646;
%Коэффициенты для альфа 0.01
St(1,3)=2.797;
St(2,3)=2.660;
St(3,3)=2.581;


INTH=zeros(3,3);
INTD=zeros(3,3);
%Интервальная оцека среднего
for j=[1 2 3]
for i=[1 2 3]
[INTH(i,j), INTD(i,j)]=interval(D(i,1),M(i,1),St(i,j),N(i,1));
end
end

%Интервальная оцека дисперсии
DISPH=zeros(3,3);
DISPD=zeros(3,3);
for j=[1 2 3]
for i=[1 2 3]
[DISPH(i,j), DISPD(i,j)]=dispersia(S(i,1),KSIH(i,j),KSID(i,j),N(i,1));
end
end


%СВЕДЕНИЕ ЗНАЧЕНИЙ В ТАБЛИЦУ
T_p90=table(N,INTD(:,1),INTH(:,1),DISPD(:,1),DISPH(:,1));
T_p95=table(N,INTD(:,2),INTH(:,2),DISPD(:,2),DISPH(:,2));
T_p99=table(N,INTD(:,3),INTH(:,3),DISPD(:,3),DISPH(:,3));

%Гистограммы
figure('Name','Гистограммы','NumberTitle','off');
subplot(2,2,1);
h=0:8;
p2=2.^(-h-1);
bar(h,p2);
title('Теоретические данные');
subplot(2,2,2);
histogram(r50,'Normalization','Probability');
title('Пр. результаты для 50 элементов');
set (gca, 'Ylim', [0 0.6]);
subplot(2,2,3);
histogram(r200,'Normalization','Probability');
title('Пр. результаты для 200 элементов');
set (gca, 'Ylim', [0 0.6]);
subplot(2,2,4);
histogram(r1000,'Normalization','Probability');
title('Пр. результаты для 1000 элементов');
set (gca, 'Ylim', [0 0.6]);
clear h p2