octave のバックアップ(No.6)

• バックアップ一覧
• 差分 を表示
• 現在との差分 を表示
• ソース を表示
• octave へ行く。
  • 1 (2014-12-31 (水) 11:15:24)
  • 2 (2014-12-31 (水) 11:16:11)
  • 3 (2015-02-14 (土) 17:23:38)
  • 4 (2015-08-25 (火) 05:20:25)
  • 5 (2015-08-25 (火) 21:10:59)
  • 6 (2015-08-26 (水) 04:22:54)
  • 7 (2015-11-05 (木) 11:52:12)
  • 8 (2015-11-10 (火) 22:33:42)
  • 9 (2015-11-11 (水) 02:18:35)

ツール?

Ubuntu12.04の問題 †

• octaveの制御パッケージはUbuntuにデフォルトインストールされない
  • sudo aptitude install octave-control
• octaveの最適化パッケージはデフォルトインストールされない
  • sudo aptitude install octave-optim
• xlim, axisなど、窓を限定するコマンドが一切使えない
  • plotで枠を書いてやるしかない

ライブラリ †

• octaveスクリプト

ファイル書き出し †

• dlmwrite('filename', val, ' ');

最適化 †

• leasqr

global verbose; verbose = true;

%Load Data
mx = [0:0.1:10]';
vy = mx.*mx;

%Set Function
function vY = myfun( mX, vP )
    a=vP(1); b=vP(2); c=vP(3);
    vY = a*mX(:,1).*mX(:,1) + b*mX(:,1) + c;
endfunction

%Inital Values
vp0 = [0.9; 0; 0];

%NLS
opt = optimset('Display', 'iter');
[vyhat vp] = leasqr( mx, vy, vp0, @myfun );
ssr = (vy-vyhat)'*(vy-vyhat);

disp('========== RESULTS ==========');
vp
ssr

%Output Results
myfile = fopen('output_nls.txt', 'w');
fprintf(myfile, '========== RESULTS ==========\n');
fprintf(myfile, 'vp = %12.7f\n', vp);
fprintf(myfile, 'SSR = %12.7f\n', ssr);
fclose(myfile);

• y-f(x)を最小化するやつ

clear
#global verbose; verbose = true;

%Load Data
mx = [0:0.1:10]';
global vy_tobe = mx.*mx;

%Set Function
function vY = myfun( mX, vP )
    global vy_tobe;
    a=vP(1); b=vP(2); 
    vY = a*mX + b - vy_tobe;
endfunction

%Inital Values
vp0 = [10; 0];

%NLS
opt = optimset('Display', 'iter');
[vyhat vp] = leasqr( mx, zeros(size(mx)), vp0, @myfun );

disp('========== RESULTS ==========');
vp
vyhat

plot(mx, vy_tobe);
hold on;
plot(mx, vp(1)*mx+vp(2))
hold off;

• 下値抵抗線（下方向のズレを過大評価する非線形関数を最適化）
  • 

clear
#global verbose; verbose = true;

%Load Data
mx = [0:0.1:10]';
global vy_tobe = 5*rand(size(mx)) + mx;

%Set Function
function vY = myfun( mX, vP )
    global vy_tobe;
    a=vP(1); b=vP(2);
   vY = a*mX + b - vy_tobe;
#    vY(vY<0)=10*vY(vY<0); #上
    vY(vY>0)=10*vY(vY>0); #下
endfunction

%Inital Values
vp0 = [-100; 0]; # だいぶエクストリームな初期値（0.98,0.44とかが本当の値）

%NLS
opt = optimset('Display', 'iter');
[vyhat vp] = leasqr( mx, zeros(size(mx)), vp0, @myfun );

disp('========== RESULTS ==========');
vp
vyhat

plot(mx, vy_tobe);
hold on;
plot(mx, vp(1)*mx+vp(2))
hold off;