22

Mx=20;
sigma=4,47;
N=10000;
c=2*pi;
a1=rand(1,N);
% сгенерировать a1 размерностью N
a2=rand(1,N);
% сгенерировать a2 размерностью N
r=sqrt(-2*log(a1));
f=a2*c; X1=r.*sin(f);
Y=Mx+sigma*X1;
% проверка результатов

figure(1)
hist(Y,50)% графическое представление выборки Х с помощью гистограммы
title(['МО нормального распределения для N=', num2str(N)])

f=sigma*sigma*a1
figure(2)
hist(f,50)% графическое представление выборки Х с помощью гистограммы
title(['Дисперсия нормального распределения для N=', num2str(N)])



z=(a1+a2)/2
figure(3)
hist(z,50)% графическое представление выборки Х с помощью гистограммы
title(['МО равномерного распределения для N=', num2str(N)])



w=a2-a1
q=w.^2/12
figure(4)
hist(q,50)% графическое представление выборки Х с помощью гистограммы
title(['Дисперсия равномерного распределения для N=', num2str(N)])