Lezione 1: Dalle funzioni iterate al caos deterministico: quando una farfalla può fare la differenza.

Lezione 2. Modelli dinamici 2dim. Biforcazioni e caos in fisica, ecologia, economia …e salotti culturali.

Gian Italo Bischi DESP-Dipartimento di Economia, Società, PoliticaUniversità di Urbino “Carlo Bo”gian.bischi@uniurb.ithttp://www.mdef.it/gian-italo-bischi/

Orientamatica Milano, 24 gennaio, 7 febbraio 2014

Fenomeni governati da leggi matematiche = fenomeni prevedibili

Leibniz, XVII secolo (1646-1716)

Vediamo allora che ogni cosa procede in modo matematico - cioè infallibilmente - nel mondo intero, in modo che se qualcuno avesse una sufficiente capacità di conoscere a fondo le cose,

e avesse abbastanza inte\lligenza e memoria per considerare tutte le circostanze e tenerne conto, questi potrebbe essere un profeta e potrebbe vedere il futuro nel presente come in uno specchio

Gottfried Wilhelm von Leibniz1646 -1716

Effetto di piccole perturbazioni ?

•Caduta di un grave x(t) = x0 + v0xt ; y(t) = y0 + v0y t ## ½ g t2

•Moto dei pianeti•Moto di un pendolo•Moto di un fluido …

Capitalizzazione (con interesse composto) interesse i%.

sia r = i/100 C(t+1) = C(t) + r C(t) = (1+r) C(t)

C(1) = (1+r) C(0) C(2) = (1+r) C(1)= (1+r)2 C(0) C(3) = (1+r) C(1)= (1+r)3 C(0)

Soluzione generale: C(t) = C(0) (1+r)t

Effetto di piccole differenze ?

Da "Il mistero di Marie Rogêt", Edgar Allan Poe, 1842.

Per quanto riguarda l’ultima parte della supposizione, si dovrà considerare che la più insignificante differenza nei fatti delle due vicende potrebbe dar luogo ai più importanti errori di calcolo, facendo divergere radicalmente le due sequenze dei fatti; proprio come in aritmetica un errore che in sé non ha valore, alla fine, moltiplicandosi da un punto all’altro del procedimento, produce un risultato lontanissimo dal vero.”

Dal racconto “La notte dei numeri” di Italo Calvino.

Questi sono tutti i libri maestri della ditta – dice il ragioniere, - nei cent’anni della sua esistenza [...] non c’è mai stato un ragioniere come Annibale De Canis, eppure quest’uomo infallibile , questo genio, vedi, il 16 novembre 1884, ... ecco, qui c’è un errore di quattrocentodieci lire. Nessuno se n’è mai accorto, io solo lo so, e sei la prima persona a cui lo dico: tientelo per te e non lo dimenticare! E poi se anche lo andrai a dire in giro, sei un ragazzo e nessuno ti darà retta... Ma adesso sai che tutto è sbagliato. In tanti anni, quell’errore di quattrocentosedici lire sai quant’è diventato? Miliardi! Miliardi! Hanno un bel girare le macchine calcolatrici, i cervelli elettronici e tutto il resto! L’errore è al fondo, al fondo di tutti i numeri, e cresce, cresce, cresce!

Caos DeterministicoCaos Deterministico: un ossimoro: un ossimoro

deterministico : regolare, prevedibile

fenomeni ordinati e pianificabili

caos : assenza di regole, irregolarità, imprevedibilità.

Il concetto di caos deterministico spezza questa dicotomia:

modelli matematici deterministici non lineari possono generare andamenti quasi indistinguibili da processi aleatori, ed estremamente sensibili a piccole perturbazioni

Outline

• Generare caos deterministico iterando semplici funzioniGenerare caos deterministico iterando semplici funzioni

• Un poUn po’’ di storia di storia

• Le proprietà del caos deterministico e lLe proprietà del caos deterministico e l’’effetto farfallaeffetto farfalla

• Un poUn po’’ di ordine nel caos: gli attrattori di ordine nel caos: gli attrattori

• Caos deterministico nella letteratura, cinema, arte ….Caos deterministico nella letteratura, cinema, arte ….

y = f(x)

x → → y f

Concetto di funzione

Feed-back

Funzione composta (con se stessa)

Sistemi dinamici a tempo discreto

Assegnato x0, la successione degli stati (traiettoria) si

ottiene per induzione: xt+1 = f ( xt ) t = 0,1,…

fx ( t ) x ( t + 1 )

Legge di evoluzione ottenuta mediante l’iterazione (applicazione ripetuta) di una funzione, che dallo stato al tempo t permette di calcolare lo stato al tempo successivo, t+1

… si ottiene una “traiettoria” del sistema dinamico

Per induzione, ossia iterando la f ...

x(1) = f (x(0)) x(2) = f (x(1)) = f (f (x(0))) = f 2 (x(0)) … x(t) = f t (x(0))

x (0) f x (2) ... fx (t) x (t+1) ...x (1) f

0t

tx x a

Mappa lineare: f ( x ) = a x.

Evoluzione xt+1 = a xt

x1 = a x0

x2 = a x1 = a ( a x0 ) = a² x0

x3 = a x2 = a ( a² x0 ) = a³ x0

…

xn = a xn1 = a ( a n-1 ) x0 = anx0

Soluzione in forma chiusa:

| a | < 1 (contraction) 0 < a < 1 Contraction and orientation preserving.It monotonously converges to x*=0 -1 < a < 0Contraction and orientation reversingThe iteration converges to x* = 0 through oscillations

