DE10211334A1

DE10211334A1 - Verfahren zum Messen der effektiven Direktivität und/oder effektiven Quelltor-Anpassung eines systemkalibrierten vektoriellen Netzwerkanalysators

Info

Publication number: DE10211334A1
Application number: DE10211334A
Authority: DE
Inventors: Thomas Reichel; Harald Jaeger
Original assignee: Rohde and Schwarz GmbH and Co KG
Current assignee: Rohde and Schwarz GmbH and Co KG
Priority date: 2002-03-14
Filing date: 2002-03-14
Publication date: 2003-09-25
Anticipated expiration: 2022-03-15
Also published as: DE10211334B4; US7076382B2; EP1483593B1; JP4141961B2; DE50303951D1; WO2003076956A3; WO2003076956A2; JP2005526965A; EP1483593A2; US20050119848A1

Abstract

Zum Messen der effektiven Direktivität und/oder der effektiven Quelltor-Anpassung eines Meßtores (M) eines systemkalibrierten vektoriellen Netzwerkanalysators (N) wird eine am Ausgang kurzgeschlossene Präzisions-Luftleitung (L) angeschaltet und der komplexe Reflexionskoeffizient am Eingang dieser Präzisions-Luftleitung (L) an einer Folge von Meßpunkten innerhalb eines vorgegebenen Frequenzbereiches gemessen. Dabei wird für die effektive Direktivität die Folge der gemessenen komplexen Reflexionskoeffizienten einer diskreten Fourier-Transformation unterworfen und aus dem dadurch gebildeten Spektrum das Basisband ausgefiltert. Durch anschließende inverse Fourier-Rücktransformation wird die Folge von effektiven Direktivitäts-Werten gewonnen.

Description

Die Erfindung betrifft ein Verfahren zum Messen der effektiven Direktivität (residual system directivity) und/oder der effektiven Quelltor-Anpassung (residual system port impedance match) eines Meßtores eines systemkalibrierten vektoriellen Netzwerkanalysators laut Oberbegriff des Anspruchs 1.

Die hohe Genauigkeit von vektoriellen Netzwerkanalysatoren (VNA) beruht darauf, daß vor der eigentlichen Messung von Betrag und Phase des komplexen Reflexionskoeffizienten der Netzwerkanalysator durch Anschalten von Kalibrierstandards an seine Meßtore kalibriert wird. Zwischenzeitlich gibt es eine Vielzahl von unterschiedlichen Kalibriermethoden. Die Häufigsten benutzen zur Systemkalibrierung Leerlauf-, Kurzschluß- und Anpassungs-Kalibrierstandards. Durch Anschluß dieser Kalibrierstandards an die Meßtore des Netzwerkanalysators können die im Netzwerkanalysator entstehenden Fehler, die zu einer Abweichung der Meßwerte vom wahren Wert führen, ermittelt und bei der anschließenden Objektmessung zur rechnerischen Fehlerkorrektur benutzt werden. Dies ist beispielsweise aus der DE 39 12 795 A1 bekannt. Diese bisher üblichen Kalibrierverfahren sind jedoch immer noch nicht genau genug.

Um die noch verbleibende restliche Unsicherheit der Direktivität bzw. Meßtor-Anpassung zu bestimmen, wird in einer EA-Richtlinie vorgeschlagen, an dem zu vermessenden Meßtor des vorher systemkalibrierten Netzwerkanalysators eine am Ausgang definiert fehlabgeschlossene bzw. kurzgeschlossene Präzisions-Koaxial-Luftleitung anzuschließen und den Reflexionskoeffizienten am Eingang dieser Luftleitung an einer Folge von Meßpunkten innerhalb eines vorgegebenen Frequenzbereiches des Netzwerkanalysators zu messen (EA-10/12, EA Guidelines on the Evaluation of Vector Network Analyzers (VNA), European co-operation for Accreditation, May 2000). Nach dieser Richtlinie wird allerdings nur die sog. Ripple-Amplitude der den Betrag des Reflexionskoeffizienten überlagernden Oszillation ausgewertet, wobei vereinfachend angenommen wird, daß diese Ripple-Amplitude in etwa identisch ist mit der effektiven Quelltor-Anpassung, was jedoch nur zutrifft, wenn die effektive Direktivität vernachlässigt wird. Diese bekannte Verifikationsvorschrift unter Verwendung einer Präzisions- Luftleitung ist daher relativ ungenau und ermöglicht keine genaue Abschätzung der zu erwartenden Meßunsicherheit geschweige denn eine Nachkorrektur der Fehlerkorrekturterme für die Quelltoranpassung.

Es ist Aufgabe der Erfindung, ein Meßverfahren aufzuzeigen und einen Kalibrierstandard-Satz zu schaffen, mit dem die effektive Direktivität und/oder effektive Quelltor-Anpassung wesentlich genauer ermittelt werden kann und zwar so genau, daß damit dann sogar die bei der Systemkalibrierung ermittelten und im Netzwerkanalysator abgespeicherten Fehlerkorrekturwerte entsprechend nachkorrigiert werden können.

