Strahlungsfeld eines Hertzschen Dipols 

(mit Maple) 

Ein Hertzscher Dipol besteht aus einer negativen Punktladung am Koordinatenursprung und einer harmonisch bewegten positiven Punktladung auf a(t). 

1) Hertzscher Dipol 

Bahnkurve der Ladung 

a_(t) = `*`(a[0], `*`(cos(`*`(omega, `*`(t))), `*`(_k))) (1.1)
 

Ladungsdichte 

Typesetting:-mprintslash([rho(r_, t) = `*`(q, `*`(`+`(delta(`+`(r_, `-`(a_))), `-`(delta(r_)))))], [rho(r_, t) = `*`(q, `*`(`+`(Dirac(`+`(r_, `-`(a_))), `-`(Dirac(r_)))))]) (1.2)
 

Strom 

Typesetting:-mprintslash([j_(r_, t) = `*`(q, `*`(diff(a_(t), t), `*`(delta(`+`(r_, `-`(a_(t)))))))], [j_(r_, t) = `*`(q, `*`(diff(a_(t), t), `*`(Dirac(`+`(r_, `-`(a_(t)))))))]) (1.3)
 

oder mit Bahnkurve 

Typesetting:-mprintslash([j_(r_, t) = `+`(`-`(`*`(q, `*`(a[0], `*`(sin(`*`(omega, `*`(t))), `*`(omega, `*`(_k, `*`(delta(`+`(r_, `-`(`*`(a[0], `*`(cos(`*`(omega, `*`(t))), `*`(_k))))))))))))))], [j_(r... (1.4)
 

Vektorpotential 

Typesetting:-mprintslash([A_ = `+`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(Int(Int(`/`(`*`(j_(rb_, tb), `*`(delta(`+`(t, `-`(`/`(`*`(r), `*`(c))), `-`(tb))))), `*`(Norm(`+`(rb_, `-`(r_))))), rb_ = `+`(`-`(i... (1.5)
 

2) Vektorpotential in Nahfeldnäherung 

Fernfeldnäherung 

Typesetting:-mprintslash([Norm(`+`(rb_, `-`(r_))) = r], [Norm(`+`(rb_, `-`(r_))) = r]) (2.1)
 

Einsetzen in Potential 

Typesetting:-mprintslash([A_ = `+`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(Int(Int(`/`(`*`(j_(rb_, tb), `*`(delta(`+`(t, `-`(`/`(`*`(r), `*`(c))), `-`(tb))))), `*`(r)), rb_ = `+`(`-`(infinity)) .. infinity)... (2.2)
 

Einsetzen des Stromes 

Typesetting:-mprintslash([A_ = `+`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(Int(Int(`+`(`-`(`/`(`*`(q, `*`(a[0], `*`(sin(`*`(omega, `*`(tb))), `*`(omega, `*`(_k, `*`(delta(`+`(rb_, `-`(`*`(a[0], `*`(cos(`*`(... (2.3)
 

Ausführen desIntegrals 

A_ = `+`(`-`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(_k, `*`(sin(`*`(omega, `*`(`+`(t, `-`(`/`(`*`(r), `*`(c)))))))))))))), `*`(r, `*`(Pi))))) (2.4)
 

In Basis der Kugelkoordinaten 

A_ = `+`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(cos(theta), `*`(sin(`/`(`*`(omega, `*`(`+`(`-`(`*`(t, `*`(c))), r))), `*`(c))), `*`(_r)))))))), `*`(r, `*`(Pi))), `-`(`/`(`*`(`/`... (2.5)
 

3) B- und E-Feld: 

Magnetisches Feld als Rotation des Vektorpotentials 

Typesetting:-mprintslash([B_ = Curl(A_)], [B_ = Curl(A_)]) (3.1)
 

also 

B_ = `+`(`-`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(sin(theta), `*`(`+`(`*`(omega, `*`(cos(`/`(`*`(omega, `*`(`+`(`-`(`*`(t, `*`(c))), r))), `*`(c))), `*`(r))), `-`(`*`(sin(`/`(... (3.2)
 

Rotation des B ist gleich der zeitlichen Änderung des E-Feldes, nach Zeitintagration also 

Typesetting:-mprintslash([E_ = Int(Curl(B_), t)], [E_ = Int(Curl(B_), t)]) (3.3)
 

Das heißt 

E_ = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(cos(theta), `*`(_r, `*`(`+`(`-`(`*`(r, `*`(sin(`+`(`-`(`*`(omega, `*`(t))), `/`(`*`(omega, `*`(r)), `*`(c))))))), `-`(`/`(`*`...
E_ = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(cos(theta), `*`(_r, `*`(`+`(`-`(`*`(r, `*`(sin(`+`(`-`(`*`(omega, `*`(t))), `/`(`*`(omega, `*`(r)), `*`(c))))))), `-`(`/`(`*`...
(3.4)
 

4) Graphische Darstellung der Felder auf der x-Achse 

Parameter: 

Typesetting:-mprintslash([`:=`(Par, q = 1, mu[0] = 1, a[0] = 1, omega = 1, c = 1)], [q = 1, mu[0] = 1, a[0] = 1, omega = 1, c = 1]) (4.1)
 

x-z-Ebene in Kugelkoordinaten 

Typesetting:-mprintslash([`:=`(XZ_Plane, phi = 0, sin(theta) = `/`(`*`(x), `*`(r)), cos(theta) = `/`(`*`(z), `*`(r)), r = `*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2))))], [phi = 0, sin(thet... (4.2)
 

y-Komponete des B-Feldes mit Parametern in XZ_Ebene 

B[y] = `+`(`-`(`/`(`*`(`/`(1, 4), `*`(x, `*`(cos(`+`(t, `-`(`*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2))))))))), `*`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `*`(Pi)))), `-`(`/`(`*`(`/`(1, 4), ... (4.3)
 

z-Komponente des E-Felder mit Parametern in XZ_Ebene 

E[z] = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(z, 2), `*`(sin(`+`(t, `-`(`*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2))))))))), `*`(Pi, `*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), 2))))), `/`(`*`...
E[z] = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(z, 2), `*`(sin(`+`(t, `-`(`*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2))))))))), `*`(Pi, `*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), 2))))), `/`(`*`...
E[z] = `+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(z, 2), `*`(sin(`+`(t, `-`(`*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2))))))))), `*`(Pi, `*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), 2))))), `/`(`*`...
(4.4)
 

B_x und -E_z zum Vergleich von Nah- und Fernzone auf der x-Achse 

Plot_2d
 

Bemerkung:  

Die Felder wurden in der Fernfeldnäherung berechnet. Bei sehr kleinen Abständen sind die E- und B-Felder auch in dieser Näherung nicht in Phase.Man sieht in der Darstellung jedoch sehr schön, wie schnell sich die gemeinsame Phase in den Fernfeldern einstellt 

3-D-Darstellung von By und Ez auf der x-Achse 

 

[x, B[y], E[z]]
Plot
 

Das E-Feld schwingt in z-Richtung, das B-Feld in y-Richtung. 

5) E-Feldlinien als Höhenlinien des B-Feldes in der x-z-Ebene 

Definition der Feldlinien des E-Feldes: 

Typesetting:-mprintslash([`&x`(E_, dr_) = 0], [`&x`(E_, dr_) = 0]) (5.1)
 

Mit E-Feld 

Typesetting:-mprintslash([E_ = Int(Curl(B_), t)], [E_ = Int(Curl(B_), t)]) (5.2)
 

also 

Typesetting:-mprintslash([`&x`(Int(Curl(B_), t), dr_) = 0], [`&x`(Int(Curl(B_), t), dr_) = 0]) (5.3)
 

Beziehungsweise 

Typesetting:-mprintslash([Int(`&x`(dr_, Curl(B_)), t) = 0], [Int(`&x`(dr_, Curl(B_)), t) = 0]) (5.4)
 

In der x-z-Ebene ist 

B_ = `*`(B[y](x, z, t), `*`(_j)), dr_ = `+`(`*`(dx, `*`(_i)), `*`(dz, `*`(_k))) (5.5)
 

also 

Typesetting:-mprintslash([`+`(`-`(`*`(_j, `*`(Int(`+`(`*`(dz, `*`(diff(B[y](x, z, t), z))), `*`(dx, `*`(diff(B[y](x, z, t), x)))), t))))) = 0], [`+`(`-`(`*`(_j, `*`(Int(`+`(`*`(dz, `*`(diff(B[y](x, z,... (5.6)
 

und  

Int(dB_y, t) = 0 (5.7)
 

Die Feldlinien von E können in der x-z-Ebene also dargestellt werden als Höhenlininen von 

B[y, int] = Int(B[y], t) (5.8)
 

6) Animation der Feldlinien des E-Feldes in der x-z-Ebene 

B-Feld 

B_ = `+`(`-`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(sin(theta), `*`(`+`(`*`(omega, `*`(cos(`/`(`*`(omega, `*`(`+`(`-`(`*`(t, `*`(c))), r))), `*`(c))), `*`(r))), `-`(`*`(sin(`/`(... (6.1)
 

y-Komponente des B-Feldes 

B_y = `+`(`-`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(`^`(omega, 2), `*`(sin(theta), `*`(cos(phi), `*`(cos(`/`(`*`(omega, `*`(`+`(`-`(`*`(t, `*`(c))), r))), `*`(c))))))))))), `*`(Pi, `*`(c,... (6.2)
 

Zeitintegral über B-Feld 

B[y, int] = int(B_y, t) (6.3)
 

also 

B[y, int] = `+`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(sin(theta), `*`(cos(phi), `*`(sin(`+`(`-`(`*`(omega, `*`(t))), `/`(`*`(omega, `*`(r)), `*`(c)))))))))))), `*`(Pi, `*`(c, `... (6.4)
 

x-z-Ebene: 

B[y, int] = `+`(`-`(`/`(`*`(`/`(1, 4), `*`(mu[0], `*`(q, `*`(a[0], `*`(omega, `*`(x, `*`(sin(`+`(`*`(omega, `*`(t)), `-`(`/`(`*`(omega, `*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2)))), `*`(c... (6.5)
 

Plot-Funktion 

`+`(`-`(`/`(`*`(`/`(1, 4), `*`(x, `*`(sin(`+`(t, `-`(`*`(`^`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `/`(1, 2))))))))), `*`(`+`(`*`(`^`(x, 2)), `*`(`^`(z, 2))), `*`(Pi)))), `/`(`*`(`/`(1, 4), `*`(x, `*`(... (6.6)
 

 

Typesetting:-mprintslash([`:=`(XMax, 20)], [20])
Typesetting:-mprintslash([`:=`(Amplitude[Max], 0.5e-1)], [0.5e-1]) (6.7)
 

Plot_2d