将列表变成对角矩阵

Question

由于数据矩阵的 SVD，我有一个奇异值列表。 Python 输出为列表而不是对角矩阵。如果不将奇异值转换为对角矩阵，则无法组合矩阵以找到回归系数。

简单地说，这就是我所拥有的：

U, W, Vt = np.linalg.svd(Absorbance_Data)
print U.shape, W.shape, Vt.shape

输出：

(71, 71), (71,), (512, 512)

为了计算 PCR 的回归系数，我需要让矩阵 W 成为 (71, 512) 对角矩阵（允许零值）。

感谢您的帮助。

Answer 1

您首先将W 对角化，然后附加零：

W = np.diag(W)
W = np.append(W, np.zeros(71, 441), axis=1)

print(W.shape)
=> (71, 512)

这会将一个 71x441 的零数组附加到 W 的 71x71 对角化。 axis=1 按列附加零。