A..7 diag()

• diag(ベクトル) とすると、ベクトルを対角成分に埋め込んだ 正方行列を返す。

  >> diag(1:5)
  
  ans =
       1     0     0     0     0
       0     2     0     0     0
       0     0     3     0     0
       0     0     0     4     0
       0     0     0     0     5
  

• diag(行列) とすると、行列の対角成分からなるベクトルを返す。

  >> a=(1:9)'*(1:9);
  >> diag(a)
  
  ans =
       1
       4
       9
      16
      25
      36
      49
      64
      81
  

  ゆえに diag(diag(行列)) とすると、 行列の対角成分以外を0 クリアした対角行列が得られる。

  >> a=(1:9)'*(1:9);
  >> diag(diag(a))
  
  ans =
       1     0     0     0     0     0     0     0     0
       0     4     0     0     0     0     0     0     0
       0     0     9     0     0     0     0     0     0
       0     0     0    16     0     0     0     0     0
       0     0     0     0    25     0     0     0     0
       0     0     0     0     0    36     0     0     0
       0     0     0     0     0     0    49     0     0
       0     0     0     0     0     0     0    64     0
       0     0     0     0     0     0     0     0    81