next up previous
Next: Sávkettéválások. Hasonlóság, renormálás. Up: Egydimenziós leképezések Previous: Kaotikus tartomány. Ljapunov-exponens. Káosz

Teljesen kifejlett káosz. Valószínűségeloszlás. Ergodicitás. Konjugálás.Mindenütt nyújtó leképezések.

A logisztikus leképezés $r=4$ esetén koordinátatranszformációval (ún. sima konjugálás) szakaszonként lineáris alakra hozható (``háztető-leképezés''):

\begin{displaymath}y=\frac{2}{\pi}\arcsin(\sqrt{x})\end{displaymath}


\begin{displaymath}x=\sin^2\left(y\frac{\pi}{2}\right)\end{displaymath}


\begin{displaymath}y_{n+1}=\left\{\begin{array}{ll}2y_n&y_n\le 0.5 2(1-y_n)&y_n\ge 0.5\end{array}\right.\end{displaymath}

Mivel a meredekség abszolút értéke állandó, a Ljapunov-exponens nyilvánvalóan $\ln 2>0$.
Ábra: Háztető-leképezés
\begin{figure}\begin{center}
\epsfbox{roof.eps}
\end{center}
\end{figure}
Ábra 19.: Bernoulli-shift
\begin{figure}\begin{center}
\epsfbox{bern.eps}
\end{center}
\end{figure}
Bernoulli-shift.
$\displaystyle x_n$ $\textstyle =$ $\displaystyle 0.1100111100101110101...$  
$\displaystyle x_{n+1}$ $\textstyle =$ $\displaystyle 0.100111100101110101...$  

$n$ időegység elteltével az $n+1$-edik jegy határozza meg, hogy az iterált pont nagyobb vagy kisebb, mint $1/2$. Ha tehát ezt a jelen állapot ($x_0$) ismerete alapján kívánjuk megjósolni, $2^{-(n+1)}$ pontosság szükséges. Ha ez nem áll rendelkezésre, a mozgás megjósolhatatlannak, véletlenszerűnek tűnik (mert az állapot ismeretlen finom részletei határozzák meg - v.ö. ``pillangó-effektus'').

Következtetés: a kaotikus rendszerek belső instabilitása, a kezdeti feltételre való exponenciális érzékenysége gyakorlatilag meghiúsítja a fizikában szokásos előrejelzési sémát (véges pontosságú kezdőállapot $\rightarrow$ mozgásegyenletek megoldása $\rightarrow$ véges, de még kielégítő pontosságú végállapot). Hosszabb időre egy kaotikus rendszer mozgásáról csak valószínűségi kijelentések tehetők.

Valószínűségi leírás ($P(x)$ valószínűségsűrűség, $\mu(x)$ valószínűségi mérték. Frobenius-Perron-egyenlet.)

A kaotikus rendszer (jelen esetben egydimenziós leképezés) pályáját hosszú ideig követve megvizsgáljuk, hogy a fázistér valamely kiválasztott tartományában a teljes idő hányad részét tölti a rendszer. Ezt a számot tekintjük az ott-tartózkodás valószínűségének.

Valószínűségi mérték: $\mu(x)$ annak a valószínűsége, hogy az iterált pont $x$-nél nem nagyobb.

\begin{displaymath}\mu(x)=lim_{n\rightarrow \infty} \frac{1}{n}\sum_{j=1}^n \Theta(x-x_j)\end{displaymath}

Valószínűségsűrűség: $P(x)dx$ annak a valószínűsége, hogy az iterált pont $x$ és $x+dx$ között tartózkodik.

\begin{displaymath}P(x)dx =lim_{n\rightarrow \infty} \frac{1}{n}\sum_{j=1}^n \delta(x-x_j)dx\end{displaymath}


\begin{displaymath}P(x)=\frac{d\mu(x)}{dx}\end{displaymath}

A Frobenius-Perron-egyenlet:

Ábra: A Frobenius-Perron-egyenlet levezetése: a függőleges tengelyen jelzett intervallumban tartózkodás valószínűsége a vízszintes intervallumokban való tartózkodás valószínűségeinek összege.
\begin{figure}\begin{center}
\epsfbox{fig12.eps}
\end{center}
\end{figure}


\begin{displaymath}P(x)dx=P(f^{-1}_l(x))\frac{dx}{\left\vert f'\left(f^{-1}_l(x)...
...(x))\frac{dx}{\left\vert f'\left(f^{-1}_u(x)\right)\right\vert}\end{displaymath}


\begin{displaymath}P(x)=\frac{P(f^{-1}_l(x))}{\left\vert f'\left(f^{-1}_l(x)\rig...
...f^{-1}_u(x))}{\left\vert f'\left(f^{-1}_u(x)\right)\right\vert}\end{displaymath}

Példa: mind a háztető-leképezés, mind a Bernoulli-shift esetén $\vert f'(x)\vert=2$, ezért a Frobenius-Perron-egyenlet:

\begin{displaymath}P(x)=\frac{1}{2}\left(P(f^{-1}_l(x))+P(f^{-1}_u(x))\right)\end{displaymath}

Ennek megoldása $P(x)=1$ (konstans valószínűségsűrűség, minden $x$ (0 és 1 között) egyformán valószínű.

Sima konjugálás (koordinátatranszformáció) hatása a valószínűségsűrűségre:


\begin{displaymath}y_n=u(x_n)\end{displaymath}


\begin{displaymath}\tilde P(y)dy=P(x)dx\end{displaymath}


\begin{displaymath}\tilde P(y)=\frac{P(u^{-1}(y)}{\left\vert u'\left(u^{-1}(y)\right)\right\vert}\end{displaymath}

vagy


\begin{displaymath}P(x)=\tilde P(u(x))\vert u'(x)\vert\end{displaymath}

Példa: logisztikus leképezés. $u(x)=\frac{2}{\pi}\arcsin(\sqrt{x})$, $\tilde P(y)=1$. Ezért

\begin{displaymath}P_{logisztikus}(x)=\frac{1}{\pi\sqrt{x(1-x)}}\end{displaymath}

Ábra: A logisztikus leképezés stacionárius valószínűségsűrűsége $r=4$ esetén (teljesen kifejlett káosz).
\begin{figure}\begin{center}
\epsfbox{fig13.eps}
\end{center}
\end{figure}

A Ljapunov-exponens kifejezése a stacionárius valószínűségsűrűséggel:


$\displaystyle \lambda$ $\textstyle =$ $\displaystyle \lim_{n\rightarrow \infty}\frac{1}{n} \sum_{j=1}^n\ln\left\vert f'(r,x_j)\right\vert$  
  $\textstyle =$ $\displaystyle \int dx \lim_{n\rightarrow \infty}\frac{1}{n} \sum_{j=1}^n\delta(x-x_j)\ln\left\vert f'(r,x)\right\vert$  
  $\textstyle =$ $\displaystyle \int dx P(x) \ln\left\vert f'(r,x)\right\vert$  

A Ljapunov-exponens nem változik sima konjugálás esetén:


\begin{displaymath}y=u(x)\end{displaymath}


\begin{displaymath}g(y)=u(f(u^{-1}(y)))\end{displaymath}


\begin{displaymath}g'(y)=u'(f(u^{-1}(y)))f'(u^{-1}(y))\left(u^{-1}(y)\right)'=\frac{u'(f(u^{-1}(y)))}{u'(u^{-1}(y))}f'(u^{-1}(y))\end{displaymath}


\begin{displaymath}g'(y_n)=\frac{u'(x_{n+1})}{u'(x_n)}f'(x_n)\end{displaymath}


\begin{displaymath}\ln\left\vert g'(y_n)\right\vert=\ln\left\vert f'(x_n)\right\...
...ft\vert u'(x_{n+1})\right\vert-\ln\left\vert u'(x_n)\right\vert\end{displaymath}


$\displaystyle \lambda$ $\textstyle =$ $\displaystyle \lim_{n\rightarrow \infty}\frac{1}{n} \sum_{j=1}^n\ln\left\vert g'(r,y_j)\right\vert$  
  $\textstyle =$ $\displaystyle \lim_{n\rightarrow \infty}\frac{1}{n} \sum_{j=1}^n\ln\left\vert f'(r,x_j)\right\vert$  
  $\textstyle +$ $\displaystyle \lim_{n\rightarrow \infty}\frac{1}{n} \left(\ln\left\vert u'(x_{n+1})\right\vert-\ln\left\vert u'(x_0)\right\vert\right)$  
  $\textstyle =$ $\displaystyle \lim_{n\rightarrow \infty}\frac{1}{n} \sum_{j=1}^n\ln\left\vert f'(r,x_j)\right\vert$  

Feladat: mutassuk meg a Ljapunov-exponens sima konjugálásra vonatkozó invarianciáját a valószínűségsűrűséggel felírt alakból kiindulva!


next up previous
Next: Sávkettéválások. Hasonlóság, renormálás. Up: Egydimenziós leképezések Previous: Kaotikus tartomány. Ljapunov-exponens. Káosz
Bene Gyula 2004-05-10