【完整解析】2024电工杯数学建模A题论文与代码

园区微电网风光储协调优化配置

1 论文
2 代码分享
- 2.1 第三题第一问
3 数据与代码

1 论文

在这里插入图片描述

2 代码分享

2.1 第三题第一问

function anss=q3w1ObjFun(ttt,id);
%tt=ttt(1);
tt=[750,0,0,1000,600,500];
lim=ttt(1)*200;
limm=ttt(2)*500*0.9-ttt(2)*500*0.1;
t1=ttt(3)*1000;t2=ttt(4)*1000;
ttt=ttt(5:28);
ned=[275 	241 	302 
275 	253 	292 
277 	329 	307 
310 	315 	293 
310 	290 	271 
293 	270 	252 
293 	307 	283 
380 	354 	223 
375 	264 	292 
281 	315 	283 
447 	313 	287 
447 	291 	362 
447 	360 	446 
405 	369 	504 
404 	389 	455 
403 	419 	506 
268 	412 	283 
313 	291 	311 
287 	379 	418 
288 	303 	223 
284 	331 	229 
287 	306 	361 
277 	285 	302 
275 	324 	291 ];
giv=[0.0000	0	0	0.2301	0.0000	0.1464
0.0000	0	0	0.3828	0.0000	0.2175
0.0000	0	0	0.2968	0.0000	0.3959
0.0000	0	0	0.4444	0.0000	0.1831
0.0000	0	0	0.5029	0.0000	0.4716
0.0000	0	0	0.3609	0.0000	0.6215
0.0000	0	0	0.2402	0.0000	0.2946
0.0058	0	0	0.0473	0.0105	0.1214
0.3026	0	0	0.1538	0.3280	0.0250
0.6020	0	0	0.1068	0.6314	0.3023
0.7711	0	0	0.0518	0.7936	0.0196
0.8555	0	0	0.2169	0.8925	0.1224
0.8531	0	0	0.3546	0.8999	0.3335
0.7842	0	0	0.2194	0.8221	0.2653
0.6437	0	0	0.1110	0.6667	0.1220
0.4242	0	0	0.2186	0.4275	0.1633
0.0619	0	0	0.3779	0.0216	0.2645
0.0000	0	0	0.3421	0.0000	0.3408
0.0000	0	0	0.5008	0.0000	0.3183
0.0000	0	0	0.4646	0.0000	0.3299
0.0000	0	0	0.2197	0.0000	0.1703
0.0000	0	0	0.1783	0.0000	0.1655
0.0000	0	0	0.1535	0.0000	0.1897
0.0000	0	0	0.0000	0.0000	0.2323];
ned=ned*1.5;
if(id==4)yuanqv_ned=0;yuanqv_giv=zeros(24,2);for i=1:3yuanqv_ned=yuanqv_ned+ned(:,i);yuanqv_giv(:,1)=yuanqv_giv(:,1)+(tt(i*2-1)+t1)*giv(:,i*2-1);yuanqv_giv(:,2)=yuanqv_giv(:,2)+(tt(i*2)+t2)*giv(:,i*2);end
elseyuanqv_ned=ned(:,id);yuanqv_giv(:,1)=(tt(id*2-1)+t1)*giv(:,id*2-1);yuanqv_giv(:,2)=(tt(id*2)+t2)*giv(:,id*2);
end
%yuanqv_giv(:,1)=yuanqv_giv(:,1)*t1;
%yuanqv_giv(:,2)=yuanqv_giv(:,2)*t1;
dianchi=0;
anss=0;
for i=1:24if(ttt(i)*dianchi>lim)anss=anss+10000;endnedi=yuanqv_ned(i)-yuanqv_giv(i,1)-yuanqv_giv(i,2)-ttt(i)*dianchi*0.95;if(nedi+ttt(i)*dianchi*0.95>0)dianchi=dianchi-ttt(i)*dianchi;endif(nedi<=0&&nedi*0.95>=-lim)dianchi=dianchi-nedi*0.95;elseif(nedi*0.95<=-lim)dianchi=dianchi+lim;endif(dianchi>limm)dianchi=limm;endif(nedi<0)anss=anss+0.4*yuanqv_giv(i,1)+0.5*yuanqv_giv(i,2);elseanss=anss+0.4*yuanqv_giv(i,1)+0.5*yuanqv_giv(i,2)+nedi;end
end
anss=anss+lim*800/3650+limm*1800/3650+t1*2500/(5*365)+t2*3000/(5*365);