help crcv4
 ==============================================================
                  SPEKTRALNI CHARAKTERISTIKY
              IDENTIFIKACE PERIODICKYCH KOMPONENT
                       TESTY PERIODICITY
 ==============================================================

colordef none
crcv4
format compact
us_deaths=[crpath,'/brockwel.dat/deaths.dat'];
[x,xinf]=getdata(us_deaths);   % nacteni a vykresleni
n=length(x);
pause;
%---------------------------------------------------------------
% Nalezneme fourierovske spektrum
  n2=fix(n./2);
  X=fft(x)./n;   % Diskretni spektrum komplexnich Fourierovych koeficientu
  mean(x)-X(1)   % X(1) = prumer
ans =
     0
  Xs=X(2:n);     % Zbytek spektra=Harmonicke komponenty
  Xs(n2-5:n2+5)  % Jsou komplexne sdruzene kolem stredu
ans =
  -2.8567 -35.9688i
  -5.0808 - 2.0484i
  13.2940 + 9.0474i
   9.9525 -25.7728i
  -2.0399 - 6.6787i
  21.3194          
  -2.0399 + 6.6787i
   9.9525 +25.7728i
  13.2940 - 9.0474i
  -5.0808 + 2.0484i
  -2.8567 +35.9688i
pause;
%---------------------------------------------------------------
  figure;
  stem(2*abs(Xs))    % Graf amplitud: stredova symetrie je tez zretelna
  title('AMPLITUDOVE SPEKTRUM');
pause;
%---------------------------------------------------------------
  figure;
  stem(180*atan2(imag(Xs),real(Xs))./pi) % Graf fazi: stredova antisymetrie
  title('FAZOVE SPEKTRUM');
  ylabel('stupen');
pause;
%---------------------------------------------------------------
% Nalezneme vyznamne periodicke komponenty pomoci periodogramu a
% Fisherova-Siegelova testu:
  [pdg,kF,kS]=tperiod(x,figure);  % periodogram pdg a jeho graf
  kF        % indexy vyznamnych harmonickych slozek uzitim Fisherova testu
kF =
     6     1    12    30     3    18     2    24
  kS        % indexy vyznamnych harmonickych slozek uzitim Siegelova testu
kS =
     6     1    12    30     3    18     2    24
% o periodach:
  if length(kS)>length(kF)
  else
     n./kF
ans =
  Columns 1 through 6 
   12.0000   72.0000    6.0000    2.4000   24.0000    4.0000
  Columns 7 through 8 
   36.0000    3.0000
  end
% Nejvyraznejsi slozka kF(1)=kS(1)=6 (6-ta harmonicka) indikuje periodu
% 12 mesicu sezonni slozky.
%---------------------------------------------------------------------------
% PROVEDTE 'type tperiod' A PROSTUDUJTE ALGORITMUS PROCEDURY 'tperiod'.
%================================== KONEC ==================================
echo off
diarz off
??? Undefined command/function 'diarz'.

diary off