如何用Numpy广播计算Rij=Aij*Bij/Cij？解决奇异矩阵报错及方法验证

核心问题：误解了公式的运算类型 #

你要计算的是逐元素算术运算 $R_{ij} = A_{ij} \times B_{ij} / C_{ij}$，但代码里误用了矩阵乘法（@），同时尝试用矩阵逆的方式完全偏离需求，还因矩阵奇异报错。

为什么np.linalg.inv(C1)报错？ #

看你定义的C矩阵：每一行j的元素是((i+1)/(j+1))/2000，相当于第j行是第0行的1/(j+1)倍——整个矩阵的秩为1，是奇异矩阵，不存在逆矩阵，所以调用inv会直接抛出线性代数错误。而且这个操作本身就不符合你要的逐元素除法逻辑。

你的D = B1 / C1[..., None]方法为什么不对？ #

这个操作会把C1从(500,500)扩展为(500,500,1)，和B1广播后得到(500,500,1)的D，再和A1做矩阵乘法@，最终得到的是(500,1)的结果，完全不是原公式要求的(500,500)逐元素计算结果，逻辑完全错误。

正确解法：直接逐元素运算 #

Numpy数组默认支持逐元素的算术操作，因为你的A、B、C都是同维度的500x500数组，直接用*和/就能完成需求：

import numpy as np

A = [[(i+j)/2000 for i in range(500)] for j in range(500)]
B = [[(i-j)/2000 for i in range(500)] for j in range(500)]
C = [[((i+1)/(j+1))/2000 for i in range(500)] for j in range(500)]

def numpy_task3_part_A(A,B,C):
    A1 = np.array(A)
    B1 = np.array(B)
    C1 = np.array(C)
    # 逐元素运算，完全匹配公式 R_ij = A_ij * B_ij / C_ij
    Rij = A1 * B1 / C1
    return Rij

# 验证结果形状（应为(500,500)）
result = numpy_task3_part_A(A,B,C)
print(result.shape)

补充说明 #

• 你的C矩阵中i+1和j+1都是1到500，不会出现除以0的情况，所以无需额外处理除数为0的问题。
• 如果后续遇到不同维度的数组需要逐元素运算，才需要用到广播规则，但这里维度完全一致，直接运算即可。

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