从列中分块对角矩阵

Question

假设我有一个 m x n 矩阵 A 。
有没有办法创建B，一个(n x m) x n 矩阵，其“对角线”由A 的列组成？

例子：

A = [1 2;
     3 4]  

B = [1 0;
     3 0;
     0 2;
     0 4]

Answer 1

这是一种方法：

1. 使用mat2cell 将A 转换为其列的cell array；
2. 从该元胞数组生成comma-separated list，并将其用作blkdiag 的输入。

代码：

A = [1 2; 3 4];                                   %// example data
C = mat2cell(A, size(A,1), ones(1,size(A,2)));    %// step 1
B = blkdiag(C{:});                                %// step 2

这会产生

B =
     1     0
     3     0
     0     2
     0     4

Answer 2

这是一个完成此任务的简短脚本。它适用于 A 的任何维度。

A=[1 2; 3 4];
[R C] = size(A);

for i=1:C
    B( 1+R*(i-1) : R*i , i ) = A(:,i);
end