Duffing Diferenciální Rovnice

ke STAŽENÍ Mathematica Notebook

nejobecnější nuceni podobě Duffing rovnice je

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

v Závislosti na zvolené parametry, rovnice může mít řadu speciálních forem. Například bez tlumení a bez vynucení delta=gamma=0 a při znaménku plus se rovnice stane

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

(Bender a Orszag 1978, str. 547; Zwillinger 1997, str. 122). Tato rovnice může zobrazovat chaotické chování. Pro beta0 rovnice představuje „tvrdou pružinu“ a pro beta0 představuje „měkkou pružinu“.“Pokud  beta0, křivky fázového portrétu jsou uzavřeny.

Pokud místo toho jsme si beta=1, omega_0=1, resetovat hodiny tak, že phi=0 a znaménkem minus, rovnice je pak

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

Toto může být napsáno jako systém prvního řádu obyčejné diferenciální rovnice jako

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

(Wiggins 1990, str. 5), která, v nevynucených případě, snižuje k

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

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

pevné body sada provázaných diferenciálních rovnic jsou dány

 x^.=y=0,
(8)

takže y=0, a

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

dává x=0,+/-1. Pevné body jsou proto (-1,0), (0,0), a (1,0).

analýza stability pevných bodů může být bod linearizací rovnic. Rozlišování dává

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

což lze zapsat jako maticová rovnice

 =.
(14)

Zkoumání stability bod (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)

Ale delta^2=0, tak lambda_+/-^((0,0)) je reálné. Od sqrt (delta^2+4)|delta / bude vždy jeden kladný kořen, takže tento pevný bod je nestabilní. Nyní se podívejte na (+/-1, 0). Charakteristická rovnice je

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

který má kořeny

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

Pro delta0, R0, takže bod je asymptoticky stabilní. Pokud delta=0, lambda_+/-^((+/-1,0))=+/-isqrt(2), takže bod je lineárně stabilní (Wiggins 1990, str. 10). Nicméně, pokud delta v (-2sqrt(2),0), radikální dává imaginární část a reálná část je 0, takže bod je nestabilní. Pokud delta=-2sqrt(2), lambda_+/-^((+/-1,0))=sqrt(2), který má kladné reálné kořeny, takže bod je nestabilní. Pokud delta-2sqrt(2), pak |delta|sqrt(delta^2-8), takže oba kořeny jsou kladné a bod je nestabilní.

DuffingOscillatorPhasePortrait

zajímavé je, že zvláštní případ delta=0 bez vynucení,

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

může být integrován do kvadratury. Rozlišování (19) a připojování (20) dává

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

vynásobením obou stran  x^. dává

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

Ale to může být napsáno

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

takže máme invariantní pohybu h,

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

Řešení pro x^.^2 dává

 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)/(částečně)
(28)
 (partialh)/(partialx)=-x+x^3=y^.,
(29)

takže Duffing rovnice oscilátoru jsou dány Hamiltoniansystem

You might also like

Napsat komentář

Vaše e-mailová adresa nebude zveřejněna.