Discussion :

[nouriy]

Bonjour
mon problème est le suivant.
j'ai une surface rectangulaire maillée sur la quelle je suis le trajet d'une ligne d'eau par conséquent je récupère pour chaque carreau traversé les coordonnées de sorties de l'eau( après des semaine ), sauf que je n'arrive pas à tracer la ligne reliant les coordonnées entre elles .
s'il vous plait, venez à secours c'est tout ce qui me reste pour finir mon programme et être dans les temps.

[Jerome Briot]

Il nous faudrait un exemple avec des données numériques simples pour pouvoir t'aider efficacement.

[nouriy]

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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
z=1
j=1
i=1
x0=0
y0=0
pent
devr %matrices réccupérées dans un fichier txt
%surface
u=zeros(10,5)
suerf(u(:,:))
        while z<=5
             while i<10 %&& j<=5
 
 
                if (i==10  && dir== double ('S'))
 
                       break
 
                elseif  (j==5 &&  dir== double ('E'))
                        break
                end
                if ((i==1 )|| (i==2 ))|| (i==4 )|| ( i==6 )|| (i==8 ) || (i==9 ))
                    i=i+1;
                else
 
                    if (pent_gch(i,j)~=0 && devr_gch(i,j)~=0)
                              p(i,j)= sqrt( (pent_gch(i,j))^2+(devr_gch(i,j))^2) 
                                        % calcule de l'angle de la lfgne de plus haute pente
                                 tan_alpha= devr_gch_(i,j)/pent_gch(i,j)
                                         %calculons delta_y
                                  delta_y=tan_alpha*L
                                         %détermfnonsles coordonnées de sortie
 
 
                                   [x1 y1 dir]=coordonnees_sortie( x0, y0,pent(i,j), devr(i,j), delta_y, tan_alpha, L, l)
 
                          else
 
                                     [x1 y1 dir]=coordonnees_sortie5( x0, y0,pent(i,j), devr(i,j),  L, l);
                          end
 
                           x(i,j)=x1;
                           y(i,j)=y1;
                           %quiver3(x, y, v,w)
                        [i,j]= cdirection2(dir,i,j) %direction de pour la prochaine facette
                        end
              disp('on a atteint un bord')
                   z=z+1
                   j=z
                   i=1
                   x0=j-1;
                   y0=0;
                   dir= double( '0')
end
disp ('on a parcouru toute la surface')
x
y

il suivre la trajectoire' de l'eau jusqu'à sa sortie des bords avec affichage en 3D
j'ai essayé plot 3; surf (x,y), surf( x, y,u) rien ( je ne suis pas très douée en programmation )


jespere que c'est compréhensible; le programme est la dans toute sa simplicité
maillage.doc

[nouriy]

je voudrais savoir si ça aurait été possible de générer une surface avec des pas de temps différents comme sur le fichier word, parce que dans ce cas ça m'aurait facilitée les choses d'avoir une surface pareille.

tu penses que ça serait possible d'afficher les mailles une apres l'autre et en traçant la ligne 'écoulement?
comme ça le surf du début, je pourrai le retirer.( à pat me permettre de visualiser la situation il e me sert pas rand choses)

[Jerome Briot]

Commençons par une version simplifiée :

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
x = 0:11;
y = [0 3 5 7 9];
 
[x,y] = meshgrid(x,y);
 
figure(1)
clf
mesh(x,y,zeros(size(x)));
 
view(2)
set(gca, 'ydir', 'reverse')
axis equal tight
 
hold on
 
x_water = [0 1 1.5 3 3.5];
 
for n = 2:size(x_water,2)
 
   quiver(x_water(n-1), y(n-1), x_water(n)-x_water(n-1), y(n)-y(n-1),1, 'r-', 'linewidth', 2)
   pause(.5)
 
end

[nouriy]

merci beaucoup pour ton aide,
je vais m’imprégner de ton code,
cependant, je comprends ce que tout ce que ton code fait mais je ne connais pas tout ce que tu as utilisé. si tu pouvais me briffer un tout petit peu dessus, ça m'aiderai.
si j'ai une question, je peux directement te la poser sur ton blog?

[Jerome Briot]

Qu'est-ce que tu ne comprends pas exactement ?

[nouriy]

Code :

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

1
2
3
set(gca, 'ydir', 'reverse')
axis equal tight
linewidht

je comprends pas ces syntaxes et elles n'apparaissent pas dans matlab, qu'est-ce qu'elles représentent ?

[Jerome Briot]

set(gca, 'ydir', 'reverse')

