Duffing Equazione differenziale

SCARICA Mathematica Notebook

La forma forzata più generale dell’equazione di Duffing è

 x^..+ deltax^.+(betax^3+ / - omega_0^2x)=gammacos(omegat+phi).
(1)

A seconda dei parametri scelti, l’equazione può assumere una serie di forme speciali. Ad esempio, senza smorzamento e senza forzatura,  delta = gamma = 0 e prendendo il segno più, l’equazione diventa

 x^..+ omega_0^2x + betax^3=0
(2)

(Bender e Orszag 1978, p. 547; Zwillinger 1997, p.122). Questa equazione può visualizzare un comportamento caotico. Per  beta0 , l’equazione rappresenta una “molla dura” e per  beta0, rappresenta una “molla morbida.”Se  beta0 , le curve di ritratto di fase sono chiuse.

Se invece prendiamo beta=1, omega_0=1, reimpostare l’orologio in modo che phi=0, e utilizzare il segno meno, l’equazione è quindi

 x^..+ deltax^.+(x^3-x)=gammacos (omegat).
(3)

Questo può essere scritto come un sistema di equazioni differenziali ordinarie del primo ordine come

x^. = y,
(4)
s^. = x-x^3-deltay+gammacos(omegat)
(5)

(Wiggins, 1990, p. 5) che, in unforced caso, si riduce a

x^. = y
(6)
s^. = x-x^3-deltay
(7)

(Wiggins 1990, p. 6; Ott 1993, p. 3).

I punti fissi di questo insieme di equazioni differenziali accoppiate sono dati da

 x^.= y=0,
(8)

quindi y = 0, e

s^. = x-x^3-deltay
(9)
= x(1-x^2)-0
(10)

mi x=0,+/-1. I punti fissi sono quindi (-1,0), (0,0), e (1,0).

Analisi della stabilità dei punti fissi può essere punto linearizzando le equazioni. Differenziare dà

x^.. = s^.
(11)
= x-x^3-deltay
(12)
y^.. = (1-3x^2) x^.- deltay^.,
(13)

che può essere scritto come equazione della matrice

 =.
(14)

Esame della stabilità del punto (0,0):

 |0-lambda 1; 1 -delta-lambda|=lambda(lambda+delta)-1=lambda^2+lambdadelta-1=0
(15)
 lambda_+/-^((0,0))=1/2(-delta+/-sqrt(delta^2+4)).
(16)

Ma  delta ^ 2=0, quindi  lambda_+/-^((0,0)) è reale. Poiché  sqrt (delta^2+4)|delta|, ci sarà sempre una radice positiva, quindi questo punto fisso è instabile. Ora guarda (+/-1, 0). L’equazione caratteristica è

 |0-lambda 1; -2 -delta-lambda|=lambda(lambda+delta)+2=lambda^2+lambdadelta+2=0,
(17)

che ha radici

 lambda_+/-^((+/-1,0))=1/2(-delta+/-sqrt(delta^2-8)).
(18)

Per  delta0,  R0, quindi il punto è asintoticamente stabile. Se  delta=0,  lambda_+/-^((+/-1,0))=+/-isqrt (2), quindi il punto è linearmente stabile (Wiggins 1990, p. 10). Tuttavia, se  delta in (- 2sqrt(2),0), il radicale dà una parte immaginaria e la parte reale è  0, quindi il punto è instabile. Se  delta= - 2sqrt(2), lambda _ +/-^((+/-1,0))= sqrt (2), che ha una radice reale positiva, quindi il punto è instabile. Se  delta-2sqrt(2), allora |delta|sqrt (delta^2-8), quindi entrambe le radici sono positive e il punto è instabile.

DuffingOscillatorPhasePortrait

È interessante notare che il caso speciale  delta = 0 senza forzatura,

x^. = y
(19)
s^. = x-x^3,
(20)

può essere integrato da quadrature. Differenziare (19) e pluggingin (20) dà

 x^..=y^.=x-x^3.
(21)

Moltiplicando entrambi i lati per  x^.

 x^..X^.-X^.x + x^.x^3=0.
(22)

Ma questo può essere scritto

 d/(dt) (1 / 2x^.^2-1/2x^2+1 / 4x^4)=0,
(23)

quindi abbiamo un invariante di movimento h,

 h=1 / 2x^.^2-1/2x^2+1 / 4x^4.
(24)

Risolvere per  x^.^ 2

 x^.^2=((dx)/(dt))^2=2h+x^2-1/2x^4
(25)
 (dx)/(dt)=sqrt(2h+x^2-1/2x^4),
(26)

so

 t=intdt=int(dx)/(sqrt(2h+x^2-1/2x^4))
(27)

(Wiggins 1990, p. 29).

Note that the invariant of motion h satisfies

 x^.=(partialh)/(partialx^.)=(partialh)/(partialy)
(28)
 (partialh)/(partialx)=-x+x^3=-y^.,
(29)

quindi le equazioni dell’oscillatore Duffing sono date dall’Hamiltoniansistema

You might also like

Lascia un commento

Il tuo indirizzo email non sarà pubblicato.