<?xml version="1.0" encoding="UTF-8"?><!-- generator="wordpress.com" -->
<rss version="2.0"
	xmlns:content="http://purl.org/rss/1.0/modules/content/"
	xmlns:wfw="http://wellformedweb.org/CommentAPI/"
	xmlns:dc="http://purl.org/dc/elements/1.1/"
	>

<channel>
	<title>kaoticna-dinamika &amp;laquo; WordPress.com Tag Feed</title>
	<link>http://wordpress.com/tag/kaoticna-dinamika/</link>
	<description>Feed of posts on WordPress.com tagged "kaoticna-dinamika"</description>
	<pubDate>Mon, 07 Jul 2008 05:18:12 +0000</pubDate>

	<generator>http://wordpress.com/tags/</generator>
	<language>en</language>

<item>
<title><![CDATA[4.669201609102990 (III.)]]></title>
<link>http://elgrunon.wordpress.com/2007/09/14/4669201609102990-iii/</link>
<pubDate>Fri, 14 Sep 2007 14:22:55 +0000</pubDate>
<dc:creator>El Gruñón</dc:creator>
<guid>http://elgrunon.wordpress.com/2007/09/14/4669201609102990-iii/</guid>
<description><![CDATA[…nastavak članka “4.669201609102990 (II.)“…
Logističko preslikavanje je pokazalo nevjeroja]]></description>
<content:encoded><![CDATA[<p><em>…nastavak članka “<a href="http://elgrunon.wordpress.com/2007/06/25/4669201609102990-ii/">4.669201609102990 (II.)</a>“…</em></p>
<p>Logističko preslikavanje je pokazalo nevjerojatno bogatstvo dinamičkih pojava. Sustav opisan tim modelom (u ovom slučaju populacija bioloških jedinki) prolazi kroz sve faze koje dinamički sustav može manifestirati. Počinje sa stabilnim fiksnim točkama (točkasti atraktori), preko stabilnih graničnih ciklusa (koje karakterizira kvaziperiodičko titranje u dvije faze: prijelazno stanje i periodičko titranje) te na koncu završava u kaosu uz tipična svojstva koja ga obilježavaju: nepredvidljivost vremenskih nizova i preosjetljivost na početne uvjete. To još uvijek nije sve. Kako ću sada pokazati, ova mala i jednostavna jednadžba manifestira još jedno, možda i najfascinantnije svojstvo kaosa: red.</p>
<p><strong>Bifurkacijski dijagram</strong></p>
<p>Kad je <a href="http://en.wikipedia.org/wiki/Robert_May%2C_Baron_May_of_Oxford">Robertu Mayu</a> postalo jasno da je u režimu kaosa nemoguće vidjeti bilo što gledajući samo pojedine vremenske nizove, pokušao je dobiti širu sliku sustava... Kao što sam pokazao, nije toliko bitno iz koje točke krenete već je bitno koliko je parametar $latex r$ te gdje ste uz takav parametar završili (u jednoj točki, graničnom ciklusu ili kaosu). To se može sumarizirati na jednom jedinom dijagramu. Takav dijagram onda daje sliku sustava za sve $latex r$.</p>
<p>Princip je sljedeći. Na<em> x-os</em> nanosite vrijednost biotičkog potencijala. Za svaku točku na <em>x-u</em> izračunate vremenski niz vrijednosti populacije. Ako za $latex r$ dobijete jednu stabilnu vrijednost ucrtate nju, ako dobijete periodički niz ucrtate sve vrijednosti niza iznad te točke na <em>x-u</em>... Uglavnom, dobijete nešto ovakvo:</p>
<p style="text-align:center;"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/bif2.png" height="212" width="300" /></p>
<p>Dijagram se naziva <strong><em>bifurakcijski dijagram, </em></strong>a proces <em>grananja </em>fiksnih točaka ovisno o nekom parametru <strong><em>bifurkacija.</em></strong></p>
<ul>
<li>Za $latex r \leq 1$ Postoji samo jedna fiksna točka - populacija izumire.</li>
<li>Za $latex 1 &#60; r \leq 3$ populacija ima također jednu fiksnu točku - preživljavanje</li>
<li>U $latex r = 3$ dolazi do prvog pucanja: pojavljuje se atraktor perioda 2, tj. populacija sada ima dvije fiksne točke. Povećanjem $latex r$ period oscilacije postaje sve duži i duži.</li>
<li>U $latex r \approx 3,5$ period se, nakon što se rastegnuo koliko god je mogao, udvostručuje te imamo atraktor perioda 4.</li>
<li>U području kaosa za $latex r &#62; 3,6$ , <em>x</em> za dani <em>r </em>prolazi kroz praktički sve vrijednosti nekog pod intervala intervala [0, 1].</li>
</ul>
<p align="left">Oke, ništa što do sada nismo znali, zar ne? Dakle, povećanjem $latex r$ dolazi do višestrukih udvostručenja perioda atraktora. Negdje oko $latex r = 3,6 $ odjednom stvar pukne i nađemo se u kaosu. Sada bismo očekivali da stvar bude jednostavna.  Daljnjim povećanjem perioda, sustav postaje sve <em>kaotičniji</em>. Na krajnjoj vrijednosti od $latex r = 4$, <em>x </em>prelazi preko svih vrijednosti u intervalu [0,1]: cijeli interval je postao atraktor.</p>
<p align="left"><strong>Kaskade perioda 3</strong></p>
<p align="left">Neće ići. Pojavljuju se praznine u dijagramu. Što one znače? Pokušajte dobiti vremenski niz s parametrom 3,835 (duboko u kaotičnom području)  i <em>voila</em>:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/period3.png" /></p>
<p align="left">granični ciklus perioda tri, populacija stabilno titra između tri vrijednosti zauvijek:</p>
<p align="center">$latex ... \rightarrow 0,1521 \rightarrow 0,4945 \rightarrow0,95686 \rightarrow ...$</p>
<p align="left"> Kaosa kao da i nije bilo... Daljnjim povećanjem r ovo titranje prolazi sve faze udvostručenja perioda, pojavljuju se atraktori perioda 6, 12, 24,... Na kraju opet završi u kaosu. To je značenje bijelih pruga u sivim vrpcama kaosa. <strong>Sustav ne samo da prelazi iz reda u kaos, nego se duboko u režimu kaosa nalaze mali <em>prozori</em> reda.</strong> Kaos i red se izmjenjuju u neobičnoj igri. No, ni to još nisu sva iznenađenja koja nam je <em>kaotični majmun</em> pripremio...</p>
<p align="left"><strong>Fraktalna svojstva bifurkacijskog dijagrama</strong></p>
<p align="left">Da bi se dobila još bolja slika sustava nije dovoljno samo povećavati biotički potencijal. Treba ga povećavati u sve sitnijim koracima. To je pak ekvivalentno tome da uzmete komadić gornjeg dijagrama, izrežete ga i povećate. Sada umjesto da promatrate r na intervalu [0, 4] gledate na intervalu npr. [3, 3.5]... Ili možda na intervalu [3.5678943, 3.5678945]... Danas, uz pomoć računala, to je relativno jednostavan posao. U doba Maya i Feigenbauma to je bio posao od par dana/tjedana čačkanja po kalkulatoru, računanja vremenskih nizova i ucrtavanja rezultata na dijagram... Nije da nisu imali računala, ali uz to što su bila spora, pristup njima je bio vremenski ograničen i većina se znanstvenika umjesto da čeka svojih par sati u danu (ili čak tjednu) pristupa računalu, i dalje služila puno pristupačnijom metodom: kalkulatorom, olovkom i papirom.</p>
<p align="left">Ako dakle napravite dijagram za podinterval intervala vrijednosti dobijete nešto ovakvo:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/zoombif.png" height="245" width="430" /></p>
<p align="left">Nad nekim podintervalom intervala [0, 4] rekonstruirani bifurkacijski dijagram ima gotovo potpuno isti izgled. Takvo <em>zoomiranje</em> se naravno može nastaviti u nedogled jer se zasniva čisto na računanju sa sve više i više decimalnih mjesta, a ne na fizičkom povećanju slike. Riječ je samo o tome koliko jako računalo imate ili, u neka davna vremena koliko strpljenja s kalkulatorom imate...</p>
<p align="left">Izrazita samosličnost, zadržavanje oblika bez obzira na skalu (povećanje)... Naravno, bifurkacijski dijagram je fraktal. Fraktali-&#62;kaos-&#62;fraktali-&#62;kaos.... jedno ne ide bez drugoga...</p>
<p align="left"><em><strong>Početak igre velikog igrača... </strong></em></p>
<p align="left">Malo povijesti. Robert May je 1971. g. proučavao logističko preslikavanje kao mogući model ponašanja populacije bioloških jedinki. Uz njega su iste godine još trojica znanstvenika (Nicholas Metrropolis, Paul Stein i Myron Stein) radili na istoj stvari i otkrili začuđujuću kompleksnost ove male jednadžbe. Feigenbaum se s druge strane interesirao za nelinearne probleme. Do tada je već isprobao sve poznate pristupe nelinearnosti, međutim nijedan mu se nije činio zadovoljavajući. Paul Stein ga je obavijestio o otkrićima vezanima uz populacijsku jednadžbu. Feigenbaum je smatrao kako su stvari postale još gore: <em>"Ako je najjednostavnije nelinearno preslikavanje praktički neshvatljivo, kakve onda nade ima za realističnu nelinearnu dinamiku?"</em></p>
<p>Prikaz dinamičkog sustava pomoću bifurkacijskog dijagrama odigrati će ključnu ulogu u Feigenbaumovom radu... Pokušavajući nastaviti tamo gdje je May stao, Feigenbaum je prvo počeo crtati to čudesno <em>smokvino stablo... </em>Naravno, trebale su mu brojke. Nije bilo tako lako dobiti pristup računalu, a njegov kalkulator, iako ga je bilo moguće programirati, je radio prilično sporo... Trebalo je nekako ubiti vrijeme dok je čekao da se iskristalizira sljedeća točka na bifurkacijskom dijagramu... Sam dijagram i ovo<em>"ubijanje vremena" </em> biti će ključni za Feigenbaumovo otkriće...</p>
<p><em>…to be kontinjued…</em></p>
]]></content:encoded>
</item>
<item>
<title><![CDATA[4.669201609102990 (II.)]]></title>
<link>http://elgrunon.wordpress.com/2007/06/25/4669201609102990-ii/</link>
<pubDate>Mon, 25 Jun 2007 17:03:45 +0000</pubDate>
<dc:creator>El Gruñón</dc:creator>
<guid>http://elgrunon.wordpress.com/2007/06/25/4669201609102990-ii/</guid>
<description><![CDATA[&#8230;nastavak članka &#8220;4.669201609102990 (I.)&#8220;&#8230;
Vidjeli smo do sada da jednostav]]></description>
<content:encoded><![CDATA[<p><em>...nastavak članka "<a href="http://elgrunon.wordpress.com/2007/06/23/4669201609102990-i/">4.669201609102990 (I.)</a>"...</em></p>
<p>Vidjeli smo do sada da jednostavna (ali nelinearna) jednadžba zvana logističko preslikavanje može biti korištena kao grub ali prilično ilustrativan model pri proučavanju populacijske dinamike:</p>
<p align="center"> $latex x_{n+1} = r x_n (1-x_n)$</p>
<p>Kao i kod svih dinamičkih sustava, zanima nas (dugoročno) predviđanje ponašanja konkretnog sustava za različite početne uvjete. Sustav smo konkretizirali odabirom parametra rasta populacije ($latex r$, poznat još i pod imenom <em>biotički potencijal</em>)  te smo odabirom početnih uvjeta ($latex x_0$) promatrali ponašanje sustava u diskretnim vremenskim razmacima iterirajući danu jednadžbu.</p>
<p>Mogli smo opaziti da bez obzira na početne uvjete sustav uvijek završi u jednom od dva stabilna stanja: populacija ili preživi ili izumre. Naš dinamički sustav se ponaša potpuno deterministički. Sada ćemo vidjeti da stvari postaju malo složenije. Ako dopustimo da biotički potencijal naraste iznad 3, stvari se polako otimaju kontroli...</p>
<p><strong>U početku bijaše...manje od 3</strong></p>
<p>Podsjetimo se. Za bilo koju vrijednost parametra $latex r$ između 0 i 3, populacija će se stabilizirati na nekoj vrijednosti ovisnoj samo o $latex r$ (ne i o $latex x_0$). Grafički, za jedan takav sustav, to izgleda ovako:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/stabilna.png" align="middle" height="313" width="386" /></p>
<p align="left">Može se pokazati da vrijedi sljedeći teorem.</p>
<p align="left"><strong>TM1 </strong>Logističko preslikavanje ima dvije fiksne točke. (<em>fiksna točka nekog preslikavanja je točka koja se preslikava sama u sebe, tj. vrijedi </em>$latex f(x_f) = x_f$)</p>
<p align="left"><u><em>dokaz</em></u><em> </em>Uvjet da bi naša jednadžba imala fiksnu točku glasi:</p>
<p align="center">$latex r x_f (1-x_f) = x_f $</p>
<p align="left">Nakon što malo raspišemo dobivamo:</p>
<p align="center">$latex x_f (rx_f + 1 - r) = 0$</p>
<p align="left">Iz čega direktno slijede rješenja kao sljedeće dvije točke $latex x_{f1} = 0$ i $latex x_{f2} = \dfrac{r-1}{r}$</p>
<p align="right"><em>QED</em></p>
<p align="left">Vidimo da druga fiksna točka postoji samo ako je <em>r &#62; 1</em> jer u protivnom fiksna točka postaje negativan broj koji nije dio domene logističkog preslikavanja (striktno govoreći jest, ali se ne uzima jer nema fizikalnog značenja - populacija ne može biti negativna).</p>
<p align="left">Uvedimo sada pojam stabilnosti fiksne točke.<em> Da bismo odredili stabilnost neke točke preslikavanja moramo promatrati ponašanje prve derivacije danog preslikavanja u točki kojoj želimo provjeriti stabilnost. Za fiksnu točku x preslikavanja f vrijedi sljedeće:</em></p>
<ul>
<li><em><strong>&#124;f' (x)&#124; &#60; 1 - stabilna</strong> Stabilne fiksne točke signaliziraju mogućnost dugoročnog predviđanja. Fizikalno predstavljaju mjesta nečeg što je poznato kao termin stabilne ravnoteže. Analogija s kuglicom na dnu zdjele: ako malo pomaknete kuglicu ona se zakotrlja i opet vrati na dno zdjele. Analogija s tekućinom: stabilne točke se nazivaju i ponorima - tekućina utječe u njih. Kako se stabilne fiksne točke ponašaju poput magneta za sve oko njih, iz latinskog je preuzet naziv koji se u dinamici često rabi: <strong>atraktori.</strong> </em></li>
<li><em><strong>&#124;f' (x)&#124; &#62; 1 - nestabilna</strong> Nestabilne fiksne točke. Analogne su pojmu labilne ravnoteže. Dokle god se nalazite u samoj nestabilnoj točki ne idete nikamo. Međutim, minimalni pomak u okolinu te točke rezultira nepovratnim udaljavanjem od nje. Ovo se može zamisliti kao kuglica na brežuljku: može se postići da stoji na njemu ali i najmanji pomak i ona se kotrlja prema dnu - sama od sebe neće se vratiti natrag u ravnotežno stanje na vrhu brežuljka. U analogiji s fluidima, ovakve točke se nazivaju izvorima.</em></li>
</ul>
<p><font color="#999999"><em>mala napomena: Ovo nije potpuna karakterizacija fiksnih točaka. Ovim problemom se bavio i elegantno ga riješio (za dvije dimenzije) Poincare. U jednom od sljedećih članaka ću to obraditi...</em></font></p>
<p align="left">Malo <a href="http://classes.yale.edu/fractals/chaos/FixedPoints/LogisticFixedPoints/StableFPProof.html">računanja</a> nam daje sljedeće:</p>
<ul>
<li>$latex x_{f1} = 0$ je stabilna za $latex 0 \leq r &#60; 1$</li>
<li>$latex x_{f2} = \dfrac{r-1}{r}$ je stabilna za $latex 1 &#60; r &#60; 3$</li>
</ul>
<p>Stabilnost danih fiksnih točaka se manifestira upravo u činjenici da za bilo koji početni uvjet $latex x_0$ naš sustav (prije ili kasnije) završi na kraju u jednoj od tih točaka (pod uvjetom da je i biotički potencijal unutar zadanih vrijednosti).</p>
<p><strong>Geometrijska interpretacija</strong></p>
<p>Uz čisto numeričko iteriranje naše funkcije, postoji još jedan, geometrijski način da se promatra iteracija . To je metoda pomoću tzv. paučinastih dijagrama. prednost ove metode je ilustrativnost i elegancija iz čega je uvijek lakše nešto zaključiti nego kada baratate hrpom suhih brojaka.</p>
<p>Ideja je sljedeća. Znamo da je $latex x_{n+1} = f(x_n)$. To znači da moramo na horizontalnu os nanijeti točku $latex x_n$ i iz nje krenuti vertikalno dok ne <em>udarimo </em>u graf funkcije $latex y = f(x_n)$. Sada tako dobivena vrijednost mora biti ponovo nanesena na horizontalnu os kako bi se mogla upotrijebiti za novu iteraciju. Najlakši (!?) način da y-vrijednosti ponovo nanesemo na x-os je da <em>presavijemo </em>x-y ravninu oko pravca y=x. Iz ovoga su ljudi napravili sljedeći algoritam za iteriranje našeg preslikavanja:</p>
<ul>
<li>počni s $latex x_o$ na <em>x-osi</em></li>
<li>crtaj vertikalnu liniju dok ne <em>udariš </em>graf funkcije <em>f(x)</em></li>
<li>pomakni se horizontalno do pravca <em>y=x</em></li>
<li>pomakni se vertikalno (gore ili dolje) do grafa <em>f(x)</em></li>
<li>ponavljaj zadnja dva koraka generirajući nove točke</li>
</ul>
<p>Zbunjeni? Izgleda vrlo jednostavno jednom kad se nacrta:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/tockasti_atr.png" align="middle" height="321" width="323" /></p>
<p align="left">Odmah se vidi da (u okviru zadanog <em>r</em>) ovaj sustav ima jedan točkasti atraktor.</p>
<p align="left"><strong>Nešto je u broju 3: put u kaos</strong></p>
<p align="left">Počnimo  sada polako povećavati biotički potencijal iznad brojke 3. Dešava se nešto neobično. Sustav je izbačen iz ravnoteže i traži novi ravnotežni položaj. Za $latex r = 3,2$ i $latex x_0 = 0,3$ nakon početne nestabilnosti niz brojeva izgleda ovako:</p>
<p align="center">$latex ... \rightarrow 0,7995 \rightarrow 0,5130 \rightarrow 0,7995 \rightarrow 0,5130 \rightarrow ...$</p>
<p align="left">tj. sustav ulazi u stabilnu oscilaciju između dva broja ili, kako se to češće kaže sustav oscilira s periodom 2:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/period2.png" align="middle" height="313" width="386" /></p>
<p align="center">&#160;</p>
<p align="left">Što se tu zapravo desilo? Fiksne točke i dalje postoje. Ako sustav krene iz jedne fiksne točke, uopće neće doći do oscilacije. populacija se neće uopće mijenjati s vremenom. Međutim, čim malo odlutate od fiksne točke, sustav odlazi u oscilaciju.</p>
<p align="left">Nastavimo povećavati biotički potencijal. Nakon nekog vremena niz brojeva doživljava još jednu nestabilnost te izgleda ovako:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/period4a.png" align="middle" height="313" width="386" /></p>
<p align="center">&#160;</p>
<p align="left">Primjećujete pravilnost? Sustav i dalje pokušava ostati u ravnoteži. Period dva nije više bio dovoljan. Sada populacija oscilira između četiri vrijednost, odnosno s periodom četiri. Period se udvostručio, a atraktor ima četiri točke. Nastavimo li i dalje povećavati biotički potencijal, period titranja se nastavlja udvostručavati. Pojavljuju se atraktori perioda 8, 16, 32, ... Međutim, u jednom trenutku stvar pukne:</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/log_kaos.png" align="middle" height="313" width="386" /></p>
<p align="left">Period se udvostručio praktički u beskonačnost. U tom trenutku sustav za dva gotovo ista početna uvjeta može dati potpuno različite vremenske nizove (na donjoj slici početni uvjeti se razlikuju za samo 0,001):</p>
<p align="center"><img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/uspored_kaos.png" /></p>
<p align="left">Nepravilnost dobivenih vrijednosti iz determinističke jednadžbe. Osjetljivost na početne uvjete. Očito je, ušli smo u područje kaosa.</p>
<p align="left">Što se zapravo desilo? May nije mogao vjerovati. Jedan te isti sustav koji se ponašao tako pristojno jednostavno je poludio. Iz potpuno stabilnog i predvidljivog stanja, preko kvazistabilnog periodičkog, sustav je otišao u nešto što je izgledalo kao hrpa slučajnih brojeva. Vrijednosti koje je dobio nije mogao objasniti. Probao je mnogo simulacija povećavajući biotički potencijal u sve manjim koracima s uvijek istim rezultatom. Bio je zadivljen kompleksnošću koju je dobio. Nije uočio uzorak koji će kasnije postati legendarnim...</p>
<p align="left">Duboko u području kaosa oscilacije izgledaju potpuno nepravilne i nepredvidljive. Ova pojava je dugo vremena bila zanemarivana i nitko se nije ni trudio dublje je proučiti. Pogledate li gornju sliku, to i ne čudi. Tko bi pomislio da ovdje postoji ikakva pravilnost? Kako će se pokazati, sustav i dalje skakuće po "ravnotežnim" stanjima. Ona su neobična<em>, </em>teško ih je na prvi pogled otkriti... Zapravo, populacija i dalje ima atraktor. Samo je malo <em><strong>čudan</strong>. </em>Da bi ga se otkrilo bilo je potrebno računalo ili čovjek koji ima volje puno tipkati po kalkulatoru i još k tome malo varati i pokušati razmišljati unaprijed... Nekoliko godina nakon izlaska Mayevog članka, njegov rad dolazi u ruke Feigenbaumu i on shvaća kako je upravo takav primjer tražio da bi mogao započeti svoj rad...</p>
<p align="left"><em>…to be kontinjued…</em></p>
<p align="left"><em><u>Reference</u></em>:</p>
<ul>
<li><a href="http://ibiblio.org/e-notes/MSet/Logistic.htm"> Interactive Logistic map</a> - nekoliko prekrasnih java apleta koji omogućuju interaktivno istraživanje kako vremenskih nizova tako i paučinastih dijagrama logističkog preslikavanja.</li>
<li><a href="http://chaos.phy.ohiou.edu/~thomas/chaos/logistic.html">Logistic Equation: Parable of parabola</a> - Još malo matematičke pozadine logističkog preslikavanja (objašnjenja stabilnosti fiksnih točaka i graničnih ciklusa)</li>
<li><a href="http://theory2.phys.cwru.edu/~pete/java_chaos/">Razni Java Apleti</a> - vezani uz logističko preslikavanje ali i kaos općenito.</li>
<li><a href="http://www.geocities.com/chaiwahwu/chaosintro/logistic_map.html">Introduction to Nonlinear Dynamics in the Logistic Map</a></li>
<li><a href="http://cnls-www.lanl.gov/People/nbt/Book/node32.html#SECTION00330000000000000000">Graphical method explained</a></li>
<li><a href="http://nedwww.ipac.caltech.edu/level5/Sept01/May/May_contents.html">"Simple mathematical models with very complicated dynamics", Robert M. May</a></li>
</ul>
]]></content:encoded>
</item>
<item>
<title><![CDATA[4.669201609102990 (I.)]]></title>
<link>http://elgrunon.wordpress.com/2007/06/23/4669201609102990-i/</link>
<pubDate>Sat, 23 Jun 2007 08:45:27 +0000</pubDate>
<dc:creator>El Gruñón</dc:creator>
<guid>http://elgrunon.wordpress.com/2007/06/23/4669201609102990-i/</guid>
<description><![CDATA[&#8220;Volim oblake&#8230; oblake, što prolaze&#8230; ondje&#8230; divne oblake!&#8221;
Charles Bau]]></description>
<content:encoded><![CDATA[<p style="text-align:right;"><span style="font-style:italic;">"Volim oblake... oblake, što prolaze... ondje... divne oblake!"</span><br />
<span style="font-style:italic;">Charles Baudelaire</span></p>
<p>Godine 1979. jedan je nadasve neobičan čovjek otkrio nadasve neobičnu činjenicu. Stvar je bila tako malo vjerojatna da ni sam u početku nije vjerovao... U <img src="http://i126.photobucket.com/albums/p87/ElGrunon/kaoslantidablog_inlined/feigenbaumm.jpg" align="left" height="100" hspace="2" vspace="2" width="100" />pokušaju dokazivanja onog što mu se na trenutak ukazalo potrošio je mnogo kofeina i nikotina te na kraju gotovo nastradao od iscrpljenosti... Uz dokaz da nije moguće preživjeti mjesece napornog rada koristeći samo kavu, cigarete i džepni kalkulator, pružio je i kvalitativan dokaz jedne neobične matematičke pojave koja se manifestira u svima (ali baš svim!) kaotičnim sustavima. Čovjek se zvao <a href="http://en.wikipedia.org/wiki/Mitchell_Feigenbaum">Mitchell Jay Feigenbaum</a>.</p>
<p>Kao jedan od vrhunski američkih fizičara tog vremena, Feigenbaum je radio u kompleksu Los Alamosa. Riječima Iana Stewarta: <span style="font-style:italic;">"Neki od njegovih kolega ne bi se složili s tom riječju - radio - jer nitko nije baš točno znao na čemu to Feigenbaum radi. Čak ni sam Feigenbaum!" </span>Kao dio znanstvene elite tog vremena, bio je oslobođen standardne groze akademskog života: dužnosti da predaje i objavljuje radove. U to vrijeme je fizika čestica bila u zamahu, svi koji su mislili o ikakvoj karijeri u fizici, bavili su se njome. Količina magisterija i doktorata iz tog područja i u to vrijeme bila je ogromna, a Los Alamos je bio jedan od velikih svjetskih centara u kojem su se mogli vršiti eksperimenti i istraživanja tog tipa. Feigenbaum je bio dio svega toga. Često se znalo desiti da mu kolege dođu s nekim problemom koji bi im on pomogao riješiti. Oni koji se tada radili tamo kažu kako je Feigenbaum to radio usputno, pokazujući golemo znanje, ali gotovo nikakav interes za uzbudljiv subatomski svijet. Sam nije započinjao nikakva istraživanja. Ono po čemu je ubrzo nakon dolaska u Los Alamos postao poznat bile su njegove dugačke šetnje po okolnim brdima s obvezatnom cigaretom u ustima i rastresenim pogledom koji zuri u nebo. Umjesto da radi, Feigenbaum je promatrao oblake...</p>
<p>Da bih pokazao što je to zapravo otkrio uzet ću školski primjer, kako ga Stewart naziva: <span style="font-style:italic;">"veseli majmun u kavezu teorije kaosa",</span> primjer koji se nalazi u svim knjigama o kaosu. Naravno, upravo na proučavanju tog <span style="font-style:italic;">majmuna </span>Feigenbaum je došao do svog otkrića. I naravno, kao i mnogo toga u teoriji kaosa, <span style="font-style:italic;">majmun </span>u početku nije imao nikakve veze s fizikom. Zapravo, u vrijeme kad je Feigenbaumu došao u ruke, više nije imao veze s ničim osim s podučavanjem studenata prve godine studija biologije i/ili ekologije kao povijesni primjer pokušaja modeliranja rasta populacije bioloških jedinki; povijesni primjer koji se više ne koristi jer postoje bolji i razvijeniji modeli za tu svrhu.</p>
<p><strong>Logističko preslikavanje (</strong><span style="font-weight:bold;"><strong>p</strong>opulacijska jednadžba)</span></p>
<p>Zamislite da imate malu baru ili jezerce i u njemu populaciju riba koju možete izbrojati. Recimo da želite predvidjeti koliko riba će biti u toj bari sljedeće godine u isto vrijeme. Svakome tko se iole bavio biologijom biti će odmah jasno da na odgovor na to pitanje utječe nekoliko faktora poput količine dostupne hrane, postojanje i broj grabežljivaca koji se hrane vašim ribama, eventualna pojava bolesti i otpornost na istu, itd... Međutim, u najjednostavnijem modelu, populacija riba iduće godine ovisit će samo o broju riba ove godine. To možemo zapisati kao jednostavnu rekurzivnu formulu:</p>
<p style="text-align:center;">$latex x_{n+1} = F (x_n) $</p>
<p style="text-align:left;">pri čemu je:</p>
<ul>
<li>$latex x_{n+1}$ - broj riba iduće godine,</li>
<li>$latex x_n $- trenutni broj riba,</li>
<li>$latex F$ - funkcijska ovisnost $latex  x_{n+1}$ o $latex x_n$</li>
</ul>
<p>Ako sad znate broj riba u tekućoj godini i na koji način broj riba u sljedećoj godini ovisi o trenutnom broju riba ($latex F$) tada možete izračunati broj riba u idućoj godini... Ponavljajući postupak (<em>iterirajući)</em>, uvijek uvrštavajući trenutnu vrijednost broja riba da biste dobili sljedeću, možete promatrati kako se populacija mijenja kroz godine. Očito je da ćete dobiti bitno različite nizove vrijednosti za različito odabrane F $latex F$. Ili možda ipak? O tome malo kasnije... Za ekologe je sigurno najbitnije pitanje: koji $latex F$ uopće odabrati?</p>
<p>Najjednostavnije bi bilo misliti kako će svake godine populacija porasti za neki postotak. Ovo bi bio idealan teoretski slučaj sustava u kojem nema ograničavajućih faktora za populaciju (nema grabežljivaca, bolesti i hranidbeni resursi su neograničeni). U tom slučaju bi naša jednadžba poprimila ovakav (čisto linearni) oblik:</p>
<p align="center">$latex x_{n+1} = rx_n$</p>
<p align="left">pri čemu $latex r$ jednostavno određuje koliko puta populacija poraste u godinu dana. Ovakav model daje:</p>
<ul>
<li> neograničen rast populacije za $latex r \geq 1$</li>
<li>neograničen pad populacije za  $latex 0 \leq r &#60; 1$</li>
</ul>
<p>Naravno, sasvim je očito da ovo nije realističan model. Realno se dešava nešto drugo. Populacija neko vrijeme raste. Zatim zbog ograničenih resursa u ekosistemu raste sve sporije dok ne dosegne neku kritičnu točku nakon koje počinje opadati. Opada sve dok se resursi ne oporave te opet počinje rasti. Kroz dugi niz godina, može se desiti da se ustali na nekoj stabilnoj vrijednosti ili ne, što opet ovisi o sustavu. Kako ovo modelirati?<br />
1971. ekolog i matematičar <a href="http://en.wikipedia.org/wiki/Robert_May%2C_Baron_May_of_Oxford">Robert May</a> bavio se jednostavnim modelom populacijske dinamike proizašlim iz općenite populacijske jednadžbe. Njegov model bio je ovakav:</p>
<p align="center">$latex x_{n+1} = r x_n (1-x_n)$</p>
<p>    Imamo sada dva faktora koji nam određuju populaciju: $latex x_n$ i $latex 1-x_n$ te naravno faktor rasta $latex r$. S obzirom da $latex 1-x_n$ opada kada $latex x_n$ raste i obrnuto, ovim jednostavnim modelom moguće je teoretski prikazati ponašanje populacije koje je puno bliže realnoj situaciji. Ostaje još pitanje izbora faktora $latex r$. I tu sada nastaju problemi. Što se faktora rasta $latex r$ tiče, možete jednostavno početi od nule i povećavati ga, te tako za svaku dobivenu jednadžbu ubaciti par početnih x-eva i pogledati kako se ponaša populacija kroz godine.  Zapravo, dosadan numerički posao kojeg je i May probao raditi i to u vrijeme kad su računala bila alat kojem su malobrojni znanstvenici vjerovali i kojeg su još malobrojniji od njih aktivno koristili.</p>
<p>Ja imam malo više sreće pa mogu na računalu isprobavati što se događa. Probajmo... Za $latex r = 1,9$ i početnu populaciju $latex x_0 = 0,359$ dobivam niz vrijednosti:</p>
<p>$latex 0,359 \rightarrow 0,4372 \rightarrow 0,4675 \rightarrow 0,4730 \rightarrow 0,4736 \rightarrow 0,4737 \rightarrow 0,4737 \rightarrow 0,4737 ...$</p>
<p>Populacija se ustali na vrijednosti 0,4737 nakon svega nekoliko godina.  Ovo nije uvijek tako. Zapravo sistematski to izgleda ovako:</p>
<ul>
<li>$latex 0 \leq r \leq 1$ - populacija će izmurijeti bez obzira na vrijednost $latex x_0$</li>
<li>$latex  1 &#60; r \leq 2$ - populacija konvergira broju $latex \dfrac{r-1}{r}$</li>
<li>$latex  2 &#60; r \leq 3$ - populacija također konvergira broju $latex \dfrac{r-1}{r}$, međutim prije nego se ustali oscilira neko vrijeme oko te vrijednosti. Kad $latex r$ dostigne 3, konvergencija je jako spora ali populacija na kraju ipak dostiže stabilno stanje.</li>
<li>Za $latex r \geq 4$ populacija divergira u $latex \pm\infty $ za gotovo svaki početni $latex x_0$</li>
<li>$latex 3 &#60; r &#60; 4$ dešavaju se zanimljive stvari o kojima će biti riječi u nastavku. Zapravo, populacija se počinje ponašati kaotički.</li>
</ul>
<p>Ono što se Mayu desilo u jednom trenutku bila je posljedica evidentne nelinearnosti gore navedene jednadžbe. Način na koji je May postupio bio je simptomatski za cijelu znanost do druge polovice dvadesetog stoljeća. U jednom trenutku vrijednosti populacije za odabrani $latex r$ su počele izgledati potpuno nepredvidljive i nepravilne. May je odredio kad se to točno desi (za koji $latex r$) i tu je stao. Zašto? Pa naravno, jednadžba je nelinearna i u određenom trenutku nelinearnost se počinje manifestirati. To je trenutak u kojem model prestaje biti zanimljiv. Znanost teži "čistim" modelima koje je (relativno) lako kontrolirati i iz njih izvlačiti zaključke. Kad se u modelu pojavi nelinearnost, model postaje neupotrebljiv za proučavanje: linearan dio modela je zanimljiv, zar ne?</p>
<p>Dakle, May jednostavno kaže: model se ponaša pristojno do te i te mjere, a nakon toga postaje neupotrebljiv. Guranje smeća pod tepih? Kosturi u ormaru? To su neke stvari koje mi padaju na pamet. Ako nagomilate probleme pod tepihom kad tad će vam eksplodirati u lice. Vidjet ćemo da su stvari puno zanimljivije nego što je May vidio. A to će se desiti onog trenutka kad Feigenbaum nabasa na <em>veselog majmuna kaosa....</em></p>
<p><em>…to be kontinjued…</em></p>
<p><u><em> Dodaci:</em></u></p>
<p>- <a href="http://elgrunon.wordpress.com/files/2007/06/logisticko-preslikavanje.doc" title="Logističko preslikavanje - tablica">Logističko preslikavanje - tablica</a> - Ako imate Microsot Excell, možete downloadati tablicu u kojoj je moguće promatrati što se događa s populacijom riba ako mijenjate parametre r i x0 kod logističkog preslikavanja...</p>
]]></content:encoded>
</item>

</channel>
</rss>
