<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 3.2 Final//EN"> <!--Converted with LaTeX2HTML 98.1p1 release (March 2nd, 1998) originally by Nikos Drakos (nikos@cbl.leeds.ac.uk), CBLU, University of Leeds * revised and updated by: Marcus Hennecke, Ross Moore, Herb Swan * with significant contributions from: Jens Lippmann, Marek Rouchal, Martin Wilck and others --> <HTML> <HEAD> <TITLE>Period analysis</TITLE> <META NAME="description" CONTENT="Period analysis"> <META NAME="keywords" CONTENT="vol2"> <META NAME="resource-type" CONTENT="document"> <META NAME="distribution" CONTENT="global"> <META HTTP-EQUIV="Content-Type" CONTENT="text/html; charset=iso-8859-1"> <LINK REL="STYLESHEET" HREF="vol2.css"> <LINK REL="next" HREF="node247.html"> <LINK REL="previous" HREF="node245.html"> <LINK REL="up" HREF="node245.html"> <LINK REL="next" HREF="node247.html"> </HEAD> <BODY > <!--Navigation Panel--> <A NAME="tex2html4462" HREF="node247.html"> <IMG WIDTH="37" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="next" SRC="icons.gif/next_motif.gif"></A> <A NAME="tex2html4459" HREF="node245.html"> <IMG WIDTH="26" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="up" SRC="icons.gif/up_motif.gif"></A> <A NAME="tex2html4453" HREF="node245.html"> <IMG WIDTH="63" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="previous" SRC="icons.gif/previous_motif.gif"></A> <A NAME="tex2html4461" HREF="node1.html"> <IMG WIDTH="65" HEIGHT="24" ALIGN="BOTTOM" BORDER="0" ALT="contents" SRC="icons.gif/contents_motif.gif"></A> <BR> <B> Next:</B> <A NAME="tex2html4463" HREF="node247.html">Comparison of two stochastic</A> <B> Up:</B> <A NAME="tex2html4460" HREF="node245.html">Examples</A> <B> Previous:</B> <A NAME="tex2html4454" HREF="node245.html">Examples</A> <BR> <BR> <!--End of Navigation Panel--> <H2><A NAME="SECTION001761000000000000000"> Period analysis</A> </H2> <P> Let us assume that table <TT>OBSERVP.tbl</TT> contains observations of a periodic phenomenon. The observations are stored in the <TT>DOUBLE PRECISION</TT> columns <TT>:TIME</TT> and <TT>:VALUE</TT>. <P> <PRE><TT> <TT><TT>CREATE/GRAPHICS ! Create graphics window <BR><TT>SET/CONTEXT TSA ! Enable TSA package <BR><TT>NORMALIZE/TSA OBSERVP :VALUE ! Subtract mean <BR><TT>BAND/TSA OBSERVP ! Find suitable frequency band <BR><TT>SHOW/TSA ! Inspect corresponding settings <BR><TT>POWER/TSA OBSERVP POWERSPEC ! Compute power spectrum and <BR><TT> ! spectral window <BR><TT>PLOT/ROW POWERSPEC ! Display power spectrum <BR><TT>PLOT/ROW POWERSPEC 2 ! Display spectral window <BR><TT>ORT/TSA ? AOVSPEC ? ? ? ! Compute multiharmonic spectrum <BR><TT>PLOT/ROW AOVSPEC ! Inspect results <BR><TT>POWER/TSA ? LINESPEC 20 .01 501 ! Compute detail of power spec. <BR><TT>PLOT/ROW LINESPEC <BR> ! Plot it <BR><TT>WIDTH/TSA LINESPEC ! Find parameters of <BR><TT> ! the oscillation <BR><TT>SINEFIT/TSA OBSERVP CLNOBS 23 1 ! Remove one particular <BR><TT> ! Fourier component <BR><TT>POWER/TSA CLNOBS POWCLN ! Inspect periodogram of data <BR><TT> ! after removal of 1 oscillation <BR> </TT></PRE></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT></TT> <P> <BR><HR> <ADDRESS> <I>Petra Nass</I> <BR><I>1999-06-15</I> </ADDRESS> </BODY> </HTML>