Duffing differenciálegyenlet

letöltés Mathematica Notebook

a Duffing egyenlet legáltalánosabb kényszerített formája

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

a választott paraméterektől függően az egyenlet számos speciális formát ölthet. Például csillapítás és kényszerítés nélkül  delta=gamma=0 és a pluszjelet figyelembe véve az egyenlet a következő lesz

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

(Bender és Orsag 1978, 547. o.; Zwillinger 1997, 122. o.). Ez az egyenlet kaotikus viselkedést mutathat.  béta0esetén az egyenlet “kemény rugót”,  béta0 esetén pedig “puha rugót” jelent.”Ha  béta0, a fázis portré görbék zárva vannak.

ha ehelyett vesszük  beta=1,  omega_0=1, állítsa vissza az órát úgy, hogy  phi=0, és használja a mínuszjelet, akkor az egyenlet

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

ez lehet írni, mint egy rendszer elsőrendű rendes differenciálegyenletek, mint

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

(Wiggins 1990, 5. o.), amely kényszerítetlen esetben a következőkre redukálódik

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

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

a kapcsolt differenciálegyenletek ezen halmazának rögzített pontjait a

 x^.= y=0,
(8)

tehát  y=0, és

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

ad  x=0,+/-1. A rögzített pontok tehát (-1,0), (0,0), és (1,0).

a rögzített pontok stabilitásának elemzése pont lehet az egyenletek linearizálásával. Differenciáló ad

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

amely lehet írni, mint a mátrix egyenlet

 =.
(14)

a pont stabilitásának vizsgálata (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)

de delta^2=0, tehát lambda_+/-^((0,0)) valódi. Mivel  sqrt(delta^2+4)|delta|, mindig lesz egy pozitív gyökér, tehát ez a fix pont instabil. Most nézd meg (+/-1, 0). A jellemző egyenlet

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

amelynek gyökerei vannak

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

mert  delta0,  R0, tehát a pont aszimptotikusan stabil. Ha  delta=0,  lambda_+/-^((+/-1,0))=+/-isqrt (2), tehát a pont lineárisan stabil (Wiggins 1990, 10. o.). Ha azonban delta in (- 2sqrt(2),0), a gyök egy képzeletbeli részt ad, a valós rész pedig  0, tehát a pont instabil. Ha delta= - 2sqrt(2), lambda_+/-^((+/-1,0))=sqrt (2), amelynek pozitív valós gyökere van, tehát a pont instabil. Ha  delta-2sqrt(2), akkor |delta|sqrt(delta^2-8), tehát mindkét gyökér pozitív és a pont instabil.

DuffingOscillatorPhasePortrait

érdekes, hogy a különleges eset  delta=0 nem kényszerítve,

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

kvadratúrákkal integrálható. Differenciáló (19) és duggingin (20) ad

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

szorozzuk meg mindkét oldalt x^ - val. ad

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

de ezt meg lehet írni

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

tehát van egy invariáns a mozgás h,

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

megoldása x^.^2 ad

 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)

tehát a Duffing oszcillátor egyenleteit a Hamilton adjarendszer

You might also like

Vélemény, hozzászólás?

Az e-mail-címet nem tesszük közzé.