140文章解读与程序——《虑多微网电能互补与需求响应的微网双层优化模型》已提供下载资源-CSDN博客

本文链接：https://blog.csdn.net/LIANG674027206/article/details/140998269

👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆

《《《《《《《《更多资源还请持续关注本专栏》》》》》》》

论文与完整源程序_电网论文源程序的博客-CSDN博客https://blog.csdn.net/liang674027206/category_12531414.html

电网论文源程序-CSDN博客电网论文源程序擅长文章解读,论文与完整源程序,等方面的知识,电网论文源程序关注python,机器学习,计算机视觉,深度学习,神经网络,数据挖掘领域.https://blog.csdn.net/LIANG674027206?type=download

代码主要做的是考虑多微网电能互补共享的微网双层优化模型，同时优化配电网运营商的动态电价以及微网用户的能量管理策略，在上层，目标函数为配电网运营商的收益最大化，决策变量为配电网运营商的交易电价；在下层，目标函数为各个用户微网最小化运行成本，实现配电网-微网双赢；

部分代码展示：

clc;clear all;
yalmip('clear');
warning off all;             %清空环境
%微网1中光伏发电预测：
%Ppv1=[0 0 0 0 0 0 5.80 10.04 15.00 35.56 45.99 56.06 58.09 55.87 45.77 36.08 35.23 18.04 15.79 0 0 0 0 0 ];
ppv1=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1.45000000000000,2.90000000000000,4.35000000000000,5.80000000000000,6.86000000000000,7.92000000000000,8.98000000000000,10.0400000000000,11.2800000000000,12.5200000000000,13.7600000000000,15,20.1400000000000,25.2800000000000,30.4200000000000,35.5600000000000,38.1675000000000,40.7750000000000,43.3825000000000,45.9900000000000,48.5075000000000,51.0250000000000,53.5425000000000,56.0600000000000,56.5675000000000,57.0750000000000,57.5825000000000,58.0900000000000,57.5350000000000,56.9800000000000,56.4250000000000,55.8700000000000,53.3450000000000,50.8200000000000,48.2950000000000,45.7700000000000,43.3475000000000,40.9250000000000,38.5025000000000,36.0800000000000,35.8675000000000,35.6550000000000,35.4425000000000,35.2300000000000,30.9325000000000,26.6350000000000,22.3375000000000,18.0400000000000,17.4775000000000,16.9150000000000,16.3525000000000,15.7900000000000,11.8425000000000,7.89500000000000,3.94750000000000,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];
%微网1中交流负荷预测：
%pload1=[67.5500000000000,66.9350000000000,66.3200000000000,65.7050000000000,65.0900000000000,64.5375000000000,63.9850000000000,63.4325000000000,62.8800000000000,64.6825000000000,66.4850000000000,68.2875000000000,70.0900000000000,71.9550000000000,73.8200000000000,75.6850000000000,77.5500000000000,77.9350000000000,78.3200000000000,78.7050000000000,79.0900000000000,80.7575000000000,82.4250000000000,84.0925000000000,85.7600000000000,85.8400000000000,85.9200000000000,86,86.0800000000000,86.2800000000000,86.4800000000000,86.6800000000000,86.8800000000000,87.9300000000000,88.9800000000000,90.0300000000000,91.0800000000000,90.8325000000000,90.5850000000000,90.3375000000000,90.0900000000000,86.5675000000000,83.0450000000000,79.5225000000000,76,75.7700000000000,75.5400000000000,75.3100000000000,75.0800000000000,75.3075000000000,75.5350000000000,75.7625000000000,75.9900000000000,76.7325000000000,77.4750000000000,78.2175000000000,78.9600000000000,81.8000000000000,84.6400000000000,87.4800000000000,90.3200000000000,92.3525000000000,94.3850000000000,96.4175000000000,98.4500000000000,100.310000000000,102.170000000000,104.030000000000,105.890000000000,106.940000000000,107.990000000000,109.040000000000,110.090000000000,111.565000000000,113.040000000000,114.515000000000,115.990000000000,112.012500000000,108.035000000000,104.057500000000,100.080000000000,96.7000000000000,93.3200000000000,89.9400000000000,86.5600000000000,83.6925000000000,80.8250000000000,77.9575000000000,75.0900000000000,73.8250000000000,72.5600000000000,71.2950000000000,56.0715876422348,58.6795056124762,66.2685292102058,60.9894187065234];
%pload1=15.*pload1;
%微网2中光伏发电预测：
%Ppv2=[0 0 0 0 0 0 5.88 10.09 13.90 15.13 20.09 22.99 25.86 30.08 33.98 30.31 20.04 20.56 18.43 12.33 0 0 0 0];
ppv2=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1.47000000000000,2.94000000000000,4.41000000000000,5.88000000000000,6.93250000000000,7.98500000000000,9.03750000000000,10.0900000000000,11.0425000000000,11.9950000000000,12.9475000000000,13.9000000000000,14.2075000000000,14.5150000000000,14.8225000000000,15.1300000000000,16.3700000000000,17.6100000000000,18.8500000000000,20.0900000000000,20.8150000000000,21.5400000000000,22.2650000000000,22.9900000000000,23.7075000000000,24.4250000000000,25.1425000000000,25.8600000000000,26.9150000000000,27.9700000000000,29.0250000000000,30.0800000000000,31.0550000000000,32.0300000000000,33.0050000000000,33.9800000000000,33.0625000000000,32.1450000000000,31.2275000000000,30.3100000000000,27.7425000000000,25.1750000000000,22.6075000000000,20.0400000000000,20.1700000000000,20.3000000000000,20.4300000000000,20.5600000000000,20.0275000000000,19.4950000000000,18.9625000000000,18.4300000000000,16.9050000000000,15.3800000000000,13.8550000000000,12.3300000000000,9.24750000000000,6.16500000000000,3.08250000000000,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];

