Skip to content
GitLab
Menu
Projects
Groups
Snippets
Loading...
Help
Help
Support
Community forum
Keyboard shortcuts
?
Submit feedback
Sign in
Toggle navigation
Menu
Open sidebar
Seitosh
Seitosh
Commits
66f5cabc
Commit
66f5cabc
authored
Oct 02, 2020
by
thomas.forbriger
Browse files
libs/libpsdxx [TASK]: specify computed values are one-sided PSD
parent
0c6e873b
Changes
3
Hide whitespace changes
Inline
Side-by-side
src/libs/libpsdxx/README
View file @
66f5cabc
...
...
@@ -13,6 +13,12 @@ This library (libpsdxx) provides a C++ class to compute the cross power
spectral density of two time series, the power spectral density of a single
time series and the coherence of two time series.
Power spectral density computed by this library is the so-called one-sided
PSD. This is, if P(f) is the Fourier transform of the normalized
auto-correlation function, the output of the program is 2P(f). The
rms-amplitude in a frequency band (f2-f1) is A=sqrt(2*P(f)*(f2-f1)) if P(f) is
constant in the interval (f1,f2).
Documentation
=============
doxygen source code documentation
...
...
src/libs/libpsdxx/doxygen.txt
View file @
66f5cabc
...
...
@@ -59,6 +59,12 @@ spectral density of two time series, the power spectral density of a single
time series and the coherence of two time series.
\date 01.01.2019
Power spectral density computed by this library is the so-called one-sided
PSD. This is, if P(f) is the Fourier transform of the normalized
auto-correlation function, the output of the program is 2P(f).
The rms-amplitude in a frequency band (f2-f1) is A=sqrt(2*P(f)*(f2-f1))
if P(f) is constant in the interval (f1,f2).
*/
/* ----- END OF doxygen.txt ----- */
src/libs/libpsdxx/dpsdcomputer_processor.cc
View file @
66f5cabc
...
...
@@ -379,6 +379,9 @@ namespace psd {
// we have an energy spectrum so far
// adjust scaling factor to obtain signal power
//
// This factor computes the one-sided PSD, i.e. twice the Fourier
// transform of the normalized auto-correlation function.
double
scalingfactor
=
2.
*
tapscaling
/
Tw
;
if
(
this
->
Mnsegments
>
1
)
{
scalingfactor
/=
this
->
Mnsegments
;
}
...
...
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
.
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment