مهندس ببخشید، این برنامه رو بصورت زیر تغییر دادم. مخرج کسر رو یک عدد ثابت گرفتم که همون ماکزیمم سود ائتلاف بین DG های یک MG هست. طبق متن بالای رابطه (15) گفتم احتمالن باید به این صورت باشه. نظر شما چیه؟ اون تکرار 1 تا 400 رو حذف کردم مشکلی توی الگوریتم بوجود نمیاد؟
clc; close all; clear;
global Q;global k;global D;global alpha;
alpha=.3; beta=1.2; k=.38; D=6.3;
V_opt=3.4713; % Tebghe jadvale 4
efficiency=zeros(1,10);
for numsupplier=2:10
Q=ones(1,numsupplier)*.5;
lb=zeros(1,numsupplier);
ub=ones(1,numsupplier);
[x1, V_ne]=fmincon(@fff,Q,[],[],[],[],lb,ub);
% for tekrar=1:400
% for supplier=1:numsupplier
% [x2, fx2]=fminbnd(@fff,0,1);
% Q(supplier)=x2;
% end
% end
efficiency(numsupplier)=V_ne/V_opt;
end
plot(2:10,efficiency(2:10).*(2:10)), grid
ylabel('Symmetric NE efficiency');
xlabel('Number of Suppliers in Microgrid')