Discussion :

[shakes06]

Bonjour,

Voici mon script Matlab (il est aussi en pièce jointe), ce sera plus facile d'expliquer ensuite mon problème :

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
Tw=450; % température de la surface chauffante en Kelvin
R0=[2.5e-3:2e-4:5e-3]; % initialisation du rayon de la goutte en mètre
tau=[0]; % initialisation du temps caractéristique d'évaporation en seconde
R=[0]; % initialisation du rayon en m
 
figure(1)
for j=1:length(R0)
   tau(j)=(2^(11/4))*[(3*mu)/(rho_v*g)]^(1/4)*[(Lv*rho*a)/(k*(Tw-Td))]^(3/4)*(R0(j)^(1/2));
   for i=1:length(t)
      R(i,j)=R0(j)*[1-(t(i)/tau(j))]^2; % calcul du rayon
    end
    hold all
    plot(t(:),R(:,j));
    cc(j)=cellstr(num2str(R0(j))); %légende de chaque courbe tracée
end
legend(cc)

Dans ce script, on fixe la température du mur Tw et on fait varier les rayons initiaux R0. Pour chaque R0, on regarde l'évolution du rayon R en fonction du temps t et on trace cette évolution.

Mon problème est le suivant : pour chaque courbe tracée, je cherche à écrire (dans la légende) R0 = ... m

J'ai pensé à inscrire quelque part dans le script : 'R0 = %f m', num2str(R0(j)) mais ça ne fonctionne pas !
Je ne sais ni comment écrire ceci proprement, ni même où l'insérer dans le script...

Serait-il possible d'avoir un peu d'aide s'il vous plaît ?

Merci par avance.

[Invité]

Bonsoir et bienvenue sur le forum,

Il suffit de faire une simple concaténation:

Code :

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

cc(j)=cellstr(['R0 = ' num2str(R0(j))]);

Ou sans l'inclure dans la boucle:

Code :

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

cc = cellfun(@(x) ['R0 = ' num2str(x)],num2cell(R0),'UniformOutput',false)

[Jerome Briot]

La même chose en utilisant STRCAT :

Code :

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

cc = strcat({'R0 = '},num2str(R0(:)));

Sinon, sous MATLAB, les crochets ont une signification différente de celle des mathématiques. Il ne faut s'en servir que pour la concaténation de valeurs :

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
Tw = 450; % température de la surface chauffante en Kelvin
R0 = 2.5e-3:2e-4:5e-3; % initialisation du rayon de la goutte en mètre
R = zeros(length(t),length(R0)); % initialisation du rayon en m
 
figure(1)
hold on
 
for j=1:length(R0)
 
    tau = 2^(11/4)*((3*mu)/(rho_v*g))^(1/4)*((Lv*rho*a)/(k*(Tw-Td)))^(3/4)*(R0(j)^(1/2));
 
    for i=1:length(t)
        R(i,j)=R0(j)*(1-(t(i)/tau))^2; % calcul du rayon
    end
 
    plot(t(:),R(:,j));
 
end
 
cc = strcat({'R0 = '},num2str(R0(:)));
legend(cc)

[shakes06]

Merci beaucoup à vous deux pour vos réponses ultra-rapides ! Mon script fonctionne maintenant !

Du coup, il me reste une question...
Je cherche à tracer deux fonctions différentes tau1 et tau2 sur le même graphe mais il faut que les intervalles des fonctions ne se recoupent pas...

Voici le script pour être plus claire :

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
R0=2e-5;            % rayon initial de la goutte en m
Tw1=[323:1:372];     % initialisation du vecteur température de la surface en Kelvin
Tw2=[373:1:1000]
tau1=[0];            % initialisation du temps caractéristique d'évaporation en s
tau2=[0];
 
figure;
for i=1:length(Tw1)
    tau1(i)=abs([2/(3*770*((Tw1(i)-Td)^2.78))]*rho*Lv*[(3*a*(R0^2))^(1/3)])
end
 
for j=1:length(Tw2)
    tau2(j)=(2^(11/4))*[(3*mu)/(rho_v*g)]^(1/4)*[(Lv*rho*a)/(k*(Tw2(j)-Td))]^(3/4)*(R0^(1/2));
end
 
plot(Tw1,tau1,'b')
hold on
plot(Tw2,tau2,'g')
hold off

Dans ce sript :
tau1 en fonction de Tw1 ne doit être tracée que pour Tw1 compris entre 0K et 372K, et les ordonnées doivent être comprises entre 0 et 50 sec
tau2 en fonction de Tw2 tracée pour Tw2 au-delà de 373K et les ordonnées doivent être comprises entre 0 et 100 sec

J'arrive à restreindre l'axe des abscisses mais je ne vois pas comment "couper" les fonctions en ordonnée... Je vous mets en pièces jointes ce que j'obtiens, mon script et ce que je souhaite obtenir...

En vous remerciant à nouveau par avance !

[Invité]

Bonjour,

Si ce n'est que pour l'affichage, une solution serait d'utiliser la fonction fplot.
Je te laisse regarder la documentation doc fplot

[shakes06]

Merci beaucoup pour votre aide !