python-使用水平堆叠和垂直堆叠生成以下块矩阵的问题

我正在尝试生成以下由子矩阵和 组成的块矩阵A，B并且N是一个正整数。到目前为止，我的代码如下：

C_lower = B

for j in range(0,N):

    for i in range(0,N-j):

        col = np.linalg.matrix_power(A,i) @ B

        C = np.hstack(np.vstack((C_lower,col)))

但是，我的代码似乎无法正常工作，因为循环会永远持续下去。有什么建议么？

同样，我在构建以下块对角矩阵时也遇到了问题。

我尝试使用block_diagfrom scipy，但我无法重复Q多次N等于等于（即，在我的情况下为 N = 50）。我必须这样做block_diag(Q,Q,Q,Q,Q,Q,Q.......)才能得到我想要的块对角矩阵。