Diese Aufgabe wird gelöst durch die Verfahren gemäß Anspruch 1 oder Anspruch 2 und einen Kalibrierstandard-Satz nach Anspruch 10.

Nach dem erfindungsgemäßen Verfahren kann aus den gemessenen komplexen Reflexionskoeffizienten die effektive Direktivität bzw. die effektive Quelltor-Anpassung wesentlich genauer bestimmt werden, als dies nach der bekannten EA-Richtlinie möglich ist. Die ermittelten Meßwerte stehen mit einem Genauigkeitsniveau zu Verfügung, das dem der verwendeten Präzisions-Luftleitung entspricht. Dabei ist es unerheblich, ob die Impedanz der verwendeten Präzisions-Luftleitung von der Bezugsimpedanz des Netzwerkanalysators abweicht, da gemäß der Erfindung auch solche Impedanzabweichungen berücksichtigt werden können, sofern sie bekannt sind. Damit ist diese Korrektur auch zu relativ niedrigen Frequenzen hin möglich, bei denen die Impedanz der Luftleitung durch den geringer werdenden Skin-Effekt zunehmend vom Nennwert abweicht. Da beim erfindungsgemäßen Verfahren eine Impedanzabweichung der Luftleitung zulässig ist, können auch preiswertere Luftleitungen für die Messung benutzt werden, solange eine ausreichende Längshomogenität der Querschnittsabmessungen gewährleistet ist. Außerdem ist prinzipiell nur die Messung mit einer kurzgeschlossenen Präzisions-Luftleitung notwendig, während beim bekannten Verfahren zusätzlich noch unbedingt eine Messung mit einem definierten Fehlabschluß nötig ist. Die gemäß der Erfindung vorgesehene zusätzliche Messung mit einem definierten Fehlabschluß einer Luftleitung kann zur Aufdeckung von Bedienfehlern benutzt werden. Dadurch wird eine routinemäßige Überprüfung und Bestimmung des Restfehlers bei Netzwerkanalysatoren stark vereinfacht.

Ein besonderer Vorteil des erfindungsgemäßen Verfahrens ist, daß mit der erzielbaren Genauigkeit bei der Messung der effektiven Quelltor-Anpassung bzw. der effektiven Direktivität die eigentlichen Fehlerkorrekturwerte korrigiert werden können, die durch die vorhergehende Kalibrierung des Netzwerkanalysators gewonnen wurden und im Netzwerkanalysator abgespeichert sind. Auf diese Weise wird die Meßgenauigkeit eines solchen vektoriellen Netzwerkanalysators wesentlich erhöht und eine Meßgenauigkeit erreicht, die der Qualität der Präzisions- Luftleitung entspricht. Das erfindungsgemäße Verfahren kann sowohl bei Netzwerkanalysatoren mit nur einem Meßtor (Reflektometer) oder auch bei solchen mit zwei oder mehr Meßtoren angewendet werden. Bei mehreren Meßtoren erfolgt die Messung der effektiven Direktivität und der effektiven Quelltor-Anpassung an jedem der Meßtore gemäß dem erfindungsgemäßen Verfahren getrennt.

Nachdem die gemäß dem erfindungsgemäßen Verfahren in Abhängigkeit von der Frequenz gemessenen Werte für effektive Direktivität und effektive Quelltor-Anpassung im wesentlichen den bei der Systemkalibrierung des Netzwerkanalysators verwendeten Kalibrierstandards zuzurechnen sind, bietet sich an, diese gemessenen Restfehlerwerte auf einem geeigneten Datenträger beispielsweise als Meßprotokolle oder Diagramme oder als Digitalwerte auf einer Diskette abzuspeichern und dem bei der Systemkalibrierung zu verwendenden Kalibrier-Kit beizugeben, so daß der Benutzer nach dem Kalibrieren seines Netzwerkanalysators mit den Kalibrierstandards unmittelbar die ermittelten Restfehlerwerte in den Netzwerkanalysator eingeben und damit die dort abgespeicherten Kalibrierdaten entsprechend korrigieren kann.

Die Erfindung wird im folgenden anhand schematischer Zeichnungen und Diagrammen an einem mathematischen Modell näher erläutert. In der Zeichnung zeigen:

Fig. 1 den Signalflußgraphen eines VNA zur Impedanzmessung;

Fig. 2 den Signalflußgraphen einer fehlabgeschlossenen Leitung;

Fig. 3 ein Ersatzschaltbild einer fehlabgeschlossenen Luftleitung;

Fig. 4 Meßwerte für den Betrag des Reflexionskoeffizienten am Eingang der kurzgeschlossenen Luftleitung;

Fig. 5 Meßwerte für den komplexen Reflexionskoeffizienten am Eingang der kurzgeschlossenen Luftleitung in parametrischer Darstellung;

Fig. 6 die diskrete Fourier-Transformierte (DFT) des Reflexionskoeffizienten am Eingang der kurzgeschlossenen Luftleitung;

