Duffing differentialligning

Hent Mathematica Notebook

den mest generelle tvungen form af duffing ligning er

 ^ ..+ deltaks^.+ (Betak^3 + / - omega_0^2 gange)=gammacos(OmegaT+phi).
(1)

afhængigt af de valgte parametre kan ligningen tage en række specielle former. For eksempel, uden dæmpning og ingen tvang, delta=gamma=0 og ved at tage plustegnet bliver ligningen

 ^ ..+ omega_0^2 gange + Betak^3=0
(2)

(1978, s. 547; 1997, s.122). Denne ligning kan vise kaotisk adfærd. For  beta0 repræsenterer ligningen en “hård fjeder”, og for  beta0 repræsenterer den en “blød fjeder.”Hvis beta0 , er faseportrætkurverne lukket.

hvis vi i stedet tager  beta=1 , omega_0=1, Nulstil uret, så phi=0 , og brug minustegnet, ligningen er så

 ^ ..+ deltaks^.+ (s^3-s)=gammacos (omegat).
(3)

dette kan skrives som et system af første ordens almindelige differentialligninger som

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

(1990, s. 5), som i det uforcerede tilfælde reducerer til

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

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

de faste punkter i dette sæt koblede differentialligninger er givet af

 ^ .=y=0,
(8)

 y=0, og

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

giver  =0,+/-1. De faste punkter er derfor (-1,0), (0,0), og (1,0).

analyse af stabiliteten af de faste punkter kan være punkt ved linearisering af ligningerne. Differentiering giver

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

som kan skrives som matricen ligning

 =.
(14)

undersøgelse af punktets stabilitet(0,0):

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

men  delta^2=0, så  lambda_+/-^((0,0)) er ægte. Siden kvm(delta^2 + 4)|delta|, vil der altid være en positiv rod, så dette faste punkt er ustabilt. Se nu på (+/-1, 0). Den karakteristiske ligning er

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

som har rødder

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

for  delta0,  R0, så punktet er asymptotisk stabilt. Hvis delta=0,  lambda_+/-^((+/-1,0))=+/-iscrt (2), så punktet er lineært stabilt (1990, s. 10). Men hvis delta i (- 2skrt(2),0), den radikale giver en imaginær del, og den virkelige del er 0, så punktet er ustabilt. Hvis delta= - 2(2), lambda_+/-^((+/-1,0))=(2), som har en positiv reel rod, så punktet er ustabilt. Hvis delta-2skrt (2), så |delta|kvm(delta^2-8), så begge rødder er positive, og punktet er ustabilt.

DuffingOscillatorPhasePortrait

interessant, det specielle tilfælde  delta=0 uden tvang,

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

kan integreres af kvadraturer. Differentiering (19) og pluggingin (20) giver

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

multiplicere begge sider med  h^. giver

 ^ ..^ .- k^.h + h^.3=0.
(22)

men dette kan skrives

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

så vi har en uforanderlig bevægelse h,

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

løsning for s^.^2 giver

 ^ .^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) / (partialh)= - + + ^3= - y^.,
(29)

så ligningerne af duffing oscillatoren er givet af Hamiltoniansystemet

You might also like

Skriv et svar

Din e-mailadresse vil ikke blive publiceret.