Fraktionalkalkül

Fraktionalkalkül

Die Fraktionale Infinitesimalrechnung bezeichnet die Erweiterung des Ableitungsbegriffs auf nichtganzzahlige Ordnungen. Der Begriff „fraktional“ ist dabei historisch bedingt, die Ableitungen können ganz allgemein von reeller oder sogar komplexer Ordnung sein.

Inhaltsverzeichnis

Vorab: wichtige Funktionen und Integraltransformationen

Diese Funktionen und Transformationen haben meist jeweils eigene Artikel in Wikipedia. Da sie aber bei der Definition der fraktionalen Integrale elementar wichtig sind, sollen sie hier kurz als Definitionen zusammengefasst werden.

(Unvollständige) Gammafunktion

Als Verallgemeinerung der Fakultätsfunktion wird die Gammafunktion wie folgt definiert:

\Gamma(x)=\int_0^\infty t^{x-1} e^{-t} \mathrm{d}t

Für ganzzahlige Argumente ergibt sich Γ(n + 1) = n!. Im Falle der unvollständigen Gammafunktion wird nicht bis unendlich, sondern nur bis zu einem bestimmten Wert y integriert:

\gamma(x,y)=\int_0^y t^{x-1} e^{-t} \mathrm{d}t.

Betafunktion

Die Betafunktion wird definiert als

\Beta(x,y)=\int_0^1 t^{x-1} (1-t)^{y-1}\, \mathrm{d}t

wobei sie sich auch als Produkt von Gammafunktionen darstellen lässt:

\Beta(x,y)=\frac{\Gamma(x) \cdot \Gamma(y)}{\Gamma(x+y)}.

Hypergeometrische Funktion

Als Erweiterung der geometrischen Reihe wird die hypergeometrische Funktion definiert als

{}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;z)=\sum_{k=0}^\infty\prod_{i=1}^p\frac{\Gamma(k+a_i)}{\Gamma(a_i)}\prod_{j=1}^q\frac{\Gamma(b_j)}{\Gamma(k+b_j)}\frac{z^k}{k!}

Sofort einsichtig ist der Spezialfall {}_0F_0(;;z)=\sum_{k=0}^{\infty}\frac{z^k}{k!}=e^z.

Fouriertransformation

Für f\in L^1(\mathbb{R}^d) (z. B. d=1) definiert man

\tilde f(k)=(\mathcal{F}f)(k)=\int_{-\infty}^\infty e^{-2\pi ikx}f(x)\mathrm{d}x

als Fouriertransformation, und

f(x)=(\mathcal{F}^{-1}\tilde f)(x)=\int_{-\infty}^\infty e^{2\pi ikx}\tilde f(k)\mathrm{d}k

als Rücktransformation.

Man beachte, dass es verschiedene Definitionsmöglichkeiten der Fouriertransformation gibt, die sich darin unterscheiden, in welche Transformation man das Minuszeichen in der e-Funktion schreibt, oder wo der Faktor von 2π auftaucht.

Translationsoperator: (\tau_af)(x)=f(x-a) \Rightarrow (\tau_af)^\sim(k)=e^{-2\pi ika}\tilde f(k)\quad\land\quad(e^{2\pi iax}f)^\sim(k)=(\tau_a\tilde f)(k)=\tilde f(k-a.)

Streckoperator: (\delta_\lambda f)(x)=f\left(\frac{x}{\lambda}\right)\quad\forall\lambda>0\qquad\Rightarrow\qquad(\delta_\lambda f)^\sim(k)=\lambda f(\lambda k) bzw. in d Dimensionen: (\delta_\lambda f)^\sim(\vec k)=\lambda^df(\lambda\vec k).

Faltung: (f*g)(x)=\int_{-\infty}^\infty f(x-y)g(y)\mathrm{d}y\quad(=\int_{-\infty}^\infty (\tau_af)(x)g(y)\mathrm{d}y).

Daraus folgt der Faltungssatz: Für f,g\in L^1(\mathbb{R}) ist

(\mathcal{F}(f*g))(k)=(\mathcal{F}f)(k)(\mathcal{F}g)(k).

Die Fouriertransformation macht also aus der Faltung zweier Funktionen die Multiplikation ihrer Fouriertransformierten.

