我想解决这个方程:
v0 = a1*v1 + a2*v2 + a3*v3 + a4*v4 + a5*v5
v0, v1, v2, v3, v4, v5
是已知列矩阵(向量),a1, a2, a3, a4, a5
是需要计算的数字。我想知道在 numpy,sympy
或 scipy
中是否有函数可以直接计算此方程,或者我如何解决此方程。
请给我一个链接或书面示例。
我想解决这个方程:
v0 = a1*v1 + a2*v2 + a3*v3 + a4*v4 + a5*v5
v0, v1, v2, v3, v4, v5
是已知列矩阵(向量),a1, a2, a3, a4, a5
是需要计算的数字。我想知道在 numpy,sympy
或 scipy
中是否有函数可以直接计算此方程,或者我如何解决此方程。
请给我一个链接或书面示例。
v0
的每个元素都是由数组v1,v2,v3,v4,v5
中对应元素的和表示的。a1,a2,a3,s4,s5
的数量等于方程的数量,也就是向量v1,v2,v3,v4,v5
的长度。from numpy import allclose,zeros_like
from numpy.random import rand
from numpy.linalg import solve
# generate the 5 vectors as random arrays
mat = rand(5,5)
v1,v2,v3,v4,v5 = mat.T
v0 = rand(5)
x= solve(mat,v0)
#first check
assert allclose(dot(mat,x),v0)
#2nd check, which is the equation of the OP
res = zeros_like(v0)
for xj,vj in zip([v1,v2,v3,v4,v5],x):
res+= xj*vj
assert allclose(res,v0)
numpy.linalg.solve
。 - Joe Kington