Duffing Differensialligning

Last Ned Mathematica Notebook

den mest generelle tvunget form Av Duffing ligningen er

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

Avhengig av parametrene som er valgt, kan ligningen ta en rekke spesielle former. For eksempel, uten demping og ingen tvang, delta=gamma=0 og tar plustegnet, blir ligningen

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

(Bender Og Orszag 1978, s. 547; Zwillinger 1997, s. 122). Denne ligningen kan vise kaotisk oppførsel. For  beta0 representerer ligningen en «hard fjær», og for  beta0 representerer den en «myk fjær».»Hvis beta0, er faseportrettkurvene stengt.

hvis vi i stedet tar  beta=1,  omega_0=1, tilbakestill klokken slik at  phi=0 , og bruk minustegnet, er ligningen da

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

dette kan skrives som et system av førsteordens ordinære differensialligninger som

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

( Wiggins 1990, s. 5) som i det uforstyrrede tilfellet reduserer til

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

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

de faste punktene i dette settet av koblede differensialligninger er gitt ved

 x^.= y=0,
(8)

y = 0 , og

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

gir  x=0,+/-1. De faste punktene er derfor (-1,0), (0,0), og (1,0).

Analyse av stabiliteten til de faste punktene kan være punkt ved å linearisere ligningene. Differensiering gir

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

som kan skrives som matriseligningen

 =.
(14)

Undersøker stabiliteten til punktet (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)

men delta^2=0 , så  lambda_+/-^((0,0)) er ekte. Siden sqrt (delta^2+4)|delta|, vil det alltid være en positiv rot, så dette faste punktet er ustabilt. Se nå på (+/-1, 0). Den karakteristiske ligningen er

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

som har røtter

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

for delta0,  R0 , så er punktet asymptotisk stabilt. Hvis delta=0, lambda_+/-^((+/-1,0))=+/-isqrt(2), så punktet er lineært stabilt (Wiggins 1990, s. 10). Men hvis  delta i (- 2sqrt(2),0), radikalet gir en imaginær del og den virkelige delen er  0, så punktet er ustabilt. Hvis delta= - 2sqrt(2), lambda_+/-^((+/-1,0))=sqrt (2) , som har en positiv reell rot, så poenget er ustabilt. Hvis delta-2sqrt(2), så |delta|sqrt (delta^2-8), så begge røttene er positive og punktet er ustabilt.

DuffingOscillatorPhasePortrait

Interessant, det spesielle tilfellet  delta=0 uten tvang,

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

kan integreres med kvadraturer. Differensiering (19) og pluggingin (20) gir

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

Multiplisere begge sider med x^. gir

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

men dette kan skrives

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

så vi har en invariant av bevegelse h,

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

Løsning for  x^.^2 gir

 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)

så ligningene Til Duffingoscillatoren er gitt Av Hamiltoniansystemet

You might also like

Legg igjen en kommentar

Din e-postadresse vil ikke bli publisert.