米脂
错误的问题陈述看来你没有解决正确的问题。正如您所说，它不能像您使用的方法一样工作：计算已确定的满秩线性矩阵方程 ax = b 的“精确”解 x。如文档中所定义。将您的问题写为 是很诱人的Ax = b。但是，您需要一个 3x3 矩阵来测试向量b=[0,0,0]T，并找到一个非平凡的解决方案x来声明所有三个向量都不是线性独立的。但这意味着矩阵没有满秩。在您的代码中，以下表达式无法满足函数签名：np.linalg.solve(np.array([v1]), v2)并引发以下错误：LinAlgError: Last 2 dimensions of the array must be square就矩阵形状而言，有效的调用numpy.linalg.solve是：np.linalg.solve(np.array([v1, v2, v3]).T, np.zeros(v1.size))但在两个方面失败了：LinAlgError: Singular matrix首先因为矩阵不是满秩的，因此不可逆。其次，这并不能解决两个向量的问题。通用解决方案import numpy as np
from scipy import linalg相反，您想要解决有关线性独立性的更普遍的问题，可以通过评估的秩[v1, v2]T来解决：在线性代数中，矩阵 A 的秩是由其列生成（或跨越）的向量空间的维数。这对应于 A 的线性独立列的最大数量。对于 numpy，您可以使用numpy.linalg.matrix_rank来执行此操作：np.linalg.matrix_rank(np.array([v1, v2]).T) # Missing dimension: 1
np.linalg.matrix_rank(np.array([v1, v3]).T) # Full rank: 2此方法对执行SVD 分解计算出的空奇异值进行计数，该算法的实现由scipy.linalg.svd提供：U, s, V = linalg.svd(np.array([v1, v2]).T) # s = [11.18033989, 0.]
U, s, V = linalg.svd(np.array([v1, v3]).T) # s = [5.55992016, 2.84381571]或者对其执行高斯消元法（假设使用LU 分解）并检查结果。scipy 包提供了scipy.linalg.lu：p, l, u = linalg.lu(np.array([v1, v2]).T) # Null pivot: u = [[5., -10.], [0., 0.]]
p, l, u = linalg.lu(np.array([v1, v3]).T) # Full rank: u = [[5., 2.], [0., 3.]]所有这些方法都会评估矩阵列所跨越的向量空间的维数。如果矩阵是满秩的，那么你的向量是线性独立的。如果不是，则至少有两个相关向量。那么解决您的问题的一个简单方法是：def indep1(*args): A = np.array(args).T return np.linalg.matrix_rank(A) == len(args)indep1(v1, v2) # Falseindep1(v1, v3) # Trueindep1(v1, v2, v3) # Falseindep1(v1, np.zeros(v1.size)) # False特例如果您停留在 3D 向量空间中并且只需要检查两个向量，您还可以利用numpy.cross提供的叉积：np.cross(v1, v2) # Colinear: [0, 0, 0]
np.cross(v1, v3) # Not colinear: [15, 0, -5]更具体的替代方案是：def indep2(a, b): return not np.allclose(np.cross(a, b), 0.)indep2(v1, v2) # Falseindep2(v1, v3) # True# indep2(v1, v2, v3) # TypeErrorindep2(v1, np.zeros(v1.size)) # False它们是一个恒定的系统正如 所建议的@MadPhycist，我们还可以评估共线向量共享标量比例因子（例如 ）的事实a = k*b，那么测试它的非防弹方法是检查坐标之比是否恒定：def indep3(a, b): r = a/b # Ratio of coordinates q = np.isfinite(r) # Filter bad ratio (zero division) return not np.allclose(r[q], r[q][0]) # Assess all ratio are equalindep3(v1, v2) # Falseindep3(v1, v3) # True# indep3(v3, np.zeros(v1.size)) # IndexError如果它比以前的解决方案需要更少的计算并且不依赖于高水平的线性代数，那么它需要处理特定的情况，并且建议的实现是臭的（不要使用它）。格拉米安法可以安全地对两个向量实现公开的方法，@dmuir如下所示：def indep4(a, b): return not np.isclose(np.dot(a,b)*np.dot(b,a), np.dot(a,a)*np.dot(b,b))indep4(v1, v2) # Falseindep4(v1, v3) # Trueindep4(v3, v4) # Falseindep4(v3, np.zeros(v1.size)) # False