How can I make a matrix with diagonals composed of matrices?

See attached image. A is a 2x2 and B is a 2x1 (although it would be nice to generalize).

Any suggestions on how to easily create this matrix? I know I could use a for loop, but is there some easier way using matlab functions?

Answer 1

[Matt J]

I don't think you can or should avoid a for-loop, but there are ways to use looping more economically,

    C=cell(1,N+1);
    C{1}=zeros(size(B));
    C{2}=B;

    for i=3:N+1

        C{i}=A*C{i-1};

    end

    idx=toeplitz(1:N+1,ones(1,N+1));

    result=cell2mat( C(idx) )

Answer 2

[Star Strider]

The blkdiag function is your friend here. It was intended for just this purpose.