Ich hätte gerne Hilfe bei der Annäherung des Oberflächenintegrals
∫Sn − 1≥ 01N^⋅ pDS
Wo
N^
ist die Einheit normal zur Kugel an dem gegebenen Punkt,
p ∈RN
ein fester Vektor mit allen positiven Komponenten ist, und
Sn − 1≥ 0
ist der Anteil der Kugel in
RN
mit allen positiven Komponenten.
Das Problem ist auf folgende Weise aufgetreten. Angenommen, man hat ein zeitabhängiges Budgetv ( t )
auf einem Markt vonN
Artikel mit jeweils einem positiven Preis, der sich über die Zeit ändert, gegeben durch den Vektorp ( t ) = ⟨P1, … ,PN⟩
. 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 ∈NN|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 ∈RN≥ 0|x ⋅ p ( t ) ≤ v ( t ) } |
Angenommen, einige
X0
erfüllt
X0⋅ p ( t ) ≤ v ( t )
, dann alle Vektoren in Richtung von berücksichtigen
X0
die auch die Ungleichung erfüllen, erhalten wir eine Linie der Länge
v ( t ) /X0⋅ 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 )∫Sn − 1≥ 01N^⋅ pDS
Das Beste, was ich tun konnte, um dieses Integral anzunähern, war, Folgendes zu berücksichtigen:
∫Sn − 1≥ 01N^⋅ pDS=1∥ p ∥∫Sn − 1≥ 01cosθDS
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= min { ∠ (e1, p ) , … , ∠ (eN, p ) }
klar können wir nach außen fegen
θ0
ohne die Grenzen zu verlassen
Sn − 1≥ 0
, 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 ∥∫Sn − 1≥ 01cosθDS∼∫π/ 20# {eich|∠ (eich, p ) ≤ θ }N⋅1cosθ⋅2πn / 2Γ ( n / 2 )θn − 1Dθ
wo die Obergrenze von
π/ 2
war relativ willkürlich, einfach weil
# {eich|∠ (eich, p ) ≤ θ }
ist garantiert null für
θ
beliebig größer. Da wir das wissen
∠ (eich, p ) = arccos(Pich/ ∥p∥)
wir können dieses Integral zu vereinfachen
=∫π/ 20# {eich|∠ (eich, p ) ≤ θ }N⋅1cosθ⋅2πn / 2Γ ( n / 2 )θn − 1Dθ2πn / 2n Γ ( n / 2 )∑ich = 1N∫arccos(Pich/ ∥p∥)0θn − 1Sekθ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 / 2n Γ ( n / 2 )∑ich = 1N∫arccos(Pich/ ∥p∥)0θn − 1Sekθ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.
zw.
Will Fischer
jorik
jorik
Will Fischer