Podręcznik

Strona: SEZAM - System Edukacyjnych Zasobów Akademickich i Multimedialnych
Kurs: Fale i prowadnice falowe
Książka: Podręcznik
Wydrukowane przez użytkownika: Gość
Data: czwartek, 16 maja 2024, 07:35

1. Wprowadzenie

Widmo promieniowania EM
Fale elektromagnetyczne EM wykorzystywane są intensywnie do rozmaitych zastosowań, opisanych niżej. Na rys.1.1. pokazano szeroki wycinek widma fal EM ze wskazaniem, które zakresy najczęściej wykorzystywana są do wybranych zastosowań. Aby dany zakres mógł być wykorzystany musi zostać „opanowany” w sensie technologicznym. Rozumiemy pod tym opanowanie techniki generacji, wzmacniania i modulacji sygnałów w tym pasmie.
Pasmo fal elektromagnetycznych jest „atakowane” z dwóch stron, od strony fal radiowych i od strony pasm optycznych. Kolejne podpasma zostają „opanowane” technologicznie.
•    Rozwój techniki radiowej to opanowanie kolejno fal długich, średnich, krótkich i UKF.
•    Rozwój techniki radarowej to opanowanie kolejnych zakresów mikrofal, od fal decymetrowych, poprzez fale centymetrowe do milimetrowych i submilimetrowych.


 
Rys.1.1. Zakresy widma promieniowania elektromagnetycznego


W zakresie promieniowania optycznego, w zakresie podczerwieni odkryto nadzwyczajne właściwości światłowodu kwarcowego, bardzo małe tłumienie przy transmisji sygnałów. Pasmo transmisji światłowodem jest stosunkowo wąskie, obejmujące długości fali 1,2 ...1,6 \mum, ale jego użyteczność jest ogromna. Opanowanie technik światłowodowych wymagało rozwoju technik generacji, modulacji, wzmacniania i detekcji promieniowania optycznego.
Między pasmem dalekiej podczerwieni a falami submilimetrowymi (1THz...100THz) istnieje zakres słabo opanowany, „atakowany z obu stron”.
Granice pasma zwanego mikrofalowym nie są dokładnie precyzowane i przyjmowane są umownie. Zwykle przyjmujemy, że mikrofale, to zakres częstotliwości fal elektromagnetycznych, rozciągający się od 300 MHz do około 1000 GHz. Poniżej wymieniono cztery cechy charakterystyczne zakresu mikrofal.
•    Rozmiary mikrofalowych elementów i obwodów są porównywalne do długości fal.
•    Czas propagacji porównywalny lub wielokrotnie dłuższy od okresu drgań.
•    W zakresie częstotliwości mikrofalowych mamy do czynienia z efektem naskórkowości. 
•    Podstawowym pomiarem zakresu mikrofal jest pomiar mocy.


Rys.1.2. Podstawowe podpasma częstotliwości zakresu mikrofal

Na rys.1.2 pokazano podział podstawowego zakresu częstotliwości pasma mikrofalowego na podpasma, które mają swoje tradycyjne, literowe oznaczenia. Pasmo fal decymetrowych to oznaczane jest przez L, pasmo 3 cm oznaczane jest przez X, itd.
     Podstawowe zastosowania mikrofal
Fale elektromagnetyczne są intensywnie wykorzystywane do rozmaitych celów. Na czele listy zastosowań jest oczywiście telekomunikacja, z telekomunikacją satelitarną, telewizją kablową i telefonią komórkową na czele. Poniżej kolejno wymieniamy najważniejsze zastosowania techniki mikrofalowej.
•    Telekomunikacja

 

  • Telekomunikacja satelitarna,
  • Telewizja satelitarna i kablowa,
  • Radiolinie,
  • Telekomunikacja ruchoma i bezprzewodowa.

•    Radiolokacja

  • Radarowe wykrywania obiektów,
  • Kontrola ruchu lotniczego,
  • Naprowadzanie rakiet i antyrakiet,
  • Radary antykolizyjne,
  • Radionawigacja.

•    Nauka

  • Radioastronomia,
  • Geonawigacja,
  • Wzorce i miernictwo czasu,
  • Spektrometria mikrofalowa,
  • Miernictwo parametrów materiałów,
  • Akceleratory cząstek.

•    Medycyna

  • Termografia mikrofalowa,
  • Hypertermia mikrofalowa.

•    Zastosowania domowe

  • Kuchnie mikrofalowe,
  • Alarmy przeciwwłamaniowe,

•    Przemysł

  • Grzejnictwo mikrofalowe,
  • Pomiary wilgotności materiałów,
  • Pomiary odległości,
  • Czyszczenie i pokrycia powierzchni.

•    Rolnictwo

  • Uzdatnianie gleby,
  • Niszczenie szkodników,
  • Kontrola stanu wegetacji roślin,
  • Wykrywanie pożarów lasów.

Jak widać z zestawienia zakres zastosowań jest bardzo szeroki. Niektóre z zastosowań są ogólnie znane, takie jak telekomunikacja satelitarna, radar, telewizja kablowa, kuchnia mikrofalowa. Inne, jak na przykład akceleratory cząstek, czy niszczenie szkodników są mniej znane, choć ich użyteczność nie budzi wątpliwości.

1.1. Wielkości charakteryzujące pole EM

Pola elektryczne i magnetyczne kiedy zmieniają się, to powiązane są ze sobą wieloma związkami, nazwanymi wspólnie równaniami Maxwella. Równania te oparte są na wielu fizycznych eksperymentach prowadzonych przez uczonych całe dziesięciolecia. Równania Maxwella pełnią rolę fundamentu dla całej techniki mikrofalowej, tworzonej i rozwijanej przez ostatnie 100 lat. 
Równania Maxwella wiążą ze sobą najważniejsze wielkości pola elektromagnetycznego. Zajmiemy się przestrzenią wypełnioną jednorodnym, izotropowym ośrodkiem, charakteryzowanym przez:
•    Przenikalność elektryczna ośrodka  \varepsilon, w próżni równa \varepsilon \equiv \varepsilon _{0} =10^{-9}/36\pi \, [F/m] .
•    Przenikalność magnetyczna ośrodka   \mu, w próżni równa  \mu \equiv \mu _{0}=4\pi \cdot 10^{-7}\, [H/m].
•    Przewodność właściwa ośrodka   \sigma (konduktywność), jej miarą jest \sigma =0 [S/m], w próżni  . 
Pole elektromagnetyczne jest kompozycją pól elektrycznego i magnetycznego. Charakteryzują go następujące wielkości:
•    Natężenie pola elektrycznego  E(x, y, z, t) jest wektorem (ma swoją wartość i kierunek), funkcją miejsca i czasu. Natężenie pola elektrycznego mierzymy w woltach na metr [V/m].
•    Natężenie pola magnetycznego H(x, y, z, t)  jest także wektorem, mierzymy jego wartość w amperach na metr [A/m].
•    Indukcja elektryczna D = \varepsilon E  jest wektorem, w ośrodku izotropowym skierowanym w tą samą stronę, co natężenie pola elektrycznego, w próżni  D = \varepsilon_{0} E, mierzymy jej wartość w amperach razy sekunda na metr kwadratowy [As/m2], czyli kulombach na m2 [C/m2] .
•    Indukcja magnetyczna B = \mu H  jest wektorem, w ośrodku izotropowym skierowanym w tą samą stronę, co natężenie pola magnetycznego, w próżni  B = \mu_{0} H, mierzymy jej wartość w voltach razy sekunda na metr kwadratowy [Vs/m2].

1.2.  Całkowa postać równań Maxwella

Jak powiedziano wyżej, równania Maxwella opierają się odkrytych eksperymentalnie prawach fizycznych. 
Pierwsze prawo Maxwella w postaci całkowej jest zapisem prawa Faradaya. Rozważane są indukcja pola magnetycznego  B = \mu H i natężenie pola elektrycznego E  w sąsiedztwie powierzchni S otoczonej konturem C, co pokazano na rys.1.3.

Rys.1.3. Powierzchnia S o konturze C w polu EM o indukcji magnetycznej B = \mu H i natężeniu pola elektrycznego E.


 
Pierwsze prawo Maxwella wiąże ze sobą zmianę strumienia indukcji magnetycznej przenikającego powierzchnię S z polem elektrycznym   całkowanym wzdłuż zamkniętego konturu C otaczającego tą powierzchnię.

 

\oint_{C}^{}Edl=-\iint_{S}^{}\frac{\partial B}{\partial t}\cdot nds

(1-1)  

Drugie równanie Maxwella w postaci całkowej (1-2) jest zapisem prawa Ampere’a. Wiąże ono ze sobą prąd płynący przez powierzchnię S z polem magnetycznym  , całkowanym wzdłuż zamkniętego konturu C.

 

\oint_{C}^{}Hdl=-\iint_{S}^{}(J+\frac{\partial D}{\partial t})\cdot nds

(1-2)  

Prawo Gaussa zastosowane dla pola elektrycznego mówi, że strumień wektora indukcji pola elektrycznego  D wypływający z objętości V przez zamkniętą powierzchnię S równy jest zgromadzonemu w tej objętości ładunkowi. Ładunek ten otrzymujemy całkując gęstość tego ładunku   mierzoną w kulombach na metr sześcienny w objętości V równanie (1-3).

 

\iint_{S}D\cdot nds=\iiint_{V}\rho dv

(1-3)  


To samo prawo zastosowane do strumienia wektora indukcji pola magnetycznego prowadzi do równanie (1-4),

 

\iint_{S}B\cdot nds=0

(1-4)  

gdyż pole magnetyczne jest bezźródłowe.
Do powyższych równań dodawane jest zapisane w formie całkowej równanie ciągłości prądu

  

\iint_{S}J\cdot nds=-\frac{\partial }{\partial t} \iiint_{V}\rho dv

(1-5)  

Mówi ono, że prąd przewodzenia wypływający przez powierzchnię S zamkniętej objętości V równy jest szybkości zmian ładunku w tej objętości.

1.3. Różniczkowa postać równań Maxwella

Równanie Maxwella w postaci całkowej opisują poparte eksperymentami zachowanie rozmaitych wielkości opisujących pole elektromagnetyczne. Ich interpretacja fizyczna jest czytelna. Z wielu względów wygodnie jest zapisać te równania w postaci różniczkowej. Otrzymujemy wtedy zestaw równań (1-6).

  

\triangledown \times E =-\frac{\partial B}{\partial t}

\triangledown \times H =J +\frac{\partial D}{\partial t}

\triangledown \cdot J =-\frac{\partial \varrho D}{\partial t}

\triangledown \cdot D = \varrho

\triangledown \cdot B = 0

(1-6)  

Pięć powyższych równań podawanych jest tradycyjnie jako równania Maxwella. Z matematycznego punktu widzenia występują związki między nimi, gdyż ostatnie dwa mogą być wyprowadzone w oparciu o pierwsze trzy. Z opisowego punktu widzenia pozostawimy w tym miejscu wszystkie z nich.
  

1.4. Równanie Maxwella w notacji zespolonej

Przyjmiemy teraz założenie, że natężenia pól elektrycznego i magnetycznego oraz prądy zmieniają się sinusoidalnie w czasie. Pozwala to wprowadzić notację zespoloną. Zespolone wektory E, D, H, B i J związane są z wektorami rzeczywistymi zależnościami w formie (1-7). 

 

A(x,y,z,t)=Re[A(x,y,z)e^{j\omega t})]

(1-7)  

W praktyce takie w przeważającej liczbie przypadków powyższe założenia są spełnione.

Pochodne po czasie wielkości zmiennych sinusoidalnie, czasami mówimy harmonicznie, w czasie obliczamy zgodnie z regułą (1-8): 

 

\frac{\partial e^{j\omega t}}{\partial t}=j\omega e^{j\omega t}

(1-8)  

Pozwala to zapisać równania Maxwella w innej, uproszczonej formie, w której usunięto zależność występujących w nich zmiennych od czasu. W notacji zespolonej podawane są zwykle tylko cztery z nich.

  

\triangledown \times \mathbf{E} =-j\omega \mathbf{B}

\triangledown \times \mathbf{H} =\mathbf{J} +j\omega \mathbf{D}

\triangledown \cdot \mathbf{D} = \varrho

\triangledown \cdot \mathbf{B} = 0

(1-9)  

Dla ośrodka jednorodnego, w którym nie ma ładunków i prądów przewodzenia równania powyższe upraszczają się do postaci (1-10):
 

  

\triangledown \times \mathbf{E} =-j\omega \mathbf{B}

\triangledown \times \mathbf{H} =j\omega \mathbf{D}

\triangledown \cdot \mathbf{D} = 0

\triangledown \cdot \mathbf{B} = 0

(1-10)  

1.5. Rodzaje ośrodków

W poprzednim punkcie założyliśmy, że pola elektryczne i magnetyczne występują w nieograniczonym ośrodku. Dokonajmy przeglądu ośrodków ze względu na własności elektryczne, które określają postać zależności między wektorami E, D, H, B i J. W szerokiej klasie ośrodków zależności między tymi wektorami można zapisać w postaci tzw. równań materiałowych:

  

D=\varepsilon E

B=\mu H

J=\sigma E

(1-12)  

Przypomnijmy, że wartości ε, μ, σ charakteryzujące ośrodek to przenikalność elektryczna, przenikalność magnetyczna oraz konduktywność, czyli parametry materiałowe ośrodka.
Związki (1-11) wskazują na liniową zależność wektorów indukcji oraz gęstości prądu od natężeń pól. Oznacza to, że wielkości ε, μ, σ nie zależą od natężeń pól i taki ośrodek nazywamy liniowym. Jeżeli przynajmniej jeden z parametrów ośrodka zależy od natężeń pól, to zależności (1-11) nie są już liniowe i ośrodek nazywa się nieliniowym.
Ośrodek jest jednorodny gdy jego parametry nie zależą od współrzędnych punktu. W przeciwnym przypadku mówimy o ośrodku niejednorodnym.
Istnieją ośrodki, których parametry ε, μ, σ zależą od częstotliwości. Ośrodki takie nazywamy dyspersyjnymi. Równania materiałowe w formie (1-11) mają sens dla ośrodków dyspersyjnych tylko w przypadku sinusoidalnej zależności pól od czasu. Dla dowolnej zależności pól od czasu związki (1-11) są słuszne dla transformat fourierowskich wektorów E, D, H, B, J.
Jeżeli wielkości ε, μ, σ są niezależne od kierunku pól, to parametry te są skalarami, a taki ośrodek nazywamy izotropowym. Odpowiednie wektory występujące w poszczególnych równaniach (1-11) są do siebie równoległe. 
Gdy parametry ośrodka zależą od kierunku pól, to mówimy o ośrodku anizotropowym, którego własności nie mogą być opisane przez skalarne wielkości ε, μ, σ. Równania (1-11) mogą być prawdziwe dla ośrodka anizotropowego, ale wtedy parametry materiałowe są reprezentowane przez tensory. Przykładowo, gdy związek między wektorami E i D zależy od kierunku wektora E to przenikalność elektryczna jest tensorem, [ε], reprezentowanym przez 9-elementową macierz, a relacja między wektorami D i E przybiera postać

  

\begin{bmatrix} D_{x}\\ D_{y}\\ D_{z} \end{bmatrix}=\begin{bmatrix} \varepsilon _{xx} & \varepsilon _{xy} & \varepsilon _{xz}\\ \varepsilon _{yxx} & \varepsilon _{yy} & \varepsilon _{yz}\\ \varepsilon _{zx} & \varepsilon _{zy} & \varepsilon _{zz} \end{bmatrix}\begin{bmatrix} E_{x}\\ E_{y}\\ E_{z} \end{bmatrix}

