Discussion :

[MarionSCN]

Bonjour à tous !

J'apprends beaucoup de chose avec ce site, merci ! C'est à mon tour de poser une question ...

Je construit une interface graphique me permettant d'analyser graphiquement un fichier CSV assez lourd.
J'ai dors et déjà construit cette interface, elle est fonctionnelle. Voilà à quoi elle ressemble:


Comme vous pouvez le voir (si vous connaissez), il s'agit d’analyser les performances d'un bateau. J'ai en tout 13 variables (vitesse, cap, tangage, force dans la voile etc...) qui sont des vecteurs d'environ 13 000 lignes.

Voici donc deux questions qui portent sur l'affichage en graphe polaire:

1/ J'aimerai pouvoir "déplacer" ma polaire de telle sorte à être "angle au vent" (0°= face au vent) plutôt qu’être au cap GPS comme c'est le cas. Ici, mon vent vient du 330, et j'aimerai qu'il soit au niveau du zéro. J'ai ajouté une boite de dialogue pour récupérer cette information (TWA: True Wind Angle) mais je ne connais pas la commande pour translater ma polaire. Voila ce que j'écrit pour le moment:

Code :

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

1
2
3
4
 
trp= polar(a,b,'.') %a et b sont les valeurs sélectionnées dans les menus déroulants
view(90,-90) %Place le diagramme avec zéro en haut et sens horaire
set(trp)

Je pense que je n'ai pas totalement compris la commande "view" meme avec l'aide MatLab, pourrait on m'expliquer ?

2/ J'ai codé une recherche d'un maximum, mais je ne sais pas comment le superposer avec mes premiers points. Le voici:

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
 
% On récupère le maximum de la sélection du menu 2
MAXb=single(max(b));
 
% Recherche l'indice de la première occurrence du max
indicemaxb= (find(b==MAXb,1)); 
 
% On cherche l'élément de a correspondant en radian si c'est le cap
if x==6
    MAXa= deg2rad(single(a(indicemaxb)));
else
    MAXa= single(a(indicemaxb));
end
 
trpmax= polar(MAXa,MAXb,'or');
view(90,-90)
set(trpmax)

Donc pour l'instant, j'affiche mes points de manière séparée de mon point max en rouge.
Or, je voudrai pouvoir le mettre sur la même polaire pour repérer facilement le point lors de l'analyse. J'ai essayé le "plot" sans succès ...


Je vous remercie d'avance pour vos réponses !

PS: J'ai construit tout le code avec l'aide de MatLab et ce site, de manière plus précise grâce à "Jerome Briot", si tu me lis, merci beaucoup !

[Jerome Briot]

1/ J'aimerai pouvoir "déplacer" ma polaire de telle sorte à être "angle au vent" (0°= face au vent) plutôt qu’être au cap GPS comme c'est le cas. Ici, mon vent vient du 330, et j'aimerai qu'il soit au niveau du zéro.

Comme ceci ?

Code :

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

polar(a-330, b, '.');

2/ J'ai codé une recherche d'un maximum, mais je ne sais pas comment le superposer avec mes premiers points.

Comme ceci ?

Code :

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

1
2
3
polar(a,b)
hold on
plot(MAXb*cos(MAXa), MAXb*sin(MAXa), 'b*')

Tu comprendras mieux si tu ais apparaitre l'objet Axes qui "se cache" derrière :

Code :

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

1
2
3
figure
polar(...,...)
axis on

[MarionSCN]

Salut ! Merci de ta réponse !

Pour le fait de tourner la polaire, j'ai réussi en modifiant un à un les termes car "a" et "b" sont des vecteurs colonnes de 13 000 lignes. Mais du coup c'était bien ça, merci.
J'ai même mis une box pour entrer le cap relevé et modifier tout automatiquement.

Pour garder mon max sur ma polaire, ça marche aussi, j'ai juste dû modifier de cette manière:

Code :

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

1
2
3
4
5
6
7
8
 
polar(a,b,'.')
grid
view(90,-90)
hold on
plot(MAXb*cos(MAXa), MAXb*sin(MAXa), 'r*')
view(90,-90)
hold off

Sinon, je perdais l'orientation de ma polaire.

Cela "interférait" avec ma courbe cartésienne, j'ai essayé "hold off" à la fin de l’exécution de mon bouton Tracer Polaire et ça marche

Merci de ton aide !