回帰の問題などでは(X'X)^-1をよく使う。 Xが横長であるときX'Xが非常に大きな行列になるので、 逆行列を求めるのが大変で高速化が重要になる。そこで、Woodbury行列反転公式(Woodbury matrix inversion formula) (A+BCD)^-1 = A^-1 - A^-1 B(C^-1+DA^-1 B)^…

MACのMATLABで mex -O hogehoge.cpp とやると， error: stdio.h: No such file or directoryが出力され，対処に思いの外，時間がかかったので，メモ． 対処法 /Applications/MATLAB_R2011a_Student.app/bin/gccopts.shの CC='gcc-4.2' SDKROOT='/Developer/S…

MATLABでJETの画像とカラー画像を重ねた画像をよく見るのだが， どう作ればいいのか分からなかったので，調べたついでにメモ． % 準備 clear; close all; img = im2double(imread('gundam.jpg')); img = imresize(img, 0.1); gimg = rgb2gray(img); h = fspe…

MATLABでeigの速度がどの程度なのか少し気になっていたので， 適当に時間を計測してみた． % eigの時間測定 clear all; nLoop = 10; sizeM = 2.^(5:12); nSizeM = numel(sizeM); times = zeros(nSizeM, nLoop); for iSizeM = 1:nSizeM disp(sizeM(iSizeM)); …

本 ブルックスの知能ロボット論―なぜMITのロボットは前進し続けるのか?作者: ロドニーブルックス,Rodney Allen Brooks,五味隆志出版社/メーカー: オーム社発売日: 2006/01メディア: 単行本購入: 2人 クリック: 26回この商品を含むブログ (21件) を見るロドニ…