| a | > 1 (expansion)a > 1, expanding and orientation preservingIt diverges monotonically a < 1 expanding and orientation reversingIt diverges through oscillations

• Particular (bifurcation) values

a= 1 xt = x0 constant

a= 1 xt = (-1)t x0 alternating values

nn xx 1

0 1 2 3 4 5 6 n0

3

xn

7 8 9

1

2

1n nx x x0=3

x0=0.5

. . . . . . . . . .

nn x

ax 1

Nota: xn+2 = xn

Dato x0 e la regola induttiva xn+1 = f(xn) …. Come andrà a finire?

Nel 1776 Laplace scriveva :

“Lo stato attuale del sistema della natura consegue evidentemente da quello che era all’istante precedente e se noi immaginassimo un’intelligenza che a un istante dato comprendesse tutte le relazioni fra le entità di questo universo, essa potrebbe conoscere le rispettive posizioni, i moti e le disposizioni generali di tutte quelle entità in qualunque istante del futuro”

Pierre-Simon Laplace 1749-1827

x0

x1 = f (x0)

Legge di evoluzione: x (t + 1) = f ( x (t) )

x0

x1

x1

x0

x1

x2

x0

x1

x2

x3

x4

Costruzione geometrica qualitativa delle traiettorie

lineare

)()1( txtx

x0 x1 x2x3 x4

x

n

x

n

2( 1) ( ( )) ( )x t f x t x t b b = 1.7

C(t+1) = (1+r)C(t) – bC(t)2

0 (1+r)/b

logistica

Interesse composto con tassa proporzionale al quadrato !!

grado 2

grado 22 = 4

grado 23 = 8 ..

C10 = ……… grado 210 = 1024 !!!!

C(1)= (1+r)C(0) – bC(0)2

C(2)= (1+r)C(1) – bC(1)2 = (1+r)[(1+r)C(0) – bC(0)2 ][1-b((1+r)C(0) – bC(0)2)]

C(3)= (1+r)C(2) – bC(2)2

Henry Poincaré (1903)Se conoscessimo esattamente le leggi della natura

e la situazione dell’universo all’istante iniziale,

potremmo prevedere esattamente la situazione

dello stesso universo in un instante successivo.

Ma se pure accadesse che le leggi naturali non

avessero più alcun segreto per noi, anche in tal

caso potremmo conoscere la situazione iniziale

solo approssimativamente.

Se questo ci permettesse di prevedere la situazione

successiva con la stessa approssimazione, non ci occorrerebbe di più e

dovremmo dire che il fenomeno è stato previsto.

Ma non è sempre così; può accadere che piccole differenze nelle condizioni iniziali ne producano di grandissime nei fenomeni finali..

Henry Poincaré, 1854-1912

xyBzdt

dz

xzyRxdt

dy

yxdt

dx

Edward Lorenz

(May 23, 1917–April 16, 2008)

Lorenz

0

2

-2

xn

x0=0.5

0 10 20 30 40 50 60 t

0

2

-2

xn

x0=0.499

21 2n nx x

0 10 20 30 40 50 60 t0

4

-4

x1(t

)-x2

(t)

x1(0)=0.5

0 10 20 30 40 50 60 t

x2(0)=0.499

Si dice che si è in presenza di dinamiche caotiche se:

(1) Sensitività rispetto alle condizioni iniziali

generando due traiettorie da diverse condizioni iniziali, ma arbitrariamente

vicine, esse si mantengono limitate ma la distanza fra esse cresce

esponenzialmente e dopo un tempo finito diventa dello stesso ordine di

grandezza delle variabili di stato.

(2) Transitività (o mixing):

i punti della traiettoria generata partendo da una generica condizione

iniziale ricoprono densamente una zona dello spazio delle fasi

cioè ciascun punto dell’intervallo su cui si muove tale traiettoria risulta

essere punto di accumulazione dei punti della traiettoria stessa.

(3) Esistenza di infiniti cicli repulsivi

con i punti periodici densi nella regione ricoperta dalle traiettorie caotiche.

Nota: (2) e (3) implicano (1)

Stretching & Folding (Stiramento e ripiegamento)

La Geometria del Caos

0.875 =

Kneading of the dough (impastare)

x(0)

t

c

c1= f (c)

c1

c

c3

c2

0 1

1

x0

0 100

1

t50

c1

c

c3

c2

I

J

I

J

a = 3.61 3 4

0

-1.80

x

-1.74

1.75 1.79

Self-similarity

(omotetia interna)

“Appello evangelico per l’introduzione di queste equazioni alle differenze

semplici in corsi elementari di matematica, cosicchè l’intuizione degli

studenti possa essere arricchita vedendo le cose bizzarre che succedono

con semplici equazioni non lineari. [...]”.

“Io vorrei sollecitare che sia presentata [l’equazione logistica] prestonell’educazione matematica. Questa equazione può essere presentatada un punto di vista fenomenologico iterandola con una calcolatrice,o persino a mano. Il suo studio non richiede più sofisticazione di quanto non richieda un corso elementare di matematica.

Tale studio potrebbe in generale arricchire l’intuito di uno studente circai sistemi non lineari. Non solo nella ricerca, ma anche nella vita politicaed economica di ogni giorno, noi saremmo più ricchi se un numeromaggiore di persone si rendesse conto che semplici sistemi non linearinon possiedono necessariamente semplici proprietà dinamiche.”

Robert May, 1976

Supponiamo che ogni anno si riproduca una frazione r di insetti

e ne muoia una frazione m.

Nell’anno successivo la popolazione è

tNmr

tmNtrNtNtN

1

1

Questa è una legge di evoluzione lineare xt+1 = axt

r-m rappresenta il tasso netto di crescita della popolazione

Popolazione che vive in un ambiente limitato.

Si fa l’ipotesi che il tasso di mortalità m non sia costante, ma

aumenti al crescere della numerosità della popolazione

ad esempio m = sN(t),termine di mortalità per sovraffollamento (carenza di cibo ecc.)

Con questa ipotesi la legge di evoluzione diventa non lineare:

Una funzione di secondo grado

Con un cambio di variabile diventa:

2)()(1)1( tsNtNrtN

)(1)()1( txtaxtx logistica

n

nn x

xkx

2

2

1

Punti fissi: 2

2

2

k xx x k x k

x

Iterazione della funzione (iperbole)2

( )2

k xf x

x

2

1'( )

2 2

kf x

x '( ) 0f k

x

0<f’(x*)<1

-1<f’(x*)<0

f’(x*)>1 f’(x*)< -1

Stabilità

Instabilità

Stabilità locale di un punto di equilibrio

Approx lineare in un intorno di x*:

xt+1= f(xt)= f (x*) + f ’(x*)(xtx*) + o(xt x*)

Da cui:

xt+1 x* f ’(x*)(xtx*)

Xt+1 = Xt dove Xt = xtx* spostamento dall’equilibrio

x* : | f ’(x*) | <1

retta tangente

Mappa logistica:

Punti fissi q*=0

p* = (a1)/a

f’ (x) = a(1 2x).

f’(q*) = a q* stabile per 1<a <1

f’(p*) = 2 a p* stabile per 1<a<3

ttt xaxx 11

c = a/4

0 1

p*

c-1

c

q*

Stabilità strutturale, biforcazioniStabilità strutturale, biforcazioni

• Fold bifurcation:– two fixed points appear, one stable and one unstable

xxx

Bifurcation diagram

Normal form: f(x,) = + x x2

f ’ (x*) through value 1

• Transcritical bifurcation (or stability exchange):– two fixed points merge, exchanging their stability

x x x

Bifurcation diagram

Normal form:f(x,) = x + x x2

f ’(x*) through value 1

• Pitchfork bifurcation– a fixed point becomes unstable (stable) and two further fixed

points appear, both stable (unstable)

x x x

subcriticalsupercritical

Normal form:f(x,) = x + x x3

f ’ (x*) through value 1

Normal form:f(x,) = x x + x3

• Flip bifurcation (period doubling bifurcation):– the fixed point becomes unstable and a stable period 2 cycle appears,

surrounding it. It corresponds to a pitchfork bifurcation of the second iterated of the map.

x

alfa

supercritical

subcritical

Normal form:f(x,) = x + x3

f ’ (x*) through value 1

x

y

x*

a = 2

0 .5 1

p*x

0

x1

x2

x3

x

y

x*

a = 2 .5

0 .5 1p *

x0

x1

x2

x3

x4

x5

x

y

a= 3 .1

0 10 .5

F x = f x( ) ( )2

F x = f x( ) ( )2

f x( )

f x( )

x*

x0

x*

a = 2

a = 3.1

a = 2.5

Let C = {c1, c2, …, ck} be a k-cycle of xt+1 = f(xt)

i.e. cic1 , i=2,…,k ; f(ci) = f(ci+1) , i=1, …, k-1, and f(ck)=c1

In other words:

C = {c1, f(c1),f 2(c1), …, f k-1(c1)} and f k(c1) = c1

Then c1 is a fixed point of f k. Indeed, any cj, j=1,…,k, is a fixed point of f k .

By the chain rule it is easy to compute the multiplier of C:

C = Df k (ci) = f ′ (c1) ∙ f ′ (c2) ∙… ∙ f ′ (ck) =

C is stable if |C| < 1

k

jjcf

1

)('

What we said for the fixed points of f , on their stability and local bifurcations etc., can be applied to k-cycles, seen as fixed points of f k

In particular:

A couple of k-cycles (one stable and one unstable) can be created by a fold bifurcation of f k

A k cycle can give rise to a 2k-ycle via a flip bifurcation of f k

Sharkovsky Theorem (1964).

If a k-cycle exists for f : II, then at least a p-cycle exists for each number p that follows k in the following total ranking of natural numbers:

3, 5, 7, 9, …, 3∙2, 5∙2, 7∙2, …, 3∙22, 5∙22, …, ….24, 23, 22, 2, 1

Li & Yorke Theorem (1975): Period 3 implies chaos

If f: II has a 3-cycle then:An uncountable set of points S I exists that does not include any

cycle and has the following properties:i) For any p, q S, pq,

