Numerische Berechnung von Nk e -bs ∫αβ e-s t - c √[t2+2bt] T(M,b;t) Jν(Nkr [t2 + 2 bt + 1]1/2) dt :

T(M,b;t) = [2M(t+b)]/[M(t+b)+(t2+2bt)½]


Autor: Klaus Rottbrand ©
Stand: 28 NOV 2005



b = (N2-1)½/N
a = Nkr
M = Dichteverhältnis
N = Brechungsindex >1
s = Nk |z0|: Laplace-Parameter (z0 ist z-Ort des Dipols)
|z|: Ort (r,z) Zylinderkoordinaten
c = Nk|z|
Es sind Ordnungen ν ≤ 100 zugelassen.


Untere Int.-Grenze α =  
Obere Int.-Grenze β =  
Bessel-Ordnung ν =  
Laplace-Parameter |z0| =  
Entfernung-Radius r =  
Wellenzahl k =  
Brechungsindex N >1 =  
Dichte-Faktor M =  
Orts-Parameter |z| =  
Teilintervall-Anzahln = 


Beschreibung:

Dieses JavaScript berechnet Laplace-transformierte, welche dem inhomogenen Anteil des transmittierten Feldes beim Sommerfeldproblem entsprechen (Raumwellen). Bei Feldberechnungen ist die Ordnung ν=0 und die untere Integrationsgrenze α=0 zu nehmen, sowie der Brechungsindex N > 1 zu wählen. Es handelt sich um einen akustischen Dipol. Das Druckfeld ist gleich dem Integral

∫Nk∞ T(M,N,k;x) e-|z0|√[x2-k2] -|z|√[x2-N2kk2] J0(rx) x dx/√[x2-k2]

T(M,N,k;x) = 2M √[x2-k2]/ (M √[x2-k2] + √[x2-N2k2])



Literatur
Milton Abramowitz und Irene A. Stegun: Handbook of Mathematical Functions, Dover (101972), Formel 25.4.20 auf Seite 887: Integration in N Teilintervallen mit Newton-Cotes Formel vom geschlossenen Typ mit n=10 (Teile) genommen und anschließend über die Teilintervalle aufsummiert wird.
Hosted by www.Geocities.ws

1