sacpz2resp.py - Convert a Sac Pole Zero file into a Amplitude/Phase Response
pz2resp.py - Convert a Sac Pole Zero file into a Resp formatted file
sacpz2resp.py --help
Option -f filename
Sac Pole Zero file should be in an output format as output from rdseed
Keys:
ZEROS nzeros
zero [ real, imaginary ]
...
POLES npoles
pole [ real, imaginary ]
...
CONSTANT constant value
Option -u units
dis - Displacement [ Zeros are not changed ]
vel - Velocity [ One Zero is removed ]
acc - Acceleration [ Two Zeros are removed ]
Output on stdout
Frequency (Hz) Amplitude (meters) Phase (radians)
Internals:
freq = frequency Hz
z_i Zeros
p_i Poles
min_freq, max_freq Minimum and Maximum Frequencies
npts Number of Points to Evaluate
delta Spacing in seconds of the time series
(max_freq - min_freq)
delta_freq = ------------------------
(npts * delta)
w = 2 * pi * f Angular Frequency, Radians
(w-z_0) (w-z_1) ... (w-z_n)
R(w) = C * -----------------------------
(w-p_0) (w-p_1) ... (w-p_n)