Внедрение системы Пекоры-Кэрролла Лоренца в Python?
У меня есть следующая система, которая представляет собой две системы Лоренца, соединенные вместе:
def coupled_dynamics_Lorenz(x,t,sigma=10.0, beta=8.0/3.0, rho=28.0, k=0.0):
xdot = np.zeros(6)
xdot[0] = sigma*(x[1] - x[0]) + k*(x[3] - x[0])
xdot[1] = x[0]*(rho - x[2]) - x[1] + k*(x[4] - x[1])
xdot[2] = x[0]*x[1] - beta*x[2] + k*(x[5] - x[2])
xdot[3] = sigma*(x[4] - x[3]) + k*(x[0] - x[3])
xdot[4] = x[3]*(rho - x[5]) - x[4] + k*(x[1] - x[4])
xdot[5] = x[3]*x[4] - beta*x[5] + k*(x[2] - x[5])
return xdot
Мне нужно реализовать настройку синхронизации реакции привода Пекора-Кэрролла с использованием этой системы тремя способами: x, y и z. Уравнения Пекоры-Кэрролла для системы Лоренца (x-движение):
def master_slave(x, t, simga=10.0, beta=8.0/3.0, rho=28.0, k=0.0):
xdot = np.zeros(5)
xdot[0] = sigma*(x[1] - x[0])
xdot[1] = x[0]*(rho - x[2]) -x[1]
xdot[2] = x[0]*x[1] - beta*x[2]
xdot[3] = x[0]*(rho - x[4]) - x[3]
xdot[4] = x[0]*x[3] - beta*x[4]
return xdot
Однако мне, очевидно, нужно соединить их, используя постоянное значение k для силы связи.
Как я могу это сделать?