Doppelintegral ∫01∫01(xy)s−log(xy)√dxdy∫01∫01(xy)s−log⁡(xy)dxdy\int\limits_0^1\!\!\int\limits_0^1\frac {(xy)^s}{\sqrt{-\log(xy)}}\,dx\,dy

Ist es möglich, eine geschlossene Form des folgenden Integrals zu erhalten?

ICH = 0 1 0 1 ( X j ) S Protokoll ( X j ) D X D j ( S > 0 ) .
Mein Versuch : Ich habe versucht, Variablen von kartesischen Koordinaten in Polarkoordinaten zu ändern:
X = R cos ( θ ) j = R Sünde ( θ ) .
Ich habe den Jacobianer berechnet:
J = | D ( X , j ) D ( R , θ ) | = | cos ( θ ) R cos ( θ ) ( R Sünde ( θ ) ) Sünde ( θ ) | = R .
Ab hier stecke ich fest.

Sie integrieren über ein Quadrat, daher haben die Integrationsgrenzen keine Rotationssymmetrie um den Ursprung, daher wird die Verwendung von Polarkoordinaten nicht empfohlen. Die Funktion hat auch keinerlei Rotationssymmetrie. Ich denke also nicht, dass Polarkoordinaten ein guter Weg sind, hier fortzufahren.
Das Problem ist, dass das innere Integral einen ziemlich erschreckenden Ausdruck liefert, der Fehlerfunktionen enthält. Ich habe im Moment keine Möglichkeit, das zu berechnen
Ein CAS kehrt zurück π / ( 2 ( 1 + S ) 3 / 2 ) für S > 1 . Also muss es jetzt von Hand gemacht werden :)
Ja, Mathematika scheint damit fertig zu werden :)

Antworten (4)

Hier ist ein Ansatz.

Wir haben z S > 0 , σ > 0 ,

(1) 0 1 0 1 X S 1 j σ 1 Protokoll ( X j ) D X D j = π S σ + σ S .

Nachweisen. Das mag man beobachten

(2) 1 A = 2 π 0 e A T 2 D T , A > 0.
Annehmen S > 0 , σ > 0 Und σ S . Dann erhalten wir
[ 0 , 1 ] × [ 0 , 1 ] X S 1 j σ 1 Protokoll ( X j ) D X D j = [ 0 , 1 ] × [ 0 , 1 ] X S 1 j σ 1 ( 1 Protokoll ( X j ) ) D X D j = [ 0 , 1 ] × [ 0 , 1 ] X S 1 j σ 1 ( 2 π 0 e T 2 Protokoll ( X j ) D T ) D X D j = 2 π 0 [ 0 , 1 ] × [ 0 , 1 ] X S 1 j σ 1 × X T 2 j T 2 D X D j D T = 2 π 0 ( 0 1 X T 2 + S 1 D X ) × ( 0 1 j T 2 + σ 1 D j ) D T = 2 π 0 1 ( T 2 + S ) ( T 2 + σ ) D T = 2 π ( σ S ) ( 0 1 ( T 2 + S ) D T 0 1 ( T 2 + σ ) D T ) = 2 π ( σ S ) ( π 2 S π 2 σ ) = π S σ + σ S .
Das gewünschte Integral erhalten wir durch Setzen S := S + 1 In ( 1 ) Und σ := S + 1 geben

(3) 0 1 0 1 ( X j ) S Protokoll ( X j ) D X D j = π 2 ( 1 + S ) 3 / 2 .

(+1) Ja, du musst nicht all dieses Fehlerfunktions-Zeug durchgehen, das ich ohne Mathematik total vermasselt habe :)

u = X j v = j X D u D v = ( u , v ) ( X , j ) D X D j = 2 j X D X D j = 2 v D X D j D u D v 2 v = D X D j

0 1 u 1 / u u S Protokoll u D v D u 2 v = 0 1 u S Protokoll u 1 2 ( Protokoll 1 u Protokoll u ) D u (1) = 0 1 u S Protokoll u D u

w = Protokoll u e w 2 = u 2 w e w 2 D w = D u

