Discussion :

[Clem44120]

Bonjour à tous,
J'ai un peu de mal à créer une boucle à deux variables sur matlab.
Les variables ini et fim sont les itérations en début et en fin de simulation. La variable tempo est le temps de simulation en secondes.
Dans l'extrait du code ci-après j'aimerais calculer la distance que mon mobile a parcouru toutes les deux minutes. Le long de la courbe que j'obtiens j'ai tracé une droite sur laquelle je mets un point avec le temps parcouru.
La réponse sera peut être évidente mais je n'arrive vraiment pas à trouver.. Matlab me dit que "d(0) must be a positive integar or logical"
Dans l'extrait du code c'est un peu brouillon, la partie avec les length est présente car j'ai tenté de définir d(i,j)
Merci d'avance

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
 
int=120;
disc=10;
for i=ini:disc:fim
    for j=0:int:tempo
   M = [cos(PSI(i)+pi/2) -sin(PSI(i)+pi/2)
         sin(PSI(i)+pi/2) cos(PSI(i)+pi/2)];
   p_int(i,1) = -Y(i)+Yini;
   p_int(i,2) = X(i)-Xini;
   desenho = M*navio';
   fill((desenho(1,:)+p_int(i,1)),(desenho(2,:)+p_int(i,2)),[0.5 0.5 0.5]);
   k=[0 0];
   l=[0 head_reach];
   q=plot(k,l);
   set(q,'Color','red','LineWidth',4);
 
 d=zeros(length(fim-ini),length(tempo));
 
 d(j+1)=d(j) + sqrt(X(i+1) - X(i)).^2;
 
 text(0,d(j),'2 min    ','HorizontalAlignment','right');
 scatter(0,d(j),'filled','LineWidth',2,'MarkerFaceColor',[0 0 1]);
 
 
 
    end 
end

[Jerome Briot]

Comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
int=120;
disc=10;
 
iterations = ini:disc:fim;
temps = 0:int:tempo;
 
for i = 1:numel(iterations)
    for j = 1:numel(temps)
        %...
    end
end

Il faut ensuite modifier le code en conséquence

[Clem44120]

Merci.
Je vaie essayer cela alors. Le seul probleme ne sera-t-il pas avec mon graphe ? Lorsque j'utilise la fonctionn fill j'appelle le i allant de ini a fim. Or desormais il ira de 1 a numel(iterations)..

[Clem44120]

J'ai essayé comme cela :
disc correspond au temps des itérations. J'ai une itérations toutes les 1/10 de s
Ce que je veux c'est calculer la distance parcourue toutes les 2 minutes, c'est à dire toutes les 1200 itérations..
Je pense que mes distances calculées ne sont pas bonnes du tout ..

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
f1 = figure;
hold on;
 
int=1200;
disc=10;
 
temps = ini:1200:fim;
 
for i= ini:disc:fim;
  for j = 1:numel(temps)
   M = [cos(PSI(i)+pi/2) -sin(PSI(i)+pi/2)
         sin(PSI(i)+pi/2) cos(PSI(i)+pi/2)];
   p_int(i,1) = -Y(i)+Yini;
   p_int(i,2) = X(i)-Xini;
   desenho = M*navio';
   fill((desenho(1,:)+p_int(i,1)),(desenho(2,:)+p_int(i,2)),[0.5 0.5 0.5]);
   k=[0 0];
   l=[0 head_reach];
   q=plot(k,l);
   set(q,'Color','red','LineWidth',4);
 
 
 
 d(j+1)=d(j) + sqrt(X(i+1) - X(i)).^2;
 
 text(0,d(j),'2 min    ','HorizontalAlignment','right');
 
 
 
       end
   end

[Clem44120]

Si ça peut aider, j'arrivais à obtenir ce que je voulais en faisant comme cela :
Mais je désire maintenant automatiser le programme car je devais, selon le temps de simulations, rajouter le nombre de d suffisants...

Code :

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

