Discussion :

[Invité]

Bonjour,

J'aimerais résoudre une équation a 2 inconnues (oui comme ça c'est difficile) mais en rajoutant une condition qui est la minimisation d'un rapport entre ces 2 inconnus.

En gros je voudrais trouver L et H solution de : L*H/(L+H)= 30, donc là il peut y avoir un paquet de solutions mais je veux retenir juste celle qui minimise L*H et la si j'ai un peu de chance il n'y en a qu'une.

Le problème, c'est que je ne voit pas comment le programmer, je sais même pas par quoi commencer ... je n'ai jamais joué avec les résolution d'équation dans matlab, je sais pas de quoi il est capable.

[Invité]

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
function [Lopt,hopt]=optimisation_section(Dh,maxL,maxh,type)
 
h=linspace(3,maxh,100);
L=zeros(1,length(h));
air=zeros(1,length(h));
 
for i=1:length(h)
    L(i)=Dh*h(i)/(2*h(i)-Dh);
    if L(i)>maxL ||L(i)<0 , L(i)=NaN; end
    air(i)=L(i)*h(i);
end
 
save air
idx=find(air == min(air));
 
Lopt=L(idx);
hopt=h(idx);

Bon j'ai fais ça, ça marche mais bon le résultat est limité à la résolution que je donne au vecteur. Il y a surement plus intelligent chez matlab

[chris05]

ou simplement :

Code :

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

1
2
3
4
5
6
7
8
9
 
H=[0:100000];  %a modifier
L=(30*H)./(H-30); %calcul de L tel que LH/(L+H)=30
prod=L.*H; %calcul de L*H
 
[m,n]=find(prod(:)==min(prod(:)));
 
disp('H=');H(n)
disp('L=');L(n)

Il faut néanmoins définir un intervalle pour H (ici entre 0 et 100000).

[Invité]

C'est pareil sans la boucle nan ?

[Klinnmas]

si tu as la toolbox optimisation, il te suffit de créer une fonction (f(L,H) = L*H/ (L+H), par exemple) et utiliser une des fonctions qui trouvent un minimum sous contraintes (par exemple fmincon). Avec cette fonction, tu lui mets des bornes (lb et ub), un point initial et la contrainte non linéaire "f" dans l'emplacement nonlcon à l'aide de @, et piouf , normalement ça marche.
Tu peux aller voir l'aide, c'est pas mal fait pour cette fonction.

Si tu as des problèmes d'utilisation, n'hésite pas !!

bon courage.

[Invité]

effectivement j'ai cette toolbox, je vais regarder