Permet d'inverser le sens de l'axe des ordonnées. Le point (0,0) se trouve donc en haut à gauche (plutôt qu'en bas à droite).

Voir "axes properties" dans la documentation :

Code :

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

docsearch axes AND properties

axis equal tight

« equal » permet d'uniformiser l'échelle entre l'axe des abscisses et celui des ordonnées.
« tight » permet d’ajuster au plus près les limites des axes.

Voir la documentation de la fonction axis :

Code :

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

doc axis

linewidht

La propriété LineWidth spécifie l'épaisseur de la ligne tracée.
Voir « Lineseries properties » dans la documentation :

Code :

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

docsearch lineseries AND properties

[nouriy]

merci pour ces infos
par contre j'ai remarqué que au début de l'affichage, l'origine est en haut à gauche, des que ça passe au second rand l'origine change de position et se place en haut à droite.

dis, tu vois un problème à cette syntaxe:

elseif (( dir== double ('NE') )||( dir== double ('SE') )||( dir== double ('SO')) || (dir== double ('NO') ) )?
??? Operands to the || and && operators must be convertible to logical scalar values.

[Jerome Briot]

par contre j'ai remarqué que au début de l'affichage, l'origine est en haut à gauche, des que ça passe au second rand l'origine change de position et se place en haut à droite.

Dans mon exemple, l'origine ne bouge pas.

??? Operands to the || and && operators must be convertible to logical scalar values.

Voir la FAQ : Operands to the || and && operators must be convertible to logical scalar values.

D'autre part, pour comparer des chaines de caractères, on utilise la fonction strcmp

Pour finir, renomme la variable dir car c'est le nom d'une fonction MATLAB.

[nouriy]

merci, franchement tu me sauves la vie

[nouriy]

par rapport au tracé de la lignes j'ai un problème qui est survenu
en tant normal ma superficie fait L=2092 et l=20

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
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
z=1% me permet d’itérer i
j=1 %me permet d’itérer la longueur suivant la colonne
i=1%me permet d’itérer la largeur suivant la ligne
x0=0
y0=0
pent
devr %matrices récupérées dans un fichier txt
while z<=5
    while i<10 %&& j<=5
 
 
        if (i==10  && dir== double ('S'))
 
            break
 
        elseif  (j==5 &&  dir== double ('E'))
            break
        end
        if ((i==1 )|| (i==2 ))|| (i==4 )|| ( i==6 )|| (i==8 ) || (i==9 ))
            i=i+1;
        else
 
            if (pent_gch(i,j)~=0 && devr_gch(i,j)~=0)
                p(i,j)= sqrt( (pent_gch(i,j))^2+(devr_gch(i,j))^2)
                % calcule de l'angle de la ligne de plus haute pente
                tan_alpha= devr_gch_(i,j)/pent_gch(i,j)
                %calculons delta_y
                delta_y=tan_alpha*L
                %déterminons les coordonnées de sortie
 
 
                [x1 y1 dir]=coordonnees_sortie( x0, y0,pent(i,j), devr(i,j), delta_y, tan_alpha, L, l)
 
            else
 
                [x1 y1 dir]=coordonnees_sortie5( x0, y0,pent(i,j), devr(i,j),  L, l);
            end
 
            x(i,j)=x1;
            y(i,j)=y1;
            %quiver3(x, y, v,w)
            [i,j]= cdirection2(dir,i,j) %direction de pour la prochaine facette
        end
        disp('on a atteint un bord')
        z=z+1
        j=z
        i=1
        x0=j-1;
        y0=0;
        dir= double( '0')
    end
    disp ('on a parcouru toute la surface')
    x
    y

mon problème est le suivant; puisque je ne veux qu'afficher qu'une partie de ma route par exemple 1077-1082,
je n'arrive plus à bien initialiser mes paramètres et en utilisant le code que tu m'as proposer pour l'affichage j'ai le schéma qui commence toujours en (0,0)
je propose ceci comme initialisation:

Code :

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

1
2
3
4
5
6
z=1
i=1
x(1,1)=1077
y(1,1)=0
while z<=5
    ........

[Jerome Briot]

Le plus simple est d'afficher toute la grille et de régler les limites de l'axe des ordonnées avec ylim.

En reprenant mon premier code :

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
x = 0:11;
y = [0 3 5 7 9];
 
[x,y] = meshgrid(x,y);
 
figure(1)
clf
mesh(x,y,zeros(size(x)));
 
view(2)
set(gca, 'ydir', 'reverse')
axis equal tight
 
hold on
 
x_water = [0 1 1.5 3 3.5];
 
for n = 2:size(x_water,2)
 
   quiver(x_water(n-1), y(n-1), x_water(n)-x_water(n-1), y(n)-y(n-1),1, 'r-', 'linewidth', 2)
   ylim([y(n-1) y(n)]) 
   pause(.5)
 
end

[nouriy]

Le plus simple est d'afficher toute la grille et de régler les limites de l'axe des ordonnées avec ylim.

ylim([y(n-1) y(n)])

j'ai essayé le code, je me suis rendu compte que avec la pause de 0.5 l'affichage était successif entre "y(n) et y(n-1)"


le problème c'est que ce sont les x que je veux figer mais quand j'essaie la même syntaxe j'ai des résultats d'erreurs
Tant bien que ça marcherai le problème c'est que au fur et à mesure que n variera j'aurais des séquences par contre ce que je veux c'est plutôt figer une séquence particulière

[nouriy]

clf

le clf est il necessaire? Parce que quand je l'utilise j'ai des segments détachés et quand il n'est pas là j'ai le graphe complet.
En effet pourquoi l'utilises tu?

[Jerome Briot]

C'est juste pour nettoyer la fenêtre avant de tracer.

Tu n'en a normalement pas besoin.

[nouriy]

J'ai réussi à trouver une méthode,
j'ai encore une question, le problème c'est que là je suis obligée de tout calculer d'abord
donc mon programme perd en efficacité. alors que pour une intervention rapide le mieux ça serait de cibler la zone des le départ, sauf que les coordonnées d'entrées restent la difficulté.
Qu'en penses tu?