我有一个8x8矩阵,例如A=rand(8,8)。我需要做的是在对角线上子集所有2x2矩阵。这意味着我需要保存矩阵A(1:2,1:2),A(3:4,3:4),A(5:6,5:6),A(7:8,7:8)。rand(8,8); for i = 1:4 end
对于小型AA矩阵和小型AA子矩阵来说,这很好,但是,使用像上面这样的for循环(在不可行的情况下)可以
我有一个漂亮的三角n x n矩阵,比如L (对于下三角),我想要解一个类似的系统其中B和X是n x k矩阵(即:我想要解一个多个右手边的三角形线性系统)。此外,我有我的三角矩阵存储的包装格式,即我只存储下三角部分。我正在使用BLAS和LAPACK,但我已经意识到没有解决我的问题的具体方法。虽然有许多功能可以解决类似的问题::采用密集格式的三角<
我将尽量不用我的解释使事情变得太复杂,但我对如何最好地用从另一个包中获得的现有相关值填充没有重复值的三角相关矩阵感到困惑。这涉及到从文本文件列表中提取特定值。M[ct,] <- as.numeric(ht) #input this value into the appropriate location}
这显然没有给出我想要的三角化输出-我知道在矩阵中输入变量'ht‘存在错误,但不确定如何解决这个问题。理想情况下,文件1_1的相关值应该插入行1,列1,文件1_2应该插入行2