io.github.kensuke1984.kibrary.datacorrection

クラス SourceTimeFunction

java.lang.Object

io.github.kensuke1984.kibrary.datacorrection.SourceTimeFunction

直系の既知のサブクラス:

SourceTimeFunctionByStackedPeaks

public class SourceTimeFunction
extends java.lang.Object

Source time function.

You have to multiply
Source time function: stf[0], .. stf[NP-1]
on
Waveform in frequency domain: U[1].. U[NP], respectively. See convolve(Complex[])

フィールドのサマリー

フィールド

修飾子とタイプ	フィールドと説明
protected int	np The number of steps in frequency domain.
protected double	samplingHz
protected org.apache.commons.math3.complex.Complex[]	sourceTimeFunction The length is NP
protected double	tlen timeLength [s].

コンストラクタのサマリー

コンストラクタ

修飾子	コンストラクタと説明
protected	SourceTimeFunction(int np, double tlen, double samplingHz)

メソッドのサマリー

修飾子とタイプ	メソッドと説明
static SourceTimeFunction	boxcarSourceTimeFunction(int np, double tlen, double samplingHz, double halfDuration) Boxcar source time function The width is determined by the half duration τ.
protected static boolean	checkValues(int np, double tlen, double samplingHz)
static SourceTimeFunction	computeSourceTimeFunction(int np, double tlen, double samplingHz, double[] obs, double[] syn) Source time function is computed simply by division.
org.apache.commons.math3.complex.Complex[]	convolve(org.apache.commons.math3.complex.Complex[] data) Operates convolution for data in frequency domain.
double[]	convolve(double[] data) Operates convolution for data in time domain.
SACData	convolve(SACData sacData)
int	getNp()
double	getSamplingHz()
org.apache.commons.math3.complex.Complex[]	getSourceTimeFunctionInFrequencyDomain()
Trace	getSourceTimeFunctionInTimeDomain() x axis: time [s], y axis: amplitude After considering that conjugate F[i] = F[N-i],
double	getTlen()
static SourceTimeFunction	readSourceTimeFunction(java.nio.file.Path sourcePath)
static SourceTimeFunction	smoothedRampSourceTimeFunction(int np, double tlen, double samplingHz, double halfDuration) Smoothed ramp source time function The width is determined by the half duration τ.
static SourceTimeFunction	triangleSourceTimeFunction(int np, double tlen, double samplingHz, double halfDuration) Triangle source time function The width is determined by the half duration τ.
void	writeSourceTimeFunction(java.nio.file.Path outPath, java.nio.file.OpenOption... options)

クラスから継承されたメソッド java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, toString, wait, wait, wait

フィールドの詳細

np

protected final int np

The number of steps in frequency domain. It must be a power of 2.

tlen

protected final double tlen

timeLength [s]. It must be a tenth of powers of 2

samplingHz

protected final double samplingHz

sourceTimeFunction

protected org.apache.commons.math3.complex.Complex[] sourceTimeFunction

The length is NP

コンストラクタの詳細

SourceTimeFunction

protected SourceTimeFunction(int np,
                             double tlen,
                             double samplingHz)

パラメータ:

np - must be a power of 2

tlen - [s] must be a tenth of powers of 2

samplingHz - 20 preferred (now must)

メソッドの詳細

triangleSourceTimeFunction

public static SourceTimeFunction triangleSourceTimeFunction(int np,
                                                            double tlen,
                                                            double samplingHz,
                                                            double halfDuration)

Triangle source time function

The width is determined by the half duration τ.
f(t) = 1/τ² t + 1/τ (-τ ≤ t ≤ 0), -1/τ ² t + 1/τ (0 ≤ t ≤ τ), 0 (t < -τ, τ < t)
Source time function F(ω) = (2-2cos(2πωτ)) /(2πωτ)²

パラメータ:

np - the number of steps in frequency domain

tlen - [s] time length

samplingHz - [Hz]

halfDuration - [s] of the source

戻り値:

SourceTimeFunction

boxcarSourceTimeFunction

public static SourceTimeFunction boxcarSourceTimeFunction(int np,
                                                          double tlen,
                                                          double samplingHz,
                                                          double halfDuration)

Boxcar source time function

The width is determined by the half duration τ.
f(t) = 1/(2×τ) (-τ ≤ t ≤ τ), 0 (t < -τ, τ < t)
Source time function F(ω) = sin(2πωτ)/(2πωτ);

パラメータ:

np - the number of steps in frequency domain

tlen - [s] time length

samplingHz - [Hz]

halfDuration - [s] of the source

戻り値:

SourceTimeFunction

smoothedRampSourceTimeFunction

public static SourceTimeFunction smoothedRampSourceTimeFunction(int np,
                                                                double tlen,
                                                                double samplingHz,
                                                                double halfDuration)

Smoothed ramp source time function

The width is determined by the half duration τ.
f(t) = (1-tanh²(2t/τ))/τ (-τ ≤ t ≤ τ), 0 (t < -τ, τ < t)
Source time function F(ω) = (π² ωτ/2)/sinh(π²ωτ/2)

パラメータ:

np - the number of steps in frequency domain

tlen - [s] time length

samplingHz - [Hz]

halfDuration - [s] of the source

戻り値:

SourceTimeFunction

checkValues

protected static boolean checkValues(int np,
                                     double tlen,
                                     double samplingHz)

readSourceTimeFunction

public static SourceTimeFunction readSourceTimeFunction(java.nio.file.Path sourcePath)
                                                 throws java.io.IOException

例外:

java.io.IOException

computeSourceTimeFunction

public static SourceTimeFunction computeSourceTimeFunction(int np,
                                                           double tlen,
                                                           double samplingHz,
                                                           double[] obs,
                                                           double[] syn)

Source time function is computed simply by division.

パラメータ:

obs - waveform of observed

syn - waveform of syn

np - steps of frequency [should be same as synthetics]

tlen - [s] length of waveform [should be same as synthetics]

samplingHz - [Hz]

戻り値:

Source time function F(obs)/F(syn) in frequency domain

getNp

public int getNp()

getTlen

public double getTlen()

getSamplingHz

public double getSamplingHz()

getSourceTimeFunctionInFrequencyDomain

public org.apache.commons.math3.complex.Complex[] getSourceTimeFunctionInFrequencyDomain()

戻り値:

source time function in frequency domain. the length is np

writeSourceTimeFunction

public void writeSourceTimeFunction(java.nio.file.Path outPath,
                                    java.nio.file.OpenOption... options)
                             throws java.io.IOException

パラメータ:

outPath - Path for a file.

options - for writing the file

例外:

java.io.IOException - if the source time function is not computed, then an error occurs

convolve

public final double[] convolve(double[] data)

Operates convolution for data in time domain.

パラメータ:

data - to be convolved in time domain. The data is convolved after FFTed.

戻り値:

convolute data in time domain

convolve

public final org.apache.commons.math3.complex.Complex[] convolve(org.apache.commons.math3.complex.Complex[] data)

Operates convolution for data in frequency domain.

パラメータ:

data - to be convolved in frequency domain. The length must be np + 1

戻り値:

convolute data in frequency domain

convolve

public final SACData convolve(SACData sacData)

パラメータ:

sacData - to convolute with this.

戻り値:

convoluted SACData

getSourceTimeFunctionInTimeDomain

public Trace getSourceTimeFunctionInTimeDomain()

x axis: time [s], y axis: amplitude After considering that conjugate F[i] = F[N-i],

戻り値:

trace of Source time function in time domain