在做数学运算的时候经常碰到这样的需求:
给定矩阵A,希望得到一个对角线方向上增广的矩阵
B =
| A 0 0 |
| 0 A 0 |
| 0 0 A …|
查了很久,没有看到numpy原生的函数可以实现,只能自己写。
朴素的方法就是写循环,不停stack加零,然后再赋值。不过这样做就少了一些python代码的优雅。
一种可能的写法:
def diag_stack(A : np.ndarray, n : int=2):
return np.tile(A, (n,n)) * \
np.repeat(np.repeat(np.eye(n), A.shape[0], axis=0), A.shape[1], axis=1).astype(int)
本质上其实是np.tile和np.repeat的拼接
测试一下:
a = np.array([[1,2,3],
[4,5,6]])
print(diag_stack(a, 3))
输出:
[[1 2 3 0 0 0 0 0 0]
[4 5 6 0 0 0 0 0 0]
[0 0 0 1 2 3 0 0 0]
[0 0 0 4 5 6 0 0 0]
[0 0 0 0 0 0 1 2 3]
[0 0 0 0 0 0 4 5 6]]
优点是代码简单。缺点是算力消耗不一定是最优。
如有更简便的方法,欢迎分享。