Matlab|计及需求响应和电能交互的多主体综合能源系统主从博弈优化调度策略

目录

主要内容

部分代码

结果一览

下载链接


主要内容

程序建立了多主体综合能源模型,采用双层模型进行求解,上层用自适应粒子群算法求解出各能源售价和需求响应补偿价格;下层采用混合整数规划算法求解出三个园区、配电网、储能电站、集中型风电场间的最优调度策略,下层程序采用matlab+cplex求解(也可替换为gurobi),程序注释清楚,方便学习!注意:程序不是完全复现,部分参考文献内容(基本模型部分、目标函数和约束条件),未按照主从博弈机制进行循环求解。

部分代码

复制代码
%参数
Ns=2;%场景数为2
ps=[0.3 0.7];%每种场景概率
Nerss=1;%储能电站个数
Nbt=1;%电化学储能电站个数
PIbt=2.3;%电化学储能单位时间折旧成本
NG=6;%常规机组数量
NDG=1;%DG数量
kc=2.5;%弃风惩罚
socmin=0.1;
socmax=0.9;
Pech=-50;%充电功率限制
Pedis=50;%放电功率限制
Se=200;%储能容量限制
Se0=100;
Pwmin=-50;%蓄水电站最小功率
Pwmax=50;%蓄水电站最大功率
Vs0=100;
Vpmin=0;%蓄水容量最小值
Vpmax=200;%蓄水容量最大值
Ca=100;Cb=125;Cc=150;Cd=150;%IDR补偿成本系数
mp=[1 150 50 0.0375 20 372.5 72 2;%常规机组参数
    2 60  20  0.175 17.5 352.3 48 2;
    3 60  15  0.625 10 316.5  30 2;
    4 50  10  0.0834 32.5 329.2 30 2;
    5 40  10  0.25  30  276.4 18 2;
    6 45  12  0.25 30 232.2  24 2];
%pload=[390 382 354 340 328 409 460 460 472 511 465 458 442 456 457 471 475 503 512 469 445 438 410 362];%负荷
​
idg=[1 2 5 8 11 13];
iw=2;
ie=2;
is=8;
%变量
%场景1变量
thetaone=sdpvar(30,T,'full');
Pgone=sdpvar(NG,T,'full');
Pwone=sdpvar(1,T,'full');
upwone=binvar(1,T,'full');
Perssone=sdpvar(1,T,'full');
uerssone=binvar(1,T,'full');
Pdgone=sdpvar(1,T,'full');
Pidrbone=sdpvar(30,T,'full');
Pidrcone=sdpvar(30,T,'full');
plossone=sdpvar(30,T,'full');
%场景2变量
thetatwo=sdpvar(30,T,'full');
Pgtwo=sdpvar(NG,T,'full');
Pwtwo=sdpvar(1,T,'full');
upwtwo=binvar(1,T,'full');

结果一览

下载链接

相关推荐
xrgs_shz42 分钟前
MATLAB读入不同类型图像并显示图像和相关信息
图像处理·计算机视觉·matlab
晶台光耦3 小时前
光耦——光伏储能的隐形动力,引领绿色能源迈向新纪元
能源·光耦·光耦应用·光耦选型·光伏储能
追梦的小猴子3 小时前
智慧监测数据集成平台:破孤岛,促互联,提数据价值强化履约,保质量,降成本为水务、城市、能源等领域发展添动力
智慧城市·能源
火山口车神丶4 小时前
某车企ASW面试笔试题
c++·matlab
今天吃饺子4 小时前
2024年SCI一区最新改进优化算法——四参数自适应生长优化器,MATLAB代码免费获取...
开发语言·算法·matlab
Evand J16 小时前
集合卡尔曼滤波(Ensemble Kalman Filter),用于二维滤波(模拟平面上的目标跟踪),MATLAB代码
matlab·平面·目标跟踪
lucy1530275107918 小时前
【青牛科技】带 ALC 双通道前置放大器电路D3308
科技·能源·安防·开关电源·白色家电
我是博博啦21 小时前
matlab例题
人工智能·算法·matlab
2402_871321951 天前
MATLAB方程组
gpt·学习·线性代数·算法·matlab
安科瑞刘鸿鹏1 天前
光伏运维监测系统在“双碳”环境中的应用
运维·服务器·网络·物联网·能源