如何将块转换为块对角矩阵（NumPy）
我在NumPy中有三个相同大小的正方形矩阵。我想将它们组合成块对角矩阵。
例：
a1 = np.array([[1,1,1],[1,1,1],[1,1,1]])
a2 = np.array([[2,2,2],[2,2,2],[2,2,2]])
a3 = np.array([[3,3,3],[3,3,3],[3,3,3]])
r = np.array([[1,1,1,0,0,0,0,0,0],[1,1,1,0,0,0,0,0,0],[1,1,1,0,0,0,0,0,0],[0,0,0,2,2,2,0,0,0],[0,0,0,2,2,2,0,0,0],[0,0,0,2,2,2,0,0,0],[0,0,0,0,0,0,3,3,3],[0,0,0,0,0,0,3,3,3],[0,0,0,0,0,0,3,3,3]])
做这个的最好方式是什么？
问题答案：
scipy.linalg具有block_diag函数以自动执行此操作
>>> a1 = np.array([[1,1,1],[1,1,1],[1,1,1]])
>>> a2 = np.array([[2,2,2],[2,2,2],[2,2,2]])
>>> a3 = np.array([[3,3,3],[3,3,3],[3,3,3]])
>>> import scipy.linalg
>>> scipy.linalg.block_diag(a1, a2, a3)
array([[1, 1, 1, 0, 0, 0, 0, 0, 0],
[1, 1, 1, 0, 0, 0, 0, 0, 0],
[1, 1, 1, 0, 0, 0, 0, 0, 0],
[0, 0, 0, 2, 2, 2, 0, 0, 0],
[0, 0, 0, 2, 2, 2, 0, 0, 0],
[0, 0, 0, 2, 2, 2, 0, 0, 0],
[0, 0, 0, 0, 0, 0, 3, 3, 3],
[0, 0, 0, 0, 0, 0, 3, 3, 3],
[0, 0, 0, 0, 0, 0, 3, 3, 3]])
>>> r = np.array([[1,1,1,0,0,0,0,0,0],[1,1,1,0,0,0,0,0,0],[1,1,1,0,0,0,0,0,0], [0,0,0,2,2,2,0,0,0],[0,0,0,2,2,2,0,0,0],[0,0,0,2,2,2,0,0,0],[0,0,0,0,0,0,3,3,3],[0,0,0,0,0,0,3,3,3],[0,0,0,0,0,0,3,3,3]])
>>> (scipy.linalg.block_diag(a1, a2, a3) == r).all()
True
-
我试图乘以两个块对称矩阵（矩阵大小矩阵大小）。我想执行块矩阵乘法（将一个矩阵分成多个块大小矩阵，并将相应的块相乘）。我已经写了一些代码，但想改进它，并存储主对角线以上的块，但我没有任何想法。如果可能的话，你们能帮忙吗？
-
我有一个编码器BCH的输出矩阵（3,63），但这个矩阵是伽罗瓦域，我需要将这个伽罗瓦域转换为矩阵二进制，因为matlab将伽罗瓦域中的元素视为字符串，我需要将这些值视为二进制数。 我需要将代码列与000010进行比较，。。。对于开关情况或if，但代码矩阵行是伽罗瓦场格式。我的问题是，遵循matlab错误是开关表达式必须是标量或字符向量。
-
我正在尝试将时间序列对象转换为矩阵。 此代码生成错误如下： 但是 出生时间序列 <代码>数据。矩阵（birthstimeseries）。它生成一维数据。如何将此数据帧的行数和列数保留到矩阵中。 我想要这样的结果作为一个矩阵。
-
这得到了我想要的，但可能没有很好地扩展？ 产量
-
特殊矩阵——三对角矩阵(Tridiagonal Matrix) 注：压缩存储的矩阵可以分为特殊矩阵和稀疏矩阵。对于那些具有相同元素或零元素在矩阵中分布具有一定规律的矩阵，被称之为特殊矩阵。对于那些零元素数据远远多于非零元素数目，并且非零元素的分布没有规律的矩阵称之为稀疏矩阵。 1. 三对角矩阵的概念 三对角矩阵就是对角线、邻近对角线的上下次对角线上有元素，其他位置均为0的矩阵。 三对角矩阵是一种特
-
我目前正在将原始矩阵转换为光栅以使用焦距函数，然后我想将光栅转换回矩阵。但是，当我尝试将光栅函数用作时，出现了一条错误消息。矩阵（）。即使有这个非常简单的例子： 以下是我得到的： 如果（！is.null（names（x）））列表（names（x），，则数组（x，c（length（x），1L）中出错： “dimnames”[1]的长度不等于数组范围 我正在使用RstuIO、R版本3.4.0和、和库。