72文章解读与程序——SCI基于合作型Stackerlberg博弈的考虑差别定价和风险管理的微网运行策略》已提供下载资源

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

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

论文与完整源程序_电网论文源程序的博客-CSDN博客icon-default.png?t=N7T8https://blog.csdn.net/liang674027206/category_12531414.html

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

论文可参考:

A cooperative Stackelberg game based energy management considering price discrimination and risk assessment

代码主要做的是一个基于主从博弈的考虑差别定价和风险管理的微网动态定价与调度策略,构建了双层能源管理框架,上层为零售商的动态定价模型,目标是社会福利最大化;下层是多个产消者的合作博弈模型,优化各产消者的能量管理策略,各产消者之间可以进行P2P交易。同时,采用纳什谈判法对多个产消者的合作剩余进行公平分配,还考虑了运行风险,采用条件风险价值(CVaR)随机规划方法来描述零售商的预期损失.

部分代码展示:

%% [SCI文章复现]A cooperative Stackelberg game based energy management considering price discrimination and risk assessment
%[中译]基于合作型Stackerlberg博弈的考虑差别定价和风险管理的微网运行策略
%International Journal of Electrical Power and Energy Systems,SCI二区
%Highlights:合作型Stackerlberg博弈,纳什谈判,差别定价
%此函数用于计算不包含产消者合作交易时的成本C_Non,设置交易量为零,即P_trading_max=0

clc
clear
close all

%% 模型参数设定
%产消者/零售商从主网购电价格 元/MW
u_Db=1e3*[0.4,0.4,0.4,0.4,0.4,0.4,0.79,0.79,0.79,1.2,1.2,1.2,1.2,1.2,0.79,0.79,0.79,1.2,1.2,1.2,0.79,0.79,0.4,0.4];
%产消者/零售商向主网售电价格 元/MW
u_Ds=1e3*[0.35,0.35,0.35,0.35,0.35,0.35,0.68,0.68,0.68,1.12,1.12,1.12,1.12,1.12,0.68,0.68,0.68,1.12,1.12,1.12,0.79,0.79,0.35,0.35];
%零售商与产消者的交易价格上下限
u_Pbmax=1e3*[0.7,0.7,0.7,0.7,0.7,0.7,1.1,1.1,1.1,1.5,1.5,1.5,1.5,1.5,1,1,1,1.5,1.5,1.5,1.1,1.1,0.7,0.7];%购价上限
u_Pbmin=u_Pbmax-0.5*1e3*ones(1,24);%购价下限
u_Psmax=u_Ds;%售价上限
u_Psmin=u_Psmax-0.35*1e3*ones(1,24);%售价下限
%产消者1-3  电负荷 MW
P_load_1=[6.62295082,5.770491803,5.442622951,5.31147541,5.37704918,5.573770492,6.295081967,6.491803279,7.213114754,7.803278689,8.131147541,8.131147541,7.93442623,7.278688525,7.016393443,7.016393443,7.147540984,8.262295082,9.442622951,9.37704918,9.37704918,7.93442623,6.819672131,5.901639344];
P_load_2=[3.344262295,3.016393443,2.754098361,2.754098361,2.754098361,2.885245902,3.147540984,3.344262295,3.639344262,3.93442623,4,4.131147541,4,3.737704918,3.475409836,3.606557377,3.606557377,4.131147541,4.721311475,4.655737705,4.721311475,4,3.409836066,3.016393443];
P_load_3=[11.60655738,10.16393443,9.442622951,9.245901639,9.114754098,9.639344262,10.75409836,11.3442623,12.45901639,13.50819672,14.10772834,14.16393443,13.63934426,12.72131148,12.19672131,12.32786885,12.59016393,14.29508197,16.59016393,16.45901639,16.26229508,13.7704918,12.13114754,10.55737705];
%产消者1-3  导入10个场景的出力和概率
Sw=10; %场景数量
load P_Gen.mat  %产消者1风电出力    P_Gen_1  维度:10*24     P_Gen_2    P_Gen_3 
%产消者1-3风电场景概率
pai_1=0.1*ones(1,10);pai_2=0.1*ones(1,10);pai_3=0.1*ones(1,10);
%其它固定参数
C_E=80; %储能充放成本
P_Pbmax=15; %最大购电量
P_Psmax=15; %最大售电量
Cap=10; %最大储能容量MW
P_Ecmax=3; %充放能功率上限
P_Edmax=3; %充放能功率上限
SOCmin=0.2; %最小存储量百分比 单位%
SOCmax=0.85; %最大容量百分比
SOCini=0.33; %初始容量百分比
SOCexp=0.85; %末段容量百分比
M=1E8; %大M法
beta=0.1; %厌恶风险系数
%% 决策变量初始化
delta=sdpvar(1,3);
eta_1=sdpvar(Sw,1); %产消者1的风险调度辅助变量
eta_2=sdpvar(Sw,1); %产消者2的风险调度辅助变量
eta_3=sdpvar(Sw,1); %产消者3的风险调度辅助变量
P_Ps_1=sdpvar(Sw,24); %零售商向产消者1售能量
P_Ps_2=sdpvar(Sw,24); %零售商向产消者2售能量
P_Ps_3=sdpvar(Sw,24); %零售商向产消者3售能量
P_Pb_1=sdpvar(Sw,24); %零售商从产消者1购能量
P_Pb_2=sdpvar(Sw,24); %零售商从产消者2购能量
P_Pb_3=sdpvar(Sw,24); %零售商从产消者3购能量
u_Ps=sdpvar(3,24); %零售商向产消者购能价格
u_Pb=sdpvar(3,24); %零售商从产消者购能价格
P_trading_1=sdpvar(Sw,24); %产消者1合作博弈交易量
P_trading_2=sdpvar(Sw,24); %产消者2合作博弈交易量
P_trading_3=sdpvar(Sw,24); %产消者3合作博弈交易量
SOC_1=sdpvar(Sw,24); %产消者1储能容量状态 单位%
SOC_2=sdpvar(Sw,24); %产消者2储能容量状态 单位%
SOC_3=sdpvar(Sw,24); %产消者3储能容量状态 单位%
P_Ec_1=sdpvar(Sw,24); %储能充电
P_Ec_2=sdpvar(Sw,24); %储能充电
P_Ec_3=sdpvar(Sw,24); %储能充电
P_Ed_1=sdpvar(Sw,24); %储能放电
P_Ed_2=sdpvar(Sw,24); %储能放电
P_Ed_3=sdpvar(Sw,24); %储能放电
Uabs_1=binvar(Sw,24); %储能充放电状态,0-1变量
Uabs_2=binvar(Sw,24); %储能充放电状态,0-1变量
Uabs_3=binvar(Sw,24); %储能充放电状态,0-1变量
Urelea_1=binvar(Sw,24); %储能充放电状态,0-1变量
Urelea_2=binvar(Sw,24); %储能充放电状态,0-1变量
Urelea_3=binvar(Sw,24); %储能充放电状态,0-1变量
%定义KKT条件中的拉格朗日乘子
lamda_pro_1=sdpvar(Sw,24);
lamda_pro_2=sdpvar(Sw,24);
lamda_pro_3=sdpvar(Sw,24);
lamda_trading=sdpvar(Sw,24);
lamda_Pb_1=sdpvar(Sw,24);
lamda_Pb_2=sdpvar(Sw,24);
lamda_Pb_3=sdpvar(Sw,24);
lamda_Ps_1=sdpvar(Sw,24);
lamda_Ps_2=sdpvar(Sw,24);
lamda_Ps_3=sdpvar(Sw,24);
lamda_Ec_1=sdpvar(Sw,24);
lamda_Ec_2=sdpvar(Sw,24);
lamda_Ec_3=sdpvar(Sw,24);
lamda_Ed_1=sdpvar(Sw,24);
lamda_Ed_2=sdpvar(Sw,24);
lamda_Ed_3=sdpvar(Sw,24);
lamda_SOCmin_1=sdpvar(Sw,24);
lamda_SOCmin_2=sdpvar(Sw,24);
lamda_SOCmin_3=sdpvar(Sw,24);
lamda_SOCmax_1=sdpvar(Sw,24);
lamda_SOCmax_2=sdpvar(Sw,24);
lamda_SOCmax_3=sdpvar(Sw,24);
lamda_SOC1_1=sdpvar(Sw,24);
lamda_SOC1_2=sdpvar(Sw,24);
lamda_SOC1_3=sdpvar(Sw,24);%公式17和18共用
lamda_SOC2_1=sdpvar(Sw,1);
lamda_SOC2_2=sdpvar(Sw,1);

