Discussion :

[ABN84]

Bonjour,
Voici un code que j'ai fait:

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
 
for i=2:length(t)
   [...]
    A(i,201)=Ac2((Y(4,: )-y(i)).^2+(Y(3,: )-y(i-1)).^2+...
                      10*(Y(4,: )-Y(1,: )).^2+...                 
                      (Y(2,: )-y(i-2)).^2+(Y(1,: )-y(i-3)).^2==...
                  min((Y(4,: )-y(i)).^2+(Y(3,: )-y(i-1)).^2+...
                      10*(Y(4,: )-Y(1,: )).^2+...
                      (Y(2,: )-y(i-2)).^2+(Y(1,: )-y(i-3)).^2));
    B(i,201)=Bc2((Y(4,: )-y(i)).^2+(Y(3,: )-y(i-1)).^2+...
                      10*(Y(4,: )-Y(1,: )).^2+...                 
                      (Y(2,: )-y(i-2)).^2+(Y(1,: )-y(i-3)).^2==...
                  min((Y(4,: )-y(i)).^2+(Y(3,: )-y(i-1)).^2+...
                      10*(Y(4,: )-Y(1,: )).^2+...                 
                      (Y(2,: )-y(i-2)).^2+(Y(1,: )-y(i-3)).^2));
end

l'idée est d'affecter à A(i,201), la valeur qui se trouve à la position qui correspond au minimum de (Y(4,: )-y(i)).^2+(Y(3,: )-y(i-1)).^2+10*(Y(4,: )-Y(1,: )).^2+(Y(2,: )-y(i-2)).^2+(Y(1,: )-y(i-3)).^2
Le soucis, ces deux lignes consomment 68% tu temps de tout mon algo, y a-t-il une autre façon d'ecrire ça pour que ça soit optimisé?

J'espère que mon explication st claire.

merci

[Jerome Briot]

Essaie ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
for i=2:length(t)
 
   [pasbesoin, idx] = min((Y(4,: )-y(i)).^2+(Y(3,: )-y(i-1)).^2+...
                      10*(Y(4,: )-Y(1,: )).^2+...
                      (Y(2,: )-y(i-2)).^2+(Y(1,: )-y(i-3)).^2);
 
    A(i,201) = Ac2(idx);
    B(i,201) = Bc2(idx);
 
end

[ABN84]

Merci DUT, ça a bien allégé l'algo