Integral von Heaviside(sin(x))

[wob]

Hallo,

es ist gerade zu warm für mich... Ich stehe gerade vor dem Problem, dass ich folgendes integrieren möchte:

$-\int_a^b \theta(-\sin(c_1(x - c_2))) \ \mathrm dx$

wobei $\theta$ die Heaviside-Funktion ist und $c_n$ Konstanten (aus $\mathbb R^+$) sind, es gilt $b>a$ sowie $a, b \in \mathbb R$, also keine Überraschungen mit den Wertebereichen und nix aus $\mathbb C$.

Wie macht ich das am geschicktesten?

Ohne das Heaviside ists kein Problem. Ich möchte aber nur die negativen Flächen zählen. Meine jetzige Idee ist, zuerst zu ermitteln, wie viele "vollständige" Zyklen ich habe und $b$ dementsprechend verringere und dann den Rest per Fallunterscheidung löse (also z.B. liegt $a$ vor dem 0-Bereich? Dann von $a$ bis zum Anfang des folgenden 0-Bereichs Integrieren. Und dann bei $b$ ebenfalls schauen, wohin ich die Grenzen legen muss.

(Das ganze will ich mit numpy durchführen und für ca. 1 Mio a/b-Paare ausrechnen - aber eigentlich hatte ich mir irgendwie einen klugen Weg erhofft, das Integral "einfach" auszurechnen)

[*john 0]

Die Heaviside-Funktion ist nur für Argumente größer Null als eins definiert sonst als null. D.h. man kann das Integral abschnittsweise als

$-\int_{l_1}^{l_2}1\;\mathrm{d}x = -x\biggr|_{l_1}^{l_2} = -(l_2-l_1)$

berechnen. Das einzige was man nun beachten muss, sind die Integralgrenzen, d.h. für welche Werte aus dem Intervall $[a,b]$ ist das Integral 0 und wann nicht. Die $\sin$-Funktion ist im Intervall $[0,\pi]$ größer 0 und für $]\pi,2\pi[$ kleiner Null. Durch das Minuszeichen in der Heaviside-Funktion vertauschen diese Wertebereiche gerade, d.h. die Heaviside-Funktion ist eins, wenn $\pi < c_1(x-c_2) < 2\pi$. Da $c_1 > 0$

$\pi/c_1 < x-c_2 < 2\pi/c_1$

und

$(\pi/c_1) + c_2 < x < (2\pi/c_1) + c_2$

Das sind nun die Grenzen in denen das obige Integral definiert ist, sonst ist es Null. Du musst nun schauen, wie $a$ und $b$ sich dazu verhalten.