(1-12)  

Z powyższego równania wynika, że jedna składowa wektora pola elektrycznego wywołuje, w ogólności, trzy składowe wektora indukcji elektrycznej.
Struktury krystaliczne, np. półprzewodniki, oraz zjonizowane gazy są przykładami anizotropowych dielektryków.
Analogiczną do (1-12) relację można określić dla związku między wektorami B i H, w którym przenikalność magnetyczna jest tensorem. W technice mikrofalowej stosuje się ferryty, które są anizotropowymi materiałami magnetycznymi.

1.6. Dielektryki i magnetyki

Materiały dielektryczne powszechnie występują w technice mikrofalowej w konstrukcji tak linii transmisyjnych jak i elementów oraz podzespołów. Rozważmy zachowanie izotropowego dielektryka w zewnętrznym polu elektrycznym o sinusoidalnej zależności od czasu, stosując wektory zespolone.
Pod wpływem pola elektrycznego cząsteczki materii ustawiają się zgodnie z prawem Coulomba tak, że wywołują własne pole elektryczne skierowane przeciwnie do pola zewnętrznego. Zjawisko to określamy mianem elektrycznej polaryzacji (lub krótko polaryzacją) ośrodka i opisujemy przez wektor polaryzacji Pe. Wektory Pe, E i D powiązane są ze sobą zależnością

  

\mathbf{D}=_{0}\mathbf{E}+\mathbf{P}_{e}

(1-13)  

W ośrodku liniowym, polaryzacja elektryczna jest liniowo związana z zewnętrznym polem elektrycznym jako

  

\mathbf{P}_{e}=\varepsilon _{0}\chi_{e}\mathbf{E}

(1-14)  

przy czym χe, która może być zespolona, nazywa się podatnością elektryczną. Tak więc:

  

\mathbf{D}=\varepsilon _{0}\mathbf{E}+\mathbf{P_{e}}=\varepsilon _{0}(1+\chi _{e})\mathbf{E}=\varepsilon \mathbf{E}

(1-15)  

gdzie:

  

\varepsilon =\varepsilon _{0}(1+\chi _{e})=\varepsilon" -j\varepsilon""

(1-16)  

jest zespoloną przenikalnością elektryczną ośrodka.
Zespolony zapis dla przenikalności elektrycznej (1-16) może być stosowany tylko w przypadku pól o sinusoidalnej zależności od czasu. Przenikalność elektryczną znormalizowaną względem przenikalności elektrycznej próżni nazywamy zespoloną względną przenikalnością elektryczną:

  

\varepsilon _{w}=\frac{\varepsilon }{\varepsilon 0}=1+\chi _{e}=\varepsilon_{w} "-j\varepsilon_{w} ""

(1-17)  

W przypadku statycznym lub przy wolnych zmianach pola wektor polaryzacji elektrycznej jest w przybliżeniu w fazie z wektorem pola elektrycznego. Dla pól szybkozmiennych opóźnienie wektora Pe względem wektora E nie jest pomijalnie małe, wymienione wektory mają różne fazy. Część urojona ε, która opisuje straty w ośrodku (grzanie) wywołane tłumieniem wibracji dipoli i tym samym opóźnienie wektora polaryzacji względem wektora pola elektrycznego, musi być ujemna (ε” jest dodatnie). 
Źródłem strat w dielektryku może być również niezerowa konduktywność ośrodka i wtedy istnieje w ośrodku wektor gęstości prądu przewodzenia opisany zależnością:

  

\mathbf{J}=\sigma \mathbf{E},

(1-18)  

która jest prawem Ohma z punktu widzenia pola elektromagnetycznego.
Całkowita gęstość prądu (przewodzenia i przesunięcia) w dielektryku, która występuje po prawej stronie drugiego równania Maxwella (1-9), wynosi

  

