131

Analiza statystyczna pomiarów 131

function [mUab,sUab,mUbc,sUbc,mUca,sUca,Un]=ustat(pomiary)

% funkcja [mUab,sUab,mUbc,sUbc,mUca,sUca,Un]=uastat(pomiary) % - czyta pomiary napiec z opliku upomiary.m % - wykreśla przebiegi czasowe napięcia % - oblicza średnie i odchylenia standardowe if nargin < 1

pomiary='upomiary';% domyślna nazwa pliku z pomiarami

end

[Uab,Ubc,Uca,Un,t]=feval(pomiary); x= [ ] ;

x= [Uab Ubc Uca] ;

% Analiza pomiarów napięcia Uab plot(t,Uab);grid on;

title(' Przebieg czasowy zmian napięcia miedzyfazowego') xlabel(' czas, min'); ylabel(' Uab, pu' ) ; fprintf('\n\n Dalej? Nacisnij ENTER'); pause;

% otwarcie pliku na wyniki fd=fopen('ustatout.m','wt¹);

fprintf(fd,    '\n    Wartość średnia    Uab_sr    =    %5.3f,

fprintf(fd,    '\n    Odch. standard.    Uab_s    =    %5.3f,

fprintf]    '\n    Wartość średnia    Uab_sr    =    %5.3f,

fprintf(    ¹\n    Odch. standard.    Uab_s    =    %5.3f,

hist(Uab);

title(' Histogram napięcia Uab¹) ; xlabel('Napięcie Uab,pu'); ylabel(¹Czestosc występowania'); grid on;

% dane o przedziałach [n,usr]=hist(Uab); np=length(n); suman=sum(n) fprintf(fd,...

'\n Nr. przedziału Środek fprintf(...

' \n Nr. przedziału Środek for i=l:np

fprintf(fd, '\n %8d %15.8f %15d %8.3f¹ , i, usr(i),n(i),npu(i)); fprintf('\n %8d %15.8f %15d %8.3f',i,usr(i),n(i),npu(i)); end

fprintf('\n\n Dalej? Nacisnij ENTER'); pause; m=mUab; s=sUab; xmin=m-3*s; xmax=m+3*s; d=6*s;

npu=n./suman; przedziału Liczebność przedziału Liczebność

pu

pu

pu

pu

,mUab); ,sUab); ,mUab); ,sUab);

Czestosc ' )

Czestosc ' )

xx=xmin:0.01*d:xmax;

yy=rngauss(xx,m,s);

plot(xx,yy,'r-',usr,npu,'b:');

legend('krzywa Gaussa','częstości z histogramu'); grid on;

xlabel('Uab'); ylabel('Czestosc ' ) ;

title('Rozkład normalny i Częstości z histogramu');