草庐IT

tridiagonal

全部标签

使用scipy.sparse的tridiagonal块矩阵

我需要使用ScipySparse构建块Tridiagonal矩阵来帮助。我的意思是平方矩阵B,我需要创建[[BI000][IBI00][0IBI0][00IBI][000IB]]现在,我希望通过编程性完成此操作,因为矩阵的大小可能会有所不同。谢谢!看答案解决了!我只是将scipy.sparse.bmat与列表综合一起使用。A=sparse.bmat([[Bifi==jelsenp.eye(n)ifabs(i-j)==1elseNoneforiinrange(n)]forjinrange(n)],format='bsr')在哪里B是一个nxn矩阵。