Fig. 7 verlängerte und gespiegelte Punktfolge;

Fig. 8 Ausfilterung des heruntergemischten Trägers mit Tiefpaß-Übertragungsfunktion;

Fig. 9 die Beträge des Reflexionskoeffizienten und der Trägerkomponente A im Vergleich;

Fig. 10 den Betrag des Seitenbandsignals;

Fig. 11 das Spektrum des Seitenbandsignals mit Tiefpaß- Übertragungsfunktion;

Fig. 12 die Komponente B_n in Polarkoordinaten-Darstellung (Parameter: Frequenz in GHz). Das Teilbild rechts erhält man bei stark reduzierter Bandbreite des Tiefpaßfilters (1/40);

Fig. 13 das Spektrum des Terms S_ac/A² mit Tiefpaß- Übertragungsfunktion;

Fig. 14 die Komponente C_n/A_n ² in Polarkoordinaten- Darstellung (Parameter: Frequenz in GHz). Das Teilbild rechts erhält man bei reduzierter Bandbreite des Tiefpaß-Filters (1/5);

Fig. 15 die vektorielle Differenz der Meßergebnisse nach Fig. 14;

Fig. 16 die Ripple-Amplitude des Meßergebnisses nach Fig. 4 im Vergleich mit den Beträgen der Meßergebnisse für "residual system directivity" und "residual system port impedance match" nach der erfindungsgemäßen Methode (f/GHz = n/10);

Fig. 17 einen Ausschnitt aus dem Spektrum des heruntergemischten Signals nach Fig. 8, dargestellt für unterschiedliche Grade der Extrapolation. I: o. Extrapolation, o. Spiegelung, 400 Punkte. II: 1600 Punkte entsprechend Fig. 7 (8). III: entsprechend Fig. 18;

Fig. 18 die stark verlängerte und bis zu lokalen Extrema geführte Punktfolge mit Meßpunkten aus Fig. 1 und

Fig. 19 den prinzipieller Aufbau eines Netzwerkanalysators

Fig. 19 zeigt den prinzipiellen Aufbau eines vektoriellen Netzwerkanalysators N und zwar der Einfachheit halber eines Netzwerkanalysators mit nur einem Meßtor M (Reflektometer). Ein eingebauter Generator G speist das Meßtor M (äußeres Meßtor bei Vergleich mit Fig. 1) mit einem Meßsignal veränderbarer Frequenz und zwar beispielsweise in äquidistanten Frequenzschritten zwischen 100 MHz bis 40 GHz. In die Verbindungsleitung zwischen Generator G und Meßtor M sind zwei Richtkoppler R1 und R2 eingefügt. Durch den einen R1 wird die Amplitude der hinlaufenden Welle a und durch den anderen R2 die Amplitude der rücklaufenden Welle b gemessen. Zum Kalibrieren eines derartigen Netzwerkanalysators werden am Meßtor M nacheinander verschiedene Kalibrierstandards wie Leerlauf, Kurzschluß und Anpassung angeschaltet und damit Kalibriermessungen durchgeführt. Die so ermittelten Fehlerwerte werden gedanklich in einem Fehlerzweitor Z zusammengefaßt und beispielsweise in einem Mikroprozessor X im Netzwerkanalysator als Fehlerkorrekturwerte abgespeichert.

Zur Abschätzung der Genauigkeit, mit der ein derartiger Netzwerkanalysator systemkalibriert ist, wird gemäß der Erfindung am Meßtor M eine koaxiale Präzisions-Luftleitung L einer vorgegebenen Mindestlänge angeschaltet, die an ihrem vom Meßtor M abgewandten Ausgang kurzgeschlossen ist. Mit dieser Luftleitung kann nun im Frequenzbereich des Generators G an einer Folge von äquidistanten Meßpunkten der komplexe Reflexionskoeffizient gemessen werden und daraus dann im Sinne des nachfolgenden mathematischen Modells die effektive Direktivität und/oder die effektive Quelltor- Anpassung bestimmt werden. Damit können dann die im Mikroprozessor X abgespeicherten Fehlerkorrekturwerte des Netzwerkanalysators ihrerseits korrigiert werden.

