Discussion :

[amineyamane]

bonjour;

Je suis amine, étudiant en maitrise en sciences appliquées
je travail actuellement sur l'optimisation d'énergie dans un véhicule à pile à combustible épaulée par une batterie. J'ai formulé le problème, et j'essaye de l'implémenter sur l'algorithme génétique en utilisant matlab. Mon problème est que l'algorithme ne me satisfait pas les contraintes non-linéaires.
merci d'avance
salutations

amine

[progfou]

Salut !

Ah là là, encore une question mal formulée...

Si tu nous disais quel problème tu as précisément ?

[amineyamane]

En fait j’ai deux sources de puissance (une pile et une batterie) qui doivent satisfaire un profile de puissance donné.
Mon objectif est de minimiser la consommation à partir de la pile (minimiser le courant de la pile Ifc) en maintenant la charge (soc) de la batterie entre deux bornes inf et sup (connues) durant un profile donne, en plus que la somme des puissances de ou vers la batterie Pb pour tout le profile demandé doit être égale à zéro.

Voila l’algorithme d’optimisation :

function y= fonction_objectif(x)

Pdem=dlmread('C:\profiles de vitesse\pdemhw.txt'); %lecture du profile
tf=length(Pdem);
for i=1:tf-1

y=y+(Mh2*Nst*x(i))/(2*F);
end
********************************************

function [c ceq]=contraintes(x)

for i=1:tf-1
Pbp(i)=Pdem(i)-(Nst*pile(x(i))*x(i));
[Vbp(i) Ibp(i) T_batp(i+1) Rintp(i+1) T_extp(i+1) C_maxp(i+1) Pmaxp(i) Ibat_maxp(i) vocp(i+1)]=batterie(socp(i),T_batp(i),Pbp(i),T_extp(i));
socp(i+1)=socp(i)-Ibp(i)/(C_maxp(i)*3600);


c=[

socp(i+1)-0.8;
0.4-socp(i+1);

Pbp(i)-Pmaxp(i);
Ibat_maxp(i)-Ibp(i);



];



ceq = [

sum(Pbp)==0;
Pbp(i)==Pdem(i)-(Nst*pile(x(i))*x(i));
];
end


quand j'appelle l'algo génétique de cette façon :

ObjectiveFunction = @fonction_objectif;

nvars = tf-1;% Number of variables
LB(1:tf-1)=5;
UB(1:tf-1)=150;
options = gaoptimset('MutationFcn',@mutationadaptfeasible,'PopulationSize',150);

[x,fval] = ga(ObjectiveFunction,nvars,[],[],[],[],LB,UB,...
ConstraintFunction,options);
j'obtiens toujour x(1:tf-1)= 5 sans satifaisant les contraintes non-lineaire

[Garulfo]

Avant de donner du code, tu dois donner les algorithmes, c'est-à-dire la conception, de ton problème.

Si celui-ci est correct alors c'est que ton erreur vient du code. Et dans ce cas c'est le forum Matlab et non Algo qui est pertinent. Si ton algo a une erreur, ça ne sert à rien de le coder.

Donc : quelle est ta conception détaillée svp ?

[amineyamane]

salut,

le problème à optimiser est comme suite :

on a 2 sources d'energie : une pile (source principale Pfc) et une batterie (source auxiliaire Pbat). les 2 sources doivent participer pour satisfaire la demande.
la demande est une série de valeur formant un profile de puissance connue apriori.donc on a à chaque instant t:
Pdem(t) = Pfc(t) + Pbat(t)
l'objectif est de minimiser la quantité totale consommée a partir de la pile durant un profile de puissance donné. Cette quantité est calculée par :

J(Ifc)=M.Nst.Ifc/2.F

M,NST,F sont des constantes.

en cherche à minimiser cette fonction tout en gardant la charge de la deuxieme sources(batterie) entre 2 bornes inf et sup en plus la somme des puissances de ou vers la batterie doit etre egale à zéro à la fin du profile + d'autre contrainte sur les limites de chaque grandeurs (courants max puissances max...).

l'algorithme complet est comme suite :

Objective:
minimize(somme (J(Ifc))) i=1....length(profile)

sujet à


Pdem(t) = Pfc(t) + Pbat(t)

SOC(t + 1)= SOC(t) -Ib/Cmax
SOC(final) = SOC(start)
SOC(t) < SOCmax
SOC(t) > SOCmin
Pfc(t)<Pfcmax
Pfc(t)>Pfcmin
Pbat(t)<Pbatmax
Pbat(t)>Pbatmin

donc si la longueur du profile demander est tf --> j'ai tf-1 variables
je veux que l'algo génétique trouve la combinaison optimale entre mes variables qui minimise la dite quantité et qui satisfait mes contraintes

si vous avez besoin plus de détails on peut les discuter.

merci d'avance
amine

[Alp]

A quel endroit précisément de la conception de tes algos génétiques rencontres-tu un problème ?

[amineyamane]

salut,

le probleme c au niveau des contraintes, en fait l'algorithme ne satisfait aucune contrainte.

merci

[kindyroot]

salut,

le probleme c au niveau des contraintes, en fait l'algorithme ne satisfait aucune contrainte.

merci

Ta fonction "fitness" prend-elle en compte toutes tes contraintes? normalement tu devrais faire une moyenne des indices de satisfaction de toutes les contraintes et classer tes individus selon la valeur de cette fonction, tu peux même rajouter des coefficients pour privilègier la satisfaction de certaines contraintes, ... parle nous un peu de ta fonction fitness!

[amineyamane]

bonjours;
Mon objectif est de minimiser une quantité proportionnelle à une variable:

f=min(somme(constante*x))

j'ai une suite des valeurs représentants des puissances demandées (Pdem), que je dois satisfaire l'une après l'autre.
la langueur de x est égale aux nombres de puissances demandée.

la minimisation de la fonction f doit respecter les contraintes suivantes :

pour i de 1 à nbpuissances

soc(x(i))<0.8
soc(x(i))>0.4
Pdem(i)=V*x(i)+Pb(i)
fin

somme(Pb)=0

merci
amine