Я хотел бы получить трехмерное поле потока системы дифференциальных уравнений в R. Система находится в коде.
Я видел, что phaseR поддерживает только одно и два измерения.
Lorenz<-function(t, state, parameters) {
with(as.list(c(state, parameters)),{
# rate of change
dX <- a*X + Y*Z
dY <- b * (Y-Z)
dZ <- -X*Y + c*Y - Z
# return the rate of change+ list(c(dX, dY, dZ))
}) # end with(as.list ...
}
Я ожидаю, что поле потока будет похоже на поле, полученное функцией flowfield()
из пакета 'phaseR' или phaseArrows()
из 'pplane', но в 3D, поскольку в системах есть 3 переменные.