CalcCorner: Calculate Corner Frequencies

In TDD: Time-Domain Deconvolution of Seismometer Response

Description

Inputs a continuous instrument response list and returns a vector of its cutoff frequencies (defined here as the -3dB point). Optionally plots results.

Usage

CalcCorners(PZ, f = 1:1000000/1000, PLOT = FALSE)

Arguments

PZ	Continuous instrument response list (from GetPZ, for example)
f	Vector of frequency values to test (Hz)
PLOT	Logical: plot response spectrum and mark corner frequencies?

Value

Vector of corner frequencies in response spectrum.

Author(s)

Jake Anderson

Examples

# Response of CMG-40T
PZ = GetPZ(12)[[1]]
CalcCorners(PZ, PLOT = TRUE)

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

[1]  0.0335 60.1255

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