4

我有一个上三角矩阵（没有对角线），由下式给出：

M = [0 3 2 2 0 0; 0 0 8 6 3 2; 0 0 0 3 2 1; 0 0 0 0 2 1; 0 0 0 0 0 0]

生成的矩阵应如下所示：

R = [0 0 0 0 0 0; 0 2 0 0 0 0; 2 3 1 0 0 0; 2 6 2 1 0 0; 3 8 3 2 0 0]

由于我找不到一个简单的解释来描述我的目标，我试图用一张图片来形象化它：

在此处输入图像描述

我已经尝试了很多不同的组合rot90，transpose等等flipud，但是我找不到合适的转换来给我矩阵R

编辑：

矩阵的行M并不总是像上面的示例那样排序。对于另一个矩阵M_2：

M_2 = [0 2 3 1 0 0; 0 0 3 6 3 9; 0 0 0 1 2 4; 0 0 0 0 2 6; 0 0 0 0 0 0]

生成的矩阵R_2需要如下：

R_2 = [0 0 0 0 0 0; 0 9 0 0 0 0; 1 3 4 0 0 0; 3 6 2 6 0 0; 2 3 1 2 0 0]

再次显示下面的可视化：

在此处输入图像描述

4

3 回答 3

5

编辑： 受@Dan 评论提示的启发，可以进一步简化为

R = reshape(rot90(M), size(M));

原答案：

这应该是一个简单的方法来做到这一点

F = rot90(M);
R = F(reshape(1:numel(M), size(M)))

返回

R =
     0     0     0     0     0     0
     0     2     0     0     0     0
     2     3     1     0     0     0
     2     6     2     1     0     0
     3     8     3     2     0     0

这个想法是，当你旋转矩阵时，你会得到

>> F = rot90(M)
F =
     0     2     1     1     0
     0     3     2     2     0
     2     6     3     0     0
     2     8     0     0     0
     3     0     0     0     0
     0     0     0     0     0

这是一个 6 x 5 矩阵。如果您考虑 F 上的线性索引，则相应的索引是

>> reshape(1:30, size(F))
     1     7    13    19    25
     2     8    14    20    26
     3     9    15    21    27
     4    10    16    22    28
     5    11    17    23    29
     6    12    18    24    30

其中元素 6, 11, 12, 16, 17, 18 和 ... 现在为零，如果您将其重塑为 5 x 6 矩阵，您将得到

>> reshape(1:30, size(M))
     1     6    11    16    21    26
     2     7    12    17    22    27
     3     8    13    18    23    28
     4     9    14    19    24    29
     5    10    15    20    25    30

现在那些对应于零值的元素在最上面，正是我们想要的。因此，通过将此索引数组传递给F我们得到所需的R.

于 2013-08-20T11:20:05.117 回答
4

不依赖顺序（只需旋转彩色条并将它们推到底部）。

M(1,3)第一个解决方案：请注意，如果“数据”值之间存在零（例如，如果在给出的示例中） ，则它不起作用0。如果可能有零，请参阅下面的第二种解决方案：

[nRows nCols]= size(M);
R = [flipud(M(:,2:nCols).') zeros(nRows,1)];
[~, rowSubIndex] = sort(~~R);
index = sub2ind([nRows nCols],rowSubIndex,repmat(1:nCols,nRows,1));
R = R(index);

第二种解决方案：即使数据中有零也有效：

[nRows nCols]= size(M);
S = [flipud(M(:,2:nCols).') zeros(nRows,1)];
mask = 1 + fliplr(tril(NaN*ones(nRows, nCols)));
S = S .* mask;
[~, rowSubIndex] = sort(~isnan(S));
index = sub2ind([nRows nCols],rowSubIndex,repmat(1:nCols,nRows,1));
R = S(index);
R(isnan(R)) = 0;
于 2013-08-20T10:39:26.700 回答
1

备用选项，使用循环：

[nRows nCols]= size(M);
R = zeros(nRows,nCols);
for n = 1:nRows
  R((n+1):nCols,n)=fliplr(M(n,(n+1):nCols));
end
于 2013-08-20T11:05:54.667 回答