%微网2中交流负荷预测：
%Load2=[99.08 101.22 93.04 102.09 101.66 100.33 108.34 82.00 85.98 90.21 115.09 130.77 161.90 143.78 158.21 143.90 125.97 120.09 114.41 115.90 140.98 135.88 108.09 100.09];
%pload2=[99.0800000000000,99.6150000000000,100.150000000000,100.685000000000,101.220000000000,99.1750000000000,97.1300000000000,95.0850000000000,93.0400000000000,95.3025000000000,97.5650000000000,99.8275000000000,102.090000000000,101.982500000000,101.875000000000,101.767500000000,101.660000000000,101.327500000000,100.995000000000,100.662500000000,100.330000000000,102.332500000000,104.335000000000,106.337500000000,108.340000000000,101.755000000000,95.1700000000000,88.5850000000000,82,82.9950000000000,83.9900000000000,84.9850000000000,85.9800000000000,87.0375000000000,88.0950000000000,89.1525000000000,90.2100000000000,96.4300000000000,102.650000000000,108.870000000000,115.090000000000,119.010000000000,122.930000000000,126.850000000000,130.770000000000,138.552500000000,146.335000000000,154.117500000000,161.900000000000,157.370000000000,152.840000000000,148.310000000000,143.780000000000,147.387500000000,150.995000000000,154.602500000000,158.210000000000,154.632500000000,151.055000000000,147.477500000000,143.900000000000,139.417500000000,134.935000000000,130.452500000000,125.970000000000,124.500000000000,123.030000000000,121.560000000000,120.090000000000,118.670000000000,117.250000000000,115.830000000000,114.410000000000,114.782500000000,115.155000000000,115.527500000000,115.900000000000,122.170000000000,128.440000000000,134.710000000000,140.980000000000,139.705000000000,138.430000000000,137.155000000000,135.880000000000,128.932500000000,121.985000000000,115.037500000000,108.090000000000,106.090000000000,104.090000000000,102.090000000000,98.2296843488321,107.740185810594,107.193887566695,98.0609888837085];
%Ppv3=[0,0,0,0,0,0,11.6,20.08,30,71.12,91.98,112.12,116.18,111.74,91.54,72.16,70.46,36.08,31.58,0,0,0,0,0];
ppv3=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,2.90000000000000,5.80000000000000,8.70000000000000,11.6000000000000,13.7200000000000,15.8400000000000,17.9600000000000,20.0800000000000,22.5600000000000,25.0400000000000,27.5200000000000,30,40.2800000000000,50.5600000000000,60.8400000000000,71.1200000000000,76.3350000000000,81.5500000000000,86.7650000000000,91.9800000000000,97.0150000000000,102.050000000000,107.085000000000,112.120000000000,113.135000000000,114.150000000000,115.165000000000,116.180000000000,115.070000000000,113.960000000000,112.850000000000,111.740000000000,106.690000000000,101.640000000000,96.5900000000000,91.5400000000000,86.6950000000000,81.8500000000000,77.0050000000000,72.1600000000000,71.7350000000000,71.3100000000000,70.8850000000000,70.4600000000000,61.8650000000000,53.2700000000000,44.6750000000000,36.0800000000000,34.9550000000000,33.8300000000000,32.7050000000000,31.5800000000000,23.6850000000000,15.7900000000000,7.89500000000000,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0];
%Load3=[109.88 ,108.05 ,128.48 ,188.85 ,109.63 ,147.82 ,199.62 ,138.86 ,89.60 ,124.17 ,199.38 ,171.14 ,265.81 ,175.12 ,177.90 ,229.92 ,189.70 ,157.82 ,172.39 ,125.40 ,150.48 ,217.24 ,187.51 ,136.36];
%pload3=[109.880000000000,109.422500000000,108.965000000000,108.507500000000,108.050000000000,113.157500000000,118.265000000000,123.372500000000,128.480000000000,143.572500000000,158.665000000000,173.757500000000,188.850000000000,169.045000000000,149.240000000000,129.435000000000,109.630000000000,119.177500000000,128.725000000000,138.272500000000,147.820000000000,160.770000000000,173.720000000000,186.670000000000,199.620000000000,184.430000000000,169.240000000000,154.050000000000,138.860000000000,126.545000000000,114.230000000000,101.915000000000,89.6000000000000,98.2425000000000,106.885000000000,115.527500000000,124.170000000000,142.972500000000,161.775000000000,180.577500000000,199.380000000000,192.320000000000,185.260000000000,178.200000000000,171.140000000000,194.807500000000,218.475000000000,242.142500000000,265.810000000000,243.137500000000,220.465000000000,197.792500000000,175.120000000000,175.815000000000,176.510000000000,177.205000000000,177.900000000000,190.905000000000,203.910000000000,216.915000000000,229.920000000000,219.865000000000,209.810000000000,199.755000000000,189.700000000000,181.730000000000,173.760000000000,165.790000000000,157.820000000000,161.462500000000,165.105000000000,168.747500000000,172.390000000000,160.642500000000,148.895000000000,137.147500000000,125.400000000000,131.670000000000,137.940000000000,144.210000000000,150.480000000000,167.170000000000,183.860000000000,200.550000000000,217.240000000000,209.807500000000,202.375000000000,194.942500000000,187.510000000000,174.722500000000,161.935000000000,149.147500000000,145.284338727498,122.862021788933,135.532009916422,116.740515720376];
jjg=[0.0520000000000000,0.0520000000000000,0.0520000000000000,0.0520000000000000,0.0560000000000000,0.0560000000000000,0.0560000000000000,0.0560000000000000,0.0620000000000000,0.0620000000000000,0.0620000000000000,0.0620000000000000,0.0700000000000000,0.0700000000000000,0.0700000000000000,0.0700000000000000,0.0800000000000000,0.0800000000000000,0.0800000000000000,0.0800000000000000,0.0840000000000000,0.0840000000000000,0.0840000000000000,0.0840000000000000,0.0920000000000000,0.0920000000000000,0.0920000000000000,0.0920000000000000,0.102000000000000,0.102000000000000,0.102000000000000,0.102000000000000,0.106000000000000,0.106000000000000,0.106000000000000,0.106000000000000,0.0930000000000000,0.0930000000000000,0.0930000000000000,0.0930000000000000,0.0940000000000000,0.0940000000000000,0.0940000000000000,0.0940000000000000,0.0800000000000000,0.0800000000000000,0.0800000000000000,0.0800000000000000,0.0780000000000000,0.0780000000000000,0.0780000000000000,0.0780000000000000,0.0820000000000000,0.0820000000000000,0.0820000000000000,0.0820000000000000,0.0740000000000000,0.0740000000000000,0.0740000000000000,0.0740000000000000,0.0810000000000000,0.0810000000000000,0.0810000000000000,0.0810000000000000,0.0930000000000000,0.0930000000000000,0.0930000000000000,0.0930000000000000,0.166000000000000,0.166000000000000,0.166000000000000,0.166000000000000,0.170000000000000,0.170000000000000,0.170000000000000,0.170000000000000,0.144000000000000,0.144000000000000,0.144000000000000,0.144000000000000,0.0930000000000000,0.0930000000000000,0.0930000000000000,0.0930000000000000,0.0950000000000000,0.0950000000000000,0.0950000000000000,0.0950000000000000,0.0780000000000000,0.0780000000000000,0.0780000000000000,0.0780000000000000,0.0580000000000000,0.0580000000000000,0.0580000000000000,0.0580000000000000];

