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
