L6p1 SPD

%Se efectueaza o selectie de volum n =20 pt cercetarea unei var aleatoare
%x. val obtinute prin masuratori fiind:
% 0,047;
% 0,228;
% 0,374;
% 0,442;
% 0,594
% 0,625;
% 0,648;
% 0,887;
% 0.933;
% 1,044;
% 1,089;
% 1,176;
% 1,234;
% 1,271;
% 1,351;
% 1,377;
% 1,430;
% 1,529;
% 1,568;
% 1,792;

% Sa se testeze cu nivelul de semnificatie alfa = 0,05 ipoteza ca media
% variabilei aleatoare X este egala cu 1;
% H0; - ipoteza m=1
% H1; -ipoteza alternativa
% Testul bilateral m !=1; testul unilateral (m>1)||(m<1)
%alfa reprezinta probabilitatea de a comite o eroare de tipul 1;-respinge
%ipoteza nula , ea fiind adevarata;
% alte valori pe care alfa le poate lua: 0,1; 0,05; 0,01;
% t= var aleatoare n-1 grade de libertate rep student

%test bilateral privind media repartitiei normale cu dispersie necunoscuta
% Ipoteza nula H0: m=1
% Ipoteza alternativa H1: m!=1

%o ipoteza este acceptata inseamna ca datele expreinmentale nu contrazic
%ipoteza data. nu contrazic nici ca ipoteza ar putea fi egala cu 1 si nici
%faptul ca media ar fi 1.2. Pt ca sa distingem dintre cele 2 ipoteze
%trebuie sa marim volumul esantionului
clc;
clear all;

n=20;
x=[0.047 0.228 0.374 0.442 0.594 0.625 0.648  0.887  0.933 1.044  1.089 1.176 1.234 1.271 1.351 1.377  1.430  1.529  1.568  1.792];
m=1.2;
alfa=0.01;

mx=sum(x)/n

s=sqrt(sum((x-mx).^2)/(n-1))
t=(mx- m)/(s/sqrt(n))
%tinv -returneaza val functiei inverse, functiei de rep student
ta=tinv(1-alfa/2,n-1)

if(abs(t)<=ta)
    disp('se accepta ipoteza nula H0: m=1');
else
    disp('se accepta ipoteza alternativa H1: m!=1');
end;