Я хочу решить линейное уравнение с тремя или более переменными. Есть ли хорошая библиотека в python для этого?
Есть ли модуль python для решения линейных уравнений?
Ответ 1
См. http://sympy.org/ и http://numpy.scipy.org/.
В частности, http://docs.scipy.org/doc/numpy/reference/routines.linalg.html
И http://docs.sympy.org/0.7.0/tutorial.html#algebra, http://docs.sympy.org/dev/modules/solvers/solvers.html
Изменить: добавлена ссылка solvers из комментария.
Ответ 2
Да, в очень популярном пакете NumPy есть функция для этого. Их пример:
Решаем систему уравнений
3 * x0 + x1 = 9
иx0 + 2 * x1 = 8
:>>> import numpy as np >>> a = np.array([[3,1], [1,2]]) >>> b = np.array([9,8]) >>> x = np.linalg.solve(a, b) >>> x array([ 2., 3.])
https://docs.scipy.org/doc/numpy-1.13.0/reference/generated/numpy.linalg.solve.html
Ответ 3
Вы можете использовать метод наименьших квадратов в python для решения системы уравнений например, для решения уравнений 3x + 4y = 7 и 5x + 6y = 8
>>> import numpy
>>> a=[[3,4],[5,6]]
>>> b=[7,8]
>>> numpy.linalg.lstsq(a,b)
(array([-5. , 5.5]), array([], dtype=float64), 2, array([ 9.27110906, 0.21572392]))
Ответ 4
Используя пример @Jeremy:
from sympy import *
x0, x1 = symbols(['x0', 'x1'])
sol = solve([3 * x0 + x1 - 9, x0 + 2 * x1 - 8], [x0, x1])
print(sol)
выход:
{x0: 2, x1: 3}
Используя пример @004 с немного другой нотацией:
from sympy import *
x, y = symbols(['x', 'y'])
system = [
Eq(3*x + 4*y, 7),
Eq(5*x + 6*y, 8)
]
soln = solve(system, [x, y])
print(soln)
{x: -5, y: 11/2}
Примечание: иногда можно увидеть следующие обозначения для символов: x, y = symbols('x, y')
, которые кажутся менее питоническими.
Ответ 5
Вы также можете добавить ограничения (верхняя и нижняя границы) на x с помощью lsq_linear: scipy.optimize.lsq_linear
Ответ 6
Вы можете написать простую функцию, которая решает систему линейных уравнений.
def solve(equations):
#the constants of a system of linear equations are stored in a list for each equation in the system
"""
for example the system below:
2x+9y-3z+7w+8=0
7x-2y+6z-1w-10=0
-8x-3y+2z+5w+4=0
0x+2y+z+w+0=0
is expressed as the list:
[[2,9,-3,7,8],[7,-2,6,-1,-10],[-8,-3,2,5,4],[0,2,1,1,0]]
"""
lists=[] # I failed to name it meaningfully
for eq in range(len(equations)):
#print "equations 1", equations
#find an equation whose first element is not zero and call it index
index=0
for i in range(len(equations)):
if equations[i][0]<>0:
index=i;
break;
#print "index "+str(eq)+": ",index
#for the equation[index] calc the lists next itam as follows
lists.append([-1.0*i/equations[index][0] for i in equations[index][1:]])
#print "list"+str(eq)+": ", lists[-1]
#remve equation[index] and modify the others
equations.pop(index)
for i in equations:
for j in range(len(lists[-1])):
i[j+1]+=i[0]*lists[-1][j]
i.pop(0)
lists.reverse()
answers=[lists[0][0]]
for i in range(1,len(lists)):
tmpans=lists[i][-1]
for j in range(len(lists[i])-1):
tmpans+=lists[i][j]*answers[-1-j]
answers.append(tmpans)
answers.reverse()
return answers