matlab、対角行列の作り方

質問

matlabで、次のようにすると、
--------------------------------------------------------------
>> N = 5;
h = 2*pi/N;
column = [0 .5*(-1).^(1:N-1).*cot((1:N-1)*h/2)]';
D = toeplitz(column,column([1 N:-1:2]));
>> D

D =

0 -0.6882 0.1625 0.1625 -0.6882
-0.6882 0 -0.6882 0.1625 0.1625
0.1625 -0.6882 0 -0.6882 0.1625
0.1625 0.1625 -0.6882 0 -0.6882
-0.6882 0.1625 0.1625 -0.6882 0
--------------------------------------------------------------
5*5のDという行列が出来ます。

今回、
X =

D 0 0 0
0 D 0 0
0 0 D 0
0 0 0 D

のように20*20の対角行列を作りたいのですが、作り方を教えて頂けないでしょうか？

回答No.1

Kules · ベストアンサー

blkdiagを使えばできそうですね。
http://www.mathworks.co.jp/help/ja_JP/techdoc/ref/blkdiag.html

参考になれば幸いです。