set parameter for dipole search [dipole_parm] = vb_set_dipole_parm_test 2007-7-12 made by M.Sato Copyright (C) 2011, ATR All Rights Reserved. License : New BSD License(see VBMEG_LICENSE.txt)
0001 function [dipole_parm] = vb_set_dipole_parm_test 0002 % set parameter for dipole search 0003 % [dipole_parm] = vb_set_dipole_parm_test 0004 % 0005 % 2007-7-12 made by M.Sato 0006 % 0007 % Copyright (C) 2011, ATR All Rights Reserved. 0008 % License : New BSD License(see VBMEG_LICENSE.txt) 0009 0010 % --- Input MEG file(relative path from proj_root) 0011 dipole_parm.megfile = 'move_artifact.meg.mat'; 0012 0013 % --- Output file(relative path from proj_root) 0014 dipole_parm.dipolefile = []; 0015 0016 % set time window & trial number to estimate 0017 dipole_parm.twin_meg = []; % Time window for analysis 0018 dipole_parm.trial_id = []; % trial number to estimate 0019 0020 % Error function for dipole search 0021 % = 0 : 'vb_dipole_error': for external current (Biot-Savart Eq.) 0022 % = 1 : 'vb_dipole_error_sarvas': for cortical current (Sarvas Eq.) 0023 % = 2 : 'vb_dipole_error_magdipole': for magnetic dipole 0024 dipole_parm.errfunc = 0; 0025 0026 % initial dipole position [Ndipole x 3] 0027 dipole_parm.Vinit = [0.1 0.0 -1.0]; 0028 %dipole_parm.Vinit = [0.1 0.0 -1.0; 0.0 0.0 -1.0]; 0029 %dipole_parm.Vinit = [0.0 0.0 -0.5; 0.1 0.0 -1.0; 0.0 0.0 -1.0]; 0030 %dipole_parm.Vinit = [0 0 0]; 0031 0032 dipole_parm.MaxIter = 5000; % max iteration number in optimization 0033 dipole_parm.TolX = 1.e-8; % Xの終了許容値 0034 dipole_parm.TolFun = 1.e-10; % 関数値の終了許容値 0035 dipole_parm.Display = 'off'; % = 'iter' or 'off' 0036 dipole_parm.mode = 1; 0037 % = 0; % Use 'fminsearch' in Matlab 0038 % = 1; % Use 'fminunc' in Optimization Toolbox 0039 % = 2; % Use 'fminunc1'