Beim nachfolgenden mathematischen Modell werden folgende Symbole verwendet:
a_e, b_e: Prädiktionskoeffizienten
A_(n): "Träger"-Signal innerhalb des aufgenommenen Satzes von Messwerten (n-ter Messpunkt)
A'_(n): Heruntergemischtes "Träger"-Signal (n-ter Messpunkt)
α: Dämpfungskonstante der Referenz-Luftleitung
B_(n): "Basisband"-Signal innerhalb des aufgenommenen Satzes von Messwerten (n-ter Messpunkt)
β: Phasenkonstante der Referenz-Luftleitung
cfft(x): Diskrete (komplexe) Fourier-Transformierte
C_(n): Signal bei der doppelten Trägerfrequenz innerhalb des aufgenommenen Satzes von Meßwerten (n-ter Meßpunkt)
δ: Residual system directivity (effektive Direktivität)
Δ: Error vector magnitude (Betrag des Fehlervektors)
ΔL₁₍₂₎: Ersatzinduktivität für den Einfluß des Steckers an Tor 1(2) der Referenz-Luftleitung
ΔC₁₍₂₎: Ersatzkapazität für den Einfluß des Steckers an Tor 1(2) der Referenz-Luftleitung
ΔX₁₍₂₎: Blindwiderstand (Reaktanz) in Serie zum Stecker 1(2) der Referenz-Luftleitung
ΔY₁₍₂₎: Blindleitwert parallel zum Stecker 1(2) der Referenz- Luftleitung
e: Index
E_a(e): Anzahl der durch lineare Prädiktion hinzugefügten Punkte am Anfang (a) bzw. Ende (e) der Folge von Meßpunkten
F: "Frequenz" des "Träger-Signals
γ: Ausbreitungskonstante der Referenz- Luftleitung
Γ_a: Reflexionskoeffizient eines Meßobjekts
Γ_ac: Reflexionskoeffizient Γ_a (nach Systemfehlerkorrektur)
Γ_ac _n: n-ter Meßpunkt für Γ_ac
Γ_{ac, mix} _n: Heruntergemischte Folge (n-ter Meßpunkt)
Γ_l: Reflexionskoeffizient des Fehlabschlusses am Ausgang der Referenz-Luftleitung
κ: Proportionalitätskonstante
k: Index
k_max: Maximalwert von k
l: Gesamtlänge der Leitung zwischen der Referenzebene des VNA und der Ebene des Fehlabschlusses
l₁: Länge der Referenz-Luftleitung zwischen den Referenzebenen der beiden Stecker
l₂: Länge des Leitungsabschnitts in der physikalischen Komponente "Fehlabschluss"
µ: residual system port impedance match (effektive Quelltoranpassung)
n, n_a, n_e: Indizes
N Anzahl Meßpunkte
ν: Index
p: Index
P: Anzahl der Prädiktionskoeffizienten
r_z: Normierte Impedanzabweichung der Referenz-Luftleitung (äqivalenter Reflexionsfaktor)
s_xy: s-Parameter der Referenz-Luftleitung
S_ac _n: Seitenbandsignal (n-ter Meßpunkt)
t: Zeit
t₁₍₂₎: Äquivalente Zeitkonstanten für den Steckereinfluß an Tor 1(2) der Referenz-Luftleitung
τ: residual system tracking
ω: Kreisfrequenz
Z₀: Bezugsimpedanz des VNA
ΔZ: Impedanzabweichung der Referenz-Luftleitung

Ausgehend von Fig. 1 gilt für den unter Berücksichtigung der Fehlerterme korrigierten Reflexionskoeffizienten eines Messobjekts:

Für den eingangsseitigen Reflexionskoeffizienten einer mit Γ₁ fehlabgeschlossenen Luftleitung gilt nach Fig. 2:

Unter Vernachlässigung von Produkten der Terme δ, τ, µ, s₁₁ und s₂₂ miteinander erhält man:

Zur Bestimmung der s-Parameter der Luftleitung wird das Ersatzschaltbild in Fig. 3 herangezogen. Der homogene Teil der Luftleitung zwischen den beiden Konnektoren habe dabei die Länge l₁ und werde durch die Ausbreitungskonstante γ = α + β sowie eine Impedanzabweichung ΔZ charakterisiert. Die kurze Leitung mit der Länge l₂ ist Teil des Fehlabschlusses mit dem Reflexionskoeffizienten Γ_l, und zwar der Abschnitt zwischen der Referenzebene und der Ebene des Fehlabschlusses. Die Impedanzabweichung dieses Abschnitts wird zu Null und die Ausbreitungskonstante so groß wie für den langen homogenen Leitungsabschnitt angenommen. Etwaige Impedanzabweichungen kann man in Form von Reaktanzen dem Halbglied zwischen den beiden Leitungsabschnitten zuschlagen. Die beiden Halbglieder repräsentieren die durch die HF-Konnektoren gebildeten Störstellen.

Mit

l = l₁ + l₂ (4)

gilt in einfacher Näherung

s₁₂ = S₂₁ ∼ e^- ^γ ^l (5)

Für die Reflexionsparameter gilt:

1. Fall Γ_l = -1 (Verifikation mit ausgangsseitigem Kurzschluss)

Nach Einsetzen der Ausdrücke für s₁₁, s₂₂ und s₂₁ sei in Gl. 3 erhält man:

Mit Re (2γ l₂) << 1 gilt die vereinfachte Beziehung

Mit

ΔX₁ = ωΔL₁ (10)

ΔX₂ = ωΔL₂ (11)

ΔY₁ = ωΔC₁ (12)

ΔY₂ = ωΔC₂ (13)

gilt:

erhält man schließlich:

Γac ∼ δ + r_z + jωt₁ + [µ - r_z + jωt₁]_e ^-4 ^γ ^l - [1 + τ + j[2r_zsin(2βl₂) - 2ωt₂cos(2βl₂)]]_e ^-2 ^γ ^l (18)

Γ_ac ist eine Funktion der Frequenz, und im Falle eines linearen Sweeps,

ω = κt (19)

auch eine Funktion der Zeit. Dann kann man sich Γ_ac als Summe einer komplexen Schwingung mit dem Träger

A = -[1 + τ + j[2r_zsin(2βl₂) - 2ωt₂cos(2βl₂)]]e^-2 ^γ ^l (20)

einem Basisband-Signal

B = δ + r_z + jωt₁ (21)

und einem Signal bei der doppelten Trägerfrequenz

C = [µ - r_z + jωt₁]e^-4 ^γ ^l (22)

vorstellen, die zu dem charakteristischen Ripple von |Γ_ac| über der Frequenz führen. Wie weiter unten gezeigt wird, lassen sich die genannten spektralen Anteile von Γ_ac durch eine diskrete Fourier-Transformation (DFT) und anschließende Ausfilterung gewinnen, so daß die gesuchten Größen δ und µ wie folgt berechnet werden können:

Gl. 23 kann weiter vereinfacht werden, ohne daß sich an der Genauigkeit, mit der µ bestimmt werden kann, viel ändert:

Die Terme r_z ± jωt₁ definieren einen Fehlervektor, der nur noch von der verwendeten Luftleitung herrührt und dessen Betrag der Impedanzabweichung (incl. Steckereinfluß) entspricht. Wenn dieser Vektor, oder zumindest die einfacher zu bestimmende Größe r_z, zur Verfügung steht, können auch Luftleitungen mit größerer Impedanztoleranz verwendet werden.

Die Korrektur der Fehlerterme "system directivity" (Direktivität) D und "system port impedance match" (Quelltoranpassung) M des VNA mit den gewonnenen Werten erfolgt mit folgenden Gleichungen:

D_neu = D_alt + δ (26)

M_neu = M_alt + µ (27)

2. Fall |Γ₁| < 0,1 (Verifikation mit kleinem ausgangsseitigen Fehlabschluss)

Für diesen Fall kann zunächst Gl. 3 vereinfacht werden:

Γ_ac ∼ δ + s₁₁ + (1 + τ + s₂₂Γ_l)s₁₂s₂₁Γ_l (28)

Nach Einsetzen der s-Parameter der Luftleitung erhält man:

Das Signal enthält drei spektrale Komponenten - Basisband (unteres Seitenband), "Träger" (. . . × e^-2 ^γ ^l) und oberes Seitenband (. . . × e^-4 ^γ ^l) - wobei das obere Seitenband auf Grund seiner Kleinheit nicht auswertbar sein dürfte. Der "Träger" enthält keine interessanten Informationen, das Basisband ist identisch mit jenem, das bei kurzgeschlossener Luftleitung erhalten wird (Term B nach Gl. 21).

Im folgenden wird die Auswertung der Messergebnisse erläutert.

1. Fall Γ_l = -1 (Verifikation mit ausgangsseitigem Kurzschluss)

Nachfolgend wird die Bestimmung der Komponenten B und C/A² aus den Messergebnissen für Γ_ac beschrieben. Dabei sei vorausgesetzt, dass N äquidistante Messpunkte

vorliegen (Fig. 4 und 5):

1. Fouriertransformation mit Suche der dominierenden Komponente F (Fig. 6).
2. Nullmischung des Trägers durch Multiplikation der Γ_ac _n mit
3. Verlängerung der Punktefolge Γ_{ac, mix} _n durch Hinzufügen von E_e Punkten am Ende (n_e = N. . . N + E_e - 1), die durch Extrapolation mit einem linearen Prädiktor über die letzten P Punkte (p = N - P - 1. . . N - 1) gewonnen werden.
4. Verlängerung der Punktefolge durch Hinzufügen von E_a Punkten am Anfang (n_a = -E_a. . . -1), die dadurch gewonnen werden, dass die gespiegelte Folge

( ν = 0. . . N + E_a - 1) auf die unter 3) beschriebene Weise extrapoliert wird.
5. Die extrapolierte gespiegelte Folge von 4) mit N + E_a + E_e Punkten wird zurückgespiegelt und daran die Folge von 4) angehängt, so dass eine neue Folge Γ_{ac, mix} _n mit k = 2 (N + E_a + E_e) Punkten entsteht. Anfangs- und Endpunkt haben per Definition denselben Wert (Fig. 7).
6. Diskrete Fourier-Transformation und Gewinnung des heruntergemischten Trägers mit einem idealen Tiefpassfilter (rechteckförmige Übertragungsfunktion, keine Laufzeit, s. Fig. 8).
7. Inverse Fouriertransformation des Trägers und Ausschneiden des interessierenden Abschnitts der A'_n(n = 0 . . . N - 1).
8. Gewinnung des Trägers in der Originallage durch Heraufmischen (Fig. 9):
9. Subtraktion des Trägers von Γ_ac _n, so daß das Basisbandsignal und das Signal bei der doppelten Trägersequenz übrig bleiben (Fig. 10):
10. Gewinnung des Basisbandsignals (B_n) von S_ac (Fig. 11, 12) analog zu 3) bis 7).
11. Gewinnung des Terms C_n/A 2|n aus S_ac/A²(Fig. 13, 14, 15) analog zu 3) bis 7).

