Równanie różniczkowe duffinga

Pobierz zeszyt Mathematica

najbardziej ogólną formą równania matematycznego jest

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

w zależności od wybranych parametrów równanie może przybierać różne formy specjalne. Na przykład, bez tłumienia i bez wymuszania,  delta=gamma=0 i biorąc znak plus, równanie staje się

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

(Bender i Orszag 1978, s. 547; Zwillinger 1997, s. 122). To równanie może wykazywać chaotyczne zachowanie. Dla beta0 równanie oznacza „twardą sprężynę”, a dla beta0 oznacza ” miękką sprężynę.”Jeśli  beta0, krzywe portret fazy są zamknięte.

jeśli zamiast tego weźmiemy  beta=1,  omega_0 = 1, zresetujemy zegar tak, aby  phi=0 i użyjemy znaku minus, równanie jest wtedy

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

można to zapisać jako układ równań różniczkowych zwyczajnych pierwszego rzędu jako

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

(Wiggins 1990, s. 5), który w niewymuszonym przypadku zmniejsza się do

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

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

PUNKTY STAŁE tego zbioru sprzężonych równań różniczkowych są podane przez

 x^.= y=0,
(8)

so  y = 0, oraz

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

daje  x=0,+/-1. PUNKTY STAŁE są więc(-1,0), (0,0), oraz (1,0).

Analiza stabilności punktów stałych może być punktowa poprzez linearyzację równań. Różnicowanie daje

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

które można zapisać jako równanie macierzowe

 =.
(14)

badanie stabilności punktu (0,0):

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

ale  delta^2=0, więc  lambda_+/-^((0,0)) jest prawdziwe. Ponieważ sqrt (delta^2+4)|delta|, zawsze będzie jeden dodatni pierwiastek, więc ten stały punkt jest niestabilny. Teraz spójrz na (+/-1, 0). Równanie charakterystyczne to

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

która ma korzenie

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

dla delta0 ,  R0, więc punkt jest asymptotycznie stabilny. If  delta=0,  lambda_+/-^((+/-1,0))=+/-isqrt (2), więc punkt jest liniowo stabilny (Wiggins 1990, str. 10). Jeśli jednak delta w (- 2sqrt(2),0), pierwiastek daje część urojoną, a część rzeczywista jest  0, więc punkt jest niestabilny. If  delta=-2sqrt(2), lambda_+/-^((+/-1,0))=sqrt (2), który ma dodatni pierwiastek rzeczywisty, więc punkt jest niestabilny. Jeśli delta-2sqrt(2) , to|delta|sqrt(delta^2-8) , więc oba pierwiastki są dodatnie, a punkt jest niestabilny.

DuffingOscillatorPhasePortrait

co ciekawe szczególny przypadek  delta=0 bez wymuszania,

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

może być zintegrowany za pomocą kwadratów. Różnicowanie (19) i pluggingin (20) daje

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

mnożenie obu stron przez  x^. daje

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

ale to można napisać

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

mamy więc niezmiennik ruchu h,

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

rozwiązywanie dla  x^.^2 daje

 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)

więc równania oscylatora Duffinga są podane przez Hamiltonianów

You might also like

Dodaj komentarz

Twój adres e-mail nie zostanie opublikowany.