pload=load('pload.txt');%读取负荷数据
pload1=pload(:,1);  %读取微网1负荷
pload1=pload1';%转置后备用
pload2=pload(:,2);%读取微网2负荷
pload2=pload2';%转置后备用
pload3=pload(:,3);%读取微网3负荷
pload3=pload3';%转置后备用
ppv1=5*ppv1;%取微网1光伏数据
ppv2=3*ppv2;%取微网2光伏数据
ppv3=3*ppv3;%取微网3光伏数据



%微网1中光伏发电预测：
%Ppv1=[0 0 0 0 0 0 5.80 10.04 15.00 35.56 45.99 56.06 58.09 55.87 45.77 36.08 35.23 18.04 15.79 0 0 0 0 0 ];
%ppv1=[0 0 0 0 0 0 5.80 10.04 15.00 35.56 45.99 56.06 58.09 55.87 45.77 36.08 35.23 18.04 15.79 0 0 0 0 0 ];
%微网1中交流负荷预测：
%pload1=[67.55 65.09 62.88 70.09 77.55 79.09 85.76 86.08 86.88 91.08 90.09 76.00 75.08 75.99 78.96 90.32 98.45 105.89 110.09 115.99 100.08 86.56 75.09 70.03];
%微网2中光伏发电预测：
%Ppv2=[0 0 0 0 0 0 5.88 10.09 13.90 15.13 20.09 22.99 25.86 30.08 33.98 30.31 20.04 20.56 18.43 12.33 0 0 0 0];
%ppv2=[0 0 0 0 0 0 5.88 10.09 13.90 15.13 20.09 22.99 25.86 30.08 33.98 30.31 20.04 20.56 18.43 12.33 0 0 0 0];