(ii) For any q S and any periodic point p I

0)()(limmin0)()(limmax

qfpfandqfpf nn

n

nn

n

0)()(limmax

qfpf nn

n

The trajectories starting from an i.c. in S (scrambled set) are chaotic, i.e. they have the 3 properies that characterize deterministic chaos

Remark: it may occur mes(S) = 0 (invisible chaos)

Altro esempio: una mappa bimodale

xn+1 = f (xn) = a (5 xn) (xn 1)2 con a > 0

f ’ (x) = - a (x1)2 + 2a(5 x)(x 1) = a(x 1)(11 3x)

Metodo di Newton per la ricerca degli zeri di una funzione

Partiamo da x0 “sufficientemente vicino” alla radice cercata. Calcolando l’intersezione fra la tangente in (x0, f(x0) e l’asse delle ascisse si ottiene x1 = x0 - f(x0) / f '(x0) E così via

xn+1 = G(xn) = xn )('

)(

n

n

xf

xf

f (x) = 3x2 1 = 0 x* = 1/3.

Esempio:

n

n

x

x

6

13 2 xn+1 = G(xn) = xn =

n

n

x

x

6

13 2

Nota: se f(x*) = 0 allora G(x*)=x* e G’(x*)=0

f(x) = x3 3x2 +1 = 0

nn

nn

xx

xx

63

132

23

nn

nn

xx

xx

63

1322

23

xn+1 = G(xn) = xn

=

=

Un altro tipo di complessità: in presenza di più attrattori occorre studiare la struttura dei bacini di attrazione

Continuous and increasing maps•The only invariant sets are the fixed points. •When many fixed points exist they are alternatingly stable and unstable: the unstable fixed points are the boundaries that separate the basins of the stable ones.• Starting from an initial condition where the graph is above the diagonal, i.e. f(x0)>x0, the trajectory is increasing, whereas if f(x0)<x0 the trajectory is decreasing

p*

q*

r*

p*

q*

r*

f(x) = a arctan (x-1)

a = 3

a = 1

a = 0.5basinboundary

fold bifurcation

a = 0.5a = 0.2

Continuous and decreasing mapsThe only possible invariant sets are one fixed point and cycles of period 2, being f 2=f °f increasingThe periodic points of the 2-cycles are located at opposite sides with respect to the unique fixed point, the unstable ones being boundaries of the basins of the stable ones. If the fixed point is stable and no cycles exist, then it is globally stable.

f(x) = – ax3 + 1

a = 0.7

Z2

Z0

c

c-1

p

q

p

q

r

q-1

Nononvertible maps. Several preimages

x

y

Z0

Z2

0 1c-1

Noninvertible map: x’ = f (x) = a x (1– x)

= 1/2

c=a/4

Se x’ < a/4 allora

a

xaaxfx

2

'4

2

1)'(1

11

211

21

11 ,)'()'()'( xxxfxfxf

