Über ein Elektrostatikintegral und einen Deltafunktionskern

Ich habe Probleme mit einem Integral und hätte gerne einige Hinweise, wie man es "nimmt":

3 γ A 2 D 3 R 4 π ( R 2 + γ 2 C 2 ( R u ) 2 + A 2 ) 5 2
Hier u , A Und γ sind Konstanten, und der Integrand konvergiert gegen ein Dirac-Delta δ ( R ) als A 0 . Das Integral muss gleich 1 sein.

Antworten (2)

Wählen Sie zunächst eine Richtung für u entlang der z-Achse. Dann ist das Integral

ICH = 1 ( X 2 + j 2 + A z 2 + B ) 5 / 2 D X D j D z

z neu skalieren um A um A loszuwerden und die Rotationsinvarianz wiederherzustellen.

ICH = 1 A 1 ( X 2 + j 2 + z 2 + B ) 2.5 D X D j D z

Jetzt finden Sie die B-Abhängigkeit sofort aus der Umskalierung von xy und z durch B (oder aus der Dimensionsanalyse - B hat Längeneinheiten zum Quadrat):

ICH = 1 A B 1 ( R 2 + 1 ) 2.5 D 3 R

Das einzige Unbestimmte ist der transzendente Faktor, der nur eine Zahl ist. Sie werten es aus, indem Sie es radial tun und eine Reihe von Substitutionen durchführen:

A B ICH = 4 π 0 R 2 ( R 2 + 1 ) 2.5 D R

Erste u = R 2 + 1 gibt

A B ICH = 2 π 1 u 1 ( u ) 2.5 D u

Dann v = 1 u macht es,

A B ICH = 2 π 0 1 1 v D v = 4 π 3

So

ICH = 4 π 3 A B

Danke schön! Das Integral ist gleich eins:
u = R 2 + 1 R = u 1 .
Ja, das ist wahrscheinlich sauberer.
@EmilioPisanty: Es ist nicht nur sauberer, dein Weg funktioniert nicht --- die sphärische Integration wird wahrscheinlich keine Funktion mit elementarer Stammfunktion ohne die Neuskalierung des z-Schritts erzeugen.
Mathematica hat keine Probleme mit der Integration. Ich kann die Schritte für die folgen θ Integration und obwohl ich das nicht versucht habe R Integral Die Stammfunktion ist kompliziert, aber elementar.
@EmilioPisanty: Ich bin überrascht --- die Winkelintegration ohne Neuskalierung von Z führt zu komplizierten Quadratwurzeln im Nenner, und ich verstehe nicht, warum die Stammfunktion existieren sollte. Haben Sie versucht, mithilfe von Mathematik eine Stammfunktion für das r-Integral zu erhalten oder nur das Integral über Theta und r zu berechnen? Mathematica führt ein bestimmtes Integral unter Verwendung von hypergeometrischen Funktionen in Zwischenschritten aus und reduziert dann mit hypergeometrischen Sonderwerten und Identitäten, sodass es nicht klar ist, dass die Stammfunktion elementar ist, da Mathematica die richtige Antwort erhält.
@RonMaimon siehe Bearbeiten in meiner Antwort.
Sehr schöne Antwort.

Stellen Sie die ein z Achse entlang der Richtung von u und sphärische Koordinaten verwenden, was Ihr Integral auf etwas wie reduziert

0 D R 0 π D θ 0 2 π D ϕ R 2 Sünde ( θ ) ( A 2 + R 2 ( 1 + γ 2 C 2 cos 2 ( θ ) ) ) 5 / 2 .
Mach das ϕ erst ganzzahlig und dann θ integral, Transformation zu u = cos ( θ ) . Achten Sie darauf, bei Bedarf absolute Werte für die Wurzeln zu verwenden. Danach die R Integral sollte hart, aber machbar sein.


BEARBEITEN, um einige Diskussionen von Kommentaren zu nehmen. Ich gab

Integrate[Sin[[Theta]]/(a^2+r^2 (1+[Gamma]^2 Cos[[Theta]]^2))^(5/2),[Theta]]/.[Theta ]->[Pi]

zu Mathematica zu bekommen

3 ( A 2 + R 2 ) + 2 R 2 γ 2 3 ( A 2 + R 2 ) 2 ( A 2 + R 2 + R 2 γ 2 ) 3 / 2 ,
und dann die radiale Integration durch

Integriere[r^2 (3 (a^2+r^2)+2 r^2 [Gamma]^2)/(3 (a^2+r^2)^2 (a^2+r^2+ r^2 [Gamma]^2)^(3/2)),r]

gibt

R 3 3 A 2 ( A 2 + R 2 ) A 2 + R 2 ( 1 + γ 2 )
für die Stammfunktion. Ich stimme zu, dass die Wurzeln nichtelementare Stammfunktionen vermuten lassen, aber es ist nur eine Wurzel, so dass elliptische Integrale ausfallen. Sobald man das Antiderivativ hat, ist es natürlich Routine zu überprüfen, ob es so differenziert, wie es sollte.

Danke schön. Aber nach einigen Ersetzungen erhielt ich
C γ u R ( A 2 + R 2 ) 2 D T ( 1 + T 2 ) 5 2 ,
Wo
T = T R u γ C A 2 + R 2 .
Aber nach der Integration bekam ich null. Können Sie mir sagen, wo ich den Fehler gemacht habe?
Zum einen integrieren Sie wahrscheinlich T aus wo Sie integrieren sollten R = 0 . Zweitens kann man das nicht ausklammern R -abhängiger Term aus dem Integral!
Verzeihung. ich meinte
T = R v γ u C A 2 + R 2 ,
Wo
v = C Ö S ( θ ) .
ist das der R integral oder die θ Integral Sie versuchen zu tun?
Wenn du das machst θ integral, dann sind Sie auf dem richtigen Weg, aber Sie müssen mit den Grenzen der Integration vorsichtig sein. Einstellung T = bräunen ( η ) sollte das Integral auf machbare trigonometrische Integrale reduzieren. Aber achten Sie auf die Integrationsgrenzen und auf die Parität des Integranden!
Ja, ich meinte Theta-Integral. Nach der Integration bekam ich
1 12 ( 3 S ich N ( ε ) + S ich N ( 3 ε ) ) .
Nach dem Ersetzen der Rahmenintegration erhielt ich Werte wie
S ich N ( 3 R γ u C ( A 2 + R 2 ) ) .
Was kann ich damit machen? Weiter integrieren? Oder habe ich einen Fehler gemacht?
nicht verwenden Sünde ( 3 ϵ ) sondern halte dich lieber dran Sünde ( ϵ ) . Danach verwenden Sie die Identitäten Sünde ( ϵ ) = T / 1 + T 2 Und cos ( ϵ ) = 1 / 1 + T 2 . (Denken Sie an Dreiecke!)
Sie haben also empfohlen, den Faktor 3 um die Gleichheit zu reduzieren
S ich N ( 3 ε ) = 3 S ich N ( ε ) 4 S ich N 3 ( ε ) ,
und danach verwenden gleich
S ich N ( A R C T G ( T ) ) = 1 1 + C T G 2 ( A R C T G ( T ) ) = T T 2 + 1 ?
Das ist nicht der richtige Weg.
Aber was ist der richtige Weg?
Du siehst für mich auf dem richtigen Weg aus. Nach der Winkelintegration sollte etwas von der Form übrig bleiben
0 R 2 ( 3 A 2 + R 2 ( 3 + 2 γ 2 ) ) 3 ( A 2 + R 2 ) 2 ( A 2 + R 2 ( 1 + γ 2 ) ) 3 / 2 D R