VF_coherenceVD_coherenceVE_coherence
FunktionKohärenzfunktion zweier Signale
Syntax C/C++#include <VFstd.h>
float VF_coherence( fVector Coher, ui specsiz, fVector X, fVector Y, ui xsiz, fVector Win );
C++ VecObj#include <OptiVec.h>
T vector<T>::coherence( const vector<T>& X, const vector<T>& Y, const vector<T>& Win );
Pascal/Delphiuses VFstd;
function VF_coherence( Coher:fVector; specsiz:UIntSize; X,Y:fVector; xsiz:UIntSize; Win:fVector ): Single;
CUDA-Funktion C/C++#include <cudaVFstd.h>
int cudaVF_coherence( float *h_coherfNyq, fVector d_Coher, ui specsiz, fVector d_X, fVector d_Y, ui xsiz, fVector d_Win );
int cusdVF_coherence( float *d_coherfNyq, fVector d_Coher, ui specsiz, fVector d_X, fVector d_Y, ui xsiz, fVector d_Win );
float VFcu_coherence( fVector h_Coher, ui specsiz, fVector h_X, fVector h_Y, ui xsiz, fVector h_Win );
CUDA-Funktion Pascal/Delphiuses VFstd;
function cudaVF_coherence( var h_coherfNyq:Single; d_Coher:fVector; specsiz:UIntSize; d_X, d_Y:fVector; xsiz:UIntSize; d_Win:fVector ): IntBool;
function cusdVF_coherence( d_coherfNyq:PSingle; d_Coher:fVector; specsiz:UIntSize; d_X, d_Y:fVector; xsiz:UIntSize; d_Win:fVector ): IntBool;
function VFcu_coherence( h_Coher:fVector; specsiz:UIntSize; h_X, h_Y:fVector; xsiz:UIntSize; h_Win:fVector ): Single;
BeschreibungDie Kohärenzfunktion der beiden Datensätze X und Y wird ermittelt. Sie ist definiert als das Betragsquadrat des Kreuzleistungsspektrums von X und Y dividiert durch das Produkt der beiden Autoleistungsspektren von X und Y:
CXY = SXY2 / (SXX * SYY).
Da das Ergebnis für negative und positive Frequenzen identisch ist, wird hier das Ergebnis nur für f=0 und positive Frequenzen gespeichert. Um specsiz als ganzzahlige Potenz von 2 definieren zu können, werden hier jedoch nur specsiz Punkte in Spc zurückgegeben. Der letzte Wert, d.h. die Kohärenzfunktion für die Nyquist-Frequenz, wird als Rückgabewert der Funktion behandelt und kann entweder vernachlässigt (durch Aufruf der Funktion wie eine void-Funktion bzw. eine procedure) oder als letztes Element in Spc gespeichert werden durch Aufruf in der Form
Coher[specsiz] = VF_coherence( Coher, specsiz, X, X, xsiz, Win );
In diesem Fall muss Coher eine Länge von specsiz+1 haben.

xsiz muss mindestens 2*specsiz betragen und specsiz eine ganzzahlige Potenz von 2 sein. Intern wird X in xsiz / specsiz - 1 Segmente unterteilt und über die für die einzelnen Segmente berechneten Kreuz- sowie Autoleistungsspektren gemittelt, bevor aus den Mittelwerten die Kohärenzfunktion bestimmt wird.

Win ist ein Fenster, das auf die Datensegmente angewendet wird. Seine Größe ist immer 2*specsiz. Innerhalb der VectorLib-Bibliotheken sind drei Funktionen vorhanden, die passende Fenster erzeugen: VF_Welch,   VF_Parzen und VF_Hann. Ein Rechteck-Fenster wird über VF_equ1 erhalten. Von seinem Gebrauch wird hier aber abgeraten.

FehlerbehandlungWenn size nicht eine Potenz von 2 ist, meldet sich VF_FFT (worauf VF_coherence basiert) mit der Fehlermeldung "Size must be integer power of 2." und bricht das Programm ab.
RückgabewertKohärenzfunktion bei der Nyquist-Frequenz
QuerverweisVF_FFT,   VF_spectrum,   VF_xspectrum,   VF_convolve,   VF_autocorr,   VF_xcorr,   VF_filter

VectorLib Inhaltsverzeichnis  OptiVec Home