Inkonsistenz der Berechnung der Kohärenzberechnung im Größenquadrat


7

Ich muss die Magnitude-Square-Kohärenz (MSC) zwischen zwei Signalen berechnen. Bei Verwendung einer Routine, die nur eine Verjüngung (oder überhaupt keine Verjüngung) verwendet, ist mein Ergebnis jedoch immer 1, obwohl die Signale deutlich unterschiedlich sind. Dies passiert nicht, wenn ich mehr als eine Verjüngung verwende. Auf der Suche nach einer Erklärung für dieses abnormale Ergebnis habe ich eine verwirrende Eigenschaft der MSC selbst. Die Definition, die ich verwende, ist diese

γ2(ω)=|X(ω)Y(ω)¯|2(X(ω)X(ω)¯).(Y(ω)Y(ω)¯)

X und Y sind die Fourier-transformierten Signale, die von der Frequenz abhängen . Wenn Sie jedoch zwei komplexe Zahlen als Wert dieser Funktionen in einer festen Frequenz verwenden, ist das Ergebnis immer 1. Wenn Sie wissen, dass istω|z|2=zz¯

γ2=(XY¯)(XY¯)¯XX¯YY¯=(XY¯)(X¯Y)XX¯YY¯=XX¯YY¯XX¯YY¯=1

Natürlich muss es etwas geben, das ich missverstehen muss, aber ich kann nicht sehen, was es ist. Kann mir jemand erklären, was der Haken ist?

Bearbeiten: Ich werde einige Matlab-Links als vertrauenswürdige Quellen verwenden. Definition der MS-Kohärenz

http://www.mathworks.com/help/signal/ref/mscohere.html

Definition der spektralen Dichte der Querleistung

http://www.mathworks.com/help/signal/ref/cpsd.html

(Die Leistungsspektraldichte ist die " Auto- Cross" -Spektraldichte, dh die Fourier-Transformation der Autokorrelation.) Eine wichtige Eigenschaft der Fourier-Transformation der Kreuzkorrelation finden Sie auf Wikipedia unter "Eigenschaften".

Eine weitere Quelle ist das Googeln unter dem Namen "Kohärenzfunktion in der biomedizinischen Signalverarbeitung". Entschuldigung, ich habe die direkten Links hier nicht gepostet, ich habe nicht genug "Ruf"


2
Es wäre hilfreich, wenn Sie einen Verweis darauf einfügen würden, woher Sie Ihre Definition von . Kohärenzmaße basieren im Allgemeinen auf der Cauchy-Schwarz-Ungleichung, die behauptet, dass mit Gleichheit, wenn wobei konstant ist. Gleichheit ist im Kontext der Kohärenz perfekte Kohärenz. Aber Ihr ist nicht ganz das Verhältnis der beiden Seiten der Cauchy-Schwarz-Ungleichung. γ2(ω)
|x,y|2x,xy,y
x=λyλγ2(ω)
Dilip Sarwate

Beachten Sie die Erwartungsoperatoren in Ihrer Definition der Kohärenz. Sie können die und Terme in der Kreuzspektraldichte nicht trennen , um mit dem Nenner aufzuheben, da der erwartete Wert ihres Produkts genommen wird. xy
Jason R

Das stimmt. Ich habe jedoch vergessen zu erwähnen, dass ich mit stationären Zufallsprozessen arbeite, so dass die spektralen Dichten und die Korrelationsfunktionen Fourier-Transformationspaare sind. Und die Fourier-Transformation der Korrelationsfunktionen kann als Multiplikation der Fourier-Transformation der einzelnen Funktionen (mit einer der konjugierten) geschrieben werden, wie in Wikipedia
Tojur

Das Problem ist, dass es für zufällige Prozesse nicht stimmt, dass und die Fourier-Transformationen der Prozesse sind (tatsächlich die Fourier-Transformation eines zufälligen Prozesses, wie Sie es zu denken scheinen macht überhaupt keinen Sinn, außer als komplexwertiger Prozess) und dass die spektrale Kreuzleistungsdichte gleich und die spektrale Leistungsdichte ist . X(ω)Y(ω)SX,Y(ω)X(ω)Y(ω)SX(ω)|X(ω)|2=X(ω)X(ω)
Dilip Sarwate

Ich bin mir nicht sicher, ob ich Ihre Definition von "Zufallsprozess" verstehe. Bei der Signalverarbeitung ist ein stationärer Zufallsprozess "eine Sammlung von Zeitverlaufsaufzeichnungen mit statistischen Eigenschaften, die für Zeitübersetzungen unveränderlich sind" (aus "Zufallsdaten - Analyse- und Messverfahren"). Tatsächlich können Sie überprüfen, wie dieses Konzept verwendet wird, um einige Eigenschaften unter dem folgenden Link
anzugeben:

Antworten:


5

Für jeden einzelnen Datenblock (Fenster) beträgt die Kohärenz, wie Sie beobachtet haben, 1. Um die Kohärenz richtig abzuschätzen, müssen Sie die Spektren und Kreuzspektren für mehrere Fenster mitteln und dann die Kohärenz berechnen.

Die Autospektren XX und YY können auf herkömmliche Weise gemittelt werden. Für das Kreuzspektrum XY müssen Sie zuerst die reellen und imaginären Komponenten mitteln, bevor Sie XY = sqrt (XY [imag] ^ 2 + XY [real] ^ 2) berechnen.

Hilft das? Die Mittelung über 8 Fenster liefert normalerweise zuverlässige Schätzungen.


Ja, das ist die Antwort. Vielen Dank John
Tojur

Hallo John . Eine kleine Frage: Kennen Sie den Namen eines Buches, in dem ich mehr darüber lesen kann? Ich habe keinen guten gefunden, um das alleine zu studieren. Thx
Tojur
Durch die Nutzung unserer Website bestätigen Sie, dass Sie unsere Cookie-Richtlinie und Datenschutzrichtlinie gelesen und verstanden haben.
Licensed under cc by-sa 3.0 with attribution required.