\mathbf{J}_{C}=\mathbf{J}+j\omega \varepsilon \mathbf{E}=\sigma \mathbf{E}+j\omega (\varepsilon "-j\varepsilon "")\mathbf{E}= j\omega \varepsilon" \mathbf{E}+(\sigma +\omega \varepsilon "")

=j\omega [\varepsilon "-j(\varepsilon ""+\frac{\sigma }{\omega })]\mathbf{E}

(1-19)  

Z zależności (1-17) wynika, że straty wynikające z tłumienia oscylacji dipoli (ωε”) 
i straty wynikające z istnienia prądu przewodzenia (σ) są nierozróżnialne. Wielkość σ+ωε” można traktować jako zastępczą konduktywność ośrodka. Z kolei wielkość ε”+σ/ω można określić jako zastępczą część urojoną przenikalności elektrycznej.
Stratność dielektryka możemy charakteryzować przez podanie tangensa kąta stratności wyrażonego wzorem

  

\mathrm{tg} \delta _{\varepsilon }=\frac{\sigma +\omega \varepsilon ""}{\omega \varepsilon "},

(1-20)  

który jest stosunkiem składowej prądu będącego w fazie z polem elektrycznym (ten prąd wywołuje straty mocy fali elektromagnetycznej w ośrodku) do składowej prądu proporcjonalnej do ωε’, odpowiedzialnej za magazynowanie energii pola elektrycznego. 
Zwykle dielektryki charakteryzuje się podając stałą dielektryczną, która jest rzeczywistą częścią względnej przenikalności elektrycznej, oznaczaną w literaturze najczęściej jako εw (z pominięciem indeksu „’ ”, co może wprawiać w zakłopotanie) oraz tangensa kąta stratności dla określonej częstotliwości.
Dla izotropowych magnetyków można zastosować podobny opis jak dla dielektryków. Zewnętrzne pole magnetyczne porządkuje magnetyczne dipole w materiale wywołując wektor magnetyzacji ośrodka opisany wektorem Pm związanym w ośrodku liniowym z zewnętrznym polem magnetycznym zależnością

  

\mathbf{P}_{m}=\mu _{0}\chi _{m}\mathbf{H}

(1-21)  

gdzie zespolona wielkość χm nazywa się podatnością magnetyczną.
Wektor indukcji magnetycznej w izotropowym magnetyku wyraża relacja 

  

\mathbf{B}=\mu _{0}\mathbf{H}+\mathbf{P}_{m}=\mu _{0}(1+\chi _{m})\mathbf{H}=\mu \mathbf{H}

(1-22)  

gdzie

  

\mu =\mu _{0}(1+\chi _{m})=\mu "-j\mu ""

(1-23)  

jest zespoloną przenikalnością magnetyczną ośrodka.
Przenikalność magnetyczną znormalizowaną względem przenikalności magnetycznej próżni nazywamy zespolona względną przenikalnością magnetyczną:

  

\mu_{w} =\frac{\mu }{\mu _{0}}=1+\chi _{m}=\mu_{w} "-j\mu_{w} ""

(1-24)  

Analogicznie jak dla dielektryków, część urojona μ lub χm opisuje straty w magnetyku wywołane tłumieniem wibracji dipoli magnetycznych i wektor polaryzacji magnetycznej jest opóźniony względem wektora pola magnetycznego. W przypadku ośrodków magnetycznych nie ma magnetycznej konduktancji, ponieważ nie ma magnetycznego prądu przewodzenia.
Przy analizie własności ośrodków magnetycznych uwzględnia się straty magnetyczne wprowadzając, analogicznie jak poprzednio, tangens kąta stratności

  

\mathrm{tg} \delta _{m }=\frac{\mu ""}{\mu "},

(1-25)  

Tak więc ośrodki magnetyczne możemy opisać podając część rzeczywistą przenikalności magnetycznej oraz tangens kąta stratności.
Równania Maxwella w formie różniczkowej (1-9) wymagają znajomości wartości brzegowych, aby uzyskać jednoznaczne rozwiązanie w danym ośrodku. 

1.7. Warunki brzegowe dla pola EM

W wielu zagadnieniach występujących w technice mikrofalowej mamy do czynienia z układem dwóch dielektryków bądź ze strukturą dielektryk – przewodnik. Rozważmy zachowanie pola elektromagnetycznego w tych szczególnych przypadkach.
Na granicy miedzy dwoma dielektrykami zwykle nie występuje ładunek powierzchniowy oraz nie płynie prąd przewodzenia. Dokładne wyprowadzenia prowadzą do wniosku, że składowe normalne wektorów D i B oraz składowe styczne wektorów E i H są ciągłe na granicy dielektryków.    
Przyjmując, że składowa normalna jest równoległa do jednej osi w trójwymiarowym układzie współrzędnych to powyższe warunki określają zachowanie sześciu składowych pola elektromagnetycznego (składowa styczna składa się z dwóch składowych). Można wykazać, że spełnienie warunków brzegowych dla składowych stycznych na granicy dwóch ośrodków dielektrycznych pociąga za sobą automatycznie spełnienie warunków ciągłości dla składowych normalnych.
W liniach transmisyjnych powszechnie występuje granica między dielektrykiem i dobrym przewodnikiem, który często można uznać za idealny. Warunki brzegowe dla pól zmiennych w czasie panujące na granicy dielektryka i idealnego przewodnika zapisuje się następująco:

 

\mathbf{n} \times \mathbf{E} =0

\mathbf{n} \times \mathbf{H} =\mathbf{J}_s

\mathbf{n} \cdot \mathbf{D_{2}} = \varrho_s

\mathbf{n}\cdot \mathbf{B}_2 = 0

(1-26)
(1-27)
(1-28)
(1-29)
 

W idealnym przewodniku pole elektryczne musi być równe zeru, gdyż w przeciwnym przypadku wywoływałoby prąd przewodzenia o nieskończonym natężeniu. Istnienie w tym ośrodku zmiennego pola magnetycznego jest również niemożliwe, ponieważ zgodnie z pierwszym prawem Maxwella musiałoby wywołać pole elektryczne, co byłoby sprzeczne z poprzednim stwierdzeniem.

Rys.1.4. Pola i prąd powierzchniowy na granicy metalu i dielektryka  

 
W rezultacie dochodzimy do wniosku, że składowa styczna pola elektrycznego na granicy idealnego przewodnika jest równa zeru, a więc pole elektryczne musi być prostopadłe do powierzchni przewodnika. Pole elektryczne indukuje na powierzchni przewodnika powierzchniowy ładunek elektryczny o gęstości równej wartości indukcji elektrycznej. Składowa normalna indukcji magnetycznej jest zerowa na brzegu idealnego przewodnika, czyli pole magnetyczne musi być styczne do granicy przewodnika i wywołuje na jego powierzchni prąd przewodzenia o gęstości równej wartości natężenia pola magnetycznego.
 

1.8. Równania falowe w dielektryku bezstratnym

Rozważamy przestrzeń nieograniczoną wypełnioną ośrodkiem liniowym, izotropowym, niedyspersyjnym, jednorodnym i bezstratnym. Zakładamy, że w ośrodku nie ma prądów 
i ładunków co oznacza, że wszelkie źródła pól są nieskończenie daleko od rozważanego obszaru. 
Parametry ośrodka ε i μ są liczbami stałymi i układ równań Maxwella (1-6) sprowadza się do czterech następujących zależności

  

\triangledown \times E =-\mu\frac{\partial H}{\partial t}

\triangledown \times H =\varepsilon\frac{\partial E}{\partial t}

\triangledown \cdot J =-\frac{\partial \varrho D}{\partial t}

\triangledown \cdot E = 0

\triangledown \cdot H = 0

(1-30)  

W celu rozwiązania układu równań (1-30), poddajemy pierwsze równanie układu obustronnie rotacji i wstawiamy do uzyskanej zależności drugie z równań (1-30), tak więc

  

 

\triangledown \times\triangledown \times E=-\mu \frac{\partial }{\partial t}\triangledown \times H=-\mu \varepsilon \frac{\partial^2 E}{\partial t^2}

(1-31)  

Lewą stronę równania (1-31) przekształcamy zgodnie z poniższą tożsamością wektorową 

  

\triangledown \times\triangledown \times E=\triangledown (\triangledown \cdot E)-\triangledown ^{2}E

(1-32)  

i uwzględniamy to, że dywergencja pola E jest równa zeru (trzecie równanie w układzie 
(1-30)). W wyniku uzyskujemy równanie falowe dla pola elektrycznego: 

  

{\color{red} {\triangledown ^{2}E-\mu \varepsilon \frac{\partial^2 E}{\partial t^2}=0;}}

(1-33)  

Podobnie, eliminując z układu równań Maxwella wektor E, otrzymujemy równanie falowe spełniane przez wektor pola magnetycznego. Postać tego równania jest analogiczna do równania (1-33).

  

{\color{red} {\triangledown ^{2}H-\mu \varepsilon \frac{\partial^2 H}{\partial t^2}=0;}}

(1-34)  

Należy zaznaczyć, że układ równań Maxwella nie jest równoważny układowi równań falowych. Równania falowe wynikają z równań Maxwella, ale wynikanie odwrotne nie zachodzi. Z tego względu każde rozwiązanie układu równań Maxwella musi spełniać równania falowe, a wśród wektorów spełniających równania falowe mogą być takie, które nie są rozwiązaniami układu równań Maxwella.
Analiza rozwiązania powyższych równań falowych prowadzi do wniosku, że prędkość  rozchodzenia się zaburzenia w przestrzeni określona jest iloczynem \mu \varepsilon  , gdyż  v=1/\sqrt{\mu \varepsilon} . W próżni jest to prędkość światła  v=1/\sqrt{\mu_0 \varepsilon}_0 .

1.9. Równania falowe w dielektryku stratnym

Utrzymajmy w mocy wszystkie założenia dotyczące ośrodka wypełniającego przestrzeń (liniowość, izotropowość, jednorodność), ale przyjmijmy, że straty występujące w ośrodku opisuje konduktywność σ różna od zera. W tym przypadku układ równań Maxwella przyjmuje postać 

  

\triangledown \times E =-\mu\frac{\partial H}{\partial t}

\triangledown \times H =\sigma E+\varepsilon\frac{\partial E}{\partial t};

\triangledown \cdot E = 0

\triangledown \cdot H = 0

(1-35)  

Przekształcając równania Maxwella, metodami analogicznymi do stosowanych w poprzednim punkcie, można uzyskać równania falowe dla ośrodka stratnego
 

  

\triangledown^{2} E -\mu \sigma \frac{\partial E}{\partial t}-\mu \varepsilon\frac{\partial^2 E}{\partial t^2}=0;

 

\triangledown^{2} H -\mu \sigma \frac{\partial H}{\partial t}-\mu \varepsilon\frac{\partial^2 H}{\partial t^2}=0;

(1-36)  

Rozwiązanie tych równań w przypadku dowolnej zależności pól od czasu jest skomplikowane
i nie prowadzi do tak prostej interpretacji fizycznej i geometrycznej jak w przypadku ośrodków bezstratnych. Celowym jest przyjęcie, że pola tworzące fale elektromagnetyczną są cosinusoidalnymi funkcjami czasu o ustalonej pulsacji i zastosowanie rachunku zespolonego (w dziedzinie częstotliwości).
Przekształcając równania falowe (1-36) do postaci zespolonej uzyskujemy następujące równania falowe w dziedzinie zespolonej (równania Helmholtza):

  

\triangledown^{2} \mathbf{E}-\gamma ^{2}\mathbf{E}=0;

\triangledown^{2} \mathbf{H}-\gamma ^{2}\mathbf{H}=0;

(1-37)  

W równaniach pojawia się zmienna zespolona  zwana współczynnikiem propagacji, o fundamentalnym znaczeniu dla opisu zjawiska propagacji fali

  

\gamma ^{2}=j\omega \mu (\sigma +j\omega \varepsilon );

(1-38)  

Pamiętamy, że współczynnik propagacji fali w nieograniczonym ośrodku określają wielkości charakteryzujące ten ośrodek.

1.10. Pola fali płaskiej

Fala wypromieniowana przez antenę nadawczą rozchodzi się we wszystkich kierunkach, choć oczywiście pewne kierunki są uprzywilejowane, o czym mówi charakterystyka anteny. Z punktu widzenia anteny odbiorczej oddalonej odpowiednio daleko, o wiele kilometrów, fala dochodząca do niej jest falą płaską. Wartości chwilowe wektorów pól elektrycznego i magnetycznego tej fali są takie same w każdym punkcie płaszczyzny prostopadłej do kierunku rozchodzenia się fali. Na podstawie równania falowego możemy znaleźć matematyczny opis takiej fali.
Rozważamy przypadek, w którym składowe pola są harmonicznymi funkcjami czasu. Wykorzystujemy równania Helmholtza w postaci (1-37). Możemy je zapisać w kartezjańskim układzie współrzędnych w formie (1-39):

  

\frac{\partial^2 \mathbf{E}}{\partial x^2}+\frac{\partial^2 \mathbf{E}}{\partial y^2}+\frac{\partial^2 \mathbf{E}}{\partial z^2}-\gamma^{2}\mathbf{E}=0

\frac{\partial^2 \mathbf{H}}{\partial x^2}+\frac{\partial^2 \mathbf{H}}{\partial y^2}+\frac{\partial^2 \mathbf{H}}{\partial z^2}-\gamma^{2}\mathbf{H}=0

(1-39)  

Pierwsze równanie musi być spełnione dla każdej ze składowych Ex, Ey i Ez pola elektrycznego, a drugie dla analogicznych składowych pola magnetycznego. 
Stwierdzamy, nie wdając się w matematyczne uzasadnienia, że dla fali płaskiej natężenie pola elektrycznego E i magnetycznego H są wektorami poprzecznymi ET i HT do kierunku propagacji, i mówimy, że fala płaska jest falą poprzeczną czyli typu TEM (ang. Transverse Electro- Magnetic).
Przyjmijmy, że fala rozchodzi się wzdłuż osi z. Wektory ET i HT leżą w płaszczyźnie xy i są niezależne od współrzędnych x i y, co oznacza    

  

\frac{\partial \mathbf{E}_{T}}{\partial x}=\frac{\partial \mathbf{E}_{T}}{\partial y}=0;\, \, \, \frac{\partial \mathbf{H}_{T}}{\partial x}=\frac{\partial \mathbf{H}_{T}}{\partial y}=0;

(1-40)  

Pamiętajmy, że dla omawianego przypadku

  

\mathbf{E}_{T}=\mathbf{i}_{x}E_{x}+\mathbf{i}_{y}E_{y};\, \, \, \mathbf{H}_{T}=\mathbf{i}_{x}H_{x}+\mathbf{i}_{y}H_{y};

(1-41)  

Równania (1-39) upraszczają się do

  

\frac{\partial^2 \mathbf{E}_T}{\partial z^2}-\gamma^{2}\mathbf{E}_T=0;

\frac{\partial^2 \mathbf{H}_T}{\partial z^2}-\gamma^{2}\mathbf{H}_T=0;

(1-42)  

Rozwiązania równań (1-42) mają prostą dwuczłonową postać (1-43):

  

\mathbf{E}_T(z)=\mathbf{E}^{+}_{T0} e^{-\gamma z}+\mathbf{E}^{-}_{T0} e^{\gamma z};

\mathbf{H}_T(z)=\mathbf{H}^{+}_{T0} e^{-\gamma z}+\mathbf{H}^{-}_{T0} e^{\gamma z};

(1-43)  

Obecność dwóch członów w rozwiązaniach, to dwie fale rozchodzące się wzdłuż osi z:
•    fala postępująca propagowana zgodnie z kierunkiem osi z o czym informuje czynnik e-γz;
•    fala powracająca rozchodzi się w kierunku przeciwnym do osi z co opisuje eγz.
Przyjmijmy, że mamy tylko falę postępującą. Gdy znamy jeden z wektorów fali płaskiej, pola elektrycznego albo pola magnetycznego, to drugi możemy wyznaczyć posługując się zależnościami:

  

\mathbf{E}_T=Z_{f}\mathbf{H}_T \times \mathbf{i}_z;

\mathbf{H}_T=\frac{1}{Z_{f}}\mathbf{i}_z \times \mathbf{E}_T;

(1-44)  

Współczynnik występujący w pierwszym z równań (1-44) zdefiniowany jako stosunek wartości wzajemnie prostopadłych (prostopadłych do kierunku rozchodzenia się fali) składowych pola elektrycznego i magnetycznego nazywamy impedancją falową

  

Z_{f}=\frac{E_{T}}{H_{T}}=\frac{E_{x}}{H_{y}}=-\frac{E_{y}}{H_{x}}

(1-45)  

Dla fali płaskiej impedancja falowa jest dana zależnością

  

Z_{f}=\sqrt{\frac{j\omega \mu }{\sigma +j\omega \varepsilon }};

(1-46)  

Równocześnie wielkość zespolona zdefiniowana wzorem

  

Z_{w}=\sqrt{\frac{j\omega \mu }{\sigma +j\omega \varepsilon }};

(1-47)  

ma wymiar impedancji, zależy od parametrów ośrodka oraz pulsacji i nazywa się impedancją właściwą ośrodka. Dla ośrodka bezstratnego impedancja właściwa jest rzeczywista, zależy tylko od parametrów ośrodka i dla próżni równa jest

  

Z_{w0}=\sqrt{\frac{\mu_0 }{\varepsilon }_0}= 120\pi \, \, [\Omega ]

(1-48)  

Jak wynika z równości (1-46) oraz (1-47), w przypadku fali typu TEM w nieograniczonej przestrzeni impedancja falowa jest równa impedancji właściwej ośrodka. Cecha ta nie jest własnością wszystkich fal elektromagnetycznych, a jedynie fal typu TEM.

Rys.1.5. Składowe pola elektrycznego i magnetycznego fali płaskiej w ośrodku bezstratnym. 

   
Ze zależności (1-44) wynika, że wektory natężenia pól elektrycznego i magnetycznego są do siebie prostopadłe oraz, że trójka wektorów ET, HT, iz jest prawoskrętna. Własność tę mają na ogół tylko fale w ośrodkach nieograniczonych i izotropowych.
Gdy położenie osi x zostanie tak dobrane, by pole elektryczne miało tylko składową ET=Ex, to wtedy pole magnetyczne ma jedynie składową HT=Hy (zgodnie z równaniami (1-44)). Przykład pól elektrycznego i magnetycznego takiej fali płaskiej ilustruje rys.1.5.
 

1.11. Parametry fali płaskiej

Dla fali rozchodzącej się w ośrodku izotropowym o parametrach ε, μ, σ zgodnie z kierunkiem osi z mającej tylko składową Ex, której wartość jest równa liczbie rzeczywistej E0 w płaszczyźnie z = 0, pola zapiszemy w dziedzinie zespolone następująco

  

E_{x}(z)=E_{0}e^{-\gamma z};

(1-49a)  
  

H_{y}(z)=\frac{E_{0}}{|Z_{w}|}e^{-j\varphi }e^{-\gamma z};

(1-49b)  

gdzie impedancję właściwą zapisano w postaci wykładniczej Zw=|Zw|e. Współczynnik propagacji \gamma , zwany niekiedy stałą propagacji, został opisany równaniem (1-38). W ogólnym przypadku \gamma  jest liczbą zespoloną:

  

\gamma =\alpha +j\beta ;

(1-50)  

 Część rzeczywistą współczynnika propagacji \alpha nazywamy stałą tłumienia, a część urojoną \beta jest stałą fazową i w ośrodku stratnym

  

\alpha =\mathrm {Re}\begin{Bmatrix} {\sqrt{j\omega \mu (\sigma +j\omega \varepsilon )}} \end{Bmatrix};

\beta =\mathrm {Im}\begin{Bmatrix} {\sqrt{j\omega \mu (\sigma +j\omega \varepsilon )}} \end{Bmatrix}=\mathrm {Im}\begin{Bmatrix} {j\omega \sqrt{ \mu \varepsilon}\sqrt{1+\frac{\sigma }{j\omega \varepsilon }}} \end{Bmatrix};

(1-51a)
(1-51b)
 

W próżni, gdy przewodność ośrodka σ=0, a ε=ε0 i μ=μ0, \gamma jest czysto urojone.

  

\gamma ^{2}=-\omega ^{2}\varepsilon _{0}\mu _{0};

\beta =\omega \sqrt{\varepsilon _{0}\mu _{0}};

(1-52a)
(1-52b)
 

Korzystając z definicji (1-50) zależności (1-49) przybierają postać

  

{\color{Red} {E_{x}(z)=E_{0}e^{-\alpha z}e^{-j\beta z};}}

{\color{Red} {H_{y}(z)=H_{0}e^{-\alpha z}e^{-j(\beta z+\varphi )};} }

(1-53)  

Jak wiemy przejście od notacji zespolonej do rzeczywistej jest proste:

  

 E_{x}(z,t)=\mathrm {Re}\begin{Bmatrix} E_{0}e^{-\alpha z} e^{-j\beta z} e^{-j\omega t} \end{Bmatrix}=E_{0}e^{-\alpha z} \cos (\omega t-\beta z)

 H_{y}(z,t)=\mathrm {Re}\begin{Bmatrix} H_{0}e^{-\alpha z} e^{-j(\beta z+\varphi )} e^{-j\omega t} \end{Bmatrix}=H_{0}e^{-\alpha z} \cos (\omega t-\beta z-\varphi )

(1-54a)
(1-54b)
 

Z zapisu pól podanych w (1-54) widać, że stała tłumienia wpływa jedynie na zmniejszanie się amplitudy danego pola, a stała fazowa informuje o tempie zmiany fazy fali wzdłuż osi z. Gdy fala rozchodzi się w ośrodku bezstratnym amplitudy pól nie zmniejszają i pole elektryczne jest w fazie z polem magnetycznym (φ=0).
Prędkość fazową fali płaskiej znajdujemy analizując ruch płaszczyzny stałej fazy. Dla tej płaszczyzny spełniony jest warunek (1-55):

  

\omega t-\beta z=\mathrm{const};

(1-55)  

Płaszczyzna ta porusza się z prędkością vf:

  

{\color{Red} {v_{f}=\frac{\omega }{\beta };}}

(1-56)  

Zwróćmy uwagę, że stała fazowa dla fali propagowanej w ośrodku stratnym opisana relacją 
(1-51b) nie jest liniową funkcją ω i prędkość fazowa zmienia się z częstotliwością. Mówimy wtedy o dyspersji fali, która w tym wypadku wynika ze strat ośrodka.
Dla próżni prędkość fazowa fali równa jest prędkości c światła:

  

 v_{f}=c=\frac{1}{\sqrt{\mu_0 \varepsilon_0 }};

(1-57)  

Droga, jaką fala poruszająca się prędkością opisaną zależnością (1-56) przebędzie w czasie okresu T, nazywa się długością fali \lambda .

  

\lambda =\frac{\omega T}{\beta }=\frac{2\pi }{\beta };

(1-58)  

Zależność (1-56) podaje wartość prędkości fazowej, ale pamiętajmy, że prędkość jest wielkością wektorową i jej zwrot jest zgodny z wektorem określającym kierunek rozchodzenia się fali. Należy podkreślić, że ruch płaszczyzny stałej fazy jest pojęciem matematycznym i nie oznacza przesuwania się żadnego obiektu materialnego, a tym samym przenoszenia energii. Z tego względu prędkość fazowa może przyjmować dowolne wartości dodatnie i nie podlega ograniczeniom wynikającym ze szczególnej teorii względności.
Aby za pomocą fali EM przesłać informację, trzeba przebieg o pulsacji nośnej ω zmodulować odpowiednim sygnałem. W wyniku modulacji z sygnału monochromatycznego otrzymujemy widmo częstotliwości. W najprostszym przypadku modulacji amplitudowej jedną częstotliwością otrzymujemy dwie częstotliwości prążków bocznych, różniące się od nośnej o Δω. Superpozycję kilku fal o zbliżonych częstotliwościach, Δω->0, i współczynnikach fazowych, Δβ->0, zwana jest fizyce grupą fal.
Obserwując ruch płaszczyzny stałej fazy obwiedni (np. płaszczyzny, w której superpozycja dwóch fal osiąga maksimum) można zapisać warunek (1-59):

  

\Delta \omega t-\Delta \beta z=\mathrm{const};

(1-48)  

(1-59)
Prędkość grupowa vg to prędkość poruszania się obwiedni sygnału: 

  

{\color{Red} {v_{g}=\frac{\mathrm{d} \omega }{\mathrm{d} \beta };}}

(1-60)  

Można wykazać, że w ośrodku izotropowym związek pomiędzy prędkością grupową i fazową przybiera postać:

  

 v_{g}=\frac{v_{f}}{1-\frac{\omega }{v_{f}}\frac{\mathrm{d} v_{f}}{\mathrm{d} \omega x}};

(1-61)  

Ze związku (1-61) wynika, że dla fali, której prędkość fazowa nie zależy od częstotliwości, prędkość grupowa jest równa prędkości fazowej.

1.12. Rodzaje fal EM w prowadnicach falowych

Do tej pory zajmowaliśmy się falą płaską rozchodzącą się w nieograniczonej przestrzeni, która jest ważnym rodzajem fali elektromagnetycznej w dziedzinie zastosowań telekomunikacyjnych bo występuje pomiędzy antenami nadawczą i odbiorczą.
Niemniej istotne jest rozchodzenie się fal elektromagnetycznych we wszelkiego rodzaju liniach transmisyjnych, w których fale prowadzone są w określonym kierunku. Linie te nazywamy prowadnicami falowymi. Prowadzenie fal uzyskuje się wzdłuż określonego układu przewodników lub dielektryków (tzn. w obszarach cylindrycznych, których granicę są przewodzące lub są to granice dwóch dielektryków, ewentualnie odpowiednia kombinacja wymienionych materiałów). W przypadku prowadzenia fal wzdłuż przewodników, możliwe jest rozchodzenie się energii elektromagnetycznej w liniach składających się z dwóch lub więcej przewodów, a także w rurach (najczęściej o przekroju prostokątnym albo kołowym) nie zawierających wewnątrz dodatkowych przewodników, tzw. falowodach. Natomiast przykładem prowadnicy falowej będącej układem warstw dielektrycznych jest światłowód, bez którego trudno wyobrazić sobie dzisiejszą telekomunikację.
Określenie postaci fali elektromagnetycznej w prowadnicy falowej wiąże się z poszukiwaniem rozwiązań równań Maxwella, które jest zagadnieniem matematycznym innego typu niż w przypadku fali w nieograniczonej przestrzeni. Ze względu na to, że mamy tu do czynienia z obszarem cylindrycznym o granicy przewodzącej lub w formie granicy dwóch dielektryków, musimy teraz uwzględnić warunki brzegowe. Dlatego własności fal elektromagnetycznych w prowadnicach są inne niż dla fali płaskiej.
Okazuje się, że fale w prowadnicach falowych nie muszą być falami typu TEM, tzn. mogą one mieć składowe pól elektrycznego i magnetycznego wzdłuż kierunku rozchodzenia się fali. Wprowadzić należy klasyfikację możliwych rodzajów fal nazywanych również modami.
Przyjmijmy, że fala rozchodzi się zgodnie z kierunkiem osi z i wtedy wyróżnia się następujące typy fal:


Tabela 1.1. Zestawienie właściwości typów fal.

Typ fali Właściwości pola E i H
Fala typu TEM Ez=0 – pole elektryczne leży w płaszczyźnie prostopadłej do kierunku propagacji fali, 
Hz=0 – pole magnetyczne leży w płaszczyźnie prostopadłej do kierunku propagacji fali;
Fala typu TM (zwana też E) Ez≠0,
Hz=0 – pole magnetyczne leży w płaszczyźnie prostopadłej do kierunku propagacji fali;
Fala typu TE (zwana też H) Hz≠0, 
Ez=0 – pole elektryczne leży w płaszczyźnie prostopadłej do kierunku propagacji fali;
Fala typu EH Ez≠0,
Hz≠0

Z powyższego wykazu wynika, że tylko pierwszy z wymienionych typów fal jest falą poprzeczną. Prowadnice falowe, w których mogą rozchodzić się rodzaje TEM nazywamy liniami TEM lub prowadnicami TEM. Struktura prowadnicy TEM musi zawierać co najmniej dwa przewody. Przykładem linii TEM jest linia współosiowa, tzw. kabel koncentryczny.
Fale E i H rozchodzą się w falowodach. Falowody stosuje się do prowadzenia fali elektro-magnetycznej z mniejszymi stratami niż w linii TEM (np. w transponderach satelitów telekomunikacyjnych) lub do przesyłania dużych mocy, których przesłanie nie jest możliwe linią współosiową (np. w radarach).
Fale typu EH występują między innymi w falowodach dielektrycznych i światłowodach.

2. Prowadnice TEM

Można wykazać, że zależności pól w prowadnicy TEM od zmiennej z są identyczne jak dla fali płaskiej w przestrzeni nieograniczonej. Podobnie możemy uzyskać równanie falowe określające wektory ET lub HT w postaci (1-42). Oznacza to, że wektory pól E i H są do siebie prostopadłe i prostopadłe do kierunku propagacji z.
Współczynnik propagacji fali w linii TEM jest taki sam jak w nieograniczonym ośrodku o parametrach materiału wypełniającego prowadnicę falową, czyli

  

\gamma =\sqrt{j\omega \mu (\sigma +j\omega \varepsilon )}=\alpha +j\beta ;

(1-62)  

Dodatkowo, relacje (1-44) między wektorami pól elektrycznego i magnetycznego obowiązują dla prowadnic TEM, czyli 

  

\mathbf{E}_T=Z_{f}\mathbf{H}_T \times \mathbf{i}_z;

\mathbf{H}_T=\frac{1}{Z_{f}}\mathbf{i}_z \times \mathbf{E}_T;

(1-63a)
(1-63b)
 

oraz spełniona jest poniższa relacja

  

Z_{f}=Z_{w}=\sqrt{\frac{j\mu }{\sigma +j\omega \varepsilon }}.

(1-64)  

 Zasadnicza różnica między cechami pół fali płaskiej w ośrodku nieograniczonym i pól w linii TEM jest związana z tym, że pola w ośrodku nieograniczonym nie zależą od zmiennych w płaszczyźnie prostopadłej do kierunku rozchodzenia się fali, natomiast w prowadnicy TEM, w której muszą być spełnione określone warunki brzegowe na powierzchni przewodników linii, pola na ogół zależą od tych zmiennych.
Warto zapamiętać, że prowadnicę falową charakteryzują dwa parametry: współczynnik propagacji \gamma oraz impedancja charakterystyczna Z0. Pierwszy z tych parametrów jest wielkością polową, której obliczenie wiąże się w ogólności z rozwiązaniem równań Maxwella. Drugi jest wielkością obwodową, wyznaczaną z zastosowaniem definicji (1-65):

  

Z_{0}=\frac{U}{I};

(1-65)  

w której U i I są amplitudami napięcia i prądu fali poruszającej się w jedną stronę.
Definicja Z0 jest przydatna przy analizie obwodów zawierających prowadnice falowe i elementy reprezentowane przez układy zastępcze o stałych skupionych.

2.1. Linia współosiowa

Najpopularniejszą prowadnicą w rodzinie TEM jest linia współosiowa o promieniach a i b (a>b), w której przestrzeń między przewodem wewnętrznym i zewnętrznym wypełniona jest małostratnym dielektrykiem o przenikalności względnej \varepsilon_r – rys.1.6A. 
Impedancja charakterystyczna linii współosiowej obliczana jest z zależności:

  

Z_{0}[\Omega ]=\frac{138}{\sqrt{\varepsilon _{r}}}\log \frac{b}{a};

(1-66)  

Rys.1.6. Struktura linii współosiowej. 
A) Przekrój poprzeczny linii współosiowej wypełnionej dielektrykiem. B) Linia współosiowa o okresowo rozmieszczonych podporach dielektrycznych.    

Zależność (1-66) wskazuje, że impedancja charakterystyczna linii współosiowej zależy od stosunku promieni przewodów i właściwości ośrodka wypełniającego prowadnicę. Wartość Z0 można ustalić dobierając stosunek średnic przewodów zewnętrznego i wewnętrznego, bądź przez dobór dielektryka (stała εr). Analiza prowadzi nas do wniosku, że dla tych samych materiałów tłumienie linii zależy od wartości Z0 i dla linii współosiowej powietrznej osiąga najmniejszą wartość dla Z0=77 Ω (rys.1.7). Z tego względu współosiowe kable telekomunikacyjne (m.in. sieci telewizji kablowej) mają impedancję charakterystyczną równą 75Ω. Problem minimalizacji tłumienia kabli współosiowych skłonił konstruktorów do szukania struktur pozwalających minimalizować ilość dielektryka w linii. Taką konstrukcję pokazano na rys. 1.6.B. Centryczność przewodu środkowego zapewniają okresowo rozłożone cienkie podpory.
W laboratoriach powszechnie używane są kable o impedancji Z0=50 Ω. 

Rys.1.7. Zależności tłumienia (linia czerwona) i największej mocy przenoszonej (linia niebieska) od impedancji charakterystycznej  linii współosiowej

    
Linia współosiowa, albo koncentryczna jest szeroko stosowana w systemach pomiarowych, a rozpowszechnionym w aparaturze standardem jest linia o impedancji Z0=50 Ω. Linie współosiowe pracują do 60 GHz.

2.2. Linia dwuprzewodowa

Na rys.1.8A przedstawiono strukturę innej linii TEM, a mianowicie linii dwuprzewodowej. Przewody zanurzone są w dielektryku o przenikalności \epsilon_r. Polowe wielkości charakteryzujące falę TEM dla tej prowadnicy są identyczne jak dla linii współosiowej. 
 
Rys.1.8. Przekroje poprzeczne i wymiary prowadnic TEM.  A) Linii dwuprzewodowa, B) Symetryczna linia paskowej.


