Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet Schmucker-Weidelt Lecture Notes, Aarhus, 1975 - MTNet

mtnet.dias.ie
from mtnet.dias.ie More from this publisher
04.08.2013 Views

- - . the same or from different sites. Let Z(w>, X(w), Y(w) be the Fouri-er transforms of Z(t), XCt) and Y(t). Then a linear relation of the form - - is established in which A and B represent the desired transfer func- - tions between Z on the one side and X and Y on the other side; 6Z is the uncorrelated "noise" in 2, assuming X and Y to be noise-free. As the best fitting transfer - functions will be considered those whicf produce minLmum noise < 16zI2 > in the statistical average. Here the average is to be taken either over a number of records -- or within extended frequency bands of the width which is L times greater than the ultimate spacing IL'T of individual spectral estimates. The noise: signal ratio defines the residual e(w) ,%atio of related to observed signal the . --- coherence R(w): The coherence in conjunction with the degree of freedom of the - - functions A and B. averaging procedure estab1ishe.s confidence limits for -the transfer The averaged products of Fourier transforms are denoted as S ZZ ., - = < Z Z * >: power spectrum of Z. - - S = < Z Y * >: cross spectrum between Z and Y ZY wi-th S = S* . ZY Y = In summary, the data reduction involves the following steps (a) Fourier transformation of tiine records (b) Calculations of power and cross-spectra (c) Calculation of transfer functions (d) )COlculation of confidence limits for the trans- ferfunctions. Steps Ca) and (bf can be substi.cuted by the foll.owing alternatives: (ax) : Calculate auto-correlation functions ) , .. . and cross- correlation f.urictf.ons R , . . . with T being a time lag, ZY

with 0 < T < 'rmax. R (T) = / Z(+-t) Z(t)dt Z z T R (T) = / Z(+-t) YCt)dt ZY T x (b ): Take the Fourier transforms of the correlation functions and obtain as in (b) power- and cross-spectral estimates, if the averaging is done within frequency bands of the width m. There is a formal correspondence between the -- -1. maximum lag and Af . The actual performance of the steps (a) to (d) with one or more setsof records is now described in detail. -- a. Fourier transformation. - The time series and their spectra are given, respectively to be found, at discrete values of t and w, which are equally spaced in time and frequency. Let Z e-tc. be an n .instantenuous value of ZCt) for t = tn, n = 0,1,2, . . . N with At = tntl - t . Outside of the record, extending from t to t n o N Z(t) is assumed to be zero. - Let Z be the Fourier transform of Z(t) for the frequency f = f . n1 IT1 Because of the finite length T of the rec6rd the lowest resolvable frequency and thereby the frequency spacing w i l l be given by T-I = Af = f and f has to be a multiple of Af. Because Z(t) is 1 m gj-ven at discrete instances of -time, A t apart, the hi.ghest resol- -, vable frequency, called the Ny quist frequenz, w i l l be the reci- procal of (At-2), i.e. and - I fH - - = MAf 2At

with 0 < T < 'rmax.<br />

R (T) = / Z(+-t) Z(t)dt<br />

Z z T<br />

R (T) = / Z(+-t) YCt)dt<br />

ZY T<br />

x<br />

(b ): Take the Fourier transforms of the correlation functions<br />

and obtain as in (b) power- and cross-spectral estimates,<br />

if the averaging is done within frequency bands of the<br />

width m. There is a formal correspondence between the<br />

-- -1.<br />

maximum lag and Af .<br />

The actual performance of the steps (a) to (d) with one or more<br />

setsof records is now described in detail.<br />

-- a. Fourier transformation. - The time series and their spectra are<br />

given, respectively to be found, at discrete values of t and w,<br />

which are equally spaced in time and frequency. Let Z e-tc. be an<br />

n<br />

.instantenuous value of ZCt) for t = tn, n = 0,1,2, . . . N with<br />

At = tntl - t . Outside of the record, extending from t to t<br />

n o N<br />

Z(t) is assumed to be zero.<br />

-<br />

Let Z be the Fourier transform of Z(t) for the frequency f = f .<br />

n1 IT1<br />

Because of the finite length T of the rec6rd the lowest resolvable<br />

frequency and thereby the frequency spacing w i l l be given by<br />

T-I = Af = f and f has to be a multiple of Af. Because Z(t) is<br />

1 m<br />

gj-ven at discrete instances of -time, A t apart, the hi.ghest resol-<br />

-,<br />

vable frequency, called the Ny quist frequenz, w i l l be the reci-<br />

procal of (At-2), i.e.<br />

and<br />

- I<br />

fH - - = MAf<br />

2At

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!