如何计算80×80矩阵中各(i,j)位置的(ij/√(ii*jj))值？

实现思路 #

核心步骤是提取矩阵主对角线元素、构造对应位置的分母矩阵，最后对原矩阵做逐元素除法即可完成计算。

R语言代码实现 #

# 示例矩阵（使用你提供的5x5测试矩阵）
b <- structure(c(0.232, 0.838, 0.488, 0.490, 0.799, 
                 0.491, 0.424, 0.334, 0.399, 0.550, 
                 0.835, 0.784, 0.426, 0.643, 0.311,
                 0.520, 0.670, 0.902, 0.992, 0.479, 
                 0.737, 0.640, 0.547, 0.011, 0.166), 
               dim = c(5L, 5L))

# 提取矩阵主对角线元素
diag_b <- diag(b)

# 构造每个(i,j)位置的分母：√(b[i,i]×b[j,j])
denominator <- sqrt(outer(diag_b, diag_b))

# 计算最终结果矩阵
result <- b / denominator

# 打印结果（可选）
print(result)

代码说明 #

• diag(b)：提取矩阵b的主对角线元素，得到长度与矩阵行数一致的向量。
• outer(diag_b, diag_b)：计算对角线向量的外积，生成与b同维度的矩阵，每个(i,j)位置的值为b[i,i]×b[j,j]。
• sqrt(...)：对整个外积矩阵逐元素开平方，得到每个位置的分母。
• b / denominator：R中矩阵的/运算符默认是逐元素除法，直接得到每个(i,j)位置的目标计算值。

80x80矩阵适配说明 #

上述代码无需修改逻辑即可直接用于80行80列的矩阵，R的矩阵运算会自动处理高维度场景，只要你的矩阵b是合法的80x80结构即可。

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