1
2
3
4
5
6
d60=0;
   for i=ini:ini+600
       d60 = d60 + sqrt(X(i+1) - X(i)).^2;
   end 
   text(0,d60,'1 min    ','HorizontalAlignment','right');
   text(0,d60,'    19.4 kts','HorizontalAlignment','left');

[Jerome Briot]

Je pense que tu confonds deux choses : les valeurs réelles du phénomène physique et leur stockage numérique (sous forme de tableau)

Prenons un phénomène P qui varie toutes les 250ms pendant 2s.

Les valeurs réelles (discrètes) du paramètre temporel : 0, 0.25, 0.5, 1, 1.25, 1.5, 1.75, 2
Le vecteur numérique associé : t = [0 0.25 0.5 1 1.25 1.5 1.75 2]

Supposons une relation très simple : P(t) = t²

Mathématiquement tu écris :

• P₀ = t₀*t₀
• P_0.25 = t_0.25*t_0.25
• P_0.5 = t_0.5*t_0.5
• ...

Numériquement tu écriras :

• P(1) = t(1)*t(1)
• P(2) = t(2)*t(2)
• P(3) = t(3)*t(3)
• ...

Et non pas :

• P(0) = t(0)*t(0)
• P(0.25) = t(0.25)*t(0.25)
• P(0.5) = t(0.5)*t(0.5)
• ...

[Clem44120]

Merci, je vois bien l'idée. Je ne sais pas si je fais un blocage mais impossible de le coder..
Je vais continuer de chercher.

[Clem44120]

Dans l'idée il me faudrait cela en fait je pense :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
for n = ini : 1200 : fim
  for j=1:numel(temps)
 
 dj+1 = dj + sqrt(X(n+1) - X(n)).^2;
 
 text(0,dj,'x min    ','HorizontalAlignment','right');
 
 
 
       end
   end

[Clem44120]

J'ai cru comprendre que ce n'était pas terrible mais je peux utiliser genvarname peut être ?

[Jerome Briot]

La variable tempo est le temps de simulation en secondes.
Dans l'extrait du code ci-après j'aimerais calculer la distance que mon mobile a parcouru toutes les deux minutes.

Reprenons depuis le début.

Le contenu du tableau numérique tempo :

Code :

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

tempo = 0:1:6*60; % exemple pour 6 minutes de simulation

Les temps en secondes à intervalle de 2 minutes : 0, 120, 240, 360

Ces temps correspondent aux indices [1 121 241 361] du tableau tempo :

Code :

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

1
2
3
4
5
t = tempo([1 121 241 361])
 
t =
 
     0   120   240   360

ce qui s'écrit aussi :

Code :

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

1
2
3
4
5
t = tempo(1:120:end)
 
t =
 
     0   120   240   360

Donc la boucle peut devenir :

Code :

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

1
2
3
4
5
for n = ini : 1200 : fim
    for j = 1:120:numel(tempo)
        % ...
    end
end

Et comme tu veux stocker le résultat dans un tableau d(1), d(2), d(3)... :

Code :

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

1
2
3
4
5
for n = ini : 1200 : fim
    for j = 1:120:numel(tempo)
        d((j-1)/120+1) = ...
    end
end

[Clem44120]

Super top !
Cela m'a mis sur la piste. En effet, je mélangeais un peu tout ...
Petite question supplémentaire, je cherche à utiliser la fonction num2str afin de légender mes points toutes les deux minutes. J'ai bien mes points avec la fonction scatter mais je n'arrive pas à mettre une légende "2 min" puis "4 min" en face des points.
J'ai essayé de la maniere suivante, sachant que j'ai 5 itérations (j de 1 a 5) :

Code :

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

1
2
 
text(0,d(j),'num2str(2*j)','HorizontalAlignment','right');

[Jerome Briot]

Soit :

Code :

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

text(0, d(j), [num2str(2*j) ' min'], 'HorizontalAlignment', 'right');

Soit :

Code :

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

text(0, d(j), sprintf('%d min', 2*j), 'HorizontalAlignment', 'right');

[Clem44120]

Merci pour tout !