3 графика на одном.

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

Модератор: Admin

zharenkov
Сообщения: 2
Зарегистрирован: Вс апр 01, 2012 4:17 pm

3 графика на одном.

Сообщение zharenkov » Вс апр 01, 2012 4:26 pm

Добрый день. Задача в следующем. Есть 3 графика, 1 - седловая поверхность, строю её Plot3D. И 2 графика - решения систем ДУ из трех уравнений. Следовательно строю их ParametricPlot3D.
Теперь задача - совместить это всё на одном графике, чтоб видно было траектории систем относительно плоскости.

Для наглядности приложу код

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

Plot3D[-y^2 + x^2, {x, -2, 2}, {y, -2, 2},AxesLabel->{x,y,z}]

w = NDSolve[{x'[t] == y[t],
   y'[t] == z[t],
   z'[t] == -100(z[t] + y[t]^2 - x[t]^2), x[0] == y[0] == z[0] == 1}, {x, y, z}, {t, 5}]
ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. w], {t, 0, 0.5},AxesLabel->{x,y,z}]
   

s = NDSolve[{x'[t] == y[t],
   y'[t] == z[t],
   z'[t] == -100(z[t] + y[t]^2 - x[t]^2)-x[t]y[t]+y[t]z[t], x[0] == y[0] == z[0] == 1}, {x, y, z}, {t, 5}]
ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /.s], {t, 0, 0.5},AxesLabel->{x,y,z}]

Буду признателен за ответы.

VSI
Сообщения: 305
Зарегистрирован: Вт мар 15, 2005 6:40 pm
Откуда: Мариуполь

Re: 3 графика на одном.

Сообщение VSI » Пн апр 02, 2012 7:55 am

zharenkov писал(а):Добрый день. Задача в следующем. Есть 3 графика, 1 - седловая поверхность, строю её Plot3D. И 2 графика - решения систем ДУ из трех уравнений. Следовательно строю их ParametricPlot3D.
Теперь задача - совместить это всё на одном графике, чтоб видно было траектории систем относительно плоскости.

Для наглядности приложу код

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

Plot3D[-y^2 + x^2, {x, -2, 2}, {y, -2, 2},AxesLabel->{x,y,z}]

w = NDSolve[{x'[t] == y[t],
   y'[t] == z[t],
   z'[t] == -100(z[t] + y[t]^2 - x[t]^2), x[0] == y[0] == z[0] == 1}, {x, y, z}, {t, 5}]
ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. w], {t, 0, 0.5},AxesLabel->{x,y,z}]
   

s = NDSolve[{x'[t] == y[t],
   y'[t] == z[t],
   z'[t] == -100(z[t] + y[t]^2 - x[t]^2)-x[t]y[t]+y[t]z[t], x[0] == y[0] == z[0] == 1}, {x, y, z}, {t, 5}]
ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /.s], {t, 0, 0.5},AxesLabel->{x,y,z}]

Буду признателен за ответы.

A := Plot3D[-y^2 + x^2, {x, -2, 2}, {y, -2, 2},
AxesLabel -> {x, y, z}]
w = NDSolve[{x'[t] == y[t], y'[t] == z[t],
z'[t] == -100 (z[t] + y[t]^2 - x[t]^2),
x[0] == y[0] == z[0] == 1}, {x, y, z}, {t, 5}]
W := ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. w], {t, 0, 0.5},
AxesLabel -> {x, y, z}, PlotStyle -> {Red}]
s = NDSolve[{x'[t] == y[t], y'[t] == z[t],
z'[t] == -100 (z[t] + y[t]^2 - x[t]^2) - x[t] y[t] + y[t] z[t],
x[0] == y[0] == z[0] == 1}, {x, y, z}, {t, 5}]
S := ParametricPlot3D[Evaluate[{x[t], y[t], z[t]} /. s], {t, 0, 0.5},
AxesLabel -> {x, y, z}, PlotStyle -> {Blue}]
Show[A, W, S]
Остается подобрать приемлемый масштаб...

zharenkov
Сообщения: 2
Зарегистрирован: Вс апр 01, 2012 4:17 pm

Сообщение zharenkov » Пн апр 02, 2012 10:12 am

спасибо.