Module to find the norm of a signal and to return a normalised array. The original signal is left untouched. More...
Prototypes | |
| void | LALInspiralWaveNormalise (LALStatus *status, REAL4Vector *in, REAL8 *norm, REAL8Vector psd) |
Module to find the norm of a signal and to return a normalised array. The original signal is left untouched.
Given the positive frequency Fourier components \(H_k,\) \(k=0,\ldots,n-1,\) of a vector and the noise PSD \(S_m,\) \(m=0,\ldots,n/2,\) this module first computes the norm \(H\) of the vector treating \(S_m\) as the measure: (note that in {fftw} notation, the zeroth frequency component is \(H_0,\) Nyquist is \(H_{n/2},\) \(H_k,\) \(k \ne 0,n/2,\) ( \(H_{n-k})\) is the real (imaginary) part of the \(k\)th harmonic)
\begin{equation} \label{eq_inspiralnorm} H = \sum_{k=1}^{n/2-1} \frac{H_k^2 + H^2_{n-k}}{S_k}. \end{equation}
(Note that the zeroth and Nyquist components are ignored in the computation of the norm.) It then replaces the original vector \(H_k\) with normalized vector using:
\begin{equation} \widehat H_k = \frac {H_k}{\sqrt H},\ \ k=0,\ldots n-1. \end{equation}
Definition in file LALInspiralWaveNormalise.c.
Go to the source code of this file.
| void LALInspiralWaveNormalise | ( | LALStatus * | status, |
| REAL4Vector * | in, | ||
| REAL8 * | norm, | ||
| REAL8Vector | psd | ||
| ) |
Definition at line 70 of file LALInspiralWaveNormalise.c.