difficulté a comprendre une fonction

**roubas** · 17/01/2009, 21h54

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

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

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

**Caro-Line** · 18/01/2009, 14h58

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 ?

**roubas** · 18/01/2009, 18h50

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 !!?