Я работаю с хаотическими аттракторами и проверяю некоторые непрерывные-> дискретные эквивалентности. Я произвел непрерывное моделирование системы Росслера таким образом
a = 0.432; b = 2; c = 4;
Rossler = {
x'[t] == -y[t] - z[t],
y'[t] == x[t] + a*y[t],
z'[t] == b + x[t]*z[t]-c*z[t]};
sol = NDSolve[
{Rossler, x[0] == y[0] == z[0] == 0.5},
{x, y, z}, {t,500}, MaxStepSize -> 0.001, MaxSteps -> Infinity]
. Теперь, когда я пытаюсь оценить дискретную эквивалентную систему с помощью RSolve, Mma ничего не делает, даже ошибки, просто не может ее решить.
RosslerDiscreto = {
x[n + 1] == x[n] - const1*(y[n] + z[n]),
y[n + 1] == 1 - a*const2)*y[n] + const2*x[n],
z[n + 1] == (z[n]*(1 - const3) + b*const3)/(1 - const3*x[n])}
Я хочу знать, есть ли для RSolve числовая функция, аналогичная NDSolve для DSolve. Я знаю, что могу произвести вычисления с некоторыми циклами For [], просто хочу знать, существует ли такая функция.