Courbes Matlab

[invite3b879c03]

Bonjour,

je voudrais tracer une courbe T=f(h) sous matlab avec T=T0 - a*h. Cette courbe je veux que lorsque h appartient à un intervalle spécifique je veux que le coefficient a change lorsque h est compris entre telle et telle valeur. Je veux également avoir sur le meme graphique toutes ces courbes bout à bout. Je voudrais obtenir cette courbe en fait pour être plus précis: http://www.rfi.fr/actufr/images/107/..._ozone_432.gif

Voici mon code qui reste expérimental (je débute sous Matlab):

function[T] =myode(h)
T0=288.16;
if h>=0 && h<=11000
a=-0.0065;
elseif h>11000 && h<=20000
a=0.0;
elseif h>20000 && h<=32000
a=1.0;
elseif h>32000 && h<=47000
a=2.8;
elseif h>47000 && h<=51000
a=0.0;
elseif h>51000 && h<=71000
a=-2.8;
elseif h>71000 && h<=86000
a=-2.0;
end
T=T0-a*h;
plot(T,h);
hold on
end

Merci pour votre aide.
-----

[fiatlux]

Salut

Tu as écrit la fonction, il ne reste plus qu'à l'utiliser. Par exemple en affichant la valeur de T pour des h entre 0 et 86000. Tu pourrais faire ça pour h=0, h=1, h=2, ... mais c'est inutile, ça fait beaucoup trop de points, tu peux par exemple y aller par pas de 100: h=0, h=100, h=200, ... , h=86000

figure, hold on;
for h=0:100:86000
    myode(h);
end

Par contre si tu fais ça (essaie), tu verras que tu n'obitendras pas la courbe que tu voulais, mais là ce n'est plus un problème avec Matlab, mais avec tes équations de segments de droite.

[invite3b879c03]

Ce que tu m'a donné je le met dans la fonction ?

Pourtant ca m'affiche ca dans la commande

Error using myode (line 3)
Not enough input arguments.

[fiatlux]

non non tu le mets pas dans la fonction. Tu sauves ta fonction telle quelle dans un fichier myode.m.
Puis sois tu appelles cette fonction directement de la fenetre de commande Matlab, soit tu te crées un script (un simple fichier .m) avec juste ces 4 lignes et tu l'exécute (Debug-->Run, ou simplement l'icone avec la flèche verte, ou encore plus simplement F5 si t'es sur windows).

[A voir en vidéo sur Futura]

[invite3b879c03]

Est ce qu'il existe une expression qui me récupère la dernière valeur T(h) de la courbe d'avant. Effectivement j'ai des courbes mais je ne vois a quoi cela correspond.

[fiatlux]

Oui. Histoire de coder un peu plus proprement:

altitudes = 0:100:86000; % les altitudes à afficher: 0, 100, 200, ..., 86000
T = zeros(size(altitudes)); % leur T équivalent, pour tous initialisés à zero
figure, hold on;
for i=1:numel(altitudes)
    h = altitudes(i);
    T(i) = myode(h);
end

[invite3b879c03]

Par exemple récupérer la valeur T(11000). est ce possible?

[invite3b879c03]

Oula je ne comprend plus grand chose là. Mon expression mathématique est fausse. En fait il faudrais que toutes ces courbes se rejoignent en fait. Je dois rajouter un offset dans mon expression mathématique.

[invite3b879c03]

En fait l'expression plus exact serait

Tm=Tm,b + Am,b *(h-hb)

avec b variant de 0à6 (correspondant au segment de courbe)
Je connais T0, A0, h0 et après les A1, A2....,A6 c'est les coefficient que j'ai rentré dans mes conditions IF.

C'est récurrent donc je vois pas comment faire ca.

[fiatlux]

Etant donné que ta courbe complète est un ensemble de segments, à mon avis la façon la plus simple de faire ça c'est de d'afficher segment par segment, et donc de faire une fonction qui plot 1 segment, puis de l'appeler N fois pour les N segments.

1. Créer une fonction qui prend en entrée:
1.a) les paramètres du segment d'un segment (pente, offset)
2.a) un vecteur avec les valeur de h dans ce segment
Donc par exemple, dans un fichier plotSegment.m :

function [] = plotSegment(h, Tmb, Amb, hb) %ici h est un vecteur, genre 11000:100:20000
Tm=Tmb + Amb *(h-hb);
plot(Tm,h)

2. Dans un autre fichier .m :

figure, hold on;
plotSegment(0:100:11000, Tmb1, Amb1, hb1)
plotSegment(11000:100:20000, Tmb2, Amb2, hb2)
etc.