Erneuerungstheorie

aus Wikipedia, der freien Enzyklopädie
Wechseln zu: Navigation, Suche

Die Erneuerungstheorie (engl. renewal theory) ist ein Spezialgebiet der Wahrscheinlichkeitstheorie und befasst sich mit Prozessen, die sich nach jedem Erreichen des Ausgangszustandes wieder so verhalten wie beim Start des Experiments.

Motivation[Bearbeiten]

Ein motivierendes Beispiel ist die regelmäßige Erneuerung der für den Betrieb einer Lampe verwendeten Glühbirne, die nach jedem Ausfall zu ersetzen ist. Die Lebensdauer einer Glühbirne wird durch eine Zufallsvariable beschrieben, deren Verteilung für alle in Frage kommenden Glühbirnen gleich ist, also eine bekannte, charakteristische Eigenschaft darstellt. Ferner wird angenommen, dass diese Lebensdauern voneinander unabhängig sind. Es ist nun von Interesse, wie häufig die Glühbirne durchschnittlich auszuwechseln ist, das heißt man fragt, wieviele Erneuerungen bis zu einer vorgegebenen Betriebszeit vorzunehmen sind.

Sehr ähnliche Aufgabenstellungen erhält man natürlich für allgemeinere Wartungsarbeiten oder für Bedienzeiten von Kunden, die nach einer vorgegebenen Verteilung an einer Abfertigungsstelle erscheinen und dort eine Warteschlange bilden. Hier liefert die Erneuerungstheorie Hinweise für optimale Wartungsintervalle oder optimale Personalvorhaltung an Servicestellen. Versteht man das Eintreten eines Wartungsfalls als Schadensfall, so wird sofort verständlich, dass die Erneuerungstheorie auch in der Versicherungsmathematik von Bedeutung ist.

Definition[Bearbeiten]

Bei jeder Erneuerungszeit S_n springt N um mindestens eine Einheit.

Ein Erneuerungsprozess wird durch eine Folge (T_n)_{n\in \N} von unabhängigen, identisch verteilten, nicht-negativen Zufallsgrößen mit P(T_n=0) < 1 gegeben, wobei P die Wahrscheinlichkeit auf dem zugrunde liegenden Wahrscheinlichkeitsraum sei.

S_n := T_1+\ldots + T_n

heißt n-te Erneuerungszeit, wobei zusätzlich S_0 die konstante Funktion 0 sei. Die Folge (S_n)_{n\in \N} heißt Erneuerungsfolge, ein Intervall (S_{n-1},S_n] heißt Erneuerungszyklus und hat definitionsgemäß die Länge T_n, die man daher auch Zykluszeit nennt. Schließlich setzt man für t\in [0,\infty)

N(t) := |\{n\in \N|\, S_n\le t\}|,

die Anzahl aller n, für die die n-te Erneuerungszeit den Zeitpunkt t noch nicht überschritten hat. Der so definierte stochastische Prozess (N(t))_{t\in [0,\infty)} heißt der Erneuerungsprozess.

Bemerkungen[Bearbeiten]

Interpretation der Definition[Bearbeiten]

Diese Definitionen werden anhand obigen Glühbirnenbeispiels sofort verständlich. T_n modelliert die Betriebsdauer der n-ten Glühbirne, S_n ist die durch n Glühbirnen hintereinander erbrachte Gesamtleuchtdauer, N(t) schließlich ist die Anzahl der bis zum Zeitpunkt t erforderlichen Glühbirnenwechsel. Die Bedingung P(T_n=0) < 1 stellt sicher, dass eine neu eingesetzte Glühbirne nicht mit Sicherheit sofort wieder ausfällt; nur dann ist die zeitliche Betrachtung regelmäßiger Erneuerungen sinnvoll. Ähnliche Interpretationen für Wartungsarbeiten, Serviceleistungen oder Schadensfälle sind naheliegend.

Verzögerter Erneuerungsprozess[Bearbeiten]

Eine häufig verwendete Variante ist der sogenannte verzögerte Erneuerungsprozess, bei dem die Verteilung von T_1 von der gemeinsamen Verteilung der übrigen T_n, n\ge 2, abweichen darf. Dies wird erforderlich, wenn man die Ausgangssituation nicht kennt und daher über T_1 eine andere Annahme treffen muss, oder wenn, etwa im Falle von Wartungsarbeiten, die Originalbauteile andere sind als die regelmäßig auszutauschenden Ersatzteile. Der eigentliche Erneuerungsprozess beginnt also erst nach T_1, was die Bezeichnung als verzögerten Erneuerungsprozess erklärt.

Auszahlungsprozess[Bearbeiten]

In der Regel sind mit dem Eintreten einer jeden Erneuerungszeit Auszahlungen, die im Kostenfall auch negativ sein können, verbunden. Daher betrachtetet man zu den in obiger Definition gegebenen Daten noch eine Folge (R_n)_{n\in \N} unabhängiger und identisch verteilter Zufallsvariablen, die für die Auszahlungen zum n-ten Erneuerungszeitpunkt S_n stehen. Die Gesamtauszahlung bis zum Zeitpunkt t ist dann

R(t) := \sum_{n=1}^{N(t)}R_n.

Der stochastische Prozess R(t)_{t\in [0,\infty)} heißt der zum Erneuerungsprozess gehörige Auszahlungsprozess.