Weiter gilt für f\in L^1(\mathbb R^d)

\left(\frac{\partial}{\partial x_j}f\right)^\sim(\vec k)=2\pi ik_j\tilde f(\vec k).

Laplacetransformation

Sei f\in L_{loc}^1(\mathbb{R}^+), dann ist die Laplacetransformation definiert als

\tilde f(u)=(\mathcal{L}f)(u)=\int_0^\infty f(x) e^{-ux}\ \!\mathrm{d}x.

Die Laplacefaltung wird ähnlich wie die Fourierfaltung definiert und liefert einen ähnlichen Zusammenhang:

(f\star g)(x)=\int_0^\infty f(x-y)g(y)\ \! \mathrm{d}y \Rightarrow \mathcal{L}(f\star g)(u)=(\mathcal Lf)(u)(\mathcal Lg)(u).

Geschichte

Bereits die Mathematiker Gottfried Wilhelm Leibniz und Leonhard Euler beschäftigten sich mit der Verallgemeinerung des Ableitungsbegriffes. Leibniz schildert in einem Brief an Guillaume François Antoine, Marquis de L’Hospital die Ähnlichkeit zwischen Potenzen und der Produktregel von Ableitungen:

(x+y)^n=x^ny^0+nx^{n-1}y+\dots+x^0y^n
\mathrm{d}^n(xy)=\mathrm{d}^nx\mathrm{d}^0y+n\mathrm{d}^{n-1}x\mathrm{d}y+\dots+\mathrm{d}^0x\mathrm{d}^ny

was sich scheinbar einfach auf

\mathrm{d}^\alpha(fg)=\mathrm{d}^\alpha f\mathrm{d}^0g+\frac{\alpha}{1}\mathrm{d}^{\alpha-1}f\mathrm{d}g+\frac{\alpha(\alpha-1)}{1\cdot2}\mathrm{d}^{\alpha-2}f\mathrm{d}^2g+\frac{\alpha(\alpha-1)(\alpha-2)}{1\cdot2\cdot3}\mathrm{d}^{\alpha-3}f\mathrm{d}^3g\dots

verallgemeinern lässt (wobei man im Falle von α-n negativ \mathrm{d}^{\alpha-n}=\int_{}^{n-\alpha} setzt). Jedoch treten bei solch naiver Verwendung von Symboliken Probleme auf. Als Beispiel wähle man eine Funktion f so, dass

\mathrm{d}f=f\mathrm{d}x \land \mathrm{d}^2f=f\mathrm{d}x^2
\Rightarrow\frac{\mathrm{d}^2f}{\mathrm{d}x}=\frac{\mathrm{d}(\frac{\mathrm{d}f}{\mathrm{d}x})}{\mathrm{d}x}=\frac{\mathrm{d}f}{\mathrm{d}x}=f \land \mathrm{d}x=\frac{\mathrm{d}f}{f}

Man beachte das mathematisch an sich nicht korrekte „durchmultiplizieren“ mit dx. Man denkt bei so einer Funktion direkt an die e-Funktion, die jedoch damals noch nicht explizit als solche bekannt war. Wo trifft man nun auf einen Widerspruch, wenn man dαf = fdxα betrachtet? Um das zu sehen setzt man einfach \alpha=\frac{1}{2}:

\frac{\mathrm{d}^\frac{1}{2}f}{\mathrm{d}x^\frac{1}{2}}=\frac{\mathrm{d}^\frac{1}{2}f}{(\frac{\mathrm{d}f}{f})^\frac{1}{2}}=\sqrt{\frac{f}{\mathrm{d}f}}\mathrm{d}^\frac{1}{2}f\neq f.

Somit kann Leibniz einfacher Ansatz nicht die geeignete Lösung des Problems sein.

Eulers Ansatz

Euler betrachtete ganzzahlige Ableitungen von Potenzfunktionen zm. Für diese gilt:

\frac{\mathrm{d}^nz^m}{\mathrm{d}z^n}=\frac{m!}{(m-n)!}z^{m-n}.

Er versuchte nun, diese Beziehung durch Ersetzen der Fakultäts- durch die von ihm gefundene Gammafunktion auf nichtganzzahlige Potenzen zu verallgemeinern:

\frac{\mathrm{d}^\alpha z^\beta}{\mathrm{d}z^\alpha}=\frac{\Gamma(\beta+1)}{\Gamma(\beta-\alpha+1)}z^{\beta-\alpha}.

Auch dieser Weg führt zu Widersprüchen. Wieder betrachte man die e-Funktion eλx, welche n-mal differenziert λneλx ergibt; verallgemeinert also:

\frac{\mathrm{d}^\alpha e^{\lambda x}}{\mathrm{d}x^\alpha}=\lambda^\alpha e^{\lambda x}.

Auf der anderen Seite jedoch ist die e-Funktion nur eine unendliche Potenzreihe, nämlich e^x=\sum_{n=0}^\infty\frac{x^n}{n!}.

Somit hat man zwei Möglichkeiten die α-Ableitung von ex zu berechnen:

  1. Direkt: \frac{\mathrm{d}^\alpha e^x}{\mathrm{d}x^\alpha}=e^x
  2. Indirekt über die Potenzreihendarstellung: \frac{\mathrm{d}^\alpha}{\mathrm{d}x^\alpha}\sum_{n=0}^\infty\frac{x^n}{n!}=\sum_{n=0}^\infty\frac{\mathrm{d}^\alpha}{\mathrm{d}x^\alpha}\frac{x^n}{n!}=\sum_{n=0}^\infty\frac{x^{n-\alpha}}{\Gamma(n-\alpha+1)}\neq\sum_{n=0}^\infty\frac{x^n}{\Gamma(n+1)}=e^x

Diesen Widerspruch kann man mit dem Beispiel α=-1 erklären, wenn wieder negative Exponenten der Differentialoperatoren als Integrale aufgefasst werden:

\frac{\mathrm{d}^{-1}}{\mathrm{d}x^{-1}}e^x=\int_{-\infty}^xe^t\mathrm{d}t=e^x
\frac{\mathrm{d}^{-1}}{\mathrm{d}x^{-1}}x^\beta=\int_0^xt^\beta\mathrm{d}t=\frac{x^{\beta+1}}{\beta+1}

Die unterschiedlichen unteren Grenzen verdeutlichen, dass man mit diesem Ansatz „wissen muss“, von wo bis wo man zu integrieren hat, um die korrekte Stammfunktion zu finden. Somit ist auch Eulers Ansatz, obwohl von der Idee und Ausführung her besser, nicht geeignet, den Differentialoperator korrekt auf reelle Potenzen zu verallgemeinern.

Definition fraktionaler Integraloperatoren I^\alpha_{a+} und I^\alpha_{a-}

Iterative und fraktionale Integrale

Eine Möglichkeit, fraktionale Integrale widerspruchsfrei zu definieren ergibt sich aus der Formel, die die doppelte Integration über zwei Variablen mit gleicher unterer Grenze in ein einziges Integral überführt. Diese Formel kann auf beliebig viele Integrale erweitert werden.

Ohne Beweis: \int_a^x\int_a^yf(z)\mathrm{d}z\mathrm{d}y=\int_a^x(x-y)f(y)\mathrm{d}y.

Führt man nun noch den Integraloperator Ia+ wie folgt ein,

(I_{a+}f)(x)=\int_a^xf(y)\mathrm{d}y=[F(y)]_a^x=F(x)-F(a)

wobei F(x) die Stammfunktion von f(x) ist, dann können beliebig hohe Potenzen dieses Operators dank obiger Formel von Mehrfachintegralen auf ein einziges Integral zurückgeführt werden:

(I_{a+}^nf)(x)=\int_a^x\int_a^{y_1}\int_a^{y_2}\dots\int_a^{y_{n-1}}f(y_n)\mathrm{d}y_n\mathrm{d}y_{n-1}\dots\mathrm{d}y_1=\frac{1}{(n-1)!}\int_a^x(x-y)^{n-1}f(y)\mathrm{d}y.

Im Gegensatz zu den Formeln zu Beginn kann man diesen Integraloperator relativ problemlos von ganzen Zahlen n auf reelle (bzw. komplexe) Zahlen α verallgemeinern, indem man n durch α und die Fakultät durch die Gammafunktion ersetzt und fordert, dass -\infty\le\alpha<x<\infty \land f\in L_{loc}^1(a,\infty):

