✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。
🍎个人主页：Matlab科研工作室
🍊个人信条：格物致知。
更多Matlab仿真内容点击👇
⛄ 内容介绍
以三峡水库的发电优化调度问题为例，其兴利库容165 亿m 3 ，库容系数0.04，具有不完全年调节能力。水库正常蓄水位175 m，防洪限制水位145 m，枯季消落低水位155 m。电站32台机组，总装机容量2250万kW，保证出力499万kW，下游航运及生态基流要求5 600 m 3/s。选择调度期为1年，计算时段为旬，调度期初水位为174 m，调度期末水位也控制为173 m，调度期内的逐旬入流过程为模拟入库过程。
⛄ 部分代码
function x=init(i,x,D,Q,qmin,dt,Vmin,Vmax)
%设置函数默认值
if nargin==2
D=37;
%各旬初的最低水位与最高水位
Vmin=[228,228,228,228,228,228,228,228,228,228,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,171.5,228];
Vmax=[393,393,393,393,393,393,393,393,393,393,393,393,393,393,344,228,179.12,179.12,179.12,179.12,179.12,179.12,179.12,179.12,179.12,300.2,300.2,393,393,393,393,393,393,393,393,393,393];
%各旬入库流量
Q=[5323 5351 4467 3967 4000 4382 4984 5037 4962 6194 7709 5354 6993 11091 12918 14778 15566 17060 ...
27816 32798 33851 23446 20049 14796 17779 21637 17266 11754 10443 9019 7725 7244 6161 5424 5471 5250];
%各旬的最小下泄流量
qmin=[6000 6000 6000 6000 6000 6000 6000 6000 6000 6000 6000 6000 5700 5700 5700 5700 5700 5700 ...
5700 5700 5700 5700 5700 5700 5700 5700 10000 10000 8000 8000 5700 5700 5700 5700 5700 5700];
%各旬的时段长，h
dt=[240 240 264 240 240 192 240 240 264 240 240 240 240 240 264 240 240 240 ...
240 240 264 240 240 264 240 240 240 240 240 264 240 240 240 240 240 240];
end
for j=2:1:D-1
%计算本时段末水库最大蓄水量，亿立方米
% Vmaxj=x(i,j-1)+3600*(Q(j-1))*dt(j-1)/10^8;
% %计算蓄水量上限
% if Vmaxj<Vmin(j)
% Vmaxj=Vmax(j);
% else
% Vmaxj=min(Vmaxj,Vmax(j));
% end
%随机产生本时段初蓄水量
x(i,j)=Vmin(j)+rand*(Vmax(j)-Vmin(j));
end
end
⛄ 运行结果
⛄ 参考文献
[1]钟平安, 张卫国, 张玉兰,等. 水电站发电优化调度的综合改进差分进化算法[J]. 水利学报, 2014, 45(10):1147-1155.
[2]王福岭. 基于协同差分演化算法的梯级水库短期发电优化调度研究[J]. 中国农村水利水电, 2013(6):4.