In vielen Anwendungen geht es darum, die mit diesem Auszahlungsprozess verknüpften Daten zu optimieren. In einer hier nicht näher betrachteten Variante kann R_n durch eine sich während des n-ten Erneuerungszyklus entwickelnde Funktion ersetzt werden, so dass obiges R_n die während des Zyklus kumulierte Auszahlung ist. Dadurch können die zu den Erneuerungszeiten auftretenden Sprünge vermieden werden.

Grundlagen[Bearbeiten]

Es liege ein wie oben beschriebener Erneuerungsprozess vor, F sei die Verteilungsfunktion der Zykluszeiten T_n. Die durchschnittliche Zykluszeit \mu = E(T_n) ist positiv, da sonst T_n fast sicher 0 wäre, was der Voraussetzung P(T_n=0) < 1 widerspräche. Für den Erwartungswert des Erneuerungsprozesses gilt

m(t) := E(N(t)) = \sum_{n=1}^{\infty}P(S_n \le t) = \sum_{n=1}^{\infty}F^{*n},

wobei F^{*n} die n-fache Faltung mit sich sei. Man nennt m in naheliegender Weise die Mittelwertsfunktion des Erneuerungsprozesses. Unter Verwendung der Bedingung P(T_n=0) < 1 kann man zeigen, dass m(t) endlich und daher N(t) fast überall endlich ist. Daraus ergibt sich weiter das für einen sinnvollen Erneuerungsprozess erwartete Grenzwertverhalten

S_n\to \infty fast sicher für n\to\infty
N(t)\to \infty fast sicher für t\to\infty.

Über das Wachstum von N(t) kann man eine viel genauere Aussage treffen:

 \frac{N(t)}{t} \to \frac{1}{\mu} fast sicher für t\to\infty.

Diese Aussage gilt auch unter dem Erwartungswert, das heißt

 \frac{m(t)}{t} \to \frac{1}{\mu} fast sicher für t\to\infty,

was auch als einfaches Erneuerungstheorem bekannt ist. Dies bestätigt die Intuition, dass die langfristig erwartete Anzahl der Erneuerungen pro Zeit mit dem Kehrwert der zu erwartenden Dauer zwischen zwei Erneuerungen übereinstimmt. Entsprechende Resultate hat man für den zugehörigen Auszahlungsprozess:

 \frac{R(t)}{t} \to \frac{E(R_1)}{\mu} fast sicher für t\to\infty,

das heißt im langfristigen Mittel ist die Auszahlung pro Zeit gleich der mittleren Auszahlung eines Erneuerungszyklus geteilt durch die mittlere Zykluslänge.

Der Poissonprozess als Erneuerungsprozess[Bearbeiten]

Der einfachste Fall liegt vor, wenn die T_n exponentialverteilt mit einem Parameter \lambda sind. Dann ist (N(t))_{t\in [0,\infty)} ein Poissonprozess zum Parameter \lambda, d.h. N(t) ist Poisson-verteilt zum Parameter \lambda t . In diesem Fall ist daher m(t) = E(N(t)) = \lambda t und das Erneuerungstheorem wird trivial, denn für die mittlere Zykluszeit gilt \textstyle \mu = \frac{1}{\lambda}, da dies der Erwartungswert einer Exponentialverteilung ist.

Anwendung[Bearbeiten]

Zur Veranschaulichung der oben eingeführten Begriffe betrachten wir folgende Strategie für das Auswechseln von Glühbirnen, deren zufällige Lebensdauern durch unabhängige und identisch verteilte Zufallsvariable X_n mit Verteilungsfunktion F gegeben seien. Wir wechseln spätestens nach einer noch zu bestimmenden Zeit T, was zu Kosten in Höhe von a führt, und nur dann früher, wenn die Glühbirne tatsächlich ausfällt, was neben a zusätzliche Kosten in Höhe von  b verursacht. Der Erneuerungszyklus hat daher die zufällige Länge

T_n = \min(X_n,T).

Zum n-ten Erneuerungszeitpunkt hat man dann Kosten R_n mit Erwartungswert

E(R_n) = E(R_1) = aP(X_n>T)+(a+b)P(X_n\le T) = a(1-F(T))+(a+b)F(T) = a+bF(T).

Die durchschnittliche Zykluszeit ist

\mu = E(T_n) = \int_0^\infty P(T_n > t)\,\mathrm{d}t = \int_0^T P(X_n > t)\,\mathrm{d}t = \int_0^T (1-F(t))\,\mathrm{d}t.

Langfristig entstehen daher Kosten pro Zeit in Höhe von

\lim_{t\to \infty} \frac{R(t)}{t} = \frac{E(R_1)}{\mu} = \frac{a+bF(T)}{\int_0^T (1-F(t))\,\mathrm{d}t}.

Zur Bestimmung des optimalen Wechselintervalls T bei bekannten Kosten a und b sowie bekannter Verteilungsfunktion F muss man die Minimalstelle dieses Ausdrucks in Abhängigkeit von T bestimmen. Das ist besonders einfach, wenn F eine stetige Dichte hat, denn dann sind F und das Integral als Funktion der oberen Grenze differenzierbar, das heißt es können die Optimierungsmethoden der Analysis verwendet werden.

Siehe auch[Bearbeiten]

Literatur[Bearbeiten]

  • Ming Liao: Applied Stochastic Processes, CRC Press 2013, ISBN 1-4665-8933-7, Kapitel 3: Renewal Processes