Impedancję charakterystyczną linii dwuprzewodowej określa zależność:

  

Z_{0}[\Omega ]=\frac{276}{\sqrt{\varepsilon _{r}}}\log \frac{s}{a};

(1-67)  

Linia dwuprzewodowa jest z historycznego punktu widzenia pierwszą linią długą, dla której znaleziono rozwiązanie falowe. Stosowana jest jeszcze w sieciach telewizyjnych i telefonicznych w postaci tzw. skrętki.

2.3. Symetryczna linia paskowa

 Strukturę symetrycznej linii paskowej pokazano na rys.1.8B. Impedancję charakterystyczną tej linii oblicza się ze wzoru: 

  

Z_{0}[\Omega ]=\frac{30\pi}{\sqrt{\varepsilon _{r}}}\frac{b}{w+0.441b}

(1-68)  

Symetryczna linia paskowa stosowana w konstrukcjach niektórych przyrządów, jak sprzęgacze, filtry, itp.

 

2.4. Falowód prostokątny

Zgodnie z pokazaną na rys.1.9 strukturą falowód prostokątny jest prowadnicą falową, w której nie występują dwa niezależne przewody, a więc nie może się rozchodzić fala elektromagnetyczna typu TEM. 
Mogą natomiast, przy spełnieniu pewnych warunków rozchodzić mody TE, oznaczane czasami jako mody E, lub TM, oznaczane jako H. Dla każdego z modów konfiguracja pól E i H jest inna. 
W materiale pominieto wyprowadzenie rozkładów pól elektrycznego i magnetycznego dla falowodu prostokątnego, dla kolejnych rodzin modów TEm,n i TMm,n. Skracając i konkludując  należy powiedzieć, że dla każdego modu można określić częstotliwość graniczną, poniżej której dany mod nie może zostać wzbudzony.


Rys.1.9. Struktura i wymiary 
falowodu prostokątnego 

   
Każdy z modów określony jest wskaźnikami „m” i „n”. Wartość częstotliwości granicznej zależy od wartości „m” i „n”, od rozmiarów a i b falowodu, oraz od wartości przenikalności elektrycznej materiału wypełniającego falowód. W miarę wzrostu częstotliwości wzbudzają się kolejne mody TEm,n i TMm,n. Mod o najniższej częstotliwości granicznej nazywany jest podstawowym. Modem podstawowym w falowodzie prostokątnym jest TE10. Dla niego wartość długości fali granicznej (jest to długość fali w wolnej przestrzeni dla częstotliwości granicznej) wynosi:

  

\lambda _{gTE10}=2a;

(1-69)  

W Tabeli 1.2 zestawiono wartości częstotliwości granicznych dla kilku pierwszych modów, dla falowodu skonstruowanego do pracy w pasmie 3 cm, bez wypełnienia dielektrykiem.
Tabela 1.2: Pierwsze mody falowodu na pasmo X, o wymiarach: a=2,286 cm, b=1,016 cm.

MOD TE10 TE20 TE01 TE11 TM11
Częstotliwość graniczna fg[GHz] 6,562  13,123 14,764 16,156 16,156

Rozkład pola elektrycznego i magnetycznego dla modu TE10 pokazano na rys.1.10.

Rys.1.10. Linie sił pola elektrycznego E i magnetycznego H dla modu podstawowego  

 
Pasmo pracy falowodu prostokątnego zawiera się między częstotliwością graniczą modu podstawowego i częstotliwością graniczną kolejnego modu, z pewnymi marginesami. 
Prędkości: fazowa vf i grupowa vg oraz długość fali \lambda _f są, dla tej samej f różne i różne dla różnych modów. Oznaczamy: prędkość v i długość \lambda dla fali płaskiej w wolnej przestrzeni wypełnionej ośrodkiem o \varepsilon _{r}\varepsilon _{0}\mu _{r}\mu _{0}.

  

v=\frac{c}{\sqrt{\mu _{r}\varepsilon _{r}}}; \lambda v=\frac{\lambda _{0}}{\sqrt{\mu _{r}\varepsilon _{r}}};

(1-70)  

Dla falowodów prędkość fazowa i długość fali w falowodzie opisują następujące zależności:

  

v_{f}=\frac{\omega }{\beta }=\frac{v}{\sqrt{1-(f_{gmn}/f)^{2}}};

v_{g}=v\sqrt{1-(f_{gmn}/f)^{2}}=1/\frac{\mathrm{d} \beta }{\mathrm{d} \omega }

(1-71)  

Długość fali \lambda _fw falowodzie jest większa, niż w wolnej przestrzeni i gdy częstotliwość zbliża się do częstotliwości granicznej długość fali rośnie do nieskończoności.

  

\lambda _{f}=\frac{\lambda }{\sqrt{1-(f_{gmn}/f^{2})}};

(1-72)  

Między prędkościami fazową i grupową istnieje związek (1-73):

  

v_{g}v_{f}=v^{2};

(1-73)  

Na rys.1.11 pokazano zależności vf(f/fg) i vg(f/fg). Gdy częstotliwość zbliża się do wartości granicznej, to prędkość fazowa rośnie do nieskończoności, prędkość grupowa maleje do zera i ustaje przepływ energii. Poniżej częstotliwości granicznej dany mod nie może zostać wzbudzony.

Rys.1.11. Prędkości fazowa i grupowa w falowodzie  

 
Straty mocy w ściankach metalowych powodują, że falowody wykazują stosunkowo duże tłumienie mocy sygnału. Aby je zmniejszyć falowody prostokątne wykonywane są  z miedzi, mosiądzu, aluminium, często są srebrzone i złocone.
 

2.5. Falowód cylindryczny

Falowód cylindryczny jest metalową rurą, najczęściej powietrzną, co pokazuje rys.1.12a.


 
 Rys.1.12. Falowód cylindryczny. A) Wymiary falowodu cylindrycznego. 
B) Oś z częstotliwościami kolejnych modów.


W falowodzie cylindrycznym można także wzbudzić nieskończenie wiele modów TEm,n i TMm,n.
Wartość częstotliwości granicznej fgmn związana jest z wartościami:

  • dla modów TMnm z m-tym pierwiastkiem funkcji Bessela Jn(x)=0,
  • dla modów TEnm z m-tym pierwiastkiem pochodnych tych funkcji J’n(x)=0.

Obecność funkcji Bessela wynika z rozwiązania równań Maxwella dla falowodu cylindrycznego.
Modem podstawowym falowodu cylindrycznego jest mod TE11. Długość fali odpowiadającej częstotliwości granicznej dla tego modu równa jest:

  

\lambda _{cTE11}=3.412a;

(1-74)  

Kolejny mod, który się wzbudzi, to TM01, a następnie TE21. tak więc pasmo pracy falowodu cylindrycznego jest niewielkie, co ogranicza zakres zastosowań – rys.1.12B. 
Falowody te stosowane są w konstrukcjach niektórych rezonatorów i filtrów, ze względu na ich duże dobrocie (będzie o tym mowa w jednym z dalszych wykładów).
 

2.6. Linia mikropaskowa

Rozwój technologii układów scalonych, planarnych z samej natury, zmusił konstruktorów do opracowania nowej rodziny prowadnic falowych, które można stosować zarówno w hybrydowych jak i monolitycznych układach scalonych. Najpopularniejszym rozwiązaniem jest linia mikropaskowa, której strukturę pokazano na rys.1.12A. Płaska, o odpowiednio dobranej grubości h warstwa dielektryka pokrywana jest obustronnie metalem. Warstwa metalizacji jest z jednej strony pozostawiona w całości, natomiast z drugiej strony pozostawione są tylko wąskie ścieżki metalizacji o odpowiednio dobranej szerokości w.
Linia mikropaskowa nazywana jest linią quasi-TEM, ponieważ fala EM porusza się w ośrodkach o 2 różnych prędkościach. Linia wykazują niewielka dyspersję.
 
 
Rys.1.13. Prowadnice planarne. A) Linia mikropaskowa rzeczywista struktura. B). Pasek linii mikropaskowej w dielektryku o przenikalności efektywnej.


Impedancja charakterystyczna linii mikropaskowej jest funkcją grubości warstwy h, szerokości paska metalizacji w oraz przenikalności elektrycznej \varepsilon_rdielektryka oddzielającego pasek od metalizacji „ziemi”. Do odpowiednich zależności wprowadza się efektywną przenikalność, której wartość leży między \varepsilon_rpodłoża a \varepsilon_0 powietrza. Z dobrym przybliżeniem można *eff obliczyć z podanej niżej zależności (1-75):

  

