Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://www.sao.ru/precise/Midas_doc/doc/94NOV/vol2/node258.html
Дата изменения: Fri Feb 23 14:02:13 1996 Дата индексирования: Tue Oct 2 18:41:37 2012 Кодировка: Поисковые слова: с р р с с п п р п п с с с р р р р р п п р р р п п р |
Let us assume that table OBSERVP.tbl contains observations of a periodic phenomenon. The observations are stored in the DOUBLE PRECISION columns :TIME and :VALUE.
CREATE/GRAPHICS ! Create graphics windowSET/CONTEXT TSA ! Enable TSA package
NORMALIZE/TSA OBSERVP :VALUE ! Subtract mean
BAND/TSA OBSERVP ! Find suitable frequency band
SHOW/TSA ! Inspect corresponding settings
POWER/TSA OBSERVP POWERSPEC ! Compute power spectrum and
! spectral window
PLOT/ROW POWERSPEC ! Display power spectrum
PLOT/ROW POWERSPEC 2 ! Display spectral window
AOV/TSA ? AOVSPEC ? ? ? ! Compute AOV test statistics
PLOT/ROW AOVSPEC ! Inspect results
POWER/TSA ? LINESPEC 20 .01 501 ! Compute detail of power spec.
PLOT/ROW LINESPEC ! Plot it
WIDTH/TSA LINESPEC ! Find parameters of
! the oscillation
SINEFIT/TSA OBSERVP CLNOBS 23 1 ! Remove one particular
! Fourier component
POWER/TSA CLNOBS POWCLN ! Inspect periodogram of data
! after removal of 1 oscillation