So ( 1 ) wird 0 e S w 2 w ( 2 w e w 2 D w ) = 2 0 w 2 e ( S + 1 ) w 2 D w .

Wir können dies auf den Erwartungswert einer Funktion einer Zufallsvariablen mit einer Standardnormalverteilung zurückführen:

( S + 1 ) w 2 = z 2 2 D w = D z 2 S + 1

Unser Integral wird

(1) 2 0 z 2 2 ( S + 1 ) e z 2 / 2 D z 2 S + 1 = 1 ( S + 1 ) 3 / 2 2 0 z 2 e z 2 / 2 D z .
Das letzte Integral ist
0 z 2 ( e z 2 / 2 2 π ) D z 2 π .

Dieses Integral ist der halbe Erwartungswert von a χ 1 2 zufällige Variable; somit ist es gleich 1 / 2 , also bekommen wir ( 1 / 2 ) 2 π , also der ganze ausdruck ( 1 ) Ist

1 ( S + 1 ) 3 / 2 2 1 2 2 π = π 2 ( S + 1 ) 3 / 2 .

Ich mochte Ihre Lösung (und es hat mich dazu gebracht, meine eigene nicht vollständig fertigzustellen). Kleine Sache: ist nicht das letzte Integral gleich π / 2 ?

Das zu zeigen ist unser Ziel

0 1 0 1 ( X j ) S Protokoll X j D X D j = 2 0 1 j 1 ( X j ) S Protokoll X j D X D j = π 2 ( 1 + S ) 3 / 2 .

Wir brauchen die Fehlerfunktion,

erf ( X ) = 2 π 0 X e T 2 D T .
Beginnen wir mit der Berechnung eines Grundelements in Bezug auf X . Mit
u = 1 + S Protokoll ( X j )
wir glauben, dass
( X j ) S Protokoll X j D X = π j 1 + S erf ( 1 + S Protokoll ( X j ) ) .
Grenzen einfügen, das finden wir
j 1 ( X j ) S Protokoll X j D X = π j 1 + S erf ( 1 + S Protokoll ( j ) ) + π j 1 + S erf ( 1 + S 2 Protokoll ( j ) )
Als nächstes wollen wir bzgl. integrieren j . Die Integrale sind sehr ähnlich. Wir glauben, dass ( u = A Protokoll j und zweimalige partielle Integration)
1 j erf ( 1 + S A Protokoll j ) D j = 2 A u erf ( 1 + S u ) D u = 1 A u 2 erf ( 1 + S u ) + 1 + S A π u 2 e ( 1 + S ) u 2 D u = 1 A u 2 erf ( 1 + S u ) 1 2 A π 1 + S u e ( 1 + S ) u 2 + 1 4 A ( 1 + S ) erf ( 1 + S u ) .
Das Einfügen überlasse ich Ihnen A = 1 Und A = 2 , und mit den Grenzen zu arbeiten.

Protokoll ( ) scheint mir die Umkehrfunktion einer Gaußschen Funktion zu sein.

Es kann sich also lohnen einzustellen X = e u 2 Und j = e v 2 zu bekommen:

ICH = ( 0 , 1 ) 2 ( X j ) S Protokoll ( X j ) D X D j = ( 0 , + ) 2 4 u v e ( S + 1 ) ( u 2 + v 2 ) u 2 + v 2 D u D v
und das sieht gut aus. Wenn wir zu Polarkoordinaten wechseln, haben wir:

ICH = ( 0 π / 2 4 Sünde θ cos θ D θ ) ( 0 + ρ 2 e ( S + 1 ) ρ 2 D ρ )
somit:

ICH = π 2 ( S + 1 ) 3 / 2 .

Siehe auch Eulers erste historische Integralformel für die Γ function , und beachten Sie, dass letztere nichts anderes als eine verallgemeinerte Gaußsche Funktion ist ( ersetzen Sie einfach die 2 im Exponenten mit n , und let T = X N ) . Nebenbei bemerkt, das Integral ist im Grunde die Stammfunktion halber Ordnung von ( 0 , 1 ) 2 ( X j ) S   D X   D j in Bezug auf s .
@Lucian: nette Bemerkung.