Fig. 16 zeigt das Ergebnis der so aus den gemessenen Reflexionskoeffizienten berechneten Werte für den Betrag der effektiven Direktivität und den Betrag der effektiven Quelltor-Anpassung und zwar im Vergleich zum Betrag der Ripple-Amplitude der Oszillation, die dem Betrag des Reflexionskoeffizienten überlagert ist. Fig. 16 zeigt, daß zwischen dieser bisher ausschließlich ausgewerteten Ripple- Amplitude und den erfindungsgemäß gemessenen Werten erhebliche Abweichungen bestehen, wobei die erfindungsgemäß gemessenen und getrennt ausgewerteten Werte wesentlich genauer sind.

Für die Anwendbarkeit des erfindungsgemäßen Verfahrens ist es ganz wesentlich, daß sich die "Spektren" der Komponenten B und C nicht mit denen des Trägers A überlappen, d. h. daß die Entfernungen der in B und C enhaltenen Reflexionen von der Referenzebene geringer sind als die Länge der Luftleitung. Das wiederum bedeutet, daß "residual system directivity" und "residual system port impedance match" keine Anteile des physischen Netzwerkanalysators mehr enthalten sollten, sondern nur noch Anteile, die von der Unzulänglichkeit der verwendeten Kalibrierstandards herrühren. Daher kann das erfindungsgemäße Verfahren immer erst nach einer Systemkalibrierung des VNA durchgeführt werden.

Eine ganz wesentliche Schwierigkeit bei der Gewinnung spektraler Anteile aus einem endlichen Zeitabschnitt ergibt sich dadurch, daß die Fourier-Transformation das Spektrum des periodisch wiederholten Signalausschnitts bildet. Abgesehen von der frequenzmäßigen Diskretisierung entstehen dadurch spektrale Komponenten, die im Signal eigentlich nicht vorhanden sind. Wenn nun aus diesem verfälschten Spektrum die gesuchten Anteile ausgeschnitten werden, entfallen zwangsweise auch einige der durch die periodische Wiederholung hinzugekommenen Linien, so daß der zeitliche Verlauf nach der Rücktransformation zusätzlich verfälscht wird.

Es sind zwei Methoden bekannt, um diesen Effekt zu verringern. Zum einen kann man den Signalabschnitt vor der Fourier-Transformation einer Fensterung unterwerfen, so daß das Signal nahe Null beginnt und endet und damit nur geringe spektrale Verfälschungen entstehen. Allerdings kann man dann auch nicht damit rechnen, daß Anfang und Ende nach der Rücktransformation stimmen, d. h. das Verfahren funktioniert nicht über den ganzen Signalabschnitt, z. B. nur über den Bereich 8 GHz bis 32 GHz bei einem sich über 100 MHz bis 40 GHz erstreckenden Satz von Meßwerten.

Bei der zweiten bekannten Methode wird an den vorhandenen Signalabschnitt noch einmal der gespiegelte Verlauf angehängt, was zumindest einen stetigen Verlauf und damit geringere spektrale Verfälschungen des periodisch wiederholten Signals erzwingt. Trotzdem sind auch hier nach der Rücktransformation Anfang und Ende des Signalabschnitts nur unter Einschränkungen verwendbar.

Das erfindungsgemäße Verfahren nutzt ebenfalls die Möglichkeit der Spiegelung, führt jedoch vorher eine Extrapolation des Signalabschnitts über das Ende und den Anfang (in Richtung "negativer Frequenzwerte") durch. Zur Extrapolation wird ein linearer Prädiktor verwendet, welcher zu einem stetigen und - in frequenzkontinuierlicher Betrachtungsweise - differenzierbaren Verhalten an den Endpunkten führt. Nach Verdoppelung des Signalabschnitts durch Spiegelung erfolgt die diskrete (komplexe) Fouriertransformation (DFT), Gewinnung des gewünschten Spektralanteils und Rücktransformation. Anschließend werden die ohnehin nur gering verfälschten Randabschnitte, d. h. die extrapolierten Anteile, abgeschnitten und das verbleibende Signal weiterverarbeitet.

Lineare Prädiktion gestattet es, die extrapolierten Werte als Linearkombination vorangegangener Werte zu berechnen:

Die Gewichte a_e und b_e werden dabei so bestimmt, dass die Anwendung von Gleichung (33) innerhalb des gemessenen Signalabschnitts zu möglichst kleinen Fehlern führt. Die mathematischen Verfahren zu ihrer Bestimmung sind in der einschlägigen Literatur beschrieben. Es hat sich gezeigt, daß eine Verlängerung des gemessenen Signalabschnitts auf den doppelten Wert für die Zwecke der Erfindung normalerweise ausreicht. Durch das Hinzufügen von noch mehr Punkten, und zwar möglichst bis zu einem lokalen Extremum des Betrags (Fig. 18), wird einerseits die Frequenzauflösung weiter erhöht und andererseits die Möglichkeit für eine erhöhte Selektivität des Verfahrens geschaffen (Fig. 17). Letzteres ist dann wichtig, wenn das Basisband und das Frequenzband bei der doppelten Trägerfrequenz relativ breit sind und nahe an den Träger heranreichen (von Reflexionen hervorgerufen, die relativ weit von der Referenzebene entfernt sind).

Die gesuchten spektralen Komponenten werden durch Tiefpaßfilterung und abschließende Rücktransformation gewonnen. Es genügt vollkommen, "ideale" Tiefpässe mit rechteckförmiger Übertragungsfunktion (1 im Durchlass-, 0 im Sperrbereich) zu verwenden, deren Bandbreite so gewählt ist, daß der Übergang vom Durchlaß- in den Sperrbereich in der Mitte zwischen den zu trennenden spektralen Komponenten liegt. Eine Einschränkung der Bandbreite unter diesen Wert führt einerseits zu einer (erwünschten) Verringerung des dem Meßergebnis überlagerten Rauschens, andererseits kann sie Verfälschungen des Meßergebnisses zur Folge haben, wenn dadurch relevante spektrale Anteile eliminiert werden. In der Praxis wird man einen Kompromiss finden müssen, wobei die Wahl der Bandbreite für die dargestellte Meßaufgabe als unkritisch eingestuft werden konnte (Fig. 14 und 15). Selbst bei extremer Reduzierung (Fig. 12 rechts) konnten noch vernünftige Meßergebnisse erzielt werden. Wenn man sich vor Augen hält, daß die vorangegangene Systemkalibrierung mit der OSM-Methode (Kalibrierstandards Open, Short und Match) durchgeführt und als "Match" eine ab 4 GHz Gleitlast (sliding load) verwendet wurde, ist das Ergebnis auch leicht erklärbar: Da sich die Luftleitungsabschnitte von Gleitlast und Referenz-Luftleitung im wesentlichen durch ihre charakteristische (frequenzunabhängige und reelle) Impedanz sowie die proportional zur Frequenz ansteigenden Reaktanzen der Steckverbindungen an der Referenzebene unterscheiden, ist auch rein theoretisch ein Verlauf wie in Fig. 12 (rechts) zu erwarten.

2. Fall

|Γ_l| < 0,1 (Verifikation mit kleinem ausgangsseitigen Fehlabschluß)

Die Auswertung kann genau nach dem gleichen Schema wie für die Verifikation mit Kurzschluß vorgenommen werden (Punkt 11 entfällt). Das Ergebnis für die effektive Direktivität ("residual system directivity") sollte identisch mit jenem sein, das bei der Verifikation mit Kurzschluß gewonnen wurde. An die Größe des Fehlabschlusses im Verhältnis zur effektiven Direktivität werden keine Anforderungen gestellt (im Gegensatz zur Ripple-Auswertung nach dem Stand der Technik).

Durch diskrete Fourier-Transformation (DFT) der verlängerten und gespiegelten Folge von Reflexionsmeßwerten entsteht eigentlich ein der Impulsantwort des Systems vergleichbares Zeitsignal. Wenn dafür vorstehend trotzdem der Begriff "Spektrum" verwendet wurde, dann nur deswegen, weil die damit zusammenhängenden Begriffe und Vorstellungen geläufiger sind. Im übrigen kann man sich die Reflexionskoeffizienten über Gl. (19) als Zeitsignal vorstellen, so daß der Begriff des "Spektrums" für die DFT der Reflexionskoeffizienten seine Berechtigung hat. Die verwendeten mathematischen Algorithmen sind ohnehin unabhängig von der Begrifflichkeit zu gebrauchen.

Claims

1. Verfahren zum Messen der effektiven Direktivität eines Meßtores (M) eines systemkalibrierten vektoriellen Netzwerkanalysators (N) durch Anschalten einer am Ausgang kurzgeschlossenen Präzisions-Luftleitung (L) und Messung des komplexen Reflexionskoeffizienten am Eingang dieser Präzisions-Luftleitung (L) an einer Folge von Meßpunkten innerhalb eines vorgegebenen Frequenzbereiches, dadurch gekennzeichnet,
daß die Folge der gemessenen komplexen Reflexionskoeffizienten einer diskreten Fourier- Transformation unterworfen wird,
daß aus dem dadurch gebildeten Spektrum das Basisband ausgefiltert wird und
daß durch anschließende inverse Fourier-Rückttansformation die Folge von effektiven Direktivitäts-Werten gewonnen wird (Formel 24).