\varepsilon _{eff}\cong \frac{\varepsilon_{r}+1}{2}+\frac{\varepsilon_{r}-1}{2}\frac{1}{\sqrt{1+12\frac{h}{w}}};

(1-75)  

Dla w/h ≤1  wartość impedancji charakterystycznej Z0 można obliczyć z zależności (1-76):

  

Z_{0}[\Omega ]\cong \frac{60}{\sqrt{\varepsilon _{eff}}}\ln (\frac{8h}{w}+\frac{w}{4h});

(1-76)  

Gdy w/h > 1 z lepszym przybliżeniem obliczamy wartość impedancji Z0 ze wzoru (1-77):

  

Z_{0}[\Omega ]\cong \frac{377}{\sqrt{\varepsilon _{eff}[\frac{w}{h}+1.393+0.677\ln (\frac{w}{h}+1.444)]}};

(1-77)  

Dla określonych materiałów, które wykorzystano projektując linię mikropaskową długość fali zależy od stosunku w/h, gdyż ze zmianą wartości w/h zmienia się εeff.

  

\lambda _{f}[\mathrm{cm}]]=\frac{30}{\sqrt{\varepsilon _{eff}f[\mathrm{GHz}]}};

(1-78)  

W praktyce obliczenia wymiarów linii mikropaskowych prowadzimy wykorzystując jeden z licznych komputerowych programów obliczeń.
Zakres częstotliwości pracy linii mikropaskowej jest bardzo szeroki, od prądu stałego DC do 30 GHz dla układów hybrydowych, i do 500 GHz dla układów monolitycznych.

2.7. Linie koplanarne

Niedogodnością struktury linii mikropaskowej jest konieczność wykonywania otworów w warstwie dielektryka zwierających pasek z warstwą metalizacji. Jest to szczególnie kłopotliwe w przypadku realizacji monolitycznych układów scalonych. Niedogodności tej nie ma całkowicie planarna struktura linii koplanarnej pokazana na rys.1.14A. W niektórych rozwiązaniach stosowane także dwuprzewodowe linie planarne – rys.1.14B.
Impedancje charakterystyczne Z0 zależy od \varepsilon_rpodło ża i wymiarów linii i mogą być dobierane w szerokich granicach: dla linii mikropaskowej 20...100 \Omega, dla linii koplanarnej 25...150 \Omega, a dla linii dwuprzewodowej koplanarnej 45...220 \Omega.
 
 
Rys.1.14. Prowadnice koplanarne. A). Linia koplanarna 
B). Linia paskowa koplanarna.


Dla potrzeb technologii mikrofalowych układów scalonych opracowano techniki wytwarzania planarnych rezystorów, kondensatorów, cewek indukcyjnych, a także diod i tranzystorów. W scalonych układach hybrydowych MIC elementy te montuje się na powierzchni układu, łączą je z paskami prowadnic falowych. W monolitycznych układach scalonych MMIC wszystkie elementy wykonuje się w procesach technologicznych na podłożu krzemu, albo arsenku galu.
 

2.8. Podłoża linii planarnych

Jak wspomniano wyżej parametry linii planarnej, takie jak impedancja charakterystyczna, długość fali, straty zależą od rodzaju użytego dielektryka. Zestawienie typowych dielektryków i ich najważniejsze parametry zestawiono w Tabeli 1.3.
 

Tabela 1.3. Zestawienie właściwości podłoży prowadnic planarnych

Podłoże \varepsilon_r

\mathrm{tg\delta }

Zastosowania
Al2O3 - ceramika alundowa  9,6 0,0001  MIC
SiO2 - kwarc 3,8 0,00006 MIC
Teflon  2,1 0,00015 MIC
GaAs - arsenek galu 12,5-13 0,002 MMIC
Si - krzem 11,2 0,004 MMIC

Wśród nich najbardziej popularnymi są ceramika alundowa Al2O3, płytki wykonane z kwarcu oraz teflon, niekiedy wymieszany z proszkiem alundowym.
Arsenek galu i krzem jako półprzewodniki samoistne o niewielkim poziomie domieszek nie są najlepszymi dielektrykami i wykonane na nich linie wykazują duże straty. Jednakże stosujemy je w układach monolitycznych, gdyż tylko na takim podłożu można wykonywać aktywne tranzystory.
 

3. Tłumienie prowadnic falowych

Tłumienie mocy sygnału propagowanego prowadnicą falową zależy od trzech najważniejszych czynników:

  • przewodności metalu, z którego wykonano przewody linii,
  • strat materiału dielektrycznego wypełniającego częściowo lub w całości prowadnicę,
  • strat mocy na promieniowanie.

Można więc w ogólnym przypadku zapisać współczynnik tłumienia następująco:

  

\alpha =\alpha _{m}+\alpha _{d}+\alpha _{rad};

(1-79)  

Gdzie \alpha_m reprezentują straty wywołane skończoną przewodnością metalu, \alpha_d to straty wywołane obecnością stratnego dielektryka, a \alpha_{rad} reprezentuje straty wywołane promieniowaniem energii na zewnątrz linii. 
Straty wywołane skończoną przewodnością metalu zależą od przewodności metalu, z którego wykonano falowód, lub elementy przewodzące prowadnicy (np. przewody wewnętrzny i zewnętrzny linii współosiowej). 

Rys.1.15. Wykładnicze zanikanie prądu w przewodniku jako efekt naskórkowości.    


Straty spowodowane niedoskonałością przewodnika są duże i rosną z częstotliwością, ze względu na efekt naskórkowości. Jak nam wiadomo pole elektryczne nie wnika do doskonałego przewodnika, ale do niedoskonałego wnika, na pewną, niewielką głębokość. rezultacie na powierzchni przewodnika płynie prąd, ale tylko w cienkiej warstwie o głębokości wnikania \delta_s[m] i szybko zanika w warstwach głębszych – rys.1.15. Głębokość wnikania może być obliczona ze wzoru (1-80):

  

\delta _{s}[\mathrm{m}]=\frac{1}{\sqrt{\pi f\mu _{r}\mu _{0}\sigma }};

(1-80)  

Na przykład dla miedzi Cu, dla częstotliwości f = 10 GHz, głębokość wnikania jest niewielka i wynosi \delta_s=0,66 \mum.
Do wzorów na wartość stałej tłumienia αm wprowadza się zwykle rezystancję powierzchniową przewodnika RS:

  

R_{s}[\Omega /\mathrm{kw}]=\frac{1}{\sigma _{s}}=\sqrt{\frac{\pi f\mu _{r}\mu _{0}}{\sigma }};

(1-81)  

Rezystancja powierzchniowa Rs[\Omega/kw.] rośnie dla każdego przewodnika z częstotliwością f, choć wolniej dla przewodników dobrze przewodzących (duża przewodność). Dla linii współosiowej straty związane ze skończoną przewodnością metalu opisuje zależność (1-82):

  

\alpha _{m}=\frac{R_{s}}{4\pi Z_{0}}(\frac{1}{a}+\frac{1}{b});

(1-82)  

Wartość tłumienia zależy także od stosunku a/b promieni przewodów, a więc od wartości impedancji charakterystycznej Z0. Dla powietrznej linii współosiowej tłumienie jest najmniejsze dla Z0=77Ω, co pokazuje rys.1.7.
Aby oszacować straty wywołane obecnością dielektryka przyjmiemy, że linia TEM wypełniona jest dielektrykiem o przenikalności ε:

  

\varepsilon =\varepsilon "-j\varepsilon ""=\varepsilon _{r}\varepsilon _{0}(1-j\mathrm{tg\delta });

(1-83)  

Dla linii współosiowej można oszacować straty następująco: 

  

\alpha _{d}=\frac{\pi \mathrm{tg\delta }}{\lambda _{f}}=\frac{\pi f\mathrm{tg\delta }}{v_{f}};

(1-84)  

Jak widać \alpha_d rośnie proporcjonalnie do częstotliwości f we wszystkich liniach TEM i quasi-TEM. Ogólnie można napisać:

  

\alpha _{d}=A \mathrm{tg\delta }

(1-85)  

a stała A proporcjonalności zależy od rozmiarów i rodzaju linii.
Straty na promieniowanie można w przypadku falowodów pominąć, dla linii koncentrycznej także, chyba, że w kablu koncentrycznym przewód zewnętrzny wykonany jest z plecionki metalowej. Jednakże w liniach planarnych nie może być pominięty, choć jest trudny do oszacowania.
Porównanie strat rozmaitych typów prowadnic falowych, zestawione na rys.1.16 prowadzi do przygnębiających wniosków: straty są duże. Podano je w decybelach na kilometr, aby można było porównać je ze stratami w światłowodach kwarcowych, które w najlepszych oknach transmisji promieniowania są rzędu 0,2…0,4 dB/km.
 
Rys.1.16. Tłumienie w funkcji częstotliwości dla różnych typów falowodów prostokątnych wykonanych z aluminium i srebra, oraz dla miedzianych kabli współosiowych.


W wielu obliczeniach projektowych straty prowadnic falowych straty są pomijane. Miejmy na uwadze, że można je pominąć w monolitycznych układach scalonych przy propagacji na odległość 1 mm, w układach hybrydowych przy propagacji na odległość 3 cm, w kablach współosiowych łączących aparaturę pomiarową pracującą w paśmie 1000 MHz na odległościach 1 metra. Ale w sieciach telewizji kablowej już nie można ich pominąć. W zestawieniu ze stratami światłowodu są to tłumienia bardzo duże. 
Tabela 1.4. Porównanie właściwości trzech popularnych typów prowadnic falowych.

Parametr prowadnicy Falowód prostokątny Linia współosiowa Linia mikropaskowa
Mod pracy TE10  TEM  Quasi-TEM
Pasmo pracy Średnie      Duże Duże
Dyspersja     Średnia      Nie występuje Mała 
Tłumienie Małe Średnie     Duże 
Rozmiary     Duże Średnie     Małe
Łatwość wytwarzania Średnia     Średnia Dużą
Możliwość integracji Trudno      Trudno      Łatwo

    
W Tabeli 1.4 zestawiono porównanie właściwości trzech najpopularniejszych typów prowadnic falowych. Każda z nich w pewnych zastosowaniach spisuje się bardzo dobrze, ale w innych zawodzi. Wiedza konstruktora pozwoli dokonać właściwy wybór.

4. Propagacja fal w linii długiej

Ten duży rozdział poświęcony jest opisaniu zjawisk zachodzących w linii długiej w procesie propagacji fali. Zacznijmy od uwagi o tym, jaką linię nazywamy „długą”.
Linię będziemy traktowali jako długą, gdy jej fizyczna długość będzie porównywalna z długością fali propagowanego przez nią sygnału. Tak więc dla fali o długości 100 cm (300 MHz) „długą” będzie kabel koncentryczny o fizycznej długości 10 cm, a dla fali o długości 3 mm (100 GHz) „długą” będzie połączenie między elementami układu scalonego wykonanego na arsenku galu o długości fizycznej 100 m.
Lista pojęć, z którymi zapoznamy się w tym rozdziale i których znaczenie powinniśmy zrozumieć, jest długa. Zaczniemy od wyprowadzenia równań opisujących zjawiska propagacji fali, potem opiszemy rozwiązania tych równań, fale rozchodzące się w układzie: generator-linia długa-obciążenie. Wprowadzimy pojęcia współczynnika odbicia i omówimy warunki dopasowania w rozumieniu impedancyjnym i energetycznym. Omówimy zjawisko fali stojącej i wprowadzimy pojęcie transformacji impedancji. 
Wprowadzimy dużo nowych pojęć i definicji, które będą wykorzystywane w dalszych rozdziałach. Poznanie ich i przyswojenie pozwoli zrozumieć materiał następnych jednostek.
 

4.1. Równania Linii Długiej

Linia dwuprzewodowa
Przeanalizowana zostanie prosta i często spotykana w praktyce struktura prowadnicy falowej, jaką jest linia dwuprzewodowa – rys.10.1a. Przewody tej linii są wykonane z dobrze przewodzącego metalu i „zanurzone” w materiale dielektrycznym. Żaden z tych materiałów nie jest idealnym przewodnikiem, czy też dielektrykiem. Znaczenie użytego przymiotnika „długa” zostanie wyjaśnione dalej.


 
 Rys.10.1. Dwuprzewodowa linia długa. a) Przewody zanurzone materiale dielektrycznym o przenikalności \varepsilon. b) Oznaczenia napięć i prądów.


Celem analizy jest opisanie procesu zmian napięcia i prądu wzdłuż takiego obwodu, gdyż łatwo przewidzieć, że wywołaniu przyrostu napięcia na jednym końcu opisywanej linii nie towarzyszy natychmiastowe pojawienie się identycznego przyrostu na drugim końcu.
Wprowadzamy następujące oznaczenia: 

  • u(t,z) chwilowa wartość napięcia,
  • i(t,z) chwilowa wartość prądu,

przy czym u(t,z) i i(t,z) są funkcjami czasu t i miejsca z.
Przyjmujemy, że propagacja zachodzi w jednym tylko wymiarze z, wzdłuż linii długiej.
Problem: Jak propagują się zmiany napięcia u(t,z) i prądu i(t,z) wzdłuż linii długiej?
     Obwód zastępczy odcinka linii
Rozpatrzymy elementarny czwórnik utworzony przez odcinek linii długiej o długości \Deltaz pokazany na rys.10.1b. Obwód zastępczy takiego czwórnika pokazano na rys.10.2.

Rys.10.2. Obwód zastępczy elementarnego odcinka linii długiej.  

 
W obwodzie tym wprowadzono następujące oznaczenia: 

 

  • R[\Omega/m] - rezystancja na jednostkę długości.
  • L[H/m] - indukcyjność na jednostkę długości.
  • G[S/m] - przewodność na jednostkę długości.
  • C[F/m] - pojemność na jednostkę długości.

     Równania telegrafistów
Zmienne u(t,z) i i(t,z) opisane są wyprowadzonym przez Kelvina równaniami różniczkowymi, zwanymi równaniami telegrafistów. Równania te poznamy w prostej formie, gdyż wyprowadzimy je i rozwiążemy dla prostych i najczęściej spotykanych przypadków, zgodnych z przyjętymi Założeniami 1 i 2.
Założenie 1: u(t) i i(t) są harmonicznymi funkcjami czasu - wielkości te są sinusoidalnymi funkcjami czasu o pulsacji \omega.
Wprowadzamy: 
U(z) - zespolona amplituda napięcia, I(z) - zespolona amplituda prądu.

  

u(z,t)=\mathrm{Re\begin{Bmatrix} \mathrm{U(z)}e^{j\omega t} \end{Bmatrix}};

i(z,t)=\mathrm{Re\begin{Bmatrix} \mathrm{I(z)}e^{j\omega t} \end{Bmatrix}};

(10-1)  

Wprowadzamy: 

  

Z[\mathrm{\Omega /m}]R + j\omega L;\, \, Y[\mathrm{S /m}]G + j\omega C;

(10-2)  
  

\gamma^{2}=YZ;

(10-3)  

Założenie 2: Linia jest jednorodna, Z i Y nie zmieniają się z odległością. 
Założenie 2 oznacza, że linia nie zmienia swoich wymiarów, średnica przewodów a, ich odległość b oraz przenikalność \varepsilon dielektryka otaczającego przewody pozostają stałe i niezależne od z.
Odpowiednie przekształcenia można prześledzić w dołączonym opisie. Końcowy rezultat przekształceń ma postać równań telegrafistów, albo równań linii długiej: 

  

\frac{d^{2}\mathrm{}\mathrm{U} }{ dz^{2}}-\gamma ^{2}\mathrm{U}=0;

