Home > functions > leadfield > bem > vb_bem_meg.m

vb_bem_meg

PURPOSE ^

magnetic field by using BEM method

SYNOPSIS ^

function B = vb_bem_meg(D,Sout,Vhead,Fhead,V,J,pick,Qpick)

DESCRIPTION ^

 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)

CROSS-REFERENCE INFORMATION ^

This function calls: This function is called by:

SOURCE CODE ^

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);

Generated on Tue 27-Aug-2013 11:46:04 by m2html © 2005