Integration Absolutbeträge < Mathematica < Mathe-Software < Mathe < Vorhilfe
|
| Status: |
(Frage) überfällig  | | Datum: | 12:13 Mi 05.09.2012 | | Autor: | kaila |
| Aufgabe | P3D = [mm] Piecewise[{{kx^2/(kx*Pi - Pi^2/4), kx > Pi/2}}, [/mm] 1]/kx*
[mm] Piecewise[{{ky^2/(ky*Pi - Pi^2/4), ky > Pi/2}}, [/mm] 1]/ky*
UnitTriangle[theta/kx, phi/ky]
etaeinz =
Abs[Cos[alpha]*Cos[theta] +
Cos[beta]*Cos[phi]*Sin[alpha]*Sin [theta] +
Sin[alpha]*Sin[beta]*Sin[phi]*Sin[theta]]
eta3D = Integrate[
P3D*Sin[theta]*etaeinz, {theta, -Pi/2, Pi/2}, {phi, -Pi/2, Pi/2}] |
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
Hallo!
Ich verzweifle langsam an diesem Problem.... durch die Unstetigkeit meiner Funktionen kann Mathematica nicht integrieren (gibt einfach die Rechnung nochmal aus nach sehr langem rechnen), da sich die Nullstellen mit verschiedenen alpha und beta laufend verschieben. Ohne Abs[] funktionierts natürlich! Das Problem in 2D ohne phi und beta funktionert auch!
Vielen Dank im Voraus!
lg Andre
|
|
| |
|
| Status: |
(Frage) beantwortet  | | Datum: | 12:09 Mi 12.09.2012 | | Autor: | kaila |
| Aufgabe | P3D = [mm] Piecewise[{{kx^2/(kx*Pi - Pi^2/4), kx > Pi/2}}, [/mm] 1]/kx*
[mm] Piecewise[{{ky^2/(ky*Pi - Pi^2/4), ky > Pi/2}}, [/mm] 1]/ky*
UnitTriangle[theta/kx, phi/ky]
etaeinz =
Abs[Cos[alpha]*Cos[theta] +
Cos[beta]*Cos[phi]*Sin[alpha]*Sin [theta] +
Sin[alpha]*Sin[beta]*Sin[phi]*Sin[theta]]
eta3D = Integrate[
P3D*Sin[theta]*etaeinz, {theta, -Pi/2, Pi/2}, {phi, -Pi/2, Pi/2}] |
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
Hallo!
Ich verzweifle langsam an diesem Problem.... durch die Unstetigkeit meiner Funktionen kann Mathematica nicht integrieren (gibt einfach die Rechnung nochmal aus nach sehr langem rechnen), da sich die Nullstellen mit verschiedenen alpha und beta laufend verschieben. Ohne Abs[] funktionierts natürlich!
Versuche es gerade über TrigToExp[], komme aber immer noch auf kein Ergebnis, kann mir niemand helfen???
Vielen Dank im Voraus!
lg Andre
|
|
|
| |
|
Hallo André ,
vermutlich wäre es sinnvoll, wenn du wenigstens
kurz die Herkunft des Problems skizzieren würdest.
Wenn man weiß, worum es geht, kommt vielleicht
auch die Bereitschaft, sich damit zu beschäftigen.
Und was macht z.B. die Funktion "UnitTriangle" ?
LG
Al-Chwarizmi
"UnitTriangle" habe ich gerade nachgeschaut.
Aber zum eigentlichen Problem: mit dem Absolut-
betrag ist das möglicherweise wirklich zu schwierig
für eine allgemeine Lösung in analytischer Form,
auch mit einem sehr starken CAS wie Mathematica.
Du musst dich halt vielleicht mit numerischen
Lösungen (NIntegrate) zufrieden geben !
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 10:32 Do 13.09.2012 | | Autor: | kaila |
vielen Dank für die Antwort!
P3D ist eine Wahrscheinlichkeitsverteilung und der Teil im Absolutbetrag ist eine Projektion auf eine gemeinsame Achse.
Ich werde jetzt versuchen die Rotation auf die Verteilfunktion zu verlegen, somit fallen alpha und beta aus dem Absolutbetrag raus.
danke nochmal,
André
|
|
|
|
