我想使用一个列表创建一个矩阵,该列表的元素将是对角线下矩阵的元素。

import numpy as np
x1 = np.array([0.9375, 0.75, 0.4375, 0.0, 0.9375, 0.75, 0.4375, 0.9375, 0.75, 0.9375])
x1

我想要的矩阵是
array([[ 1.    ,  0.9375,  0.75  ,  0.4375,  0.    ],
   [ 0.9375,  1.    ,  0.9375,  0.75  ,  0.4375],
   [ 0.75  ,  0.9375,  1.    ,  0.9375,  0.75  ],
   [ 0.4375,  0.75  ,  0.9375,  1.    ,  0.9375],
   [ 0.    ,  0.4375,  0.75  ,  0.9375,  1.    ]])

我以为你可以用np.tril做这个,但它给出了一个我不期望的结果。
mat = np.tril(x1, k = -1  )
print(mat)

我错过了什么?
如果这是一个微不足道的问题,我会事先道歉,但如果没有循环,我就不知道该怎么做。

最佳答案

你可以:

x = np.ones((5, 5), dtype=float)
x[np.triu_indices(5, 1)] = x1        # sets the upper triangle
x[np.triu_indices(5, 1)[::-1]] = x1  # sets the lower triangle

在最后一行中,索引被反转,因为x1是为上三角形排序的。如果感觉更直观,也可以使用x[np.tril_indices(5, -1)] = x1[::-1]

09-27 06:47