ЛР / ЛР1 / lr1

x = [0.4480617715 0.2641551216 -2.906882284 2.031714874 1.687223292 ...
3.285391929 -1.914812285 -2.263417057 6.329957276 1.987025714 ...
-1.568758646 -1.843516995 4.454997568 -2.227162324 -1.384504114 ...
-1.483698154 -4.293654795 -1.660694918 -0.729183623 -1.707611202 ...
2.222858828 -2.475120138 1.723317726 2.419616638 5.006546988 ...
2.301841621 -0.4986608221 -1.38481789 -1.840283742 -1.914075594 ...
3.28477802 -0.9601702661 5.455314411 -6.507696864 -0.2056693777 ...
-5.364625733 3.012461664 0.860533535 -0.1694316981 0.6340829966 ...
-1.37884251 4.321677013 0.9760166247 1.736793027 -0.7582590317 ...
9.297647582 -0.9034177967 2.801222425 4.618215487 -2.231671144];

%Построение вариационного ряда
x = sort(x);
%Поиск максимального и минимального элементов выборки
xmax = max(x);
xmin = min(x);
%Определим количество интервалов по формуле Стерджесса
b = 3.332;
n = 50;
r = ceil(1+fix(b*log10(n)));
%Длина интервала
stp = (xmax-xmin)/r;
%Определяем середины интервалов
centr = [];
centr(1) = xmin+(stp/2);
for i=2:1:r
centr(i) = centr(i-1)+stp;
end
%Определяем абсолютную частоту
k1 = xmin;
i = 1;
while i<=r
k2 = 0;
for j=1:n
if (x(j)>=k1) & (x(j)<=k1+stp)
k2 = k2+1;
end
end
freqn(i) = k2;
k1 = xmin+stp*i;
i = i+1;
end
% Числовые характеристики выборки:
% Выборочное среднее
m = mean(x);
% Дисперсия
D = var(x);
% Ср. кв. отклонение
SKO = std(x);
% Мода
moda = mode(x);
% Медиана
med = median(x);
% Коэффициент эксцесса
kurt = kurtosis(x);
% Коэффициент асимметрии
skew = skewness(x);



% Числовые характеристики выборки:
% Выборочное среднее
m = mean(x);
% Дисперсия
D = var(x);
% Ср. кв. отклонение
SKO = std(x);
% Мода
moda = mode(x);
% Медиана
med = median(x);
% Коэффициент эксцесса
kurt = kurtosis(x);
% Коэффициент асимметрии
skew = skewness(x);

% Вывод значений
fprintf('Максимальное значение = %f\n',xmax);
fprintf('Минимальное значение = %f\n',xmin);
fprintf('Количество интервалов = %f\n',r);
fprintf('Длина одного интервала = %f\n',stp);
fprintf('Выборочное среднее = %f\n',m);
fprintf('Выборочная дисперсия = %f\n',D);
fprintf('Ср. кв. отклонение = %f\n',SKO);
fprintf('Мода = %f\n',moda);
fprintf('Медиана = %f\n',med);
fprintf('Коэффициент эксцесса = %f\n',kurt);
fprintf('Коэффициент асимметрии = %f\n',skew);

% Построение полигона частот
figure()
plot(centr,freqn/n,'r-o')
xlabel('Интервалы');
ylabel('Относительная частота')
grid on

% Построение гистограммы
figure()
%histogram(x,r)
hist(x,r)
xlabel('Интервалы')
ylabel('Частота')
grid on

% Построение эмпирической
% функции распределения
figure()
ecdf(x)
% Подпись оси 0X
xlabel('x')
% Подпись оси 0Y
ylabel('F(x)')
% Добавление сетки на график
grid on