No, Not a full spectrum.

I mean a subroutine I might use within my own program
to input a single frequency then test the samples
given for the presence of that single frequency ?

I am sure there exists something like a
FFT which will only look for the
energy contained within a single
given frequency and not like powers of two (512) all at once.

The only routine I have found is one where you
plug in 1024 samples and get out of it
512 different discrete frequencies.

I am not a math expert, so I need a public domain
subroutine which anyone is free to use at will.

thanks for your answers,
geoff

-----Original Message----- 
From: Gary Lindgren 
Sent: Saturday, April 23, 2011 5:09 AM 
To: psnlist@.............. 
Subject: RE: Modulated Seismic 

Sure Geoffrey,
It can be done.
Just select the time period in WinSDR and go to Winquake and select
Calculate/FFT/Full and you will get a display of response vs frequency. 
It can be done.
Gary




-----Original Message-----
From: psnlist-request@.............. [mailto:psnlist-request@...............
On Behalf Of Geoffrey
Sent: Friday, April 22, 2011 8:17 PM
To: psnlist@..............
Subject: Re: Modulated Seismic

Is there some way to pick a frequency
like 0.123 Hz then process the sampled
seismic signal to see if any energy of only
that frequency is present ?

Can you point me to a complete subroutine
[complete meaning does not use a library of functions]
to obtain such an answer ?

It would be like an FFT but you give the data,
sample rate, and what SINGLE freq to look for ?


-----Original Message----- 
From: Chuck / Judy Burch 
Sent: Sunday, April 17, 2011 3:44 PM 
To: psnlist@.............. 
Subject: Modulated Seismic 


All,

Thanks to Randal for raising this interesting topic.  The Teager-Kaiser 
Algorithm and Randal's rectification scheme are efforts to find the 
instantaneous amplitude (or envelop) of a seismic record.  Think of a 
short enough piece of a seismic record -- it can be modeled by a 
waveform of a given frequency, phase and amplitude.  Throw away the 
frequency and phase information and you are left with just the amplitude 
as a function of time.  So you can think of a record (very 
approximately) as a collection of unit-amplitude waveforms from an 
appropriate frequency range that has been AM modulated by the envelop.

The FFT of the envelop, then, can be thought of as the Fourier Transform 
of the "modulating" function.

There are different ways of approximating the envelop.  Many texts on 
signal processing discuss the "analytic trace" or "analytic signal" 
method, which is mathematically rigorous.  If the original record is X, 
then the quadrature (90 degree advanced) record, Y, is the Hilbert 
Transform of X.  The envelop of X is SQRT ( X2 + Y2 ).

In the (remote) chance that anyone wants to pursue this, contact me and 
I'll show you an easy way to calculate the Hilbert Transform.


Chuck Burch
__________________________________________________________

Public Seismic Network Mailing List (PSNLIST)

To leave this list email PSNLIST-REQUEST@.............. with 
the body of the message (first line only): unsubscribe
See http://www.seismicnet.com/maillist.html for more information.
__________________________________________________________

Public Seismic Network Mailing List (PSNLIST)

To leave this list email PSNLIST-REQUEST@.............. with 
the body of the message (first line only): unsubscribe
See http://www.seismicnet.com/maillist.html for more information.

__________________________________________________________

Public Seismic Network Mailing List (PSNLIST)

To leave this list email PSNLIST-REQUEST@.............. with 
the body of the message (first line only): unsubscribe
See http://www.seismicnet.com/maillist.html for more information.
__________________________________________________________

Public Seismic Network Mailing List (PSNLIST)

To leave this list email PSNLIST-REQUEST@.............. with 
the body of the message (first line only): unsubscribe
See http://www.seismicnet.com/maillist.html for more information.