如何用Numpy创建100×100无环二元非对称共引网络矩阵？

你的代码存在几个明显问题：

• 变量名不统一：初始化的矩阵变量是array，但循环中使用了未定义的A，会触发NameError。
• 循环语法错误：for i,j in range(0,n)写法无效，range()返回单个整数序列，无法解构为两个变量。
• 赋值操作错误：使用比较运算符==而非赋值运算符=，导致无法修改矩阵元素。
• 未处理自环：共引网络通常不允许自环（对角线元素为1），原代码未对此做限制。

针对你的需求，这里提供几种实现方式：

方式一：嵌套循环实现（直观易懂） #

import numpy as np

n = 100
# 生成初始随机二元矩阵，先将对角线置0（消除自环）
array = np.random.randint(0, 2, size=(n, n))
np.fill_diagonal(array, 0)

# 遍历上三角区域，若当前位置为1，则将对称下三角位置置0
for i in range(n):
    for j in range(i + 1, n):
        if array[i][j] == 1:
            array[j][i] = 0

方式二：Numpy向量化操作（高效，适合大矩阵） #

利用Numpy的三角矩阵函数直接生成，避免循环：

import numpy as np

n = 100
# 生成上三角随机矩阵，再将对角线置0
upper_tri = np.triu(np.random.randint(0, 2, size=(n, n)))
np.fill_diagonal(upper_tri, 0)
# 下三角自动为0，直接得到符合要求的矩阵
array = upper_tri

方式三：掩码法（灵活控制随机区域） #

通过掩码精准控制保留随机值的区域：

import numpy as np

n = 100
# 创建掩码：仅上三角区域（不含对角线）为True
mask = np.triu(np.ones((n, n), dtype=bool), k=1)
# 初始化全0矩阵，仅在掩码区域填充随机0/1
array = np.zeros((n, n), dtype=int)
array[mask] = np.random.randint(0, 2, size=mask.sum())

以上代码生成的矩阵均满足：

• 对角线全为0（无自环）
• 任意i≠j时，array[i][j]与array[j][i]不同时为1（保证不对称，无双向环）

内容的提问来源于stack exchange，提问作者yang zhang