a

xaaxfx

2

'4

2

1)'(1

22

Punto di massimo c = a/4 2

1)()( 1

21

21 cfcfc

f ’ (c-1) = 0 e c = f (c-1)

Z3

Z1

cmax

p

q

cmin

Z1

z

r

Z3

Z1

cmax

p

qcmin

Z1

z

r

c-1

q-11

q-12

Dinamica dei prezzi: modello della Ragnatela

Quantità richiesta al tempo t dai consumatori

Qd = D ( p ) funzione di domanda

Quantità che viene immessa sul mercato dai produttori

Qs = S ( p ) funzione di offerta

Esempio: funzioni D e S lineari:

D(p) = a b p ; S(p) = c + d p

a, b, c, d costanti positivep

Q

SD

Equilibrio: Qd = Qs

p*

Introduciamo il tempo.

Al tempo t i consumatori decidono in base al prezzo pt osservato, ma la merce nel mercato al tempo t è stata prodotta in base a decisioni prese in un tempo precedente, perché la produzione richiede un certo lasso di tempo.

Sia t = 1 il tempo di produzione che intercorre fra decisione e immissione nel mercato

•I consumatori decidono la quantità da acquistare in base a pt

• I produttori decidono in t1 la quantità da immettere nel mercato al tempo t in base al prezzo atteso

Introduciamo aspettative statiche (o naive) :

pt = D-1 (S(pt-1 )) = f (pt-1)

Con funzioni di domanda e offerta lineari: D(p) = a b p ; S(p) = c + d p

a b pt = c d pt-1, da cui:b

cap

b

dp tt

1

atttp

1 tattt pp

)()( atttt pSpD

D(pt) = S(pt-1) diventa

a bpt = arctan ((pt-11))

da cui si ottiene il modello dinamico

pt = F(pt-1)= [a arctan ( (pt-1 1))]

La mappa F(p) è monotona decrescente

Funzione di offerta non lineare:

QS = S ( p) = arctan ((p1))

pe

S

D

D ( p ) = a - b pS ( p ) = arctan ((p - 1))

p

Q

pe

p

p1

p

p2

F F

p

00.3 0.4

5

atttp 1 = att

tp + ( pt atttp ) con 0 < 1

Aspettative adattive

Inserendo:

pt = F ( atttp ) mapping from beliefs to realizations

nell'equazione delle aspettative adattive:

atttp 1 = att

tp + ( F ( atttp ) att

tp ) = (1 atttp +

b

[a arctan ( ( att

tp 1))]

attp

)()1( attatt pFp

)( attpF

atttp

atttppe

Modelli con aspettative

In economia e nelle scienze sociali lo stato attuale consegue sì da quelli del passato, ma dipende anche dalle decisioni degli individui che lo compongono, decisioni che sono influenzate dalle aspettative che essi hanno sul futuro.

xt+1 = f ( ) oppure xt = f ( )

Le aspettative degli agenti sul futuro si riflettono sul modo in cui i sistemi evolvono: mappings from beliefs to realizations.

Gli agenti economici dei modelli devono essere dotati della capacità fare congetture sulla distribuzione di probabilità dei possibili stati futuri dell’Economia

)(1

etx

)(1

etx

Ipotesi delle aspettative razionali (Muth, 1961, Lucas, 1972)

Gli agenti economici sono in grado di prevedere correttamente il futuro dei sistemi che studiano, così come un fisico conosce le leggi della natura.

Così nasce l’agente economico rappresentativo razionale

1)(1 t

et xx

Aspettative razionali e caos deterministico. Una evidente antinomia

Se si parte da un modello con aspettative razionali e si scopre che esso genera caos deterministico, allora le previsioni non possono essere razionali (cioè perfette) per definizione di dinamiche caotiche.

Un corollario che contraddice un’ipotesi del teorema!

Benhabib, Day (1982) “A characterization of erratic dynamics in the overlapping generations model” Journal of Economic Dynamics and Control, 4, 37-55.

Boldrin, Montrucchio. (1986) “On the Indeterminacy of Capital Accumulation Paths.” Journal of Economic Theory 40: 26—39.

Grandmont, J.M. (1985) “Endogenous Competitive Business Cycles” Econometrica 53: 995—1045.

IPOTESI STANDARD

ASPETTATIVE RAZIONALI

AGENTE RAPPRESENTATIVO

RAZIONALITÀ ILLIMITATA (O PIENA)

ASPETTATIVE ADATTIVE

AGENTI ETEROGENEI

RAZIONALITÀ LIMITATA

IPOTESI DEI MERCATI EFFICIENTI

FINANZA COMPORTAMENTALE

IPOTESI ALTERNATIVE

MODELLO CON FUNZIONI LINEARI

Pt+1 = Pt + γ [DtC + Dt

F]

Comportamento dei FONDAMENTALISTI:

DtF = α (F – Pt )

Comportamento dei CHARTISTI:

α > 0

DtC = β (Pt – F) β > 0

Identifichiamo due categorie di investitori: fondamentalisti e chartisti

D = eccesso di domanda

Pt+1 = Pt + γ [β (Pt – F) + α (F – Pt)]

Cosa succede al passare del tempo al prezzo del titolo?

Occorre identificare il punto fisso o stato stazionario (P*) e successivamente verificarne la stabilità.

Il punto fisso è il valore Fondamentale.

Pt+1 = Pt = P* P* = F

Stabilità

Pt+1 = f (Pt) = (1 + γβ – γα) Pt – γβF + γαF

Modello lineare, condizione per la stabilità : 1 < (1 + γβ – γα) < 1

Pendenza POSITIVA e MINORE di 1

Pendenza NEGATIVA e MAGGIORE di -1

Pendenza POSITIVA e MAGGIORE di 1 Pendenza NEGATIVA e MINORE di -1

I chartisti sono più reattivi dei fondamentalisti (dominano il mercato)

I fondamentalisti sono troppo più reattivi dei chartisti

Cioè i fondamentalisti predominano, ma non eccessivamente

0 < (α – β) < 2 /γ

I prezzi CONVERGONO quando:

β > α

I prezzi DIVERGONO quando una delle seguenti condizioni si verifica:

1 + γ(β – α) > 1 1 + γ(β – α) < –1

(α – β) > 2 /γ

MODELLO CON CHARTISTI PRUDENTI

Comportamento dei FONDAMENTALISTI:

DtF = α (F – Pt) α > 0

Comportamento dei CHARTISTI:

β > 0 DtC = β arctan (Pt – F)

Punto fisso Pt+1 = Pt = F

1 + γ(β – α) > 1 1 + γ(β – α) < –1

Equazione di evoluzione del prezzo:

Pt+1 = f (Pt) = Pt + γ [βarctg (Pt – F)+ α (F – Pt)]

β > α(α – β) > 2 /γ

Stabilità: 1< f’(F) <1

Cosa succede quando il prezzo Fondamentale è instabile?

Indice FTSE MIB Telecom Italia

• Dinamica del prezzo:

MODELLO CON FONDAMENTALISTI PIÙ AGGRESSIVI

Comportamento dei FONDAMENTALISTI:

α > 0 Df = a(F – Pt)3

(con x = P-F)

xt+1 = f(xt) = xt + γ [βxt – αxt3] = h(xt) = (1+ γβ) xt– γα xt

3

Punti fissi:

x0* = 0 (P = F)

Gli altri due stati stazionari sono quelli che risolvono l'equazione:

1 = 1 + γβ − γα(x*)2 x1,2* = ± (P = F )

Stabilità degli equilibri

Nel punto fisso centrale (quello del prezzo fondamentale):

Quindi il punto fisso fondamentale è sempre instabile.

Negli altri punti fissi:

xd

df= 1+ γβ – 3γαx2

= 1 + γβ

= 1 – 2γβ

Stabili se 2γβ > 2 β > 1/γ

xd

df

xd

df

CAOS DETERMINISTICO E LETTERATURADal romanzo: Jurassic Park, di Michael Crichton (23 ottobre 1942, 4 novembre 2008)

Un passo tratto dalla Seconda Iterazione[…] Ian Malcom era uno dei più famosi rappresentanti di quella nuova generazione di matematici che mostravano un vivo interesse per i “meccanismi del mondo reale”. Questi studiosi,sotto molti aspetti, avevano rotto la tradizione di isolamentodei matematici.Per prima cosa si servivano continuamente del computer,cosa che i matematici tradizionali non vedevano di buon occhio.Poi lavoravano quasi esclusivamente con equazioni non lineari,nel campo emergente del cosiddetto caos.Terza cosa, sembravano voler fare di tutto il possibile affinché iloro sistemi matematici descrivessero qualcosa che di fattoesisteva nel mondo reale.

Ancora Ian Malcom, da Jurassic Park, terza iterazione.

“I computer vennero costruiti verso la fine degli anni 40, perché matematici come John Von Neumann , il massimo matematico della sua generazione, pensavano che avendo a disposizione una macchina capace di gestire contemporaneamente molte variabili, si sarebbe stati in grado di fare previsioni meteorologiche a lungo termine. […]. La teoria del caos manda all’aria tutto questo, non si può prevedere il tempo se non per pochi giorni. […] Tutto il denaro speso per previsioni meteorologiche lungo termine - circa mezzo miliardo di dollari negli ultimi decenni- è buttato via. È un’impresa vana quanto cercare di trasformare il piombo in oro. Oggi gli sforzi degli alchimisti ci fanno ridere, ma generazioni future guarderanno noi e rideranno nello stesso modo”.

Jurassic Park, terza iterazione:

“Un simile controllo è impossibile” dichiarò Ian Malcom“Invece sì” disse Hammond“Mi scusi, ma lei non sa quello che dice” ribattè Malcom“Piccolo stronzo arrogante” disse Hammond. Si alzò e uscì.“Mi spiace” disse Malcom “ma il punto è che ciò che definiamonatura è di fatto un sistema complesso, non lineare.Ci costruiamo una immagine lineare della natura e poi combiniamo pasticci.Io non sono uno di quegli ambientalisti dal cuore tenero, ma dovete capire ciò che non capite. Quante volte bisogna sbattere il muso contro l’evidenza dei fatti? Abbiamo costruito la diga di Assuan sostenendo che avrebbe rivitalizzato l’Egitto, e invece distrugge il fertile delta del Nilo,produce infestazioni da parassiti e rovina l’economia.Abbiamo costruito...

Carlo Emilio Gadda (1953) nel racconto "L’egoista" "Se una libellula vola a Tokio, innesca una catena di reazioni che raggiungono me".

Gadda (1974) Meditazione milanese

"L'ipotiposi della catena delle cause va emendata e guarita, se mai, con quella di una maglia o rete. Ogni anello o grumo o groviglio di relazioni è legato da infiniti filamenti a grumi o grovigli infiniti. Come gli gnocchi. Unti, agglutinati, filamentosi per formaggio e per salse, e uno cento ne traina, e ognuno dei cento poi mille e ognuno dei mille, milioni. Altro che le ciliegie, delle quali sogliono li esperti affermare che una tiri l’altra!"

Gadda (1957) Quer pasticciaccio brutto de via Merulana

«Il dottor Ingravallo sosteneva, fra l'altro, che le inopinate catastrofi non sono mai la conseguenza o l'effetto che dir si voglia d'un unico motivo, d'una causa al singolare: ma sono come un vortice, un punto di depressione ciclonica nella coscienza del mondo, verso cui hanno cospirato tutta una molteplicità di causali convergenti. Diceva anche nodo o groviglio, o garbuglio, o gnommero, che alla romana vuoi dire gomitolo. L'opinione che bisognasse «riformare in noi il senso della categoria di causa» quale avevamo dai filosofi, da Aristotele o da Emmanuele Kant, e sostituire alla causa le cause era in lui una opinione centrale e persistente: una fissazione, quasi» […]

Caos deterministico al cinema

2

( 1) ( ) ( ):

( 1) ( )

x t ax t y tT

y t x t b

x

y

T

T

T

T

Mappe iterate del piano

Point mapping f : Rn Rn x’ = f ( x ) x Rn

Takes a point in Rn and moves it to a new position

If S is a set of points then f(S) = { f(x) | xS}

Lineland : f : R R x’ = f(x)

f

f

A B

x x’

A’ B’f

A B A’B’f

A B A’ B’

Flatland f : R2 R2

f

S

),('

),('

2

1

yxfy

yxfx

f(x,y)

S’

(x’,y’)

2-dimensional linear maps: contractions, expansions, rotations etc.

2dim-linearx’ = ax + by +cy’ = dx + ey +f

rotation

AB

C

A’

B’C’

F’

T

Linear map T:

y

x

aa

aa

y

x

2221

1211

'

'

T is orientation preservingif det A > 0

area (F’) = |det A | area (F) |det A | < 1 (>1) contraction (expansion)

Meaning of the sign of |det A |

F

A’B’

C’

T is orientation reversingif det A < 0 F’

FA

B

C T

Henon map:nonlinear, invertible

bxy

axyxT

'

1':

2

1''

'1

:2

2

1

yb

axy

yb

xT

))),(((),( 123 yxTTTyxT

21 1,),(: axyxyxT

ybxyxT ,),(:2

xyyxT ,),(:3

transforms a line y=k into a parabola

is a linear contraction in x direction for |b|<1

is a reflection through the diagonalS

T(S)

T2(S) T3(S)

det DT = b

LC = T(LC -1)

Equivalently, we say that p’ has several rank-1 preimages

LC

LC-1

SH2

SH1

R1 R2

Z2 Z0

11T

12T

UU-1,2

U-1,1

x’

y’

y

x

T

p’

p1

p2 T.. .

Noninvertible (Many-to-One) map: Distinct points are mapped into the same point

Folding action of TLC-1

. T1-1

p’

p1

p2 T2-1. .

Unfolding action of T-1

-2

x

..

P = T(P1) = T(P2)

x

y

.1 2 3

1

2

3

-1-2-3-1

1P

2P

T

T

..

P

y

.1 2 3

1

2

3

-1-2-3-1

-2

11P

12P

11T

12T

2

( 1) ( ) ( ):

( 1) ( )

x t ax t y tT

y t x t b

byxy

byxT

byxy

byxT

''

':

''

': 1

21

1

2 inverses

Noninvertible maps: many to one

12

1

t t t

t t

x ax y

y x b

2D example T:

byxy

byxT

byxy

byxT

''

':

''

': 1

21

1

02

1

x

aDT

det DT = -2x =0 for x=0

T({x=0}) {y=b}

LC = {(x,y) | y = b }LC-1 = {(x,y) | x = 0 }

FT

LC-1

LC

F’= T(F)

bxy

yaxxT 2'

':

LC LC

LC-1LC-1

A B

B’

A’

A B

O’ C’

x x

yyCD

C’

D’

O

C

A’ B’

LC LC

LC-1LC-1

A

B

B’

A’

A

B

C’

x x

yy

C

C’

C

A’ B’

Itera…

2

( 1) ( ) ( ):

( 1) ( )

x t ax t y tT

y t x t b

miraquad

12

1

t t t

t t

x ax y

y x b

LC-1LC1

LC

LC2

LC3

LC4

LC5

T: 12

1

t t t

t t

x ax y

y x b

miraquad

miraquad

miraquad

12

1

t t t

t t

x ax y

y x b

http://paulbourke.net/fractals/

Fractals, Chaos

http://paulbourke.net/fractals/clifford/

Clifford Attractors

Definitionxn+1 = sin(a yn) + c cos(a xn)yn+1 = sin(b xn) + d cos(b yn)

where a, b, c, dAre parameters thatdefine each attractor.

a = -1.4, b = 1.6, c = 1.0, d = 0.7

a = 1.1, b = -1.0, c = 1.0, d = 1.5

a = 1.6, b = -0.6, c = -1.2, d = 1.6

a = 1.7, b = 1.7, c = 0.06, d = 1.2 a = 1.3, b = 1.7, c = 0.5, d = 1.4

a = 1.5, b = -1.8, c = 1.6, d = 0.9

Lovely renderings by Thomas Burt.

Let us choose a polynomial of the form:

f(x,y) = a x2+ b y2 + c x y + d x + e y + f

where a, b, c, d, e, f are constants discussed later and x, y are the usual coordinates in 2-space. The simplest way to turn this polynomial into a map is known as coordinate rotation : ynew = f(xold, yold)

xnew = yold

OK, say you, what about those 6 constants? A willy-nilly choice of those will not give you a strange attractor. There may be (probably is) some way to telling if a particular set of constants produces an attractor but, I am not a mathematician, which means I don't really care -- I basically write a program to pick them at random and see what comes out. Most sets diverge, some converge to a point, some converge to a boring loop, and only few produce good-looking pictures. These things are called attractors, not because they're attractive, but because they attract reasonable points Here is a short C/C++ program I wrote.

What is a strange attractor?

I, of course, do not know the formal, mathematical, definition of Chaotic Attractors, but I will do my best to correctly guess it. Strange Attractor is a collection of points such that each point is a function of another point What kind of function? Everything from polynomials to transcendentals. Like a randomly-appearing mosaic - instead of individual features appearing one after the other, dots light up and eventually compose distinct shapes.

Coupled Twisted Logistic: four parameter algebraic map of 4th degree

)()1('

)()1(':

22

11

xyxxy

yxyyxT

Chaos and symmetry

1 2 1 2

Tx y y x y

y x x y xS :' ( ) ( )

' ( ) ( )

1

1

Symmetric case:

Coutl simmetrica

http://people.mbi.ohio-state.edu/mgolubitsky/

Chaos and Symmetry. Mike Field , Martin Golubitsky,

Symmetric Icons Program

DEFDBL I, P-Q, X-Z ON ERROR GOTO errortrap DEF fnxpix (x) = nstartx + scalex * (x + scale) DEF fnypix (y) = npixely - scaley * (y + scale) GOSUB initialize GOSUB menuloops: GOSUB iterate x = xnew: y = ynew PSET (fnxpix(x), fnypix(y)) a$ = INKEY$ IF a$ = "c" THEN iterates = 1: CLS : GOSUB parameters IF a$ = "i" THEN GOSUB parametersrestart: IF a$ = "m" THEN GOSUB menu iterates = iterates + 1 GOTO loopsiterate: zzbar = x * x + y * y p = alpha * zzbar + lambda zreal = x: zimag = y FOR i = 1 TO n - 2 za = zreal * x - zimag * y zb = zimag * x + zreal * y zreal = za: zimag = zb NEXT i zn = x * zreal - y * zimag p = p + beta * zn xnew = p * x + gamma * zreal - omega * y ynew = p * y - gamma * zimag + omega * x RETURN

menu: GOSUB parameters PRINT USING "(X,y) = ##.#### ##.####"; x; y PRINT "Scale =", scale PRINT "ESC to exit program" PRINT "R for RETURN"CLS GOTO menuinitialize: CLS scale = 1! nscreen = 12: npixelx = 640: npixely = 480 nstartx = 160 SCREEN nscreen GOSUB setscreen x = .01: y = .003: n = 4: iterates = 1 lambda = -1.8: alpha = 2: beta = 0: gamma = 1!: omega = 0 RETURNinitialpoint: CLS PRINT "Enter r to reset coordinates automatically" PRINT "Enter x to INPUT coordinates"setscreen: CLS scaley = npixely / (2 * scale) scalex = (npixelx - nstartx) / (2 * scale) RETURN

Two kinds of complexity

k = 1; v1 = v2 = 0.851 ; 1= 2 =0.6 ; c1 = c2 = 3

y

x

1.5

1.500

E*

(a)

From: G.I. Bischi and M. Kopel “Multistability and path dependence in a dynamic brand competition model”, Chaos, Solitons & Fractals, vol. 18 (2003) pp.561-576

strutture complesse deibacini di attrazione nel casodi più attrattori coesistenti

http://aleph0.clarku.edu/~djoyce/newton/newton.html

Newton's method is an iterative method for finding solutions of equations.Given an estimate zn for a solution of an equation f(z) = 0, it generates the estimate zn+1, which would be exact if f  were linear. This works out to give the formula

zn+1 = zn - f(zn) / f '(zn).

Thus, beginning with an initial estimate z0, we obtain a series of estimates z0, z1, z2, ..., which, one hopes, converges to a solution of f(z) = 0.

The question of convergence turns out to be a most interesting one. Given a function f  with f(z) = 0 having more than one solution, we may ask, for an initial value z0, not only whether the series z0, z1, z2, ... converges to a solution, but, if so, which solution it converges to. In this way we arrive at the notion of a Newton basin: The Newton basin of a solution is the set of starting points z0 for which Newton's method converges to that solution.

Although you can enjoy the Newton Basin images without knowing the technical basis for them, it helps to understand it. Newton's method for finding roots of functions.Newton's method began as a method to approximate roots of functions, equivalently, solutions to equations of the form f(x)=0. Not only is the method easy to comprehend, it is a very efficient way to find the solution to the equation.

From z0 we get z1 = z0 - f(z0) / f '(z0) and so on….

The method does not work for any initial guess z0, and if several roots exist then the method may converge to one or another according to the initial condition

The method works for complex functions, too. Now, a function can have several roots. A polynomial of degree n has n roots. That's the "fundamental theorem of algebra." (Of course, you have to count multiplicities of roots) If you start close to a root, then Newton's method brings you closer to the root faster and faster. But if you start far from a root, then Newton's method may take you to some root other than the closest one. Suppose you paint the complex plane in colors so that all the points that approach a particular root are painted one color. We'll say those points lie in one Newton basin for the function. An interesting pattern emerges.

Introducing the problem: finding the basins of attraction of the roots of z2 - 1

How much harder can it be to find the basins of attraction of the roots of z3 - 1?The complex function f(z) = z3 - 1 has three roots:z1 = +1, z2 = -1/2 + i⋅(√3)/2, and z3 = -1/2 - i⋅(√3)/2.These points are equally spaced around the unit circle, so one might expect the basins of attraction of these three roots would be three 120 deg pie slices, symmetrically placed about each root.

This is not true, because the solution is very complicated, as we see here.

b=Re(a)

c=Im(a)

2

-1

1

Iterazione di Mandelbrot

f(z) = z2 a con z = x + iy e a = b+ic numeri complessi

schema iterativo sul piano complesso: 21n nz z a

2 21 1 1 1, , 2n n n n n nx y x y b x y c

equivalente a :

L’insieme di Mandelbrot è l’insieme dei valori del parametro complesso a tali che, partendo dalla condizione iniziale z0=0+i0, lo schema iterativo genera successioni limitate,

Benoit MandelbrotVarsavia, 20 novembre 1924

Cambridge, Massachusetts, 14 ottobre 2010

David Ruelle “Caso e Caos”, Bollati Boringhieri, 1992

James Gleick “Caos. La nascita di una nuova scienza”, Sansoni 1997, 3° edizione(edizione inglese : “Chaos. The amazing science of the unpredictable”)

Ian Stewart “Dio gioca a dadi ? La nuova matematica del caos” Bollati Boringhieri, 1993

Angelo Vulpiani “Determinismo e caos” La Nuova Italia Scientifica, 1994.

Douglas R. Hofstadter “Strani attrattori : schemi matematici collocati fra l’ordine e il caos” su “Le Scienze”, Febbraio 1982.

James P. Crutchfield, J. Doyne Farmer, Norman H. Packard, Robert S. Shaw “Il Caos”“Le Scienze”, Febbraio 1987

A.K. Dewdney “Alla scoperta delle strane attrattive del caos” su “Le Scienze” , Settembre 1987.

F. Di Stefano “Il caos deterministico” in La Fisica nella scuola” n.6, 1991.

Alberto Rebaglia “Il Caos e i Frattali” inserto di “Scienza e Vita”, giugno 1993.

Marco dal Bosco “Comportamenti Casuali di un sistema deterministico” in “La Fisica nella scuola” n.2, 1998.

Michele Fontana “Metti ordine nel caos” in Panorama del 21 febbraio 1988.

Pietro Greco “Il caos minaccia Newton”da “La Repubblica” del 7 febbraio 1990.

Giorgio Israel “Grande è la confusione sotto il cielo della scienza : il determinismo non è morto” da “La Repubblica” del 11 dicembre 1991.

Franco Prattico “I sacerdoti del caos”da “La Repubblica” del 30 aprile 1993.

Giampiero Borrella “Avanti caos” in Panorama del 18 luglio 1993.

Carlo Bernardini “Finiremo tutti in un grande frattale”da “La repubblica” del 3 dicembre 1986

Umberto Bottazzini “Il mondo del pressappoco” Il Sole 24 ore, 8 marzo 1987.

Omar Calabresi “Matematicamente belli” in Panorama del 17 gennaio 1988.

Pier Luigi Sacco “La finanza turbolenta si spiega coi frattali”da Il sole 24 ore del 28 febbraio 1988.

If you think you are too small to make a difference, try sleeping with a mosquitothe Dalai Lama

Caos deterministico, libero arbitrio, predeterminazione, controllo …

Anche in un mondo così rigido un piccolo evento, una minuscola azione, può provocare una rivoluzione.Questo risulta perfettamente compatibile anche nell'ambito di un sistema governato da un rigido e predeterminato modello matematico.

Contrapposizione fra:•completa e rigida regolamentazione e totale casualità degli eventi•predeterminazione e il libero arbitrio