ReWiA 2.0 html manual

POWER          obliczenie dyskretnej transformaty Fouriera, widma mocy i okna.

Zapis   POWER freqmin freqmax freqstep outsht

Parametry

   freqmin, freqmax, freqstep   
      typu Real.  Określają zakres  częstotliwości dla  którego będzie liczone
      widmo mocy  i transformata.  Wielkość widma  wynikowego będzie wynosić :
      int((freqmax -  freqmin) /  freqstep) +  1. Ponieważ  dane wejściowe  są
      rzeczywiste (bez  części urojonych) amplitudy ujemnych częstotliwości są
      identyczne jak dodatnich, więc można z powodzeniem ustawić freqmin na 0.
      Należy wtedy  pamiętać, że  wszystkie  amplitudy  (oprócz  amplitudy  na
      częstotliwości 0) należy pomnożyć przez 2.

   outsht
      typu Sheet.  Określa widmo  wynikowe.  Domyślnie  widmo  bieżące.  Widmo
      wynikowe będzie  miało 6  pól zawierających  kolejno :  amplitudę, okno,
      część rzeczywistą  amplitudy, część urojoną amplitudy, część rzeczywistą
      okna oraz  część urojoną  okna.  Wskaźniki  będą  ustawione  :  XPTR=-2,
      YPTR=0, EPTR=1, X0=freqmin, X1=freqstep.

Komenda ta wykonuje się długo nawet na szybkich komputerach (parę minut na 486
gdy widmo  wynikowe ma  długość 10000 punktów). W celu przyspieszenia obliczeń
należy używać tak dużego freqstep jak to tylko możliwe.

   Dane wejściowe nie muszą być posortowane ani też ze stałym krokiem. Komenda
POWER zakłada,  że dane  wejściowe są  podane w  okresach. Wyniki  są podane w
jednostkach 1/okres  (komenda POWER  zwraca informację  o zjawiskach krótszych
niż 1 okres). Tak więc, jeśli wprowadzimy czas w dniach to uzyskamy informację
co się dzieje w ciągu 1 dnia. Aby uruchomić przykładowe makro należy napisać :

        M "POWER.MAC"