Annäherung an Integral über Sphäre

Ich hätte gerne Hilfe bei der Annäherung des Oberflächenintegrals

S 0 N 1 1 N ^ P D S
Wo N ^ ist die Einheit normal zur Kugel an dem gegebenen Punkt, P R N ein fester Vektor mit allen positiven Komponenten ist, und S 0 N 1 ist der Anteil der Kugel in R N mit allen positiven Komponenten.

Das Problem ist auf folgende Weise aufgetreten. Angenommen, man hat ein zeitabhängiges Budget v ( T ) auf einem Markt von N Artikel mit jeweils einem positiven Preis, der sich über die Zeit ändert, gegeben durch den Vektor P ( T ) = P 1 , , P N . Eine Möglichkeit, eine "Kaufkraft" zu berücksichtigen, wäre die Anzahl der Möglichkeiten, unterschiedliche Mengen dieser Artikel zu kaufen, d.h.

PP ( T ) = # { X N N | X P ( T ) v ( T ) }
Um dies anzunähern, wollen wir es "undeskretisieren". Unter der Annahme, dass beide gleichmäßig verteilt sind, die Anzahl solcher ganzzahliger Komponenten X sollte proportional zum Volumen der Menge sein, wenn sie auf realwertige Komponenten erweitert wird, d.h.
PP ( T ) | { X R 0 N | X P ( T ) v ( T ) } |
Angenommen, einige X 0 erfüllt X 0 P ( T ) v ( T ) , dann alle Vektoren in Richtung von berücksichtigen X 0 die auch die Ungleichung erfüllen, erhalten wir eine Linie der Länge v ( T ) / X 0 P ( T ) . Somit kann dieses Volumen gut angenähert werden, indem die Längen dieser Linien über alle möglichen Richtungen integriert werden
PP ( T ) v ( T ) S 0 N 1 1 N ^ P D S
Das Beste, was ich tun konnte, um dieses Integral anzunähern, war, Folgendes zu berücksichtigen:
S 0 N 1 1 N ^ P D S = 1 P S 0 N 1 1 cos θ D S
Wo θ ist der Winkel dazwischen P und die aktuelle normal N ^ . Um dies anzunähern, habe ich mich für die Integration in Bezug auf entschieden θ auf der Kugel um den Winkel aus nach außen fegen P beginnt um P Wo θ = 0 . Wenn wir lassen θ 0 = Mindest { ( e 1 , P ) , , ( e N , P ) } klar können wir nach außen fegen θ 0 ohne die Grenzen zu verlassen S 0 N 1 , bewegt sich jedoch an diesem Winkel vorbei P wird diese Grenzen verlassen. Ich gehe jedoch davon aus (möglicherweise fälschlicherweise), dass der Beitrag bei einem bestimmten Winkel ungefähr proportional zum Anteil der Achsen innerhalb dieses Winkels ist P . Das ist
1 P S 0 N 1 1 cos θ D S 0 π / 2 # { e ich | ( e ich , P ) θ } N 1 cos θ 2 π N / 2 Γ ( N / 2 ) θ N 1 D θ
wo die Obergrenze von π / 2 war relativ willkürlich, einfach weil # { e ich | ( e ich , P ) θ } ist garantiert null für θ beliebig größer. Da wir das wissen ( e ich , P ) = arccos ( P ich / P ) wir können dieses Integral zu vereinfachen
0 π / 2 # { e ich | ( e ich , P ) θ } N 1 cos θ 2 π N / 2 Γ ( N / 2 ) θ N 1 D θ = 2 π N / 2 N Γ ( N / 2 ) ich = 1 N 0 arccos ( P ich / P ) θ N 1 Sek θ D θ
Alles in allem habe ich also das Gefühl (obwohl mir ein vollständig formaler Beweis fehlt).
PP ( T ) v ( T ) P ( T ) 2 π N / 2 N Γ ( N / 2 ) ich = 1 N 0 arccos ( P ich / P ) θ N 1 Sek θ D θ
aber das ist das Beste, was ich tun kann. Ich kann dieses endgültige Integral anscheinend auch nicht vereinfachen oder annähern. Ich könnte evtl. eine Serienerweiterung von gebrauchen Sek ( X ) und termweise integrieren.

Vermuten N = 2 Und P = ( 1 , 1 ) . Dann integrieren Sie sich
0 π / 2 1 cos T Sünde T D T .
Das hat ein Problem T = π / 4. Welchen Sinn macht das Integral in diesem Fall?
Mein Fehler, P sollte alle positiven Komponenten in diesem Problem haben, da es Preise darstellt. Wenn überhaupt 0 Dann PP wäre per Definition unendlich und hätte Probleme bei allen Näherungen.
Von v ( T ) / X 0 P ( T ) was meinen Sie v ( T ) / ( X 0 P ( T ) ) ?
Zwei Dinge: a) Sie bekommen nicht die Lautstärke von { X R 0 N | X P ( T ) v ( T ) } durch Integrieren der Längen dieser Linien. Sie berücksichtigen nicht, dass die Linie umso mehr zur Lautstärke beiträgt, je weiter sie vom Ursprung entfernt ist. ZB in drei Dimensionen sollten Sie die Volumina von Kegeln entlang dieser Linien integrieren. Der richtige Integrand wäre proportional zu R N , nicht zu R . b) Warum integrieren Sie dies in Kugelkoordinaten statt in kartesischen Koordinaten? Das Integral ist einfach in kartesischen Koordinaten.
@joriki Ah du hast recht. Wie einfach ist es in kartesischen Koordinaten? Ist so etwas wie
0 v / P 1 0 ( v P 1 X 1 ) / P 2 0 ( v P 1 X 1 P N 1 X N 1 ) / P N D X N D X 1
was hattest du vor?

Antworten (1)

Sie wollen die Lautstärke des Sets { X R 0 N | X P v } . Dies ist das Volumen von R 0 N definiert durch die Koordinatenebenen und die Ebene X P = v . Dies ist ein Simplex mit Koordinatenachsenabschnitten X ich = v / P ich . So hat es das Volumen 1 / N ! der Einheit simplex skaliert um ich v / P ich .

Dies ist ein unglaublich einfacher Ansatz. Danke!