difficulté a comprendre une fonction

Version imprimable

bonjour , alors je veux un aide de vous part pour faire marcher et comprendre cette fonction :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
function [lf] = calcule(rd, pnt, cond);
 
siz = size(rd);
if any(siz==1)
  % positions are specified as a single vector
  Ndipoles = prod(siz)/3;
elseif siz(2)==3
  % positions are specified as a Nx3 matrix -> reformat to a single vector
  Ndipoles = siz(1);
  rd = rd';
  rd = rd(:);
else
  error('incorrect specification of dipole locations');
end
 
Npnt     = size(pnt,1);
lf       = zeros(Npnt,3*Ndipoles);
s1       = size(rd);
 
if s1(1)>s1(2)
  % make sure that the dipole position is a row-vector
  rd = rd';
end
 
for i=1:Ndipoles
  r = pnt - ones(Npnt,1) * rd((1:3) + 3*(i-1));
  R = (4*pi*cond) * (sum(r' .^2 ) .^ 1.5)';
  if any(R)==0
    warning('dipole lies on boundary of volume model');
  end
  lf(:,(1:3) + 3*(i-1)) = r ./ [R R R];
end

jai essaiyer avec ces argument en entré mai ça na pas marcher :

rd=[0 0.5 0.5]
pnt=[0 0 1]
cond=[0.33 0.0252 0.33]

alors aider moi svp

18/01/2009, 14h58
Caro-Line

Où as-tu trouvé cette fonction ?
Que sont censés représenter physiquement ces arguments ?
Pourquoi as-tu testé avec ces valeurs là et qu'est-ce qui ne fonctionne pas quand tu le fais ?
18/01/2009, 18h50
roubas

Alors c’est une méthode qui calcule la méthode BEM ( Boundry element méthode ) qui calcule un potentielle électrique dans un maillage , les argument en entré ces
- le point de maillage
- le point de mesure de
- la conductivité dans ces surfaces
-
Et on a trois surface, alors aider moi svp a exécuté cette fonction !!?
18/01/2009, 18h53
Jerome Briot

Et la réponse à la première question ?

Citation:

Envoyé par Caro-Line

Où as-tu trouvé cette fonction ?
18/01/2009, 19h02
Caro-Line

Et les 3ème et 4ème questions ?

Citation:

Envoyé par Caro-Line

Pourquoi as-tu tester avec ces valeurs là et qu'est-ce qui ne fonctionne pas quand tu le fais ?
18/01/2009, 19h50
roubas

je les trouver sur internet
18/01/2009, 20h19
Jerome Briot

Citation:

Envoyé par roubas

je les trouver sur internet

Et où exactement ?
Sur quel site ?
19/01/2009, 20h24
roubas

je me rapelle pas .... !!!
19/01/2009, 20h29
Jerome Briot

Citation:

Envoyé par roubas

je me rapelle pas .... !!!

C'est bien dommage...

Bon sinon, c'est quoi le problème ?
Le code ne fonctionne pas ? Il y a un message d'erreur qui apparait ?

alors voila la function :

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
function [lf] = inf_medium(rd, pnt, cond);
 
% INF_MEDIUM_LEADFIELD calculate the infinite medium leadfield
% on positions pnt for dipole position R and conductivity cond
%       
% [lf] = inf_medium_leadfield(R, pnt, cond)
 
rd=[0 0.5 0.5;0 0 0.5]
pnt=[0 0 1; 1 1 1]
cond=[0.3 0.0025 0.3]
%----------------------------------------------------------------
siz = size(rd);
if any(siz==1)
  % positions are specified as a single vector
  Ndipoles = prod(siz)/3;
elseif siz(2)==3
  % positions are specified as a Nx3 matrix -> reformat to a single vector
  Ndipoles = siz(1);
  rd = rd';
  rd = rd(:);
  rd
else
  error('incorrect specification of dipole locations');
end
%---------------------------------------------------------------------
Npnt     = size(pnt,1);
lf       = zeros(Npnt,3*Ndipoles);
s1       = size(rd);
 
if s1(1)>s1(2)
  % assurez-vous que la position est un dipôle-vecteur ligne
  rd = rd';
end
 
for i=1:Ndipoles
  r = pnt - ones(Npnt,1) * rd((1:3) + 3*(i-1));
  R = (4*pi*cond) * (sum(r' .^2 ) .^ 1.5)';
  if any(R)==0
    warning('dipôle est situé sur la frontière du volume modèle');
  end
  lf(:,(1:3) + 3*(i-1)) = r ./ [R R R];
end

il me genair cette érreur :

Code:

1
2
3
4
5
??? Error using ==> mtimes
Inner matrix dimensions must agree.
 
Error in ==> inf_medium_leadfield at 37
  R = (4*pi*cond) * (sum(r' .^2 ) .^ 1.5)';

je ponce que la (cond) doit avoir une seul valeur et non vecteur, car la fonction doit etre integré dans une autre fonction

Code:

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
function [lf] = lead(pos, elc, vol)
 
% [lf] = lead(pos, elc, vol)
%
% with the input arguments
%   pos		position dipole (1x3 or Nx3)
%   elc		position electrodes (optional, can be empty)
%   vol		volume conductor model
%
% the volume conductor model is a structure and should have the fields
%   vol.bnd	structure array with vertices and triangles of each boundary
%   vol.cond	conductivity of all compartments
%   vol.mat 	system matrix, which can include the electrode interpolation
%   vol.tra 	transfer matrix for interpolation from vertices to electrodes
%
% the compartment boundaries are described by a structure array with
%   vol.bnd(i).pnt
%   vol.bnd(i).pnt
 
ncmp = length(vol.cond);
cond = vol.cond(vol.source);
% concatenate the vertices of all compartment boundaries in a single Nx3 matrix
  pnt = [];
  for i=1:ncmp
    pnt = [pnt; vol.bnd(i).pnt];
  end
 
  lf = inf_medium(pos, pnt, cond);
end

Aider moi SVP , j'attendrai vous réponce avec impatience !!

12/02/2009, 15h01
Caro-Line

Il y a manifestement un problème de dimensions dans ton expression.
A toi de regarder la dimension de chacun des éléments et de comprendre pourquoi cela ne fonctionne pas.

Utilise par exemple le debuggueur (help debug) pour t'arrêter sur la ligne incriminée et regarde les valeurs de tes éléments avant que ça ne plante.
12/02/2009, 15h54
Loïc B.

Apparemment, tu as raison, il faut une unique valeur dans cond.
En même tps, ça se comprend puisque il s'agit d'une conductivité de ton milieu. Si tu en définissait plusieurs, tu aurais un code bien plus compliqué, puisqu'il faudrait définir les zones correspondantes à ces conductivités.

Donc au final :
- rd, matrice de taille (nbre de dipole, 3) contenant sur chaque ligne les coordonnées 3d de tes dipôles
- pnt, matrice de taille (nbre de points de mesure,3) contenant sur chaque ligne les coordonnées 3d de ces points
- cond, scalaire, valeur unique de conductivité de ton milieu

Le résultat est une matrice (nbre de points de mesures,3) contenant la valeur du champ sous forme d'un vecteur 3d.

Bonne journée.

PS : ton code, tu l'as peut être récupéré ici : https://svnserv.cbs.mpg.de/eeglab/br...um_leadfield.m