探索云世界
基于条件风险价值CVaR的微网动态定价与调度策略（matlab代码）
2024-05-0851
版权
版权声明：
本文内容由阿里云实名注册用户自发贡献，版权归原作者所有，阿里云开发者社区不拥有其著作权，亦不承担相应法律责任。具体规则请查看《 阿里云开发者社区用户服务协议》和 《阿里云开发者社区知识产权保护指引》。如果您发现本社区中有涉嫌抄袭的内容，填写 侵权投诉表单进行举报，一经查实，本社区将立刻删除涉嫌侵权内容。
简介：基于条件风险价值CVaR的微网动态定价与调度策略（matlab代码）
1 主要内容
程序复现文章《A cooperative Stackelberg game based energy management considering price discrimination and risk assessment》，建立基于主从博弈的考虑差别定价和风险管理的微网动态定价与调度策略模型，构建了双层能源管理框架，上层为零售商的动态定价模型，目标是社会福利最大化；下层是多个产消者的合作博弈模型，优化各产消者的能量管理策略，各产消者之间可以进行P2P交易。同时，采用纳什谈判法对多个产消者的合作剩余进行公平分配，还考虑了运行风险，采用条件风险价值（CVaR）随机规划方法来描述零售商的预期损失。
该模型通过建立双层模型，出现了双线性和非线性问题，为了解决该问题，引入 kkt条件和大M法将原模型转化为等效单层模型，从而进行求解。而且程序采用三种对比算例进行分析，非常方便对照学习，算例1和2采用matlab+cplex求解，算例3采用matlab+mosek求解！
2 部分代码
%% 模型参数设定%产消者/零售商从主网购电价格 元/MWu_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];%产消者/零售商向主网售电价格 元/MWu_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 电负荷 MWP_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; %最大储能容量MWP_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条件中的拉格朗日乘子
3 程序结果
目录
相关文章
【2023高教社杯】D题 圈养湖羊的空间利用率 问题分析、数学模型及MATLAB代码
本文介绍了2023年高教社杯数学建模竞赛D题的圈养湖羊空间利用率问题，包括问题分析、数学模型建立和MATLAB代码实现，旨在优化养殖场的生产计划和空间利用效率。
【2022年华为杯数学建模】B题 方形件组批优化问题 方案及MATLAB代码实现
本文提供了2022年华为杯数学建模竞赛B题的详细方案和MATLAB代码实现，包括方形件组批优化问题和排样优化问题，以及相关数学模型的建立和求解方法。
【2023五一杯数学建模】 B题 快递需求分析问题 建模方案及MATLAB实现代码
本文介绍了2023年五一杯数学建模竞赛B题的解题方法，详细阐述了如何通过数学建模和MATLAB编程来分析快递需求、预测运输数量、优化运输成本，并估计固定和非固定需求，提供了完整的建模方案和代码实现。
车辆行驶控制运动学模型的matlab建模与仿真,仿真输出车辆动态行驶过程
该课题在MATLAB2022a中建立了车辆行驶控制运动学模型并进行仿真，展示车辆动态行驶过程。系统仿真结果包含四张图像，显示了车辆在不同时间点的位置和轨迹。核心程序定义了车辆参数和初始条件，使用ode45求解器模拟车辆运动。车辆运动学模型基于几何学，研究车辆空间位姿、速度随时间变化，假设车辆在平面运动且轮胎无滑动。运动学方程描述位置、速度和加速度关系，模型预测控制用于优化轨迹跟踪，考虑道路曲率影响，提升弯道跟踪性能。
地震波功率谱密度函数、功率谱密度曲线，反应谱转功率谱，matlab代码
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
耐震时程曲线，matlab代码，自定义反应谱与地震波，优化源代码，地震波耐震时程曲线
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
热门文章
最新文章
1
基于深度学习网络的十二生肖图像分类matlab仿真
58
2
m基于LDPC编译码的matlab误码率仿真,对比SP,MS,NMS以及OMS四种译码算法
211
3
基于直方图的图像曝光量分析FPGA实现,包含tb测试文件和MATLAB辅助验证
66
4
MATLAB热传导方程模型最小二乘法模型、线性规划对集成电路板炉温优化
42
5
基于DCT变换的彩色图像双重水印嵌入和提取算法matlab仿真
74
6
Matlab决策树、模糊C-均值聚类算法分析高校教师职称学历评分可视化
50
7
基于直方图的图像阈值计算和分割算法FPGA实现,包含tb测试文件和MATLAB辅助验证
67
8
m基于Yolov2深度学习网络的智能零售柜商品识别系统matlab仿真,带GUI界面
80
9
MATLAB改进模糊C均值聚类FCM在电子商务信用评价应用：分析淘宝网店铺数据|数据分享
47
10
MATLAB模糊C均值聚类FCM改进的推荐系统协同过滤算法分析MovieLens电影数据集
83