我想在CVXPY中生成(对角线)块矩阵(最好是稀疏的)。有些块可以是eye(m)或任何东西,但是我有一个块,它是:import cvxpy as cvximport scipyW = cvx.Variable(m,1)然后,我尝试用W_diag作为块来形成块对角线矩阵,例如:
T = scipy.sparse
因此,我正在编写一段代码,要求我使用大小分别为n×n的较小的‘正方形’矩阵J和M构造一个大型矩阵M,并重复如下:即M的维数使得M沿对角线重复'L‘次,J沿上第二对角线重复,J在下第二对角线上重复请注意,我正在研究Julia v 1.0.0,据我所知,在Julia中没有直接分配块矩阵的方法,这与Mathematica不同。我试图使用Kronecker产品来解决我的问题:
=Diagonal(ones(L)) #IDENTITY matrix of