magnetic field transfer matrix from dipole potential on boundary
Bvol = vb_bem_inverse_meg(D,Vs,Fs,Sout,X,Q)
境界面モデルのポテンシャルを計算するための
双極子ポテンシャルからセンサ磁場への変換行列計算
D : 境界面頂点係数行列
Vs : 境界三角面の座標
Fs : 境界三角面頂点インデックス
X : 磁気センサー座標 (Npick, 3)
Q : 磁場の観測方向 (Npick, 3)
Sout : 三角面の外向き面積法線
Bvol : 双極子ポテンシャルからセンサ磁場への変換行列
体積電流による境界面ポテンシャル が X に作る磁場の Q 方向射影
2004-02-06 Taku Yoshioka
2004-12-26 M. Sato modified
Copyright (C) 2011, ATR All Rights Reserved.
License : New BSD License(see VBMEG_LICENSE.txt)0001 function Bvol = vb_bem_inverse_meg(D,Vs,Fs,Sout,X,Q) 0002 % magnetic field transfer matrix from dipole potential on boundary 0003 % Bvol = vb_bem_inverse_meg(D,Vs,Fs,Sout,X,Q) 0004 % 境界面モデルのポテンシャルを計算するための 0005 % 双極子ポテンシャルからセンサ磁場への変換行列計算 0006 % 0007 % D : 境界面頂点係数行列 0008 % Vs : 境界三角面の座標 0009 % Fs : 境界三角面頂点インデックス 0010 % X : 磁気センサー座標 (Npick, 3) 0011 % Q : 磁場の観測方向 (Npick, 3) 0012 % Sout : 三角面の外向き面積法線 0013 % 0014 % Bvol : 双極子ポテンシャルからセンサ磁場への変換行列 0015 % 体積電流による境界面ポテンシャル が X に作る磁場の Q 方向射影 0016 % 0017 % 2004-02-06 Taku Yoshioka 0018 % 2004-12-26 M. Sato modified 0019 % 0020 % Copyright (C) 2011, ATR All Rights Reserved. 0021 % License : New BSD License(see VBMEG_LICENSE.txt) 0022 0023 % 境界三角形の単位電流 が X に作る磁場の Q 方向射影 0024 BS = vb_sensor_mag(Vs,Fs,Sout,X,Q); % (Nsensor x Nvertex) 0025 0026 % 双極子ポテンシャル P0 からセンサ磁場への変換行列 0027 % D * P = P0 0028 % P = inv(D) * P0 0029 % B = BS * P = ( BS * inv(D) )* P0 0030 % = Bvol * P0 0031 % Bvol = BS * inv(D) = BS/D 0032 0033 Bvol = BS/D; 0034