2. Verfahren zum Messen der effektiven Quelltor-Anpassung eines Meßtores (M) eines systemkalibrierten vektoriellen Netzwerkanalysators (N) durch Anschalten einer am Ausgang kurzgeschlossenen Präzisions-Luftleitung (L) und Messung des komplexen Reflexionskoeffizienten am Eingang der Luftleitung (L) an einer Folge von Messpunkten innerhalb eines vorgegebenen Frequenzbereiches, dadurch gekennzeichnet,
daß die Folge der gemessenen komplexen Reflexionskoeffizienten einer komplexen Fourier- Transformation unterworfen wird,
daß aus dem dadurch gebildeten Spektrum die dominierende Komponente (Träger) ausgefiltert wird und durch inverse Fourier-Rücktransformation die Größe (A) dieser dominanten Komponente bestimmt wird,
daß dann das obere Seitenband bei der doppelten Frequenz der dominanten Komponente ausgefiltert wird und wiederum durch inverse Fourier-Rücktransformation die Größe (C) dieses ausgefilterten oberen Seitenbandes bestimmt wird und
daß schließlich durch jeweilige Division der Größe (A) der dominanten Komponente (Träger) durch das Quadrat der Größe (C) des oberen ausgefilterten Seitenbandes die Folge von effektiven Quelltor-Anpassungs-Werten gewonnen wird (Formel 25).

3. Verfahren zum Messen der effektiven Direktivität und der effektiven Quelltor-Anpassung eines Meßtores (M) eines systemkalibrierten vektoriellen Netzwerkanalysators (N) durch Anschalten einer am Ausgang kurzgeschlossenen Präzisions-Luftleitung (L) und Messung des komplexen Reflexionskoeffizienten am Eingang der Luftleitung (L) an einer Folge von Meßpunkten innerhalb eines vorgegebenen Frequenzbereiches, dadurch gekennzeichnet, daß die Verfahrensschritte nach Anspruch 1 und 2 kombiniert und gleichzeitig parallel ausgeführt werden.

4. Verfahren nach einem der vorhergehenden Ansprüche, gekennzeichnet durch die Verwendung einer Präzisions-Luftleitung (L), deren Impedanz von der Systemimpedanz des Netzwerkanalysators (N) um einen bekannten Wert abweicht, wobei diese bekannte Impedanzabweichung bei der Bestimmung der Folge von Werten der effektiven Direktivität bzw. effektiven Quelltor- Anpassung als entsprechende Korrekturwerte berücksichtigt wird (Gleichungen 24 und 25).

5. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, daß vor der komplexen Fourier-Transformation der Folge der gemessenen komplexen Reflexionskoeffizienten diese Folge von komplexen Reflexionskoeffizienten durch lineare Prädikation extrapoliert wird und durch Spiegelung dieser extrapolierten Folge von komplexen Reflexionskoeffizienten die ursprüngliche Folge verlängert wird.

6. Verfahren nach Anspruch 5, dadurch gekennzeichnet, daß vor der linearen Prädikation die Folge von komplexen Reflexionskoeffizienten mit einem solchen Faktor multipliziert wird, daß die dominierende Komponente (Träger) durch Abwärtsmischung im Spektrum bei der Frequenz Null liegt (Gleichung 30).

7. Verfahren nach Anspruch 1 oder 3 bis 6, dadurch gekennzeichnet, daß an dem zu vermessenden Meßtor (M) des systemkalibrierten vektoriellen Netzwerkanalysators (N) eine am Ausgang leicht fehlabgeschlossene (|Γ_e| < 0,1) Luftleitung (L) angeschlossen wird und anschließend die gemessenen komplexen Reflexionskoeffizienten gemäß den Verfahrensschritten nach Anspruch 1 oder der Ansprüche 3 bis 6 behandelt werden.

8. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, daß mit den gewonnenen Werten der effektiven Direktivität und/oder der effektiven Quelltor-Anpassung die Systemfehler- Terme des Netzwerkanalysators (N) korrigiert werden (Formeln 26 und 27).

9. Verfahren nach einem der vorhergehenden Ansprüche 1 bis 7, dadurch gekennzeichnet, daß die gewonnene Folge von Werten der effektiven Direktivität und/oder der effektiven Quelltor-Anpassung auf einem Datenträger abgespeichert wird, der zusammen mit den zur Systemkalibrierung des Netzwerkanalysators (N) benutzten Kalibrierstandards vertreibbar ist.

10. Kalibrierstandard-Satz zur Systemkalibrierung von vektoriellen Netzwerkanalysatoren (N) mit verschiedenen vorzugsweise in Koaxialleitungstechnik ausgebildeten Kalibrierstandards für Leerlauf, Kurzschluß und Anpassung, dadurch gekennzeichnet, daß ihm ein Datenträger mit Werten der effektiven Direktivität und/oder effektiven Quelltor-Anpassung beigegeben ist, die gemäß dem Verfahren nach einem der vorhergehenden Ansprüche gewonnen sind.