Wie kann man ein Zeitdomänensignal "aufhellen"?

12

Ich versuche zu verstehen, wie man einen so genannten "Pre-Whitening" -Filter oder einfach einen "Whitening" -Filter genau umsetzt.

Ich verstehe, dass der Zweck darin besteht, ein Delta als Autokorrelationsfunktion zu verwenden, aber ich bin mir nicht sicher, wie ich das genau machen soll.

Der Kontext ist hier der folgende: Ein Signal wird an zwei verschiedenen Empfängern empfangen und deren Kreuzkorrelation wird berechnet. Die Kreuzkorrelation kann wie ein Dreieck oder eine andere gottverlassene Form aussehen. Aufgrund dessen wird es schwierig, die Spitze des Kreuzkorrelationssignals zu finden. In diesem Fall muss ich die Signale "aufhellen", bevor eine Kreuzkorrelation durchgeführt wird, sodass die Kreuzkorrelation nun mehr delta-artig ist.

Wie geht das?

Vielen Dank!

autocorrelation

— Spacey
quelle

Beachten Sie, dass Ihre als Weißmacher bezeichnete Frage im Kontext von Kommunikationssystemen im Wesentlichen die Funktion eines Equalizers erfüllt. Das klingt für mich genauso. es könnte einfach eine andere Nomenklatur sein.

— Jason R

Ja, die schlecht definierte Nomenklatur macht es umso verwirrender, was sie manchmal zu tun versuchen.

— Spacey

7

Angenommen, Sie haben die Signale und deren Kreuzkorrelationsfunktion Ihnen nicht gefällt. Sie möchten, dass impulsartig ist. Man beachte , daß in der Frequenzdomäne, $x(t)$ $y(t)$ $R_{x,y}(t)$ $R_{x,y}$

F [R_{x, y}] = S_{x, y} (f) = X (f) {Y.}^{*} (f) .

$\mathcal{F}[R_{x,y}] = S_{x,y}(f) = X(f)Y^*(f).$ So , um die Signale , die durch lineare Filter filtert ,

und

jeweils erhalten

,

, und

,

, und jetzt ist ihre Kreuzkorrelationsfunktion

g

$g$

h

$h$

\hat{x} (t) = x * g

$\hat{x}(t) = x*g$

\hat{X} (f) = X (f) G (f)

$\hat{X}(f) = X(f)G(f)$

\hat{y} = y * h

$\hat{y} = y*h$

\hat{Y} (f) = Y (f) H (f)

$\hat{Y}(f) = Y(f)H(f)$

, dessen FourierTransformation

R_{\hat{x}, \hat{y}}

$R_{\hat{x},\hat{y}}$

das

die Kreuzkorrelation ist

mit

. Noch wichtiger ist, dass Sie

und

sowählen, dass dieKreuzspektraldichte

\begin{aligned} F [R_{\hat{x}, \hat{y}}] = S_{\hat{x}, \hat{y}} (f) & = [X (f) G (f)] [Y. (f) H (f)]^{*} \\ = [X (f) {Y.}^{*} (f)] [G (f) H^{*} (f)] \\ = [X (f) {Y.}^{*} (f)] [G^{*} (f) H (f)]^{*}, \end{aligned}

$\begin{align*} \mathcal{F}[R_{\hat{x},\hat{y}}] = S_{\hat{x},\hat{y}}(f) &= [X(f)G(f)][Y(f)H(f)]^*\\ &= [X(f)Y^*(f)][G(f)H^*(f)]\\ &= [X(f)Y^*(f)][G^*(f)H(f)]^*, \end{align*}$

R_{\hat{x}, \hat{y}}

$R_{\hat{x},\hat{y}}$

R_{x, y}

$R_{x,y}$

R_{h, g}

$R_{h,g}$

g

$g$

h

$h$

von

und

ist die multiplikative Inverse der Kreuzspektraldichte

von

und

oder etwas in der Nähe davon. Wenn Sie nur ein Signal und einen Filter haben, erhalten Sie das von Hilmar gegebene Ergebnis (mit der Änderung, die durch meinen Kommentar dort gegeben wurde). In jedem Fall bleibt das Problem der Kompensation von spektralen Nullen oder allgemein von Frequenzbändern, in denen die Signale noch wenig Energie haben, bestehen.

G (f) H^{*} (f)

$G(f)H^*(f)$

g

$g$

h

$h$

X (f) Y^{*} (f)

$X(f)Y^*(f)$

x

$x$

y

$y$

— Dilip Sarwate
quelle

Danke für die Nachricht - können Sie uns die Länge erklären? Was ist zum Beispiel die Länge der Leistungsübertragungsfunktion von X, wenn x [n] die Länge N hat? (

— Spacey

Ok - ich werde die Antwort akzeptieren, aber irgendwann heute Abend eine brandneue Frage dazu schreiben, und wir können sie von dort übernehmen. Danke noch einmal.

— Spacey

7

Das Voraufhellen kann durch Filtern mit einer Übertragungsfunktion erfolgen, die ungefähr der Umkehrung des Leistungsspektrums des Signals entspricht. Angenommen, Sie haben ein Audiosignal, das ungefähr rosa ist. Um dies aufzuhellen, würden Sie einen inversen Pink-Filter anwenden (Frequenzgang steigt um 3 dB pro Oktave).

Ich bin mir jedoch nicht sicher, ob dies bei Ihrem Problem helfen wird. Das Voraufhellen verstärkt in der Regel die energiearmen Teile des Signals, die verrauscht sein können und daher das Gesamtrauschen in Ihrem System erhöhen. Wenn Sie versuchen zu bestimmen, ob zwei Signale zeitlich ausgerichtet sind (oder wie die zeitliche Ausrichtung ist), ist das Problem, das mit der Bandbreite des Signals zusammenhängt, inhärent verschwommen. Dies ist genau in der Zeitbereichsform der Autokorrelationsfunktion dargestellt.

— Hilmar
quelle

Vielen Dank für Ihre Antwort - ja das Spektrum invertiert , wie Sie wahrscheinlich hier nicht funktionieren , sagen ... die Verwendung von ‚pre-Aufhellern‘ scheint so allgegenwärtig , neige ich dazu , zu denken , dass es viele Möglichkeiten es neben dem zu tun ...?

— Spacey

2

$\bf x$ $\bf x$ $\bf x$

$\bf x$ $C_{ij} = \frac{1}{N}\sum_{{\bf x}\in Data}x_i x_j$ $N$ $i,j$ $\bf x$

Sobald Sie diese Kovarianzmatrix haben, können Sie eine Aufhellungstransformation in Form einer Matrix berechnen, um die Daten zu multiplizieren und die aufgehellte Version zu erhalten. Die Kovarianz dieser neuen aufgehellten Daten ist die Identitätsmatrix.

${\bf y} = C^{-1/2}{\bf x}$

$C^{-1/2}$ $C=LL^T$ ${\bf y}=L^{-1}{\bf x}$ , die seit $L$ Dieses untere Dreieck kann effizient durch gewöhnliche Löser berechnet werden, ohne das Inverse zu bilden.

— Roboterwanzen
quelle

0

Wenn es nur darum geht, die energiearmen Teile des Signals zu filtern, können Sie einen Tiefpassfilter verwenden? Hierzu gibt es einige Implementierungen.

Wenn dies hilfreich sein sollte: Dieser Artikel von Karjalaien et. Bei al handelt es sich um das Bleaching-Filter und die Methode der verzerrten linearen Vorhersage, die vom Filter verwendet wird.

— jcomouth
quelle