magnetic field by using BEM method B = vb_bem_meg(D,Sout,Vhead,Fhead,V,J,pick,Qpick) D : 境界面係数行列 Vhead : 境界三角面の座標 Fhead : 境界三角面頂点インデックス Sout : 三角面の外向き面積法線 V : 電流双極子座標 (m) (Ndipole,3) J : 電流双極子モーメント (A) (Ndipole,3) pick : 磁気センサー座標 (Npick, 3) Qpick : 磁気センサー方向 (Npick, 3) B : pick における磁場 （テスラ） (Npick, Ndipole) 2004-02-06 Taku Yoshioka 2004-12-26 M. Sato modified 2005-8-5 M. Sato modified Copyright (C) 2011, ATR All Rights Reserved. License : New BSD License(see VBMEG_LICENSE.txt)
0001 function B = vb_bem_meg(D,Sout,Vhead,Fhead,V,J,pick,Qpick) 0002 % magnetic field by using BEM method 0003 % B = vb_bem_meg(D,Sout,Vhead,Fhead,V,J,pick,Qpick) 0004 % D : 境界面係数行列 0005 % Vhead : 境界三角面の座標 0006 % Fhead : 境界三角面頂点インデックス 0007 % Sout : 三角面の外向き面積法線 0008 % V : 電流双極子座標 (m) (Ndipole,3) 0009 % J : 電流双極子モーメント (A) (Ndipole,3) 0010 % pick : 磁気センサー座標 (Npick, 3) 0011 % Qpick : 磁気センサー方向 (Npick, 3) 0012 % 0013 % B : pick における磁場 （テスラ） (Npick, Ndipole) 0014 % 0015 % 2004-02-06 Taku Yoshioka 0016 % 2004-12-26 M. Sato modified 0017 % 2005-8-5 M. Sato modified 0018 % 0019 % Copyright (C) 2011, ATR All Rights Reserved. 0020 % License : New BSD License(see VBMEG_LICENSE.txt) 0021 0022 0023 % 自己重み係数：基底関数を三角面上で積分した寄与を加える 0024 fprintf('--- Auto solid angle calculation\n'); 0025 D = vb_solid_auto_grk(D,Fhead,Sout); 0026 0027 0028 % Bvol : 双極子ポテンシャルからセンサ磁場への変換行列 0029 fprintf('--- Inverse Matrix\n'); 0030 Bvol = vb_bem_inverse_meg(D,Vhead,Fhead,Sout,pick,Qpick); 0031 0032 clear D 0033 0034 % 磁場計算 0035 fprintf('--- Magnetic field calculation\n'); 0036 B = vb_bem_meg_linear_grk(Bvol,Vhead,Fhead,V,J,pick,Qpick);