PZ2Resp: Calculate Instrument Response from Poles/Zeros

In TDD: Time-Domain Deconvolution of Seismometer Response

Description

Returns complex instrument response for a vector of frequencies and set of analog poles and zeros. Optionally plots the magnitude of the instrument response.

Usage

PZ2Resp(PZ, f, PLOT = TRUE)

Arguments

PZ	Analog poles and zeros
f	Vector of frequencies for which response is calculated
PLOT	Logical: whether to plot magnitude of response

Details

The response is calculated by the following equation: R(s) = PZ$Sense * PZ$Knorm * prod(s - PZ$zeros)/prod(s - PZ$poles)

where s = 2 * pi * 1i * f.

PZ requires the following elements: poles: Vector of poles np: number of poles zeros: Vector of zeros nz: Number of zeros Knorm: Normalization constant Sense: Instrument sensitivity (V/(m/s))

Value

Vector of instrument response values corresponding to the frequencies in f.

Author(s)

Jake Anderson

Examples

# Response of Guralp CMG-40T

PZ = GetPZ(12)[[1]]

f = (1:10000 - 1)/1000

PZ2Resp(PZ, f)

Example output

Loading required package: signal

Attaching package: 'signal'

The following objects are masked from 'package:stats':

    filter, poly

Loading required package: RSEIS
Loading required package: pracma

Attaching package: 'pracma'

The following objects are masked from 'package:RSEIS':

    detrend, hypot, logspace, peaks, trapz

The following objects are masked from 'package:signal':

    conv, ifft, interp1, pchip, polyval, roots


Attaching package: 'TDD'

The following object is masked from 'package:RSEIS':

    ReadInstr

TDD documentation built on May 2, 2019, 4:51 a.m.