clear all
close all
clear clc

%cwiczenie1

M=8;
A3=1.099;

N=3;
mi=10;
sigma=1;

Samples=sigma*randn(N,M)+mi;


subplot(3,1,1)

for i=1:1:N
    plot(Samples(i,:),'.');
    xlabel('numer próbki')
    ylabel('wartość')
    hold on
end



Average=zeros(M);

for m=1:1:M
    for n=1:1:N
        Average(m)=Average(m)+Samples(n,m);
    end
end

for m=1:1:M
        Average(m)=Average(m)/N;
end

subplot(3,1,2)
plot(Average(:,1),'.')
xlabel('numer próbki')
ylabel('wartość średnia')

Standard=zeros(N,M);
Dev=zeros(M);
for m=1:1:M
    for n=1:1:N
        Standard(n,m)=Samples(n,m)-Average(m);
    end
end

for m=1:1:M
    for n=1:1:N
        Dev(m)=Dev(m)+(Standard(n,m)*Standard(n,m));
    end
end

for m=1:1:M
        Dev(m)=sqrt(Dev(m));
end

subplot(3,1,3)
plot(Dev(:,1),'.')
xlabel('numer próbki')
ylabel('wartość odchylenia')

figure(2)
%ćwiczenie 2
boxplot(Samples)

%ćwiczenie 3
AveAverage=0;
for a=1:1:M
AveAverage=AveAverage+Average(a);
end
AveAverage=AveAverage/M;

S=mean(std(Samples));

GLK=AveAverage+A3*S;
DLK=AveAverage-A3*S;

%ćwiczenie 4

B3=0;
B4=2.568;

DLK2=B3*S;
GLK2=B4*S;



figure(3)
subplot(3,1,1)

for i=1:1:N
    plot(Samples(i,:),'.');
    xlabel('numer próbki')
    ylabel('wartość')
    hold on
end



Average=zeros(M);

for m=1:1:M
    for n=1:1:N
        Average(m)=Average(m)+Samples(n,m);
    end
end

for m=1:1:M
        Average(m)=Average(m)/N;
end

subplot(3,1,2)
plot(Average(:,1),'.')
xlabel('numer próbki')
ylabel('wartość średnia')
refline([0 GLK])
refline([0 DLK])

Standard=zeros(N,M);
Dev=zeros(M);
for m=1:1:M
    for n=1:1:N
        Standard(n,m)=Samples(n,m)-Average(m);
    end
end

for m=1:1:M
    for n=1:1:N
        Dev(m)=Dev(m)+(Standard(n,m)*Standard(n,m));
    end
end

for m=1:1:M
        Dev(m)=sqrt(Dev(m));
end

subplot(3,1,3)
plot(Dev(:,1),'.')
xlabel('numer próbki')
ylabel('wartość odchylenia')
refline([0 GLK2])
refline([0 DLK2])