(10-4)  
  

\frac{d^{2}\mathrm{}\mathrm{I} }{ dz^{2}}-\gamma ^{2}\mathrm{I}=0;

(10-3)  

Jak widać, zespolone amplitudy prądu U(z) i I(z) jednorodnej linii długiej związane są prostymi równaniami różniczkowymi ze stałą \gamma zwaną stałą propagacji. Stała propagacji \gamma reprezentuje parametry linii długiej, rozmiary przewodów, parametry ośrodka dielektrycznego.
Przypomnijmy jeszcze, że identyczny kształt równań uzyskujemy z równań Maxwella dla pól E i H. Równania te, opisane w wcześniej, zwane są równaniami falowymi.

4.2. Rozwiązanie równań Linii Długiej

Fale: postępująca i odbita
Równania telegrafistów są równaniami różniczkowymi. Ta postać równań różniczkowych ma znaną i prostą postać rozwiązań:
 

  

\mathrm{U}(z)=\mathrm{U}_{+}+\mathrm{U}_{-}= \mathrm{U}_{1}e^{-\gamma z}+\mathrm{U}_{2}e^{\gamma z}

(10-6a)  
  

\mathrm{I}(z)=\mathrm{I}_{+}+\mathrm{I}_{-}= \mathrm{I}_{1}e^{-\gamma z}+\mathrm{I}_{2}e^{\gamma z}

(10-6b)  

Rozwiązanie jest dwuczłonowe, składniki z indeksem „1” reprezentują falę rozchodzącą się wzdłuż osi z, składniki z indeksem „2” reprezentują falę rozchodzącą się w przeciwną stronę, niż kierunek osi z. 
Pamiętamy prostą i oczywistą interpretację rozwiązań:

  •  U1, I1 - stałe całkowania – zespolone amplitudy napięcia i prądu fali rozchodzącej się w kierunku z, jest to fala postępująca.
  • U2, I2 - stałe całkowania - zespolone amplitudy napięcia i prądu fali rozchodzącej się w kierunku przeciwnym do z, nazywamy ją falą odbitą, albo wtórną.

Pamiętamy: Dla każdego typu prowadnicy falowej, w której propagowany jest jeden mod fali, można przyjąć obwód zastępczy w postaci  linii dwuprzewodowej. W każdym takim przypadku rozwiązanie równania linii długiej mają postać (10-6) i (10-7) i ich interpretacja jest identyczna.
     Stała propagacji
Gdy mówimy o propagacji fali, to powinniśmy wyznaczyć tłumienie fali, długość fali i prędkości rozchodzenia. Wprowadzona zależnością (10-3) i występująca w rozwiązaniach (10-6) i (10-7) stała propagacji \gamma jest bardzo ważnym parametrem zjawiska propagacji fali. Stała propagacji jest wielkością zespoloną i można zapisać ją w następującej postaci: 

  

\gamma = \alpha +j\beta ;

(10-7)  

Zależność (10-6a) może być zapisana następująco:

  

\mathrm{U}(z)=\mathrm{U}_{1}e^{-\alpha z}e^{-j\beta z}+\mathrm{U}_{2}e^{\alpha z}e^{j\beta z};

(10-8)  

Interpretacja fizyczna obu składników ​​\alpha +j\beta jest oczywista:

  • Część rzeczywista \alpha stałej propagacji \gamma nazywana jest stałą tłumienia. Stała tłumienia \alpha(Np/m) decyduje o szybkości strat mocy fali biegnącej wzdłuż linii.
  • Część urojona \beta stałej propagacji \gamma nazywana jest stałą fazowa. Stała fazowa \beta(rad/m) decyduje o szybkości zmian fazy fali biegnącej wzdłuż linii, a tym samym o długości fali \lambda.

Powróćmy do zależności (10-3), aby znaleźć, jak \alpha i \beta zależą od parametrów R,G,L i C obwodu zastępczego z rys.10.2. Wracamy do podstawowych zależności opisujących stałą propagacji  \gamma w zależności od impedancji Z i admitancji Y.

  

\gamma =\sqrt{YZ}=\sqrt{(R+j\omega L)(G+j\omega C)};

(10-9)  

Zwykle spełnione są następujące warunki:  R/\omega L , gdyż w praktycznych rozwiązaniach konstrukcji linii dwuprzewodowych przewody wykonane są z dobrze przewodzącego metalu i otoczone małostratnym dielektrykiem. Wtedy:

  

\alpha \cong \frac{G}{2}\sqrt{\frac{L}{C}}+\frac{R}{2}\sqrt{\frac{C}{L}};

(10-10)  
  

\beta \cong \omega \sqrt{LC};

(10-11)  

W zależności (10-10) drugi ze składników zwykle dominuje nad pierwszym.
     Prędkości fazowa i grupowa
Gdy mówimy o prędkości propagacji fali musimy wyróżnić prędkość fazową i prędkość grupową.
Prędkość fazowa vf propagowanej fali jest prędkością z jaką przesuwa się płaszczyzna stałej fazy. Prędkość vf związana jest z wartością stałej fazowej \beta:
 

  

v_{f}=\frac{\mathrm{d} z}{\mathrm{d} t}=\frac{\omega }{\beta };\, \, \beta =\frac{\omega }{v_{f}}=\frac{2\pi }{\lambda _{f}};

(10-12)  

Prędkość grupowa vg propagowanej fali jest prędkością przepływu energii.

  

v_{g}=\frac{\mathrm{d} \omega }{\mathrm{d} \beta }

(10-13)  

Przypomnienie:
W prowadnicach falowych typu TEM prędkości fazowa i grupowa są sobie równe. W falowodach prostokątnych i cylindrycznych, w których propagowane są mody TE albo TM, prędkości fazowa i grupowa różnią się.
     Impedancja charakterystyczna
Zespolone amplitudy napięcia U(z) i prądu I(z) opisane są zależnościami (10-6) i (10-7). Określimy teraz związki między nimi. Stosunki zespolonych amplitud napięcia i prądu dla obu propagowanych fal są sobie równe z dokładnością do znaku i nazwane impedancją charakterystyczną Z0.

  

{\mathrm{Z_{0}} }=\frac{\mathrm{U_{1}}}{\mathrm{I_1}}=-\frac{\mathrm{U_{2}}}{\mathrm{I_2}};

(10-14)  

Wartość impedancji charakterystycznej jest bardzo ważnym parametrem prowadnicy falowej. Impedancja charakterystyczna Z0 jest funkcją rozmiarów prowadnicy i parametrów ośrodka.

  

{\mathrm{Z_{0}} }=\sqrt{\frac{R+j\omega L}{G+j\omega C}}\cong \sqrt{\frac{L}{C}};

(10-15)  

Dla prowadnicy bezstratnej Z0 jest rzeczywiste. Dla prowadnicy z małymi stratami przyjmuje się także, że z dobrym przybliżeniem Z0 jest rzeczywiste.

4.3. Fale w jednorodnej prowadnicy falowej

Napięcie i prąd wzdłuż linii
Powracamy do układu generator – prowadnica – obciążenie: 
 


Rys.10.3. Układ: generator, linia długa, obciążenie

Napiszemy najpierw rozwiązania równań linii długiej z nowymi oznaczeniami. Oznaczymy przez:

  • UP, IP - zespolone amplitudy napięcia i prądu fali pierwotnej, padającej.
  • UW, IW - zespolone amplitudy napięcia i prądu fali odbitej, wtórnej.

Odległość l liczona jest teraz od końca linii w stronę generatora, podczas gdy z liczona była od generatora w stronę obciążenia. Równania propagacji zapiszą się teraz następująco:

  

\mathrm{U(l)}=\mathrm{U_{P}}e^{\gamma l}+\mathrm{U_{W}}e^{-\gamma l};

\mathrm{I(l)}=\mathrm{I_{P}}e^{\gamma l}+\mathrm{I_{W}}e^{-\gamma l};

(10-16)  

 Dla bezstratnej prowadnicy falowej, gdy  \gamma =j\beta ;:

  

 

\mathrm{U(l)}=\mathrm{U_{P}}e^{j\beta l}+\mathrm{U_{W}}e^{-j\beta l};

\mathrm{I(l)}=\mathrm{I_{P}}e^{j\beta l}+\mathrm{I_{W}}e^{-j\beta l};

(10-17)  

Znamy już definicję impedancji charakterystycznej: 

  

\mathrm{Z_{0}}=\mathrm{\frac{U_{P}}{I_{P}}}=-\mathrm{\frac{U_{W}}{I_{W}}};

(10-18)  

Aby znaleźć wartości amplitud UP, IP, UW i IW musimy znaleźć kolejne związki między nimi.

     Wpływ obciążenia na falę odbitą
Obciążenie reprezentowane przez impedancję:

  

\mathrm{Z_{L}}=\mathrm{\frac{U_{L}}{I_{L}}}=\mathrm{\frac{1}{Y_{L}}};

(10-19)  

Należy zauważyć, że UL=U(l=0) oraz IL=I(l=0), czyli: 

  

\mathrm{Z_{L}}=\mathrm{\frac{U_{L}}{I_{L}}}=\mathrm{\frac{U_{P}+U_{W}}{I_{P}+I_{W}}};

(10-20)  

w rezultacie: 

  

\mathrm{\frac{U_{W}}{U_{P}}}=\mathrm{\frac{Z_{L}-Z_{0}}{Z_{L}+Z_{0}}};

(10-21)  

Wniosek: Wartość amplitudy UW napięcia fali odbitej zależy nie tylko od ZL, ale także od wartości impedancji charakterystycznej Z0. Gdy ZL=Z0, w prowadnicy nie pojawi się fala odbita; obciążenie jest dopasowane do impedancji charakterystycznej prowadnicy falowej, obciążenie jest bezodbiciowe.
W równaniach (10-16) i (10-17) na U(l) i I(l) można zastąpić UP, UW, IP i IW przez ZL, UL i IL.
W ogólnym przypadku linii ze stratami otrzymuje się następujące zależności:

  

\mathrm{U}(l)=\mathrm{U_{L}}(\cosh \gamma l+\frac{\mathrm{Z_{0}}}{Z_{L}}\sinh \gamma l);

\mathrm{I}(l)=\mathrm{I_{L}}(\cosh \gamma l+\frac{\mathrm{Z_{L}}}{Z_{0}}\sinh \gamma l);

(10-22)  

Dla linii bezstratnej funkcje hiperboliczne znikają i równania przyjmują prostszą postać:

  

\mathrm{U}(l)=\mathrm{U_{L}}(\cos \beta l+j\frac{\mathrm{Z_{0}}}{Z_{L}}\sin \beta l);

\mathrm{I}(l)=\mathrm{I_{L}}(\cos \beta l+j\frac{\mathrm{Z_{L}}}{Z_{0}}\sin \beta l);

(10-23)  

Można też zapisać je w innej formie, w pewnych przypadkach wygodniejszej w użyciu:

  

\mathrm{U}(l)=\mathrm{U_{L}}(\cos \beta l+j\mathrm{I_{L}}{\mathrm{Z_{0}}}\sin \beta l);

\mathrm{I}(l)=\mathrm{I_{L}}(\cos \beta l+j\frac{\mathrm{U_{L}}}{\mathrm{Z_{0}}}\sin \beta l);

(10-23a)  

 Współczynnik odbicia - definicja
Zdefiniowany zostanie teraz bardzo ważny parametr określający związek między falą odbitą i padającą. Współczynnik odbicia \Gamma jest miarą stosunku zespolonych amplitud fali odbitej do padającej. Definiujemy go następująco:

  

\Gamma _{L}=\mathrm{\frac{U_{-}}{U_{+}}}=\mathrm{\frac{U_{W}}{U_{P}}}e^{-2\gamma l};

(10-24)  

Na końcu linii, dla l=0, współczynnik odbicia przyjmuje wartość \GammaL:

  

\Gamma (l=0)=\Gamma _{L}=\mathrm{\frac{U_{W}}{U_{P}}}=\mathrm{\frac{Z_{L}-Z_{0}}{Z_{L}+Z_{0}}}=\mathrm{\frac{Y_{0}-Y_{L}}{Y_{0}+Y_{L}}}

(10-25)  

Współczynnik odbicia \GammaL - podobnie jak ZL lub YL - jest parametrem charakteryzującym jednowrotnik/obciążenie umieszczone na końcu linii, inaczej mówiąc, jest on zespoloną miarą niedopasowania obciążenia do impedancji charakterystycznej Z0
Współczynnik odbicia \Gamma(l) zależy od wartości \GammaL na końcu linii oraz od odległości l od końca linii. Zależność ta ma następującą postać:

  

\Gamma _{L}=\Gamma _{L}e^{-2\gamma l}\mid _{\alpha =0}=\Gamma _{L}e^{-2\beta l};

(10-26)  

Napisana wyżej zależność (10-26) nazywana jest równaniem transformacji współczynnika odbicia.

Rys.10.4. Transformacja współczynnika odbicia wzdłuż bezstratnej linii.    
Ilustracja procesu transformacji współczynnika \Gamma pokazana jest na rys.10.4. Wskaz \Gamma wiruje zgodnie ze wskazówkami zegara. Dla linii ze stratami długość wskazu \left |\Gamma \right |maleje wykładniczo z odległością, dla linii bezstratnej  \left |\Gamma \right |=const.
     Współczynnik odbicia - różne przypadki
Obciążenie umieszczone na końcu linii reprezentowane jest zwykle przez impedancję ZL=RL+jXL. Wtedy zależność (10-26) przyjmie postać (10-27):

  

\left |\Gamma _{L} \right |^{2}=\mathrm{\frac{(R_{L}-Z_{0})^{2}+X_{L}^{2}}{(R_{L}+Z_{0})^{2}+X_{L}^{2}}}

(10-27)  

Logarytmiczną miarą współczynnika odbicia są straty odbicia Lr (ang. return loss), definiowane następująco:

  

\mathrm{L_{r}}=-20\log\left | \Gamma _{L} \right |

(10-27)  

Wykorzystamy wyprowadzoną zależność (10-26) do analizy kilku charakterystycznych przypadków.
Przypadek 1: Mówimy, że umieszczony na końcu prowadnicy jednowrotnik, nazywany też obciążeniem, jest dopasowany do impedancji charakterystycznej tej prowadnicy jeżeli  \left |\Gamma \right |=0. Zgodnie z (10-27) stan dopasowania powstanie, gdy Z= Z0.
Przypadek 2: Stan pełnego odbicia mocy powstaje wtedy, gdy \left |\Gamma \right |=1  i amplitudy obu fal: padającej i odbitej są sobie równe. Pełne odbicie mocy ma miejsce, gdy obciążenie jest czystą reaktancją ZL = jXL. Wartość reaktancji XL ma wpływ na argument współczynnika odbicia, jego moduł równy jest 1.

  

\Gamma _{L} =\mathrm{\frac{(X_{L}-Z_{0})}{jX_{L}+Z_{0}}}=e^{j(\pi -2\textrm{arctg}\frac{X_{L}}{Z_{0}})}

(10-28)  

Przypadek 3: Najczęściej impedancja obciążenia obok części urojonej ma część rzeczywistą, przy czym  RL > 0. Wtedy część mocy ( \left | \mathrm{I_{L}} \right |\mathrm{R_{L}/2}) fali padającej zostaje pochłonięta przez obciążenie i amplituda fali odbitej jest zawsze mniejsza od amplitudy fali padającej, a  \left |\Gamma \right | < 1.
Przypadek 4: Gdy amplituda fali odbitej jest większa od amplitudy fali padającej, mamy do czynienia ze wzmocnieniem mocy, z obciążeniem aktywnym. W modelu impedancyjnym obciążenie takie reprezentowane jest przez impedancję z ujemną rezystancją. Gdy  \left |\Gamma \right |>>1, wtedy  RL < 0.

