Numerické modelování elektromagnetických polí se s rozvojem výpočetní techniky a neustále rostoucí výkonností počítačů stalo spolu s optimalizačními technikami nepostradatelnou složkou návrhu konstrukcí nových elektrotechnických a elektronickýchzařízení i zařízení z ostatních oblastí technické praxe. Numerické modelování je také bezesporu nedílnou součástí komplexních analýz chování časoprostorových polí, které jsou důležité pro posouzení nových požadavků na kvalitu zařízení jako je elektromagnetická kompatibilita. Složité problémy řešené v současné technické praxi nelze zvládnout ve většině případů jinými prostředky než pomocí vhodných numerických metod za použití výkonných počítačů.
.
Metoda nejmenších čtverců
integruje kvadrát zbytku přes celou oblast.
Oblast rozdělíme prvků uzly.,, 2
21
I kvadratickou formou uzlových potenciálů tvaru
022113132112
2
222
2
111 .. permitivitě předpokládáme, je
na prvku konstantní nezávisí intenzitě pole.8. integraci dostaneme výraz, který již neobsahuje souřadnice ani derivace
podle nich
( )∫=
Ω
Ωφφφφφ dz,y,x,I NU1NU ..8: okrajovým podmínkám
( )∑∑
==
−=−≈
NU
j
jj
NU
j
jj zyxNzyxN
11
,,grad,,grad εφφεD
Operátor derivace grad aplikuje známou funkci Nj, derivace tak rovnice
odstraněna. Uveďme stručně
některé metody.. Dosazením aproximace rovnice pro
indukci nebude tato pro libovolné splněna přesně, ale vykazuje zbytek res(x, φj)
(reziduum), který funkcí souřadnic uzlových potenciálů
( )NU1
NU
j
jj z,y,x,zyxN φφρεφρ . 4.. 4...
φ1, Γe1
Ε
Γn
Ε, Γn
un
un
φ0, Γe0
Ω
...,res,,graddivdiv
1
=−−=− ∑
=
D
Je třeba nalézt takové hodnoty φj, NU, které dají minimální zbytek..
V Poissonově rovnici, napsané tvaru
( ρφεε =−== graddivdivdiv ED
aproximujeme indukci
Obr.. Protože každá aproximační funkce součtem tvarových funkcí derivace součtu
je rovna součtu derivací, derivují tyto tvarové funkce.. Výpočet otevřených úloh (pole
venkovních zařízení, pole antén) vyžaduje doplňující matematický aparát, zajišťující expanzi
pole neomezené oblasti. cccccccI ++++++++= φφφφφφφ
Koeficienty cij jsou integrály aproximačních funkcí jejich derivací přes oblast Ω.
MKP vhodná výpočtu polí uzavřených oblastech. Proto postačí počítat pole
v libovolné výseči mezikruží. Podrobnosti lze nalézt monografiích MKP. Pole není potřeba počítat celé oblasti, neboť ohledem symetrii
známe tvar siločar Plochy směru siločar tvoří hranici Γn. Integrují známé funkce jejich derivace
podle souřadnic.
Potenciál aproximujeme
∑
=
=≈
NU
j
jj zyxN
1
a ),,(φφφ
z uzlových hodnot aproximačními funkcemi Nj..Modelování elektromagnetických polí 43
Význam oblasti naznačený Obr..,res. něm oblastí část mezikruží
v koaxiálním kabelu. Různé způsoby
minimalizace zbytku dají soustavy rovnic pro hledané uzlové potenciály