效果展示:

72号资源-源程序:论文可在知网下载《基于合作型Stackerlberg博弈的考虑差别定价和风险管理的微网运行策略》资源-CSDN文库icon-default.png?t=N7T8https://download.csdn.net/download/LIANG674027206/89101340👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆下载资源链接👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆👆

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

论文与完整源程序_电网论文源程序的博客-CSDN博客icon-default.png?t=N7T8https://blog.csdn.net/liang674027206/category_12531414.html

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

相关推荐

最近更新

  1. docker php8.1+nginx base 镜像 dockerfile 配置

    2024-04-13 17:42:05       94 阅读
  2. Could not load dynamic library ‘cudart64_100.dll‘

    2024-04-13 17:42:05       100 阅读
  3. 在Django里面运行非项目文件

    2024-04-13 17:42:05       82 阅读
  4. Python语言-面向对象

    2024-04-13 17:42:05       91 阅读

热门阅读

  1. 探索未来学术:ChatGPT如何引领论文写作革命

    2024-04-13 17:42:05       32 阅读
  2. jvm中堆与栈的区别详细讲解

    2024-04-13 17:42:05       39 阅读
  3. 2024-04-09 问AI: 在深度学习中,dropout层是什么?

    2024-04-13 17:42:05       37 阅读
  4. python项目练习——22、人脸识别软件

    2024-04-13 17:42:05       36 阅读
  5. 手机重启手app没了

    2024-04-13 17:42:05       39 阅读
  6. 数组名的指针用法

    2024-04-13 17:42:05       34 阅读