Discussion :

[ABN84]

bonjour,
y a-t-il moyen de se passer du for dans cette expression:

Code :

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

1
2
3
4
5
6
7
for k=2*pas+1:length(i1)-1*pas
    for l=2*pas+1:length(i1)-1*pas
        x1=x1+(j1(k,l)*j2(k+i,l+j));
        y1=y1+(j1(k,l)*j1(k,l));
        z1=z1+(j2(k,l)*j2(k,l));
    end
end

merci

[Jerome Briot]

A tester :

Code :

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

1
2
3
4
5
6
7
8
k = 2*pas+1:length(i1)-1*pas;
 
idx = sub2ind(size(j1),k,k);
y1 = sum(j1(idx).*j1(idx));
z1 = sum(j2(idx).*j2(idx));
 
idx = sub2ind(size(j1),k+i,k+j);
x1 = sum(j1(idx).*j2(idx));