Duffing differentialekvation

ladda ner Mathematica Notebook

den mest generella tvångsformen av Duffing-ekvationen är

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

beroende på de valda parametrarna kan ekvationen ta ett antal specialformer. Till exempel, utan dämpning och ingen tvingning, delta = gamma=0 och med plustecknet blir ekvationen

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

(Bender och Orszag 1978, s. 547; Zwillinger 1997, s.122). Denna ekvation kan visa kaotiskt beteende. För beta0 representerar ekvationen en ”hård fjäder” och för beta0 representerar den en ”mjuk fjäder.”Om  beta0 är fasporträttkurvorna stängda.

om vi istället tar  beta=1, omega_0=1, Återställ klockan så att phi=0 och använd minustecknet, ekvationen är då

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

detta kan skrivas som ett system av första ordningens vanliga differentialekvationer som

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

(Wiggins 1990, s. 5) som i det ofrivilliga fallet minskar till

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

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

de fasta punkterna i denna uppsättning kopplade differentialekvationer ges av

 x^.= y=0,
(8)

 y = 0, och

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

ger  x=0,+/-1. De fasta punkterna är därför (-1,0), (0,0), och (1,0).

analys av stabiliteten hos de fasta punkterna kan vara punkt genom att linjärisera ekvationerna. Differentiering ger

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

som kan skrivas som matrisekvationen

 =.
(14)

undersöka stabiliteten i punkten (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)) är verklig. Eftersom  sqrt (delta^2+4)|delta| kommer det alltid att finnas en positiv rot, så denna fasta punkt är instabil. Titta nu på (+/-1, 0). Den karakteristiska ekvationen är

 |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)

för  delta0,  R0, så är punkten asymptotiskt stabil. Om  delta=0,  lambda_+/-^((+/-1,0))=+/-isqrt (2), så punkten är linjärt stabil (Wiggins 1990, s. 10). Men om  delta i (- 2sqrt(2),0), radikalen ger en imaginär del och den verkliga delen är 0, så punkten är instabil. Om delta= - 2sqrt(2), lambda_+/-^((+/-1,0))=sqrt (2), som har en positiv verklig rot, så punkten är instabil. Om  delta-2sqrt(2), då |delta|sqrt (delta^2-8), så båda rötterna är positiva och punkten är instabil.

DuffingOscillatorPhasePortrait

intressant, det speciella fallet delta=0 utan att tvinga,

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

kan integreras med kvadraturer. Differentierande (19) och pluggingin (20) ger

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

multiplicera båda sidor med  x^. ger

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

men detta kan skrivas

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

så vi har en invariant av rörelse h,

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

lösa för  x^.^2 ger

 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å ekvationerna för Duffing-oscillatorn ges av Hamiltoniansystemet

You might also like

Lämna ett svar

Din e-postadress kommer inte publiceras.