PS.VarUntilF: Variance of a powerlaw process if integrated from until...

View source: R/PS.VarUntilF.R

PS.VarUntilFR Documentation

Variance of a powerlaw process if integrated from until frequency f

Description

Integral of PSD=f^(-beta) from f1=1/N to f2=f this equals the variance of a lowpass filtered powerlaw process WARNING: The result is not normalized

Usage

PS.VarUntilF(f, beta, N)

Arguments

f

frequency until which to integrate

beta

powerlaw slope

N

length of the timeseries

Value

non-normalized variance

Author(s)

Thomas Laepple

Examples

beta <- 1
signal <- ts(SimPowerlaw(beta,100000))
spec <- SpecMTM(signal)
v1 <- GetVarFromSpectra(spec,f=c(1/length(signal),0.5))
v2 <- GetVarFromSpectra(spec,f=c(1/length(signal),0.01))
PS.VarUntilF(0.01,beta,length(signal))/PS.VarUntilF(0.5,beta,length(signal))
v2$var/v1$var

EarthSystemDiagnostics/paleospec documentation built on Feb. 17, 2024, 4:36 p.m.