nDiki : Octave

変数への代入

'=' を使う

 x = 5

代入で最後に ';' をつけると代入結果は表示されない。

変数の中身を表示

変数名

 x

変数一覧を表示

 who

変数のクリア

'clear' を使う

 clear x

行列

[, ] を使う。1行目から順に各要素を入力していき、行と行の間は ';' で区切る

 [1 2; 3 4]

行列の加算・減算

'+', '-'

  [1 2; 3 4] + [5 6; 7 8]

行列の積

'*'

 [1 2; 3 4] * [5 6; 7 8]

転置行列

行列の後ろに' をつける

 [1 2; 3 4]'

行列の固有値を求める

eig を使う。

 eig([3 1 1; 1 2 0; 1 0 2])

固有値分解

eig を使う

 m = [1 2; 3 4]
 [v, lambda] = eig(m)
 m' * m
 v * lambda * v'

固有値は、lambda の対角成分に値の大きい順に並ぶ(絶対値順ではない)。

特異値分解 (singular value decomposition)

svd を使う

 [u, s, v] = svd([1 2; 3 4])
 u * s * v'

行列の各列ごとに2ノルム正規化

 function retval = normalize_column(M)
   [m, n] = size(M);
   for i = 1:n
      nr = norm(M(:, i));
      if (nr != 0)
        nv([1:m], i) = nr;
      else
        nv([1:m], i) = 1;
      endif
   endfor
   retval = M ./ nv;
 endfunction

 m = [1, 2; 3 4]
 a = normalizeColumn(m)

もっと良い方法があると思う。

別スクリプトを include する。

 source("script.m")

source は遅い。データ(行列)を読み込むだけなら load の方が良い。

sscanf

 [val, count] = sscanf(line, "%s")

では val は column vector にならずに、文字列が連結されてしまう。

web

• http://www.octave.org/
  • マニュアル
• http://adlib.rsch.tuis.ac.jp/...
• http://www.higuchi.ecei.tohoku.ac.jp/...
• http://www.bekkoame.ne.jp/~ponpoko/
• http://www.mlb.co.jp/.../science/octave/

関連情報

• PDL