NAME: s_trans PURPOSE: Calculate the local power spectrum for given time-series data and return the structure which contains {st:, ph: freq:, time:} where st is the S Transform. SYNTAX: results = s_trans(time-series) KEYWOARDS: \help explains all the keywords and parameters \verbose flags errors and size \samplingrate if set returns array of frequency \maxfreq \minfreq \freqsmplingrate \piecewisenumber divides the time series, and passes back array \power returns the power spectrum \abs returns the absolute value spectrum \Rremoveedge removes the edge with a 5% taper, and takes out least-sqares fit parabola Added a "mask out edges" keyword, which STs a line (ie st of edges) and thresholds the returned st matrix. The value of masked edges is the percent at which to make the threshold default = 5 %. Added an EXAMPLE keyword, will display a time series and the amplitude of the ST to the current graphics device WARNING, will call !P.multi=0 ===================== Modified hilbert transform ================================ ;;;;;;;;;;;;;;;;;;;; + NAME: hilbert PURPOSE: Return a series that has all periodic terms shifted by 90 degrees. CODE: A. Shinbori, 04/06/2012. MODIFICATIONS: ACKNOWLEDGEMENT: $LastChangedBy: jwl $ $LastChangedDate: 2014-01-22 15:54:40 -0800 (Wed, 22 Jan 2014) $ $LastChangedRevision: 13976 $ $URL $