%微网2中交流负荷预测：
%Load2=[99.08 101.22 93.04 102.09 101.66 100.33 108.34 82.00 85.98 90.21 115.09 130.77 161.90 143.78 158.21 143.90 125.97 120.09 114.41 115.90 140.98 135.88 108.09 100.09];
%pload2=[99.08 101.22 93.04 102.09 101.66 100.33 108.34 82.00 85.98 90.21 115.09 130.77 161.90 143.78 158.21 143.90 125.97 120.09 114.41 115.90 140.98 135.88 108.09 100.09];
%Ppv3=[0,0,0,0,0,0,11.6,20.08,30,71.12,91.98,112.12,116.18,111.74,91.54,72.16,70.46,36.08,31.58,0,0,0,0,0];
%ppv3=[0,0,0,0,0,0,11.6,20.08,30,71.12,91.98,112.12,116.18,111.74,91.54,72.16,70.46,36.08,31.58,0,0,0,0,0];
%Load3=[109.88 ,108.05 ,128.48 ,188.85 ,109.63 ,147.82 ,199.62 ,138.86 ,89.60 ,124.17 ,199.38 ,171.14 ,265.81 ,175.12 ,177.90 ,229.92 ,189.70 ,157.82 ,172.39 ,125.40 ,150.48 ,217.24 ,187.51 ,136.36];
%pload3=[109.88 ,108.05 ,128.48 ,188.85 ,109.63 ,147.82 ,199.62 ,138.86 ,89.60 ,124.17 ,199.38 ,171.14 ,265.81 ,175.12 ,177.90 ,229.92 ,189.70 ,157.82 ,172.39 ,125.40 ,150.48 ,217.24 ,187.51 ,136.36];
%jjg=[0.0520000000000000,0.0560000000000000,0.0620000000000000,0.0700000000000000,0.0800000000000000,0.0840000000000000,0.0920000000000000,0.102000000000000,0.106000000000000,0.0930000000000000,0.0940000000000000,0.0800000000000000,0.0780000000000000,0.0820000000000000,0.0740000000000000,0.0810000000000000,0.0930000000000000,0.166000000000000,0.170000000000000,0.144000000000000,0.0930000000000000,0.0950000000000000,0.0780000000000000,0.0580000000000000];
%价格
jg=[0.0468000000000000,0.0468000000000000,0.0572000000000000,0.0468000000000000,0.0504000000000000,0.0504000000000000,0.0504000000000000,0.0504000000000000,0.0558000000000000,0.0558000000000000,0.0558000000000000,0.0558000000000000,0.0630000000000000,0.0630000000000000,0.0630000000000000,0.0630000000000000,0.0720000000000000,0.0720000000000000,0.0720000000000000,0.0720000000000000,0.0756000000000000,0.0756000000000000,0.0756000000000000,0.0756000000000000,0.0828000000000000,0.0828000000000000,0.0828000000000000,0.0828000000000000,0.0918000000000000,0.0918000000000000,0.0918000000000000,0.0918000000000000,0.0954000000000000,0.0954000000000000,0.0954000000000000,0.0954000000000000,0.0837000000000000,0.0837000000000000,0.0837000000000000,0.0837000000000000,0.0846000000000000,0.0846000000000000,0.0846000000000000,0.0846000000000000,0.0720000000000000,0.0720000000000000,0.0720000000000000,0.0720000000000000,0.0702000000000000,0.0702000000000000,0.0702000000000000,0.0702000000000000,0.0738000000000000,0.0738000000000000,0.0738000000000000,0.0738000000000000,0.0666000000000000,0.0666000000000000,0.0666000000000000,0.0666000000000000,0.0729000000000000,0.0729000000000000,0.0729000000000000,0.0729000000000000,0.0837000000000000,0.0837000000000000,0.0837000000000000,0.102300000000000,0.182600000000000,0.182600000000000,0.182600000000000,0.182600000000000,0.187000000000000,0.187000000000000,0.161600000000001,0.187000000000000,0.158400000000000,0.129600000000000,0.129600000000000,0.129600000000000,0.0837000000000000,0.0837000000000000,0.0837000000000000,0.0837000000000000,0.0855000000000000,0.0855000000000000,0.0855000000000000,0.0855000000000000,0.0702000000000000,0.0702000000000000,0.0702000000000000,0.0702000000000000,0.0522000000000000,0.0522000000000000,0.0522000000000000,0.0522000000000000];
ytac=0.95;ytad=0.95;cn(1)=50;cn(2)=50;cn(3)=50;%与储能效率有关的系数
pchmax=10;pchmin=0;pdismax=10;pdismin=0;
jrfh=[0.5.*pload1,0.2.*pload2,0.2.*pload3];%柔性负荷
jgfh=[0.5.*pload1,0.8.*pload2,0.8.*pload3];%固定负荷
cdis=0.1;cch=0.1;
tm=96;
%jg=sdpvar(1,tm);
pch=sdpvar(1,3*tm);%储能充电
pdis=sdpvar(1,3*tm);%储能放电
bch=binvar(1,3*tm);%表征储能是否充电的0-1变量

效果展示：