程序简介
本程序参考电力系统自动化文献《含冷热电联供及储能的区域综合能源系统运行优化》，下面具体对文章和代码进行介绍。
文章使用场景分析法对可再生电源出力随机性进行建模，以揭示可再生能源的随机波动性对含有冷热电联供的区域综合能源系统的运行优化所带来的影响。在考虑系统内设备特性后，建立了包含能源转换设备、能源储存设备在内的区域综合能源系统模型，并采用Yalmip工具箱和商业软件Cplex在MATLAB中进行了仿真分析，求得了系统内各个单元最佳出力、机组组合和不同调度模式下总运行成本。
综合能源系统框架：
新能源不确定性处理
新能源随机性建模步骤：
风电预测概率结果：
文中的目标函数为运行成本最小，目标函数和约束条件可以参考论文，提供的数据很全面，文中较为详细的介绍了蓄电池的等式和不等式约束，感兴趣的同学可以看一下。
结果对比
部分程序
%%束条件% 机组出力上下限约束Constraints = [];for k = 1:HorizonConstraints = [Constraints, onoff(:,k).*Pmin <= P(:,k) <= onoff(:,k).*Pmax];% indicator will be 1 only when switched onindicator = onoff(unit,k)-onoff(unit,k-1);range = k:min(Horizon,k+minup(unit)-1);% Constraints will be redundant unless indicator = 1Constraints = [Constraints, onoff(unit,range) >= indicator];s% indicator will be 1 only when switched offindicator = onoff(unit,k-1)-onoff(unit,k);range = k:min(Horizon,k+mindown(unit)-1);% Constraints will be redundant unless indicator = 1Constraints = [Constraints, onoff(unit,range) <= 1-indicator]; end% 蓄电池约束% 蓄电池初始能量soc_initial = 7000; % ESS存储能量下限soc_min = 5000;% ESS存储能量上限soc_max = 8500;% 储热约束% 储热初始能量hsoc_initial = 2000;% 储热存储能量下限hsoc_min = 1000;% 储热存储能量上限hsoc_max = 6000;% 储热能量限制Constraints = [Constraints, hsoc_max >= -sum(P(6,1:k)) + hsoc_initial >= hsoc_min];Constraints = [Constraints, -sum(P(6,1:k)) == 0];% 功率平衡约束% 热Constraints = [Constraints, (P(1,k) + P(2,k) + P(3,k))*0.37/0.34 + P(4,k)*0.9*0.8 + P(6,k)*0.9*0.9 + P(8,k)>= he(:,k)];% CCHP燃料费用fuelcost1 = fuelcost1 + Q(1)*P(1,k)*P(1,k) + C(1)*P(1,k) + D(1);fuelcost2 = fuelcost2 + Q(2)*P(2,k)*P(2,k) + C(2)*P(2,k) + D(2);fuelcost3 = fuelcost3 + Q(3)*P(3,k)*P(3,k) + C(3)*P(3,k) + D(3)% 交换功率费用Cgrid = 0;Cgrid = Cgrid + P(7,k)*cec(k);% 储能运行费用Cess = 0; Cess = Cess + abs(P(5,k))*0.001 + abs(P(6,k))*0.0010000;% 燃气锅炉费用Cgb = 0;Cgb = Cgb + P(4,k)*0.073/(0.9*0.7);fuelcost = fuelcost1 + fuelcost2 + fuelcost3;O = fuelcost + Cgrid + Cess + Cgb;%% 目标函数Objective = 0;
部分内容源自网络，侵权联系删除！
欢迎感兴趣的小伙伴点击文末阅读原文获取完整版代码，小编会不定期更新高质量的学习资料、文章和程序代码，为您的科研加油助力！
更多创新性预测及优化程序请点击公众号首页：《预测优化》|《综合能源》链接！
