Discussion :

[pen_duick]

Bonjour,

Dans le cadre de mon stage, j'ai créé une fonction à deux variables (Longueur et rayon) qui me permet d'obtenir en sortie la masse totale d'un réservoir, la quantité de LH2 qui y est stocké ainsi que la longueur totale de ce réservoir:

Code :

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

function [Mtotal, MLH2, Ltotal] = Tank_sizingLRtest(L,r)

Je souhaite réutiliser cette fonction dans un autre programme. Dans ce nouveau programme, MLH2 sera fixé, et je souhaite calculer la longueur et le rayon optimums (avec 1<r<4 et 50<L<80) pour avoir Mtotal minimum. Je difinis également un intervalle d'incertitude sur cette valeur de MLH2 (+ ou - 100 kg)

Voici le script que j'ai commencé à écrire:

Code :

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

1
2
3
4
5
6
7
8
9
10
MLH2=52296;
for r=1:4
    for L=50:80  
            [Mtotal, MLH2, Ltotal] = Tank_sizingLRtest(L,r)
        if MLH2-100<MLH2<MLH2+100       
            disp(L) 
            disp(r)
        end  
    end
end

Etant novice sur ce logiciel, je me sens un peu perdu et un peu d'aide ne serait pas de refus.

Merci par avance,

[Jerome Briot]

Stocke les valeurs de chaque itération dans des tableaux :

Code :

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

1
2
3
4
5
for r = 1:4
    for L = 50:80  
        [Mtotal(r,L-49), MLH2(r,L-49), Ltotal(r,L-49)] = Tank_sizingLRtest(L,r);
    end
end

Tu pourras visualiser la surface de réponse de la masse totale :

Code :

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

1
2
figure
surf(50:80, 1:4, Mtotal)

Ensuite, tu crées un masque binaire avec la valeur de MLH2 fixée :

Code :

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

1
2
MLH2_fixed = 52296;
idx = MLH2>(MLH2_fixed-100) & MLH2<(MLH2_fixed+100);

Tu mets toutes les valeurs de Mtotal qui ne correspondent pas à ce masque à +Inf :

Code :

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

Mtotal(~idx(:)) = Inf;

Tu cherches le minimum :

Code :

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

[Mmin, k] = min(Mtotal(:));

Tu récupères les valeurs de r et de L correspondantes :

Code :

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

1
2
 
[Lmin,rmin] = ind2sub(size(Mtotal), k)

[pen_duick]

Merci pour avoir répondu aussi rapidement. J'ai pu tester votre solution et le programme ne semble pas donner les bons résultats. En effet, les valeurs de Lmin et rmin obtenues ne sont pas celles qui correspondent à Mtotal minimum.

J'ai regardé en détail le code, celui-ci détecte parfaitement Motal minimum pour MLH2 compris dans l'intervalle souhaité. Le problème réside dans la récupération des valeurs de r et L correspondantes.

Le problème se situerait donc à la dernière ligne du code ?

[pen_duick]

Je pense avoir trouvé l'erreur, voici ma proposition :

Code :

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

1
2
3
[rmin,Lmin] = ind2sub(size(Mtotal), k);
Lmini=Lmin+49
rmini=rmin

Inversion de rmin et Lmin dans la matrice ainsi que rajout des 49 soustraits précédemment ? Ca semble fonctionner..

[pen_duick]

Serait-il possible également d'ajouter un message indiquant, lorsque cela est nécessaire, que les bornes de L et de r initialement choisies ne permettent pas de trouver une masse MLH2 dans l'intervalle choisi ?

[Jerome Briot]

Serait-il possible également d'ajouter un message indiquant, lorsque cela est nécessaire, que les bornes de L et de r initialement choisies ne permettent pas de trouver une masse MLH2 dans l'intervalle choisi ?

Il suffit de tester si il y a au moins un 1 logique dans le vecteur idx :

Code :

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

1
2
3
4
5
6
idx = MLH2>(MLH2_fixed-100) & MLH2<(MLH2_fixed+100);
if any(idx(:))
   % recherche du minimum
else
   fprintf('Pas de masse  MLH2 dans l intervalle choisi');
end

[pen_duick]

Merci, le programme marche bien maintenant. Cependant, je rencontre un autre problème:

En effet, je souhaiterais que r varie de 1 à 2 par pas de 0.1.

J'ai donc modifié ma boucle for en conséquence, mais une erreur apparaît puisque les indices r et L ne peuvent seulement être des réels positifs, ou des entiers..

[Jerome Briot]

Deux possibilités

Code :

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

1
2
3
4
5
6
7
8
9
kr = 0;
for r = 1:0.1:2
    kr = kr+1;
    kL = 0;
    for L = 50:80  
        kL = kL+1;
        [Mtotal(kr,kL), MLH2(kr,kL), Ltotal(kr,kL)] = Tank_sizingLRtest(L,r);
    end
end

ou

Code :

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

1
2
3
4
5
6
7
8
r = 1:0.1:2;
L = 50:80;
 
for i = 1:numel(r)
    for j = 1:numel(L) 
        [Mtotal(i,j), MLH2(i,j), Ltotal(i,j)] = Tank_sizingLRtest(L(j),r(i));
    end
end

A toi de gérer les indices dans le reste du code…

[pen_duick]

Je vous remercie beaucoup pour la rapidité et la précision de vos réponses. Je mets ce sujet en "résolu".