4.4. Fala stojąca

 Napięcie i prąd wzdłuż linii – wykres wskazowy

W tym punkcie wyprowadzimy odpowiednie formuły opisujące rozkład napięcia wzdłuż linii długiej. Powracamy do układu z rys.10.3 z generatorem, prowadnicą i obciążeniem. Wykorzystamy zależność (10-24) opisującą współczynniki odbicia \Gamma(l) aby określić wartości amplitud napięcia i prądu na linii. Prowadzi to do zależności (10-29):

  

\mathrm{U}(l)=\mathrm{U_{+}}[1+\Gamma (l)]=\mathrm{U_{P}}e^{j\beta l}(1+\Gamma _{L}e^{-j2\beta l})

\mathrm{I}(l)=\mathrm{I_{+}}[1-\Gamma (l)]=\frac{\mathrm{U_{P}}}{\mathrm{Z_{0}}}e^{j\beta l}(1-\Gamma _{L}e^{-j2\beta l})

(10-29)  

Zauważmy, że na końcu linii napięcie UL jest proporcjonalne do (1+\GammaL) a prąd IL jest proporcjonalny do (1-\GammaL). Wskazy napięcia UL i prądu IL pokazane są na rys. 10.5A.


 
 
Rys.10.5. Ilustracja zmian napięcia i prądu wzdłuż linii. 
A) Napięcie i prąd na końcu linii. B) napięcie i prąd w odległości l od końca.


W odległości l od końca wartości napięcia i prądu zmieniają się, ponieważ wskaz +\GammaL zmienił położenie – rys. 10.5B.
Kąt fazowy \phi _{L} między UL i prądu IL zależy od impedancji obciążenia: 

  

\phi _{L}=\arctan \mathrm{\frac{X_{L}}{R_{L}}}

(10-30)  

Gdy impedancja obciążenia jest rzeczywista prąd i napięcie są w fazie.
     Napięcie i prąd wzdłuż linii – minima, maksima
Moduł napięcia \left | U(l) \right | można wyznaczyć wychodząc z zależności (10-29). Dla linii bezstratnej otrzymujemy zależność (10-33), w której \Psi_L jest argumentem współczynnika odbicia \Gamma:

  

\left | U(l) \right |=\left | U_P \right |\sqrt{1+\left | \Gamma \right |^{2}+2\left | \Gamma \right |\cos (2\beta l-\Psi_L )}

(10-33)  

Przykład przebiegu \left | U(l) \right | pokazano na rys.10.6.

Rys.10.6. Moduł U(l) napięcia wzdłuż linii długiej.

Ponieważ przyjęto założenie bezstratności linii, to wszystkie maksymalne i minimalne wartości napięcia są sobie równe. 
Wnioski: Napięcie \left | U(l) \right | określone wzdłuż linii długiej jest okresową funkcją odległości o okresie równym połowie długości fali \lambda/2, co oznacza, że:

  • odległość między kolejnymi maksimami, lub minimami równa jest \lambda/2,
  • odległość między maksimum a minimum równa jest \lambda/4.

     Czysta fala stojąca
W przypadku, gdy \left | \Gamma \right |=1 amplitudy fali padającej i odbitej są sobie równe i mamy do czynienia z czystą falą stojącą. Moduł napięcia wzdłuż linii zapisze się dla przypadku \Psi_k=0 następująco: (10-34)

  

\left | \mathrm{U(l)} \right |=2\left | \mathrm{U_P} \right |\left | \cos \beta l \right |;

(10-34)  

Moduł prądu wzdłuż linii opisuje zależność (10-35): 

  

\left | \mathrm{I(l)} \right |=2\frac{\mathrm{U_{P}}}{\mathrm{Z_{0}}}\left | \sin \beta l \right |

(10-35)  


Na rys. 10.7 pokazano przebiegi modułów U(l) i I(l) dla czystej fali stojącej.
 


Rys.10 .7. Napięcie i prąd wzdłuż linii dla czystej fali stojącej.

Jak widać wartości napięć i prądów okresowo osiągają wartości maksymalne i spadają do zera, przy czym maksymalnej wartości napięcia towarzyszy zero wartości prądu i na odwrót. Kolejne zera oddalone są o pół fali.
     Współczynnik fali stojącej
Ważnym parametrem opisującym rozkład napięcia na linii i tym samym stan dopasowania obciążenia do impedancji charakterystycznej Z0 jest współczynnik fali stojącej. Zgodnie z definicją współczynnik fali stojącej \varrho jest stosunkiem maksymalnej i minimalnej wartości modułu napięcia na linii. 

  

\varrho =\frac{\left | \mathrm{U}(l) \right |_{max}}{\left | \mathrm{U}(l) \right |_{min}}=\frac{1+\left | \Gamma \right |}{1-\left | \Gamma \right |}\geq 1

(10-36)  

Współczynnik fali stojącej, często oznaczany jako WFS, jest liczbą rzeczywistą, co oznacza, iż daje nam tylko jedną informację o jednowrotniku/obciążeniu. Pamiętajmy, że współczynnik odbicia \Gamma daje – jako liczba zespolona – dwie informacje o obciążeniu. Między tymi wielkościami istnieje prosty i oczywisty związek:

  

\left | \Gamma \right |=\frac{\varrho -1}{\varrho +1};

(10-37)  

Graficzna ilustracja powyższej zależności pokazana jest na rys. 10.8.

Omówimy dwa charakterystyczne przypadki obciążenia linii:
Przypadek 1: Na końcu prowadnicy umieszczono rezystancję RL > Z0 . Współczynnik fali stojącej dla takiego obciążenia obliczamy prosto jako  \varrho =RL/Z0.
Przypadek 2: Na końcu prowadnicy umieszczono rezystancję  RL < Z0. Współczynnik fali stojącej dla takiego obciążenia obliczamy prosto jako    \varrho =RL/Z0..

4.5. Przepływ mocy

Postawienie problemu 
Kolejny raz wracamy do rys.10.3 i rozważymy zjawiska zachodzące w prostym układzie generator-linia długa-obciążenie. Układ ten powtórnie pokazano na rys.10.9, jednakże z użyciem nieco innych oznaczeń elementów.
Celem rozważań jest określenie mocy występujących w tym prostym układzie. Wyznaczymy:

  • moce fal pierwotnej i odbitej, 
  • moc wydzieloną w obciążeniu, 
  • maksymalną moc, którą może dostarczyć generator,
  • warunek, przy którym to może nastąpić.

Rys.10.9. Układ: generator, prowadnica i obciążenie.

    
Rozważania będą prowadzone przy następujących oznaczeniach i założeniach:

  • generator reprezentowany parametrami źródła EG i ZG,
  • prowadnica falowa jest jednorodna i bezstratna, opisana przez: Z0 i \beta l=2\pi l/\lambda
  • obciążenie/jednowrotnik  charakteryzowany jest przez ZL, YL bądź \Gamma_L

Ponadto przyjmiemy, że w prowadnicy rozchodzą się fale o amplitudach Uw i Up.
     Generator mikrofalowy
Generator umieszczony w obwodzie z rys.10.9 jest generatorem idealnym. Pobudza prowadnicę falową monoczęstotliwościowym sygnałem i jest „odseparowany od obciążenia”, co oznacza, że zmiany impedancji obciążenia nie zmieniają warunków generacji, m.in. jego częstotliwości.

Rys.10.10. Generator mikrofalowy. A) Szeregowy obwód zastępczy generatora, B) Obwód równoległy, C) Graf przepływu sygnału.    

Na rys.10.10 pokazano trzy różne sposoby opisywania generatora:

  • Szeregowy obwód zastępczy składa się z: idealnego źródła napięciowego EG o impedancji wewnętrznej ZG.
  • Równoległy obwód zastępczy to źródło prądowe o wydajności IG i admitancji wewnętrznej YG, przy czym wartości IG i YG związane są z EG i ZG następująco:
  

\mathrm{I_{G}}=\mathrm{\frac{E_{G}}{Z_{G}}};\, \, \mathrm{Y_{G}}=\mathrm{\frac{1}{Z_{G}}};

(10-38)  
  •  Zachowanie się generatora mikrofalowego można także opisać zależnością (10-39), wiążącą amplitudy aG i bG fal wypływającej i powracającej do generatora:
  

\mathrm{a_{G}}=\mathrm{E+b_{G}\Gamma _{G}};

(10-39)  
  

\mathrm{E=\frac{E_{G} \sqrt{Z_{0}}}{Z_{G}+Z_{0}}}=\mathrm{\frac{I_{G} \sqrt{Y_{0}}}{Y_{G}+Y_{0}}}

(10-40)  
  

\mathrm{\Gamma _{G}=\frac{Z_{G} -Z_{0}}{Z_{G}+Z_{0}}}=\mathrm{\frac{Y_{0} -Y_{G}}{Y_{G}+Y_{0}}};

(10-41)  

Zależności (10-40) i (10-41) pokazują, że wszystkie trzy formy prezentacji generatora są sobie równoważne, a elementy tych form są ze sobą związane prostymi zależnościami. Miejmy na uwadze, że niezależnie od sposobu prezentacji właściwości generatora dochodzimy do takich samych wyników końcowych. 
     Moce fal
Jako punkt wyjścia przyjmiemy warunki dopasowanego obciążenia:

  

\mathrm{Z_{L}=Z_{0}};\, \, \mathrm{\Gamma _{L}=0};\, \, \mathrm{U _{w}=0};

(10-42)  

W obwodzie płynie fala pierwotna do obciążenia i nie ma fali odbitej. Napięcie na zaciskach obciążenia równe jest: 

  

\mathrm{U_{L}=U_{p}(1+\Gamma _{L})= U_{p}};

(10-43)  

Moc PL wydzielona w obciążeniu: 

  

\mathrm{P_{L}=\frac{\left | U_{p} \right |^{2}}{2Z_{0}}= \frac{\left | U_{p} \right |^{2}Y_{0}}{2}};

(10-44)  

Moc wydzielona w obciążeniu jest mocą niesioną przez falę pierwotną, nie ma fali odbitej, czyli moc fali pierwotnej opisana jest następującym wzorem P+:

  

\mathrm{P_{+}=\frac{\left | U_{p} \right |^{2}}{2Z_{0}}= \frac{\left | a_{L} \right |^{2}}{2}};

(10-45a)  

Przez analogię moc fali odbitej P-

  

\mathrm{P_{-}=\frac{\left | U_{w} \right |^{2}}{2Z_{0}}= \frac{\left | b_{L} \right |^{2}}{2}};

(10-45b)  

Zauważamy, że w zależnościach (10-44) i (10-45) wprowadziliśmy nowe oznaczenia a i b zwane zespolonymi i znormalizowanymi amplitudami napięć fal i definiowane następująco:

  

\mathrm{a=\frac{U_{p}}{\sqrt{Z_{0}}}};\, \, \mathrm{b=\frac{U_{w}}{\sqrt{Z_{0}}}};

(10-46)  

W wielu przypadkach, które zostaną opisane później, wielkości te zostaną wykorzystane do opisu dwuwrotników mikrofalowych.
     Moc wydzielona w obciążeniu
Do rozważań w tym punkcie przyjmiemy warunki (10-47). 

  

\mathrm{Z_{L}\neq Z_{0}} ; \Gamma _{\mathrm{G}}=0;

(10-47)  

Obciążenie jest niedopasowane i część mocy P+ niesionej przez falę pierwotną/padającą  zostaje odbita i jako moc P- wędruje w stronę generatora. Oznaczając przez PL moc wydzieloną w jednowrotniku można napisać oczywisty bilans mocy:

  

\mathrm{P_{L}=P_{+}-P_{-}};

(10-48)  

Można teraz połączyć ze sobą moce: padającą i wydzieloną w obciążeniu ze współczynnikiem odbicia. Otrzymujemy:

  

\mathrm{P_{L}=P_{+}(1-\frac{\left | U_{w} \right |^{2}}{\left | U_{p} \right |^{2}})=P_{+}(1-\left | \Gamma _{L} \right |^{2})};

(10-49)  

Jak widać argument współczynnika odbicia nie ma wpływu na bilans mocy. Do powyższej zależności można dopisać dwie kolejne:

  

\mathrm{\frac{P_{-}}{P_{+}}=\left | \Gamma _L \right |^{2}=\left | \Gamma \right |^{2}};

(10-50)  
  

\mathrm{\frac{P_{-}}{P_{+}}=\frac{4}{\varrho +1/\varrho +2}};

(10-51)  

 Stosunek mocy P- odbitej do padającej P+ jest zależny tylko od modułu współczynnika odbicia, co oznacza, że znajomość współczynnika fali stojącej WFS pozwala określić stosunki wszystkich trzech mocy.
     Moce w warunkach obustronnego niedopasowania
W tym punkcie przyjmiemy założenie (10-52), że generator i obciążenie są niedopasowane do impedancji charakterystycznej Z0 prowadnicy falowej.

  

\mathrm{\Gamma _{L}\neq 0};\, \, \mathrm{\Gamma _{G}\neq 0};

(10-52)  

 Jest to przypadek ogólny i często spotykany. Amplitudy a i b fal pierwotnej i odbitej oznaczone na rys.10.9 powiązane są ze sobą następującymi zależnościami:

  

\mathrm{a_{G }=E+\Gamma _{G}b_{G}}; \, \, \mathrm{b_{L}=\Gamma _{L}a_L}

(10-53)  
  

\mathrm{a_{L }=a _{G}}e^{j\beta l}; \, \, \mathrm{b_{G }=b _{L}}e^{-j\beta l};

(10-54)  

Wykorzystując te zależności oraz zależność (10-39) opisującą amplitudę fali na wyjściu generatora można skonstruować graf przepływu sygnału pokazany na rys.10.11. 

Rys.10.11. Graf przepływu sygnału w obwodzie z rys.10.9.  

 
Na podstawie grafu znajdujemy amplitudę aL fali padającej na obciążenie:

  

\mathrm{a_L}=\frac{\mathrm{E}e^{-j\beta l}}{1-\mathrm{\Gamma _{G}\Gamma _{L}}e^{-j2\beta l}};

(10-54)  

Jak widać z otrzymanej zależności odległość między generatorem a obciążeniem wpływa w istotny sposób na wartość amplitudy \mathrm{\left | a_L \right |} fali, jaka ustali się na skutek odbić od obciążenia i generatora. Moc P+ niesiona przez falę zmieni się w jeszcze szerszych granicach, gdyż z kwadratem amplitudy \mathrm{\left | a_L \right |}, co widać z zależności (10-55), w której \mathrm{P_{G0}=\left |E \right |^{2}/2}  jest mocą niesioną przez falę pierwotną w warunkach dopasowania:

  

\mathrm{P_{+}=\frac{P_{G0}}{\left | 1-\Gamma _{G} \Gamma _{L}e^{j2\beta l}\right |^{2}}};

(10-55)  

W zależności powyższej P+ jest > lub < od PG0; stosunek maksymalnej do minimalnej mocy może zmieniać się w szerokich granicach:

  

\mathrm{\frac{P_{+max}}{P_{+min}}}=(\frac{1+\left | \mathrm{\Gamma _{G}\Gamma _{L}} \right |}{1-\left | \mathrm{\Gamma _{G}\Gamma _{L}} \right |});

