Bistable

Z Skrypty dla studentów Ekonofizyki UPGOW

(Różnice między wersjami)
(Obligacja ze stałym kuponem)
(Oscylator bistabliny)
Linia 3: Linia 3:
<math>m \ddot x = -\frac{dU(x)}{dx}-\gamma x  ,</math>  
<math>m \ddot x = -\frac{dU(x)}{dx}-\gamma x  ,</math>  
 +
 +
równoważnie:
 +
 +
<math>m \dot x = v</math>
 +
<math>-\frac{dU(x)}{dx}-\gamma x  ,</math>
 +
 +
jako potencjal weżmy funkcję z dwoma minimami np.:
 +
 +
<math>x.^4-4*x.^2,</math>
 +
 +
<source lang="matlab">
 +
fplot(@(x) x.^4-4*x.^2,[-2.1,2.1],200)
 +
</source>
który możemy zaimplementować jako funkcję w matlabie:
który możemy zaimplementować jako funkcję w matlabie:

Wersja z 09:07, 27 paź 2010

Oscylator bistabliny

Rozważmy oscylator nieliniowy:

\(m \ddot x = -\frac{dU(x)}{dx}-\gamma x ,\)

równoważnie:

\(m \dot x = v\) \(-\frac{dU(x)}{dx}-\gamma x ,\)

jako potencjal weżmy funkcję z dwoma minimami np.:

\(x.^4-4*x.^2,\)

fplot(@(x) x.^4-4*x.^2,[-2.1,2.1],200)

który możemy zaimplementować jako funkcję w matlabie:

function dx = ODEbistable(X,T)
    global gama;
    dx = zeros(2,1);
    dx(1) = X(2);
    dx(2) = -gama*X(2)-4*X(1).^3+8*X(1);
    return
end
Baseny przyciągania w nieliniowym oscylatorze.