ДУ (Фазовая траектория)

Форум пользователей пакета Mathematica

Модератор: Admin

denverus
Сообщения: 5
Зарегистрирован: Вс ноя 07, 2010 12:35 pm

ДУ (Фазовая траектория)

Сообщение denverus » Вс ноя 07, 2010 12:38 pm

Имеется система уравнений.

x'=2x+y^2-1
y'=sin(x)-y^2+1

Нужно найти фазовую траекторию, построить и исследовать её.

В интернете находил информацию, но у самого не получается это сделать :cry: .

denverus
Сообщения: 5
Зарегистрирован: Вс ноя 07, 2010 12:35 pm

ДУ

Сообщение denverus » Вт ноя 09, 2010 2:30 pm

Ошибки(( почему не могу понять.
На примере похожий пример работал.

Gordon
Сообщения: 300
Зарегистрирован: Вт сен 08, 2009 8:16 pm

Сообщение Gordon » Вт ноя 09, 2010 2:48 pm

Код: Выделить всё

S0 = 1;
sol = NDSolve[{2 x[t] + y[t]^2 - 1 == x'[t],
   Sin[x[t]] - y[t]^2 + 1 == y'[t], x[0] == S0, y[0] == S0}, {x[t],
   y[t]}, {t, 0, 2}]
ParametricPlot[Evaluate[{x[t], y[t]} /. sol[[1]]], {t, 0, 2},
 AspectRatio -> 1]

Gordon
Сообщения: 300
Зарегистрирован: Вт сен 08, 2009 8:16 pm

Сообщение Gordon » Вт ноя 09, 2010 4:08 pm

Ваш файл почти рабочий, ошибка в еу была.
При отрицательных x0, y0 ~ < -0.7 система не решается

denverus
Сообщения: 5
Зарегистрирован: Вс ноя 07, 2010 12:35 pm

ДУ

Сообщение denverus » Вт ноя 09, 2010 4:18 pm

Получилось спасибо.

А почему не работает векторное поле?

VectorPlot[{ex[x,y],ey[x,y]},{x,0,2},{y,0,2},VectorPoints->15]

P.S. Пример брал с Mathematica 5, реализация идет на Mathematica 7

denverus
Сообщения: 5
Зарегистрирован: Вс ноя 07, 2010 12:35 pm

ДУ

Сообщение denverus » Чт ноя 11, 2010 11:20 am

Без Sin пашет как час,с Sin выдаёт ошибку((

denverus
Сообщения: 5
Зарегистрирован: Вс ноя 07, 2010 12:35 pm

Сообщение denverus » Вт ноя 23, 2010 11:13 pm

Спасибо огромное за помощь !!!