МП-2 / eterver

function [] = TerVer( A )

format rational;

n1=length(A);

%сортировка

X=sort(A);

disp('Вариационный ряд:');

disp(X);

fprintf('Размах выборки = x(%g)-x(1) = %g\n\n',n1,X(n1)-X(1));

i=1;

%статистический ряд

disp('Статистический ряд:');

for xx=-50:50

K=find(X==xx);

if length(K)~=0

B(1,i)=xx;

B(2,i)=length(K);

i=i+1;

end

end

disp(B);

%гистограмма и полином(группа)

figure(1);

n2=length(B);

d=(B(1,n2)-B(1,1))/7;

Z(1)=B(1)+d/2;

for i=1:6

Z(i+1)=Z(i)+d;

end

disp('Середины разрядов:');

disp(Z);

for i=1:7

K=find(Z(i)-d/2<=X & X<Z(i)+d/2);

N(i)=length(K);

end

N(7)=N(7)+B(2,n2);

for i=1:7

N1(i)=N(i)/d;

end

disp('Высота:');

disp(N1);

bar(Z,N1,1,'g');

for i=1:6

line([Z(i) Z(i+1)],[N1(i) N1(i+1)],'Color','red')

end

title('Гистограмма и полином')

grid on

%эмпирическая функция (групп)

figure(2)

for i=1:7

F1(i)=0;

for j=1:7

if (j<=i)

F1(i)=F1(i)+N(j)/n1;

end

end

end

F1(8)=1;

Z(8)=Z(7)+3;

stairs(Z,F1);

grid on

title('Эмперическая функция(групп.)')

%эмпирическая функция (не групп)

figure(3)

for i=1:n2

F2(i)=0;

B2(i)=B(1,i);

for j=1:n2

if (j<=i)

F2(i)=F2(i)+B(2,j)/n1;

end

end

end

B2(n2+1)=B2(n2)+3;

F2(n2+1)=1;

stairs(B2,F2);

grid on

title('Эмперическая функция(не групп.)')

end