(I_{a+}^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_a^x(x-y)^{\alpha-1}f(y)\mathrm{d}y

Dies wird rechtsseitiges fraktionales Riemann-Liouville-Integral genannt. Analog dazu kann durch

(I_{b-}^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_x^b(y-x)^{\alpha-1}f(y)\mathrm{d}y

das linksseitige Äquivalent definiert werden.

Rückführung fraktionaler Integrale auf Faltungen

Definiert man die Distribution x_{a+}^\alpha=\begin{cases}x^\alpha & \forall x>a \\ 0 & \forall x\leq a\end{cases}, kann das fraktionale Integral auf eine Laplacefaltung zurückgeführt werden:

\mathcal L(I_{a+}^\alpha f)(u)=(\mathcal Lf)(u)(\mathcal L\frac{x_{a+}^\alpha}{\Gamma(\alpha)})(u),

da \mathcal L(x_{a+}^\alpha)(u)=\Gamma(\alpha-1)u^{-\alpha-1}

Fraktionale Weylintegrale

Lässt man in den Gleichungen oben a bzw. b betragsmäßig gegen unendlich gehen erhält man die sogenannten Weylintegrale und die entsprechenden partiellen Integraloperatoren

(I_+^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_{-\infty}^x(x-y)^{\alpha-1}f(y)\mathrm{d}y
(I_-^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_x^\infty(y-x)^{\alpha-1}f(y)\mathrm{d}y

für -\infty<x<\infty und der Definitionsmenge \{f\mathcal|I_\pm^\alpha f<\infty\} (z. B. f\in L^p(\mathbb{R}), 1\leq p\leq\frac{1}{\alpha}).

Fraktionale Weylintegrale und Faltungen

Auch fraktionale Weylintegrale lassen sich auf Faltungen zurückführen. Allerdings sind dies Fourierfaltungen, da Weylintegrale eine unendliche untere beziehungsweise obere Grenze haben.

(I_+^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_{-\infty}^x(x-y)^{\alpha-1}f(y)\mathrm{d}y

was durch xy = t überführt werden kann in

\frac{1}{\Gamma(\alpha)}\int_0^\infty t^{\alpha-1}f(x-t)\mathrm{d}t=\frac{1}{\Gamma(\alpha)}\int_{-\infty}^\infty t_+^{\alpha-1}f(y)\mathrm{d}y mit t_+^\alpha=\begin{cases}t^\alpha & \forall t>0 \\0 & \forall t\leq0 \end{cases}
\Rightarrow (I_+^\alpha f)(x)=(f*\frac{x_+^{\alpha-1}}{\Gamma(\alpha)})(x)

Daher ergibt die Fouriertransformation für 0 < α < 1

(I_+^\alpha f)^\sim(k)=(\mathcal FI_+^\alpha f)(k)=\left(\mathcal F\frac{x_+^{\alpha-1}}{\Gamma(\alpha)}\right)(k)(\mathcal Ff)(k)=(2\pi ik)^{-\alpha}\tilde f(k)

Man sieht also, dass der fraktionale Riemann-Liouville-Integraloperator durch die Laplacefaltung, der fraktionale Weyl-Integraloperator entsprechend durch die Fourierfaltung diagonalisiert wird.

Beispiele

\underline{f(x)=x^\beta}:

(I_{0+}^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_0^x(x-y)^{\alpha-1}y^\beta\mathrm{d}y=\frac{x^{\alpha-1}}{\Gamma(\alpha)}\int_0^x\left(1-\frac{y}{x}\right)^{\alpha-1}y^\beta\mathrm{d}y

Substituiere z(y)=y/x \rightarrow dy=xdz

(I_{0+}^\alpha f)(x)=\frac{x^{\alpha+\beta}}{\Gamma(\alpha)}\int_0^1(1-z)^{\alpha-1}z^\beta\mathrm{d}y=\frac{\Beta(\beta+1,\alpha)x^{\alpha-\beta}}{\Gamma(\alpha)}=\frac{\Gamma(\beta+1)}{\Gamma(\alpha+\beta+1)}x^{\alpha+\beta}

Im Spezialfall α=1 wird daraus

(I_{0+}^1f)(x)=\frac{\Gamma(\beta+1)}{\Gamma(\beta+2)}x^{1+\beta}=\frac{\Gamma(\beta+1)}{(\beta+1)\Gamma(\beta+1)}x^{\beta+1}=\frac{1}{\beta+1}x^{\beta+1}
\underline{f(x)=e^{ax}} (mit der Substitution z=a(x-y):
(I_+^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_{-\infty}^x(x-y)^{\alpha-1}e^{ay}\mathrm{d}y=-\frac{1}{\Gamma(\alpha)}\int_\infty^0(\frac{z}{a})^{\alpha-1}e^{ax-z}\frac{1}{a}\mathrm{d}z=\frac{a^{-\alpha}e^{ax}}{\Gamma(\alpha)}\int_0^\infty z^{\alpha-1}e^{-z}\mathrm{d}z=a^{-\alpha}e^{ax}

Man erkennt also, dass man auch bei diesem Integraloperator, ähnlich wie bei Eulers Ansatz, „wissen muss“, von wo bis wo man zu integrieren hat, um die eigentliche Stammfunktion einer Funktion zu erhalten, jedoch steckt dies in der Operatordefinition explizit drin. Somit muss die untere Grenze so gewählt werden, dass in F(x) − F(a) (siehe ganz oben, Def. von Ia+) das F(a) verschwindet und man F(x) (bzw. das fraktionale Äquivalent dazu) erhält. So haben wir in diesem zweiten Beispiel eax von -∞ bis x integriert, wohlwissend, dass eax für a -> -∞ gegen 0 geht. Daher integrieren wir diese Funktion einfach noch einmal, diesmal jedoch mit unterer Grenze 0 (und der Substitution z=x-y):

(I_{0+}^\alpha e^{ax})(x)=\frac{1}{\Gamma(\alpha)}\int_0^x(x-y)^{\alpha-1}e^{ay}\mathrm{d}y=\frac{1}{\Gamma(\alpha)}\int_0^xz^{\alpha-1}e^{a(x-z)}\mathrm{d}z=\frac{e^{ax}}{\Gamma(\alpha)}\int_0^xz^{\alpha-1}e^{-z}\mathrm{d}z

Substituiert man hier nun noch az mit t, dann ergibt sich:

=\frac{e^{ax}}{\Gamma(\alpha)a^{\alpha-1}}\int_0^{ax}t^{\alpha-1}e^{-t}\frac{1}{a}\mathrm{d}t=\frac{\gamma(\alpha,ax)}{\Gamma(\alpha)}a^{-\alpha}e^{ax}

\underline{f(x)=(x+c)^\beta}:

(I_{0+}^\alpha f)(x)=\frac{1}{\Gamma(\alpha)}\int_0^x(x-y)^{\alpha-1}(y+c)^\beta\mathrm{d}y=\frac{c^\beta x^{\alpha-1}}{\Gamma(\alpha)}\int_0^x(1-\frac{y}{x})^{\alpha-1}(1+\frac{y}{c})^\beta\mathrm{d}y

Substitution von z mit y/x führt auf

(I_{0+}^\alpha f)(x)=\frac{c^\beta x^\alpha}{\Gamma(\alpha)}\int_0^1(1-z)^{\alpha-1}(1+\frac{xz}{c})^\beta\mathrm{d}z

Man vergleiche dies mit \int_0^1y^{a-1}(1-y)^{c-a-1}(1-zy)^{-b}\mathrm{d}y=\frac{\Gamma(c-a)\Gamma(a)}{\Gamma(c)}{}_2F_1(a,b;c;z). Man sieht, dass man einfach a=1, b=-β, c=α+1 und z=-x/c setzen muss, um das obige Integral zu erhalten. Also ist

(I_{0+}^\alpha f)(x)=\frac{x^\alpha c^\beta}{\Gamma(\alpha+1)}\frac{\Gamma(\alpha+1-1)\Gamma(1)}{\Gamma(\alpha+1)}{}_2F_1(1,-\beta;\alpha+1;-\frac{x}{c})=\frac{c^\beta}{\Gamma(\alpha+1)}{}_2F_1(1,-\beta;\alpha+1;-\frac{x}{c})x^\alpha

Integration von hypergeometrischen Funktionen

Da sich mit hypergeometrischen Funktionen sehr viele andere Funktionen darstellen lassen bietet es sich an, hier eine Formel zu deren Integration darzustellen.

f(x)={}_pF_q(a_1,\dots,a_p;b_1,\dots,b_q;x)
(I_{0+}^\alpha f)(x)=\frac{x^\alpha}{\Gamma(\alpha+1)}{}_{p+1}F_{q+1}(a_1,\dots,a_p,1;b_1,\dots,b_q,1+\alpha;x)

Dualität der + und - Operatoren

Ganz allgemein gilt

\int_a^bf(x)(I_{a+}^\alpha g)(x)\mathrm{d}x=\int_a^b(I_{b-}^\alpha f)(x)g(x)\mathrm{d}x

die beiden Riemann-Liouville- und die beiden Weylintegrale sind also jeweils dual zueinander. So kann man in Integralen das fraktionale Integral von einer Funktion auf die leichter zu integrierende verschieben.

Definition allgemeiner frakt. Integraloperatoren Iα und {I^\star}^\alpha

Erste Verallgemeinerung der fraktionalen Integration

Durch den Ansatz

(I^\alpha f)(x)=\frac{1}{C(\alpha)}\int_{-\infty}^\infty \vert x-y\vert^{\alpha-1}f(y)\mathrm dy

soll versucht werden, einen allgemeineren Integraloperator zu definieren. Die Betragsstriche statt einfach runde Klammern deuten bereits an, dass für diesen eine Art Kugelsymmetrie vorausgesetzt wird. C(α) soll so bestimmt werden, dass die Additivität der Ordnung (IαIβ=Iα+β) weiterhin gilt. Man kann schon vermuten, dass dieser Operator einfach eine Linearkombination aus den bereits bekannten Weylintegraloperatoren ist, was man auch beweisen kann:

(I^\alpha f)(x)=\frac{1}{C(\alpha)}\left[\int_{-\infty}^x(x-y)^{\alpha-1}f(y)\mathrm dy+\int_x^\infty(y-x)^{\alpha-1}f(y)\mathrm dy\right]=\frac{\Gamma(\alpha)}{C(\alpha)}\left[(I_+^\alpha f)(x)+(I_-^\alpha f)(x)\right]

Also ist

I^\alpha=\frac{\Gamma(\alpha)}{C(\alpha)}(I_+^\alpha+I_-^\alpha)

was man auch auf höhere Dimensionen verallgemeinern kann:

(I^\alpha f)(\vec x)=\frac{1}{C(d,\alpha)}\int_{\mathbb R^d}\vert \vec x-\vec y\vert^{\alpha-d}f(\vec y)\mathrm d\vec y

Nun ist die Frage, wie C(d,α) bestimmt werden kann. Wenn man die Wahl so treffen möchte, dass (I^\alpha f)^\sim\left(\vec k\right)=\left|2\pi \vec k\right|^{-\alpha}\tilde f\left(\vec k\right) gilt, dann ergibt sich nach eingehendem Studium der Fouriertransformation für C(d,α):

C(d,\alpha)=\frac{2^\alpha\pi^\frac{d}{2}\Gamma(\frac{\alpha}{2})}{\Gamma(\frac{d-\alpha}{2})}.

Unter Ausnutzung der Formeln \Gamma(2z)=\frac{2^{2z-1}}{\sqrt{\pi}}\Gamma(z)\Gamma\left(z+\frac{1}{2}\right) und \Gamma\left(\frac{1-z}{2}\right)\Gamma\left(\frac{1+z}{2}\right)=\frac{\pi}{\cos(z\frac{\pi}{2})} ergibt sich somit im eindimensionalen Fall:

(I^\alpha f)(x)=\frac{1}{2\Gamma(\alpha)\cos(\alpha\frac{\pi}{2})}\int_{-\infty}^\infty\vert x-y\vert^{\alpha-1}f(y)\mathrm dy.

Dies wird fraktionales Riesz-Feller-Integral genannt.

Weitere Verallgemeinerung der fraktionalen Integration

Die Formel I^\alpha=\frac{\Gamma(\alpha)}{C(\alpha)}(I_+^\alpha+I_-^\alpha) lässt den Schluss zu, dass man weitere derart allgemeinen Integraloperatoren durch

I^\alpha=\frac{\Gamma(\alpha)}{C(\alpha)}(c_+I_+^\alpha+c_-I_-^\alpha)

definieren kann, was das Riesz-Feller-Integral zum Spezialfall c+=c-=1 macht. Z. B. ergibt sich für c+=1 und c-=-1

({I^\star}^\alpha f)(x)=\frac{\Gamma(\alpha)}{C(\alpha)}((I_+^\alpha f)(x)-(I_-^\alpha f)(x))=\frac{1}{2\Gamma(\alpha)\sin(\alpha\frac{\pi}{2})}\int_{-\infty}^\infty\operatorname{sign}(x-y)\vert x-y\vert^{\alpha-1}f(y)\mathrm dy

Diese beiden Operatoren sind verknüpft durch die Hilberttransformation:

{I^\star}^\alpha=\mathcal HI^\alpha \land I^\alpha=\mathcal H{I^\star}^\alpha

Feller hat für Integrale der Form

(I_\delta^\alpha f)(x)=\frac{1}{\Gamma(\alpha)\sin(\alpha\pi)}\int_{-\infty}^\infty\sin\left(\alpha\left(\frac{\pi}{2}+\delta\operatorname{sign}(y-x)\right)\right)\vert x-y\vert^{\alpha-1}f(y)\mathrm dy

bewiesen, dass die Additivität der Ordnung gilt. Diese Integrale lassen sich ebenfalls als Linearkombination der obigen Form darstellen, dazu muss man nur c_\pm=\frac{\sin(\alpha\frac{\pi}{2}\mp\delta\alpha)}{\sin(\alpha\pi)} wählen.

Beispiele fraktionaler Integralgleichungen in der Physik

Tautochronenproblem

Problem in zwei Dimensionen: ein Massepunkt fällt unter Einfluss der Schwerkraft entlang einer festen, aber unbekannten Bahn y=h(x) von der Höhe y0 auf Höhe y1; die Zeit die er dafür benötig wird angegeben mit T_{y_0}(y_1) = Zeit des Falles von festem y0 zu variablem y1. Die Frage ist nun: Lässt sich aus Kenntnis des Fallzeiten T_{y_0}(y)\;\forall\;y\in[y_1,y_0] allein bereits h(x) bestimmen?

Wir setzen v(y) gleich dem Betrag der Momentangeschwindigkeit, dann ergibt sich für die Zeitdauer des Falles von P0 auf P1: T_{y_0}(y_1)=\int_{P_0(x_0,y_0)}^{P_1(x_1,y_1)}\frac{\mathrm ds(y)}{v(y)} mit s(y) gleich dem zurückgelegten Wert als Funktion der Höhe. Wenn nun y=h(x) invertierbar ist, dann ist x=h-1(y)=Φ(y) und das Bogenlängendifferential \mathrm ds=\sqrt{1+\Phi'(y)^2}\mathrm dy mit der Bogenlänge s=\int_{y_0}^{y_1}\sqrt{1+\Phi'^2}\mathrm dy

Aus dem Energiesatz mg(y_0-y)=\frac{m}{2}v(y)^2 folgt v(y)=\sqrt{2g(y_0-y)}. Einsetzen in die Gleichung für T ergibt

T_{y_0}(y_1)=\int_{y_0}^{y_1}\sqrt{\frac{1+\Phi'(y)^2}{2g(y_0-y)}}\mathrm dy

Definiert man nun f(y)=\sqrt{\frac{\pi}{2g}}\sqrt{1+\Phi'(y)^2} (und bedenkt, dass \sqrt{\pi}=\Gamma\left(\frac{1}{2}\right) ist), dann ergibt sich

T_{y_0}(y_1)=\frac{1}{\Gamma(\frac{1}{2})}\int_{y_0}^{y_1}\frac{f(y)}{(y-y_0)^\frac{1}{2}}\mathrm dy=\left(I_{y_0+}^\frac{1}{2}f\right)(y_1)

Literatur

  • R. Herrmann: Fraktionale Infinitesimalrechnung. Eine Einführung für Physiker. BoD, Norderstedt 2008, ISBN 978-3837059588

Wikimedia Foundation.

Игры ⚽ Поможем написать курсовую

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”