104文章解读与程序——西南交通大学学报《基于氢储能的热电联供型微电网优化调度方法》已提供下载资源

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

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

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

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

针对质子交换膜燃料电池和电解槽的热电联供特性，为避免氢能系统的热能浪费并进一步提高氢能系统的效率，搭建了一种考虑氢能系统的热电联供型光伏/风机/燃料电池/蓄电池/电锅炉/燃气锅炉微电网系统，提出一种包括日前调度与实时优化的两阶段优化调度方法. 所建系统考虑了电氢转换时的余热回收，将氢能系统作为热电氢耦合设备，实现了电、热、氢能的协调利用与相互转换，有效提高了能量利用率. 在第一阶段调度中，根据日前的风光发电出力及负荷需求预测，以微电网整体运行成本最小为目标，采用混合整数线性规划方法实现日前最优全局调度；在第二阶段调度中，根据超短期预测结果，使用模型预测控制嵌入混合整数二次规划算法，减小预测误差带来的经济性影响.

部分代码展示：

%% 初始化导入数据
clc;
clear;
close all;
tic
%% 24小时数据
pbat=[     0.0000    1.9205         0   -0.0000         0         0   -0.0000   11.9166         0         0         0    0.0000   22.2494    1.8446    0   -0.0000   -0.0000    0.0000         0         0         0         0         0  -30.7242];
phss=[  1.5300   -7.6450   -7.6839    3.6000    5.9000    9.9000   11.5750   -7.6666   14.7000   10.8000   14.6500   12.7750  -13.0994   17.8554   20.5000   24.9250   27.2500   24.5000    5.7750   -1.8177   -7.6500   -5.6250   -9.5225   39.0492 ];
qtst=[   -6.3976    0.0000         0  -13.0500   -0.0000   -0.0000    0.0000    2.4499    0.6250    1.2500    2.6125    3.1437   10.1394    2.1138  1.6750    0.0000    0.0000   -0.0000         0         0   -0.0000    0.0000    0.0000    0.0000];
soc=[0.5,0.500000000000035,0.505761456256818,0.505761456256808,0.505761456256809,0.505761456256812,0.505761456256825,0.505761456256826,0.541511165559164,0.541511165559166,0.541511165559180,0.541511165559193,0.541511165559207,0.608259296788454,0.613793189405949,0.613793189405962,0.613793189405968,0.613793189405973,0.613793189406014,0.613793189406028,0.613793189406033,0.613793189406033,0.613793189406033,0.613793189406046];
sohc=[0.5,0.507598999999959,0.354190037037037,0.200000000000000,0.217880000000002,0.247183333333322,0.296353333333322,0.353842499999992,0.200000000000000,0.273010000000000,0.326650000000000,0.399411666666667,0.462860833333317,0.200000000000000,0.288681666666613,0.390498333333280,0.514292499999960,0.649634166666641,0.771317499999974,0.800000000000000,0.763525666339920,0.610015666339919,0.497140666339870,0.306055833006481];
hoc=[0.5,0.476305000009099,0.476305000009105,0.476305000009114,0.427971666675758,0.427971666675769,0.427971666675772,0.427971666675785,0.435321405047890,0.437196405047914,0.440946405047927,0.448783905047930,0.458215155047941,0.488633473154265,0.494974999999908,0.499999999999921,0.499999999999925,0.499999999999941,0.499999999999945,0.499999999999945,0.499999999999958,0.499999999999947,0.499999999999961,0.499999999999974];
Qbl=[  17.2199   10.9500   10.9500   10.9500   23.6250   16.2250   14.5063   10.9500   10.9500   10.9500   10.9500   10.9500   10.9500   10.9500 0.9500   11.2688   12.6875   17.0750   23.1563   23.3407   18.1150   19.5375   10.9500   15.2377];
Qeb=[0    6.5695    6.8345         0    0.0000         0         0         0         0         0         0         0         0         0   0         0         0   -0.0000   -0.0000    0.6234         0         0    4.4797         0];
peb=Qeb/0.9;
%% 96数组 
A=[pbat;phss;peb;qtst;soc;sohc;hoc;Qbl;Qeb];
B=zeros(9,96);
for i1=1:9
    for j1=1:24
        for k1=1:4
            B(i1,(j1-1)*4+k1)=A(i1,j1);
        end
    end
end
pbata=[B(1,:) 0 0 0 0];
phssa=[B(2,:)  0 0 0 0];
peba=[B(3,:)   0 0 0 0];
Qtsta=[B(4,:)  0 0 0 0];
soca=[B(5,:)  0.5 0.5 0.5 0.5];
sohca=[B(6,:)  0.5 0.5 0.5 0.5];
hoca=[B(7,:)  0.5 0.5 0.5 0.5];
Qbla=[B(8,:)   14.6 14.6 14.6 14.6];
Qeba=[B(8,:)   0 0  0   0 ];
% RES及LOAD实测出力 96组数据
ppv=[0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0.944356427378540,2.24688011934705,3.67596375164204,5,6.11217021720757,7.40481046425189,9.39504547917028,12.6000000000000,17.2819627037912,22.6838780236454,27.7938543316769,31.6000000000000,33.4443539676276,34.0846774411666,34.6326621941223,36.2000000000000,39.5624964256983,44.1524122116882,49.0661218918340,53.4000000000000,56.4744103295793,58.5056737120806,59.9341002385416,61.2000000000000,62.6336122559847,64.1248929399896,65.4537271539997,66.4000000000000,66.8223906464819,66.8947545279611,66.8697411454598,67,67.4580751580878,68.0960889481659,68.6860582641611,69,68.8515587211668,68.2208896793751,67.1297757978959,65.6000000000000,63.6825649572449,61.5453523343336,59.3854635442555,57.4000000000000,55.7181814498535,54.1977009832905,52.6283700250822,50.8000000000000,48.5728342433409,46.0888437325045,43.5604313554158,41.2000000000000,39.1404815767828,37.1969240866917,35.1049045532547,32.6000000000000,29.4746144495280,25.7484599207289,21.4980754315653,16.8000000000000,11.8079356251054