Решение системы нелинейных уравнений с питоном

sympy ответ @ Russ, это может быть легко достигнуто в sympy . Например:

 In [1]: import sympy as sp In [2]: x, y, z = sp.symbols('x, y, z') In [3]: rho, sigma, beta = sp.symbols('rho, sigma, beta') In [4]: f1 = sigma * (y - x) In [5]: f2 = x * (rho - z) - y In [6]: f3 = x * y - beta * z In [7]: sp.solvers.solve((f1, f2, f3), (x, y, z)) Out[7]: [(0, 0, 0), (-sqrt(beta*rho - beta), -sqrt(beta*(rho - 1)), rho - 1), (sqrt(beta*rho - beta), sqrt(beta*(rho - 1)), rho - 1)] 

где выходной формат – 3 возможных набора возможных значений для (x, y, z) .

Предупреждение. Я разработчик мудреца, поэтому я не могу быть нейтральным.

Я не знаю, как это сделать на чистом Python, но я бы порекомендовал систему Sage , интерфейс которой находится в Python (на самом деле это командная строка – специально сконфигурированный IPython) и позволяет это делать:

 +--------------------------------------------------------------------+ | Sage Version 5.10, Release Date: 2013-06-17 | | Type "notebook()" for the browser-based notebook interface. | | Type "help()" for help. | +--------------------------------------------------------------------+ sage: var("sigma yx rho beta z") (sigma, y, x, rho, beta, z) sage: sys = [sigma*(yx), x*(rho-z)-y, x*y-beta*z] sage: solve(sys, x, y, z) [[x == sqrt(beta*rho - beta), y == (beta*rho - beta)/(sqrt(rho - 1)*sqrt(beta)), z == rho - 1], [x == -sqrt(beta*rho - beta), y == -(beta*rho - beta)/(sqrt(rho - 1)*sqrt(beta)), z == rho - 1], [x == 0, y == 0, z == 0]] 

Обычно это проще использовать следующим образом:

 sage: solve(sys, x, y, z, solution_dict=True) [{z: rho - 1, x: sqrt(beta*rho - beta), y: (beta*rho - beta)/(sqrt(rho - 1)*sqrt(beta))}, {z: rho - 1, x: -sqrt(beta*rho - beta), y: -(beta*rho - beta)/(sqrt(rho - 1)*sqrt(beta))}, {z: 0, x: 0, y: 0}] 

Главным недостатком является то, что Sage – это полное распространение математического программного обеспечения, которое отправляет собственный интерпретатор Python (вместе с огромным множеством других вещей, написанных на многих языках, включая C / C ++, Cython, lisp, fortran) и, как правило, трудно установить, если вы хотите использовать собственный интерпретатор.

Хорошей новостью для вашей проблемы является то, что Scipy уже поставляется с шалфеем.

SymPy может помочь; Я не знаю, насколько это хорошо при решении нелинейных уравнений: http://scipy-lectures.github.io/advanced/sympy.html#id23

Вы должны иметь возможность выполнять код, как (следуя приведенным выше примерам):

 from sympy import * x = Symbol('x') y = Symbol('y') z = Symbol('z') beta = Symbol('beta') rho = Symbol('rho') sigma = Symbol('sigma') solve([sigma*(yx), x*(rho-z)-y, x*y-beta*z], [x, y, z]) 

Я не тестировал, если он работает (у меня нет этой возможности для этой машины).