(10-56)  

 Im silniejsze jest obustronne niedopasowanie, im większą wartość ma moduł   w tym szerszych granicach zmienia się moc niesiona przez falę padającą na obciążenie i tym samym moc PL wydzielona w obciążeniu. Moc tą obliczamy opierając się na (10-57):

  

\mathrm{P_L=\frac{1}{2}\left ( \left | a_L \right |^{2}- \left | b_L \right |^{2}\right )=\frac{1}{2}\left | a_L \right |^{2}(1-\left | \Gamma _L \right |^{2})};

(10-57)  

Po prostych przekształceniach otrzymujemy:

  

\mathrm{P_L=P_{G0}\frac{1-\left | \Gamma _{\mathrm{L}} \right |^{2}}{\left | 1- \Gamma _{\mathrm{G}}\Gamma _{\mathrm{L}}e^{-j2\beta l}\right |}};

(10-58)  

 Zależność powyższa jest podstawą do wprowadzenia pojęcia dopasowania energetycznego.
     Dopasowanie energetyczne i moc dysponowana generatora
Rozważymy następujący problem: bezstratna linia długa zasilana jest przez generator niedopasowany, dla którego \mathrm{\left | \Gamma _G \right |>0}. Jak dobrać warunki obciążenia generatora, to znaczy jak dobrać \mathrm{\left | \Gamma _L \right |}i długość linii, aby w obciążeniu wydzieliła się maksymalna moc?
Rozwiązanie problemu można znaleźć analizując zależność (10-58). Moc w PL wydzielona w obciążeniu jest maksymalna, gdy spełniony jest następujący warunek:

  

\mathrm{ \Gamma _L }e^{-j2}={\mathrm{\Gamma _G}^{*}};

(10-33)  

(10-59)
Warunek ten nazywamy warunkiem dopasowania energetycznego. Oznacza on, że współczynnik odbicia „widziany” przez generator w jego wrotach wyjściowych powinien być równy sprzężonej wartości jego własnego współczynnika odbicia.
Problem ten był także rozważany w obwodach o stałych skupionych. Zapisany formułą impedancyjną ma znaną postać:

  

\mathrm{Z_{L}^{"}=Z_{G}^{*}};

(10-60)  

gdzie ZL jest impedancją obciążenia widzianą przez generator.
W warunkach dopasowania energetycznego moc PGA wydzielona w umieszczonym na końcu prowadnicy jednowrotniku jest maksymalna i nazywana mocą dysponowaną generatora (ang. available power). 

  

\mathrm{P_{GA}=\frac{ P_{G0} }{1-\left | \Gamma _G \right |^{2}}};

(10-61)  

Zależność (10-61) zapisana w formule impedancyjnej/admitancyjnej ma postać:

  

\mathrm{P_{GA}=\frac{ \left |E_{G} \right |^{2} }{4Re(Z_{G})}=\frac{ \left |I_{G} \right |^{2} }{4Re(Y_{G})}};

(10-62)  

Znamy teraz dopasowanie dwojakiego rodzaju: 

  • dopasowanie impedancji obciążenia ZL do impedancji charakterystycznej Z0 prowadnicy falowej, co jest równoznaczne warunkowi bezodbiciowości, 
  • dopasowanie impedancji obciążenia ZL do impedancji wewnętrznej generatora ZG, co jest warunkiem dopasowania energetycznego.

Powinniśmy umieć odróżniać opisane warunki dopasowania i analizując kolejne problemy właściwie je interpretować. Warunkom dopasowania poświęcone będą oddzielne wykłady, gdyż umiejętność dopasowania odgrywa w technice mikrofalowej dużą rolę. 

4.6. Transformacja impedancji

Postawienie problemu 
Odpowiemy teraz na pytanie, jak zmieni się impedancja ZL przez dodanie odcinka prowadnicy falowej o odpowiedniej długości l i przez dobór jej impedancji charakterystycznej Z0 – rys.10.12. Rozwiązanie tego problemu oznacza, że impedancję ZL i odcinek prowadnicy l,Z0 zastąpimy teraz impedancją Z(l) o takiej wartości, że rozkłady prądów i napięć na lewo od płaszczyzny l nie ulegną zmianie.


 Rys.10.12. Odcinek prowadnicy falowej zakończony impedancją ZL

Aby rozwiązać postawiony problem należy wyznaczyć wartości napięcia U(l) i prądu i I(l) w płaszczyźnie odległej o l od końca. Jeśli to się uda zrobić, to odcinek prowadnicy o długości l i impedancji charakterystycznej Z0 oraz impedancję ZL można zastąpić impedancją Z(l), równą:

  

\mathrm{Z}(l)\equiv \mathrm{Z}(l, \mathrm{Z_0}, \mathrm{Z_L})=\frac{\mathrm{U}(l)}{\mathrm{I}(l)};

(10-63)  

 Równanie transformacji impedancji – linia bezstratna
Aby rozwiązać postawiony problem wykorzystamy znane związki (10-29) między współczynnikiem odbicia Γ(l) a napięciem i prądem, otrzymamy wtedy: 

  

\mathrm{Z}(l))=\mathrm{Z_0}\frac{1+\Gamma (l)}{1-\Gamma (l)};

(10-64)  

Przyjmiemy dalej, że znamy wartość współczynnika odbicia na końcu linii ΓL(ZL), a linia jest bezstratna, to znaczy stała propagacji zapisze się jako \gamma =j\beta . Wtedy: 

  

\mathrm{Z}(l)=\mathrm{Z_0}\frac{1+\Gamma_L e^{-j2\beta l}}{1-\Gamma_Le^{-j2\beta l}};

(10-65)  

Teraz należy wykorzystać znaną z teorii liczb zespolonych tożsamość ejx=cosx+jsinx. Po przekształceniach otrzymujemy równanie transformacji impedancji z tangensami:

  

\mathrm{Z}(l)=\mathrm{Z_0}\frac{\mathrm{Z_L}+j\mathrm{Z_0} \mathrm{tg}\beta l}{\mathrm{Z_0}+j\mathrm{Z_L} \mathrm{tg}\beta l};

(10-66)  

Odwrotność wyrażenia (10-66) daje admitancję Y(l) i po niewielkich przekształceniach otrzymuje się podobne wyrażenie na transformację admitancji Y(l) wzdłuż bezstratnej prowadnicy falowej o impedancji charakterystycznej Z0 i długości l.

  

\mathrm{Y}(l)=\mathrm{Y_0}\frac{\mathrm{Y_L}+j\mathrm{Y_0} \mathrm{tg}\beta l}{\mathrm{Y_0}+j\mathrm{Y_L} \mathrm{tg}\beta l};

(10-67)  

Analizując wyrażenie (10-66) prowadzi do kilku wniosków:

  • Impedancja Z(l) jest funkcją aż 3 zmiennych: ZL, Z0, \betaI.
  • Impedancja Z(I) jest okresową funkcją odległości, Z(I) = Z(I+\lambda/2), a okresem jest pół fali \lambda/2.

 
Rys.10.13. Ilustracja transformacji impedancji. A) Transformacja na płaszczyźnie współczynnika odbicia. B) Transformacja i na płaszczyźnie impedancji.


Na rys.10.13 pokazano ilustrację graficzną procesu transformacji impedancji przy odsuwaniu się obciążenia w stronę generatora (patrz rys.10.9). Na rys.10.13A widzimy klasyczną transformację współczynnika odbicia Γ(l). Jak przy tej transformacji zmienia się impedancja z(l) pokazuje ry.10.13B. Okrąg transformacji umieszczono na płaszczyźnie impedancji. Zaznaczono charakterystyczne punkty tego okręgu: rMAX i rMIN.
Zależność (10-66) wskazuje na bardzo interesujące właściwości linii długiej, umożliwiające komponowanie żądanych parametrów obwodów. Zapoznamy się z nimi w kolejnych punktach.
 

4.7. Przykłady zastosowań

Transformacja impedancji – szczególne przypadki
Przypadek 1: Linia długa jest zakończona impedancją ZL =Z0. W takim przypadku, zgodnie z (10-66), Z(l)=Z0
Wniosek: W każdym punkcie linii impedancja ma tą samą wartość. 
Przypadek 2: Obliczymy impedancję w odległości równej wielokrotności pół fali l=n\lambda/2 od obciążenia. Łatwo zauważyć, że  Z(l=n\lambda/2) = ZL, impedancja okresowo przyjmuje taką wartość, jaką ma no końcu linii.
Wniosek: Linia o długości n\lambda/2 jest - z punktu widzenia transformacji impedancji - przezroczysta.
Przypadek 3: Obliczymy impedancję w odległości równej ćwierć fali l=\lambda/4 od obciążenia. 

  

\mathrm{Z}(l=\lambda /4)=\frac{\mathrm{Z_L}}{\mathrm{Z_L}};\, \, \mathrm{Z}(l=\lambda /4)\mathrm{Z_L}=\mathrm{Z_0}^{2};

(10-68)  

Linia o długości l=(2n-1)\lambda/4 ma specjalne właściwości i dlatego nazywana jest transformatorem ćwierćfalowym.  Linia.
Wnioski:

  • Transformator ćwierćfalowy jest inwerterem impedancji. Zamienia on duże (małe) wartości rezystancji na rezystancje małe (duże).
  • Transformator ćwierćfalowy zamienia impedancje obciążenia o charakterze indukcyjnym (pojemnościowym) na impedancje wejściowe pojemnościowe (indukcyjne).
  • Jeśli obciążeniem linii jest obwód rezonansu szeregowego, to impedancja wejściowa zachowuje się jak dla obwodu rezonansu równoległego, i vice versa.

Przypadek 4: W ogólnym przypadku obciążenia linii impedancją ZL=RL+jXL, gdy RL>0, to współczynnik odbicia równy jest wtedy\left | \Gamma \right |=\left | \Gamma \exp (j\psi _{L}) \right |, przy czym \left | \Gamma_L \right |(patrz rys.10.5). W miarę odsuwania się od obciążenia zmienia się Arg{\Gamma}. Gdy odsuniemy się na odległość l1, dla której spełniony jest warunek (10-69):

  

\psi _{L}-2\beta l_{1}=2n\pi ;

(10-69)  

to napięcie U(l1) i prąd I(l1) są w fazie. Oznacza to, że impedancja Z(l1) jest czysto rzeczywista i równa: 

  

\mathrm{Z}(l_1))=\mathrm{Z_0}\frac{1+\Gamma (l)}{1-\Gamma (l)}=\varrho \mathrm{Z_0};

(10-70)  

gdzie \varrho jest współczynnikiem fali stojącej na linii.
Podobnie, gdy odsuniemy się na odległość l2, dla której spełniony jest warunek (10-71):

  

\psi _{L}-2\beta l_{2}=(2n+1)\pi ;

(10-71)  

sytuacja powtarza się i także wtedy napięcie U(l1) i prąd I(l1) są w fazie, a więc:

  

\mathrm{Z}(l_2))=\mathrm{Z_0}\frac{1+\Gamma (l)}{1-\Gamma (l)}= \frac{\mathrm{Z_0}}{\varrho};

(10-72)  

Oba miejsca l1 i l2 oddalone są od siebie o ćwierć długości fali \lambda/4. Oba te przypadki mogą być wykorzystane przy projektowaniu obwodów dopasowujących. 
     Linia zwarta na końcu
Rozważymy efekty zachodzące w linii długiej zwartej na końcu. Oznacza to, że:  \mathrm{Z_L}=0\, i\, \Gamma _{\mathrm{L}}=-1 . Zgodnie z zależnością (10-66) impedancja wejściowa linii zwartej na końcu jest w każdym miejscu czystą reaktancją:

  

\mathrm{Z}(l))=j\mathrm{X}(l)=j\mathrm{Z_0}\mathrm{tg}\beta l;

(10-73)  

Rozkład prądu i napięcia dla linii zwartej na końcu pokazano na rys.10.14. 
Prąd I(l) i napięcie U(l) są przesunięte w fazie o \pi/2, a kolejne zera napięcia lub prądu odległe są od siebie o \lambda/2.

  

\mathrm{I}(l)=\mathrm{I_L}\cos \beta l;

\mathrm{U}(l)=j\mathrm{I_LZ_0}\sin \beta l;

(10-74)  

Kąt fazowy między prądem I(l) i napięciem U(l) jest cały czas równy 900, jednakże co ćwierć fali zmienia się jego znak. Dlatego X(l) ma dla pewnych zakresów l charakter indukcyjny, dla innych pojemnościowy, co pokazano na rys.10.14B.
 
 


Rys.10.14. Linia długa zwarta na końcu. 
A) Rozkład prądu i napięcia dla linii.
B) Reaktancja wejściowa linii zwartej na końcu.  

 
Zastępcze wartości indukcyjności Leq i pojemności Ceq znajdujemy ze wzorów:

  

L_{eq}=\frac{\mathrm{Z_0 tg}\beta l}{\omega };

(10-75)  
  

C_{eq}=\frac{1}{\omega\mathrm{Z_0 tg}\beta l };

(10-76)  

Dla małych długości l, gdy \beta<0.5 (l<0.08\lambda) to tg\betal~\betal i indukcyjność Leq tego odcinka zapisze się następująco (vf jest prędkością fazową fali):

  

\omega L_{eq}\cong =\mathrm{Z_0}\beta l=2\pi \frac{\mathrm{Z_0}l}{\lambda _{f}}=\omega \mathrm{Z_0}\frac{l}{v_f};

(10-77)  

Dla zakresów częstotliwości w sąsiedztwie \beta=(2n-1)\pi/2 (nieparzysta liczba ćwiartek fali) linia zwarta na końcu zachowuje się jak obwód rezonansu równoległego.    
Dla zakresów częstotliwości w sąsiedztwie \beta=n\pi (wielokrotność połowy fali) linia zwarta na końcu zachowuje się jak obwód rezonansu szeregowego.
     Linia rozwarta na końcu
Impedancja wejściowa linii rozwartej na końcu zapisuje się zależnością: 

  

\mathrm{Z}(l)=j\mathrm{X}(l)=-j\mathrm{Z_0}\mathrm{ctg}\beta l;

(10-78)  

Charakter zmian prądu I(I) i napięcia U(I jest taki, jak na rys.10.14A, z tą różnicą, że na końcu linii rozwartej I(I=0)=0. Prąd i napięcie w każdym miejscu linii przesunięte są w fazie o \pi /10. 
Charakter zmian impedancji jak dla linii zwartej, tylko przesunięty o \lambda/4. W zależności od l linia raz jest pojemnością, raz indukcyjnością rys.10.15. 

 


Rys.10.15. Reaktancja wejściowa linii rozwartej na końcu.    


Podobnie jak w przypadku linii zwartej, linia rozwarta może realizować pojemności i indukcyjności, zależnie od odległości od rozwarcia.
Dla małych długości l, gdy l<0.08\lambda i słuszne jest przybliżenie tg\beta l~\beta l, odcinek linii rozwartej można zastąpić równoważną pojemnością Ceq.

  

C_{eq}\cong \frac{\omega l}{\mathrm{Z_0}v_{f}}=\frac{2\pi l}{\mathrm{Z_0}\lambda _f};

(10-79)  

Dla zakresów częstotliwości w sąsiedztwie \beta l = (2n-1)\pi/2 (nieparzysta liczba ćwiartek fali) linia rozwarta na końcu zachowuje się jak obwód rezonansu szeregowego.
Dla zakresów częstotliwości w sąsiedztwie \beta l = n\pi (wielokrotność połowy fali) linia rozwarta na końcu zachowuje się jak obwód rezonansu równoległego.