如何生成满足特定条件的随机二进制可逆矩阵？

要生成满足你要求的n阶矩阵，我们可以通过构造特定结构的0-1矩阵来保证三个条件同时成立，以下是具体方法和实现：

核心思路：利用三角矩阵的特性 #

主对角线全为1的下三角/上三角0-1矩阵天生满足所有要求：

1. 元素仅为0或1；
2. 行列式等于主对角线元素的乘积（即1），必然可逆；
3. 逆矩阵的每行元素之和一定非负（推导见下文）。

这类矩阵的数量足够多：下三角版本有$2^{n(n-1)/2}$个不同矩阵，完全覆盖你的需求。

具体实现代码 #

生成下三角型符合条件矩阵 #

import numpy as np

def generate_lower_triangular_valid_matrix(n):
    # 生成下三角0-1矩阵，主对角线强制设为1
    A = np.tril(np.random.randint(0, 2, (n, n)), k=0)
    np.fill_diagonal(A, 1)
    # 验证（理论上可省略，下三角主对角线全1必然可逆）
    inv_A = np.linalg.inv(A)
    row_sums = inv_A.sum(axis=1)
    # 浮点误差容忍，确保行和非负
    assert np.all(row_sums >= -1e-9), "逆矩阵行和不符合要求"
    return A

生成上三角型符合条件矩阵 #

def generate_upper_triangular_valid_matrix(n):
    # 生成上三角0-1矩阵，主对角线强制设为1
    A = np.triu(np.random.randint(0, 2, (n, n)), k=0)
    np.fill_diagonal(A, 1)
    inv_A = np.linalg.inv(A)
    row_sums = inv_A.sum(axis=1)
    assert np.all(row_sums >= -1e-9), "逆矩阵行和不符合要求"
    return A

为什么这类矩阵满足条件？ #

以下三角矩阵为例：

• 可逆性：下三角矩阵的行列式等于主对角线元素的乘积，这里主对角线全为1，行列式为1，必然可逆。
• 逆矩阵行和非负：设逆矩阵第i行的和为$r_i$，由$A \cdot A^{-1} = I$可得：
  $$r_i + \sum_{k=1}^{i-1} A[i][k] \cdot r_k = 1$$
  用归纳法可证：
  • 第1行的$r_1=1 \geq 0$；
  • 假设前i-1行的$r_k \geq 0$，则$\sum_{k=1}^{i-1} A[i][k] \cdot r_k \geq 0$，因此$r_i = 1 - \text{非负数} \geq 0$。

上三角矩阵的推导逻辑完全一致。

扩展方案：对角占优0-1矩阵 #

如果你需要更多样的矩阵结构，可以生成主对角线全1、每行非主对角线1的个数≤1的不可约0-1矩阵：

• 这类矩阵是弱对角占优且不可约，因此可逆；
• 其逆矩阵的行和必然非负（可通过线性方程组$A\mathbf{x}=\mathbf{1}$存在非负解推导）。

示例代码（n≥3）：

def generate_diagonally_dominant_matrix(n):
    A = np.eye(n, dtype=int)
    # 每行随机选一个非主对角线位置设为1（保证不可约）
    for i in range(n):
        j = np.random.choice([x for x in range(n) if x != i])
        A[i][j] = 1
    inv_A = np.linalg.inv(A)
    row_sums = inv_A.sum(axis=1)
    assert np.all(row_sums >= -1e-9), "逆矩阵行和不符合要求"
    return A

内容的提问来源于stack exchange，提问作者Erel Segal-Halevi