Welche Verbindung besteht zwischen dem analogen Signal-Rausch-Verhältnis und dem Signal-Rausch-Verhältnis in der IQ-Ebene in einem Quadratur-Demodulationssystem?


9

Wir möchten die quantitative Beziehung zwischen analogem Rauschen in der Nähe der LO-Frequenz und der Statistik der Punkte berechnen, die in der IQ-Ebene nach der IQ-Demodulation gefunden wurden. Um die Frage vollständig zu verstehen, geben wir zunächst eine detaillierte Beschreibung des IQ-Demodulationssystems.

IQ-Demodulationssystem

Ein IQ-Mischer nimmt Signale mit hoher Frequenz auf und bringt sie auf eine niedrigere Frequenz, damit sie leichter verarbeitet werden können. 1 zeigt eine schematische Darstellung eines IQ-Mischers. Das lokale Oszillatorsignal (LO) ) wird verwendet, um das HF-Signal auf eine niedrigere Frequenz herunterzumischen.cos(Ωt

Komplettes IQ-Demodulationssystem Abbildung 1: Vollständige Signalverarbeitungskette. Das Mikrowellenfrequenzsignal (und das Rauschen) gelangen über den HF-Anschluss in den IQ-Mischer. Dieses Signal wird mit einem lokalen Oszillator (LO) gemischt, um in Zwischenfrequenzsignale und umzuwandeln . Die Zwischenfrequenzsignale werden dann gefiltert, um die verbleibende Hochfrequenzkomponente (siehe Text) zu entfernen, und digital abgetastet. Die Erfassung der Amplitude und Phase jeder Frequenzkomponente erfolgt über eine diskrete Fourier-Transformation in digitaler Logik.IQ

Kohärentes Signal - Gleichstromfall

Angenommen, im eingehenden HF-Signal wären . Dann wären die und Signale Wir leiten diese Signale durch Tiefpassfilter, um die Terme zu entfernen , was Wie wir sehen können, sind DC undMcos(Ωt+ϕ)IQ

I(t)=M2cos(ϕ)+M2cos(2Ωt+ϕ)Q(t)=M2sin(ϕ)M2sin(2Ωt+ϕ).
2Ω
IF(t)=M2cos(ϕ)QF(t)=M2sin(ϕ).
IQSpannungen können als kartesische Koordinaten betrachtet werden, die die Amplitude und Phase des ursprünglichen Signals angeben. Daher hat der Mischer seine Aufgabe erfüllt, die Amplitude und Phase eines Hochfrequenzsignals zu ermitteln, indem nur Niederfrequenzmessungen durchgeführt wurden.

Kohärentes Signal - Wechselstromfall

In der Praxis demodulieren wir das HF-Signal normalerweise nicht in Gleichstrom. Dafür gibt es mehrere Gründe:

  1. Die spektrale Rauschdichte steigt bei niedrigen Frequenzen fast immer stark an.

  2. Wenn wir gleichzeitig die Amplitude und Phase mehrerer sinusförmiger Komponenten bei verschiedenen Frequenzen messen wollen, können wir im analogen Teil des Systems nicht direkt auf Gleichstrom demodulieren.

Als Beispiel für # 2 könnten wir Um die Amplitude und Phase beider Frequenzkomponenten zu ermitteln, müssen wir eine etwas komplexere Signalverarbeitung verwenden. Die und sind in diesem Fall Um beide Amplituden und beide Phasen zu finden, müssen wir im Wesentlichen eine Fourier-Transformation durchführen. Dazu digitalisieren wir die Wellenformen und geben nach

RF(t)=M1cos([Ω+ω1]t+ϕ1)+M2cos([Ω+ω2]t+ϕ2).
IFQF
IF(t)=M12cos(ω1t+ϕ1)+M22cos(ω2t+ϕ2)QF(t)=M12sin(ω1t+ϕ1)M22sin(ω2t+ϕ2).
In=M12cos(ω1nδt+ϕ1)+M22cos(ω2nδt+ϕ2)Qn=M12sin(ω1nδt+ϕ1)M22sin(ω2nδt+ϕ2)
wobei das digitale Abtastintervall ist. Dann konstruieren wir in der digitalen Logik die komplexe Reihe die durch . Für die oben geschriebenen Signale ist dies Wenn wir jetzt in der digitalen Logik die Summe berechnen δtznznIn+iQn
zn=M12exp(i[ω1nδt+ϕ1])+M22exp(i[ω2nδt+ϕ2]).
Z(ωk)=1Nn=0N1zneiωknδt
Wir stellen die Amplitude und Phase für die Komponente bei der Frequenz . Wenn wir zum Beispiel berechnen würden, würden wir .ωkZ(ω1)(M1/2)exp(iϕ1)

Lärm

In der Praxis kommt das Signal immer mit Rauschen. Das Rauschen bewirkt, dass eine Zufallsvariable anstelle eines deterministischen Werts ist. Mit anderen Worten, für jedes zufällig und wird für jede Durchführung des Experiments unterschiedlich sein.Z(ω) ω Z(ω)

Wir können aus der Intuition erraten, dass in Gegenwart von Rauschen eine zirkularsymmetrische Verteilung in der IQ-Ebene hat, deren Mittelwert dem deterministischen Wert . Die Frage ist, wie genau ist die statistische Verteilung von bei Vorhandensein von Rauschen?Z(ω)(M/2)exp(iϕ)Z

Antworten:


9

Da jeder Schritt in der Verarbeitungskette linear ist, betrachten wir einen Fall mit nur Rauschen und ohne kohärentes Signal. Bezeichne das Rauschen . Die und Signale sind Wir drücken die Wirkung des Filters als Faltung mit der Zeitantwortfunktion , und ähnlich für . Da der Filter kausal ist, ist für . Die Abtastung wählt einfach den Wert von undξ(t)IQ

 I(t)=ξ(t)cos(Ωt)Q(t)=ξ(t)sin(Ωt).
h
IF(t)=dtξ(t)cos(Ωt)h(tt)
QFh(t)=0t<0IFQF zu den Zeiten , und ähnlich für . Nach der oben für den digitalen Teil der Verarbeitungskette beschriebenen Konstruktion haben wir Unser Problem ist daher, die Statistik dieses Ausdrucks zu berechnen.{nδt}
In=dtξ(t)cos(Ωt)h(nδtt)
Qn
Z(ω)=1Nn=0N1dtξ(t)eiΩth(nδtt)eiωnδt.

Das Ändern der Variablen erzeugt In diesem Stadium können wir eine Überprüfung der geistigen Gesundheit durchführen, indem wir den Durchschnittswert von berechnen . Denken Sie daran, dies ist ein Ensemble- Durchschnitt. Mit anderen Worten, wir berechnen den Durchschnittswert von den wir finden würden, indem wir viele Fälle von demoduliertem Rauschen in IQ-Punkte umwandeln und dann den Mittelwert aller dieser Punkte nehmen. In jedem Fall ist das Ergebnis nδttt

Z(ω)=1Nn=0N1dtξ(nδtt)eiΩ(nδtt)h(t)eiωnδt.
Z(ω)Z(ω)
Z(ω)=1Nn=0N1dtξ(nδtt)0eiΩ(nδtt)h(t)eiωnδt=0.
Dies ist sinnvoll, da wir erwarten, dass das Rauschen den Durchschnittswert des demodulierten IQ-Punkts nicht ändert, sondern nur eine Zufälligkeit hinzufügt, die um den deterministischen Wert zentriert ist.

Ich weiß nicht, wie ich die Statistik von direkt berechnen soll , daher verfolgen wir einen alternativen Ansatz, indem wir stattdessen das mittlere Quadrat von berechnen . Nach dem zentralen Grenzwertsatz sollten der Real- und Imaginärteil von mindestens annähernd guassisch verteilt (und, wie wir noch betonen werden, unkorreliert sein), sodass das Finden des mittleren quadratischen Moduls von uns tatsächlich alles sagt, was wir wissen müssen.Z(ω)Z(ω)ZZ

Wir fahren fort, indem wir direkt konstruieren und den statistischen Durchschnitt nehmen (der statistische Durchschnitt wird mit ). |Z(ω)|2

|Z(ω)|2=dtdt1N2n,m=0N1eiΩ(tt)h(t)h(t)ξ(nδtt)ξ(mδtt)ei(Ω+ω)(nm)δt.()
ξ(t)ξ(τ)ξ(0)Sξüber die folgende Gleichung: Verwendung dieser Formel für ergibt
ξ(τ)ξ(0)=12dω2πSξ(ω)eiωτ.
ξ(nδtt)ξ(mδtt)
|Z(ω)|2=12dtdtdω2π1N2n,m=0N1eiΩ(tt)h(t)h(t)Sξ(ω)eiω((nm)δt(tt))ei(Ω+ω)(nm)δt=12dω2π|h(ωΩ)|2Sξ(ω)|1Nn=0N1ei(Ω+ωω)nδt|2=12Ndω2π|h(ωΩ)|2Sξ(ω)1N(sin([Ω+ωω]δtN/2)sin([Ω+ωω]δt/2))2Nth order Fejer kernel=12Ndω2π|h(ωΩ)|2Sξ(ω)FN([Ω+ωω]δt/2)
Dabei ist der Fejer-Kernel der Reihenfolge . Wenn wir die Variablen wir Bisher waren die Ergebnisse genau und genaue Ergebnisse können durch numerische Auswertung der Integrale gefunden werden. Wir treffen nun eine Reihe relativ schwacher Annahmen, um zu einer praktischen Formel zu gelangen. Der Fejer-Kernel hat ein Gewicht, das nahe . Deshalb integrieren wir überFNNthΩωω
|Z(ω)|2=12Ndω2π|h(ω)|2Sξ(Ωω)FN([ω+ω]δt/2).
FN(x)x=0Sξnur für Frequenzen in der Nähe von und so können wir in diesem Integral als Konstante , was Wir können hier bereits sehen, dass die Rauschstatistik des demodulierten IQ-Punkts nur von der HF-Spektraldichte nahe der LO-Frequenz abhängt. Das macht Sinn; Der IQ-Mischer ist so ausgelegt, dass er den Signalinhalt in der Nähe der LO-Frequenz auf eine niedrigere ZF reduziert, wo er verarbeitet werden kann. Die Anti-Aliasing-Filter entfernen alle Frequenzkomponenten, die zu weit vom LO entfernt sind.ΩSξS(Ωω)Sξ(Ω)
|Z(ω)|2=12NSξ(Ω)dω2π|h(ω)|2FN([ω+ω]δt/2).

Die erste Null von tritt bei , und der größte Teil des Gewichts ist in den ersten paar Lappen enthalten. Die ersten Nullen befinden sich daher bei Dies bedeutet, dass das Integral über von Frequenzen in einem Bereich dominiert wird, der durch die Abtastfrequenz geteilt durch . In den meisten praktischen Anwendungen ist dieser Bereich so klein, dass über diesen Bereich ungefähr konstant ist. Wenn dies der Fall ist, können wir durch ersetzen (beachten Sie, dass ) FN(x)x=2π/N

ωnull2π=ω2π±1Nδt.
ωNh(ω)h(ω)h(ω)h(ω)=h(ω)
|Z(ω)|2=12NSξ(Ω)|h(ω)|2dω2πFN([ω+ω]δt/2N)1/δt=Sξ(Ω)2T|h(ω)|2
wobei die Gesamtmesszeit ist.TNδt

Signal-Rausch-Verhältnis

Es ist hinreichend bekannt, dass, wenn eine Zufallsvariable Gaußsche und unabhängig verteilte Real- und Imaginärteile und einen durchschnittlichen quadratischen Modul , die Verteilungen der Real- und Imaginärteile dieser Variablen die Standardabweichung . Nehmen wir daher unser Ergebnis für , unsere Beobachtung, dass der Real- und Imaginärteil von Gauß-verteilt sind, und die Tatsache, dass sie nicht korreliert sind, wir wissen, dass die Standardabweichungen der Verteilungen des Real- und Imaginärteils ZRR/2[a]|Z(ω)|2Z[b]

σ=Sξ(Ω)|h(ω)|2/4T.
Wie zu Beginn diskutiert, wird ein Signal in der IQ-Ebene zu . Natürlich haben wir dort den Effekt des Filters ignoriert, der einfach darin besteht, die Amplitude zu skalieren, um Angenommen, wir verwenden, wie in Abbildung 2 dargestellt, das IQ-Demodulationssystem, um zwischen zwei oder mehr Signalen mit jeweils unterschiedlicher Phase, jedoch mit derselben Amplitude . Aufgrund des Rauschens führt jede der möglichen Amplituden / Phasen zu einer Punktwolke in der IQ-Ebene mit dem radialen Abstand vom Ursprung. Der Abstand zwischen zwei Wolkenzentren beträgtMcos([Ω+ω]t+ϕ)(M/2)eiϕ
Z(ω)=M|h(ω)|2eiϕ.
MM|h(ω)|/2g(M/2)|h(ω)|Dabei ist ein geometrischer Faktor, der von den Phasen der Wolken abhängt. Wenn der Bogenwinkel zwischen zwei Wolken und das Zentrum jeder Wolke gleich weit vom Ursprung entfernt ist, dann ist . Wenn zum Beispiel die beiden Phasen dann ist . Geometrisch liegt dies daran, dass der Abstand zwischen den Wolkenzentren doppelt so groß ist wie der Abstand einer der Wolken vom Ursprung.gθg=2sin(θ/2)±π/2g=2sin(π/2)=2

Das Signal-Rausch-Verhältnis (SNR) ist wobei die eingehende analoge Leistung ist. Beachten Sie, dass das SNR nicht von abhängt . Um sich an dieses Ergebnis zu erinnern, ist zu beachten, dass die Rauschleistung die spektrale Dichte multipliziert mit einer Bandbreite . Wenn wir , sehen wir, dass unser Ergebnis nur besagt, dass das SNR in der IQ-Ebene genau gleich dem analogen SNR multipliziert mit dem geometrischen Faktor .

SNRseparation22×(cloud std deviation)2=(gM|h(ω)|/2)22Sξ(Ω)|h(ω)|2/4T=(gM)2T2Sξ(Ω)=g2PTSξ(Ω).
PM2/2hBB=1/Tg2

Geben Sie hier die Bildbeschreibung ein

Abbildung 2: Zwei IQ-Wolken. Der Abstand zwischen den Wolkenzentren ist proportional zu ihrer radialen Größe , jedoch skaliert mit einem geometrischen Faktor . Projiziert auf die Linie, die ihre Zentren verbindet, wird jede Wolke zu einer Gaußschen Verteilung mit der Breite .MgSξ(Ω)|h(ω)|2/4T

[a] : Schlagen Sie die Chi-Quadrat-Verteilung nach .

[b] : Wir können sehen, dass der Real- und Imaginärteil von tatsächlich nicht korreliert sind, indem wir das Äquivalent von Gleichung schreiben, jedoch für . Wenn wir dies tun, würden wir feststellen, dass die Summe, die sich im Fall von in den Fejer-Kernel verwandelte , (zumindest ungefähr) auf Null gehen würde, da dies ungefähr die Überlappung von Sinus und Cosinus wäre. die orthogonal sind.Z()ZZ|Z|2

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.