power22.m

来自「发电厂竞价」· M 代码 · 共 20 行

pd=[1022 963 978 993 1037 1066 1095 1124 1138 1154 1182 1197 1242 1271 1314 1372 1314 1271 1285 1290 1372 1459 1459 1081];
p0 = [0.1; 0.1;0.1;0.1;0.1;0.1;0.1;0.1;0.1;0.1];% Make a starting guess at solution
aeq=[1 1 1 1 1 1 1 1 1 1];
lb=[30;20;25;50;50;15;120;75;250;75];
ub=[100;80;120;150;150;60;320;200;520;280];
for i=1:length(pd);
beq=pd(i);
[p,fval,maxfval,exitflag]  = fminimax(@myfun,p0,[],[],aeq,beq,lb,ub);
DLMWRITE('PowerMarket1.doc',i,'-append',  'precision', '%.3f', 'newline', 'pc');
DLMWRITE('PowerMarket1.doc','p','-append',  'precision', '%.3f', 'newline', 'pc');
DLMWRITE('PowerMarket1.doc',p,'-append',  'precision', '%.3f', 'newline', 'pc');
DLMWRITE('PowerMarket1.doc','C','-append',  'precision', '%.3f', 'newline', 'pc');
DLMWRITE('PowerMarket1.doc',maxfval,'-append',  'precision', '%.3f', 'newline', 'pc');
DLMWRITE('PowerMarket1.doc',' ','-append',  'precision', '%.3f', 'newline', 'pc');
%disp(i);
%disp('p');
%disp(p);
%disp('Com');
%disp(maxfval);
end;