MATLAB|【免费】高比例可再生能源电力系统的调峰成本量化与分摊模型

   

目录

   

主要内容   

  部分代码   

  结果一览   

下载链接

主要内容   

程序复现文献《高比例可再生能源电力系统的调峰成本量化与分摊模型》，从净负荷波动的角度出发，建立了调峰成本的量化与分摊模型，构造了无调峰需求的替代场景，将负荷和可再生能源出力曲线分别转换为无波动的均值线。其次，建立了含深度调峰和抽水蓄能的调度优化模型，用于计算不同场景下的调峰成本，并将有无调峰需求两种场景下的系统调峰成本之差作为单一主体导致的边际调峰成本，然后采用 Shapley值计算不同主体导致的调峰成本。最后，根据成本的引发程度分摊调峰成本。程序采用matlab+cplex编写，程序已经深度优化，求解速度很快，注释清晰，代码免费分享给大家学习参考！

  部分代码   

for t=1:24C=[C, sum(sum(P300_Diot(:,:,t)))+sum(sum(P600_Diot(:,:,t)))+sum(P200_it(:,t))+P_PV(t)+P_Wind(t)+sum(P_PuG(:,t)-P_PuS(:,t))+P_guding == L_jun(t)+D_Eload_FW(t)+D_Eload_JM(t)+D_Eload_ZZ(t),P_Wind(t) == 0;     P_PV(t) == 0;    ];                       
end                                 
%%  公式（8-9）系统正负旋转备用约束
u200Git = binvar(20,24,'full');   %uit 表示机组运行状态的0-1变量,1开机，0关机
u300Git = binvar(18,24,'full');  
u600Git = binvar(18,24,'full');  
R_Ureq = (Eload*0.05+E_Wind*0.05+E_PV*0.05);    %正旋转备用
R_Dreq = (Eload*0.05+E_Wind*0.05+E_PV*0.05);    %负旋转备用
for t=1:24         C=[C,   200*sum(u200Git(:,t))+300*sum(u300Git(:,t))+600*sum(u600Git(:,t))-sum(sum(P300_Diot(:,:,t)))-sum(sum(P600_Diot(:,:,t)))-sum(P200_it(:,t)) >= R_Ureq(t),   %正旋转备用约束sum(sum(P300_Diot(:,:,t)))+sum(sum(P600_Diot(:,:,t)))+sum(P200_it(:,t))-200*0.5*sum(u200Git(:,t))-300*0.35*sum(u300Git(:,t))-600*0.35*sum(u600Git(:,t)) >= R_Dreq(t),  %负旋转备用约束     ];      
end   
%%  公式（10-12）火电机组功率大小约束，这里换一种方法去写
uP300_Diot = binvar(18,4,24,'full');
uP600_Diot = binvar(18,4,24,'full');
for t= 1:24   for i = 1:18 C = [C,  sum(uP300_Diot(i,:,t))<= u300Git(i,t),  sum(uP600_Diot(i,:,t))<= u600Git(i,t),  u200Git(i,t)*200*0.5<=P200_it(i,t),P200_it(i,t)<=u200Git(i,t)*200, uP300_Diot(i,1,t)*300*0.5<=P300_Diot(i,1,t),P300_Diot(i,1,t)<=uP300_Diot(i,1,t)*300,uP300_Diot(i,2,t)*300*0.45<=P300_Diot(i,2,t),P300_Diot(i,2,t)<=uP300_Diot(i,2,t)*300*0.5,uP300_Diot(i,3,t)*300*0.4<=P300_Diot(i,3,t),P300_Diot(i,3,t)<=uP300_Diot(i,3,t)*300*0.45,uP300_Diot(i,4,t)*300*0.35<=P300_Diot(i,4,t),P300_Diot(i,4,t)<=uP300_Diot(i,4,t)*300*0.4,uP600_Diot(i,1,t)*600*0.5<=P600_Diot(i,1,t),P600_Diot(i,1,t)<=uP600_Diot(i,1,t)*600,uP600_Diot(i,2,t)*600*0.45<=P600_Diot(i,2,t),P600_Diot(i,2,t)<=uP600_Diot(i,2,t)*600*0.5,uP600_Diot(i,3,t)*600*0.4<=P600_Diot(i,3,t),P600_Diot(i,3,t)<=uP600_Diot(i,3,t)*600*0.45,uP600_Diot(i,4,t)*600*0.35<=P600_Diot(i,4,t),P600_Diot(i,4,t)<=uP600_Diot(i,4,t)*600*0.4,];endfor i = 19:20C = [C,  u200Git(i,t)*200*0.5<=P200_it(i,t),P200_it(i,t)<=u200Git(i,t)*200, ];end
end
%公式（11-12）不用再写了，已经隐含在上式（10）里了。
%%   公式（13-14）
%约束的含义是，持续运行需要受到上下爬坡约束，但启停不受爬坡约束。火电机组爬坡率为2%/min,小时调度体现不出来
for  i= 1:18for t = 2:24  C = [C,sum(P300_Diot(i,:,t))-sum(P300_Diot(i,:,t-1)) <= 0.3*300+(1-u300Git(i,t-1))*300,sum(P300_Diot(i,:,t-1))-sum(P300_Diot(i,:,t)) <= 0.3*300+(1-u300Git(i,t))*300,sum(P600_Diot(i,:,t))-sum(P600_Diot(i,:,t-1)) <= 0.3*600+(1-u600Git(i,t-1))*600,sum(P600_Diot(i,:,t-1))-sum(P600_Diot(i,:,t)) <= 0.3*600+(1-u600Git(i,t))*600,               P200_it(i,t)-P200_it(i,t-1) <= 0.3*200+(1-u200Git(i,t-1))*200,P200_it(i,t-1)-P200_it(i,t) <= 0.3*200+(1-u200Git(i,t))*200,];end       for i = 19:20C = [C, P200_it(i,t)-P200_it(i,t-1) <= 0.3*200+(1-u200Git(i,t-1))*200,P200_it(i,t-1)-P200_it(i,t) <= 0.3*200+(1-u200Git(i,t))*200,];end
end       
%%   公式（15-16）   火电机组的启停状态标识位置
C = [C,y200_it(:,2:24)-z200_it(:,2:24)==u200Git(:,2:24)-u200Git(:,1:23),  %公式(15)y300_it(:,2:24)-z300_it(:,2:24)==u300Git(:,2:24)-u300Git(:,1:23),  y600_it(:,2:24)-z600_it(:,2:24)==u600Git(:,2:24)-u600Git(:,1:23),  y200_it(:,1)==0,y300_it(:,1)==0,y600_it(:,1)==0,z200_it(:,1)==0,z300_it(:,1)==0,z600_it(:,1)==0,y200_it+z200_it<=1,%公式(16)y300_it+z300_it<=1,y600_it+z600_it<=1,];   
%%   公式（17-18）   最短运行时间，最短停机时间限值，取20h和10h
%需要注意的是，调度的维度只有24，因此，编写代码的时候，需要避免调用越界24
Ton = 20;%最短开机时间          
Toff = 10;%最短关机时间          
for t=1:(25-Ton)   %公式(17)for i=1:18C = [C,sum(u200Git(i,t:(t+Ton-1))) >= Ton*y200_it(i,t)];C = [C,sum(u300Git(i,t:(t+Ton-1))) >= Ton*y300_it(i,t)];C = [C,sum(u600Git(i,t:(t+Ton-1))) >= Ton*y600_it(i,t)];endfor i=19:20   C = [C,sum(u200Git(i,t:(t+Ton-1))) >=Ton*y200_it(i,t)];end
end
for t=1:(25-Toff)    %公式(18)for i=1:18   C = [C,sum(1-u200Git(i,t:(t+Toff-1))) >= Toff*z200_it(i,t),sum(1-u300Git(i,t:(t+Toff-1))) >= Toff*z300_it(i,t),sum(1-u600Git(i,t:(t+Toff-1))) >= Toff*z600_it(i,t)];endfor i=19:20    C = [C,sum(1-u200Git(i,t:(t+Toff-1))) >=Toff*z200_it(i,t)];end  
end   
%%  公式（19）  可再生能源出力约束
C = [C, 0<=P_PV,P_PV<=E_PV,0<=P_PVcur,P_PVcur<=E_PV,P_PV+P_PVcur==E_PV,0<=P_Wind,P_Wind<=E_Wind,0<=P_Windcur,P_Windcur<=E_Wind,P_Wind+P_Windcur==E_Wind,];    
%%  公式（20-24）  抽水蓄能电站功率约束
U_PuG = binvar(4,24,'full');
U_PuS = binvar(4,24,'full');
u_PuG = binvar(1,24,'full');
u_PuS = binvar(1,24,'full');
C = [C, U_PuG*300*0.3<=P_PuG,P_PuG<=U_PuG*300,   %公式（20）U_PuS*300*0.3<=P_PuS,P_PuS<=U_PuS*300,   %公式（21）u_PuG+u_PuS<=1,   %公式（22）];  
for t=1:24C = [C, U_PuG(:,t)<=u_PuG(1,t),   %公式（23）U_PuS(:,t)<=u_PuS(1,t),   %公式（24）];  
end    
​

  结果一览   

下载链接

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/681735.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

在Ubuntu22.04上部署FoooCUS2.1

在Ubuntu22.04上部署FoooCUS2.1

Fooocus 是一款基于 Gradio的图像生成软件&#xff0c;Fooocus 是对 Stable Diffusion 和 Midjourney 设计的重新思考&#xff1a; 1、从 Stable Diffusion 学习&#xff0c;该软件是离线的、开源的和免费的。 2、从 Midjourney 中学到&#xff0c;不需要手动调整&#xff0c;…
阅读更多...
优先级队列(堆)_PriorityQueue

优先级队列(堆)_PriorityQueue

前言 想要看如何使用可以通过目录跳转到 PriorityQueue的使用 优先级队列 概念 队列是一种先进先出(FIFO)的数据结构&#xff0c;但有些情况下&#xff0c;操作的数据可能带有优先级&#xff0c;一般出队 列时&#xff0c;可能需要优先级高的元素先出队列&#xff0c;该中场…
阅读更多...
前端秘法引言(配置vscode, 以及html的基础)

前端秘法引言(配置vscode, 以及html的基础)

目录 一.配置环境vscode 二.配置插件 三.vscode的实用小技巧 四.标题段落换行标签 五.格式化标签 一.配置环境vscode vscode官网https://code.visualstudio.com/ 点击右上角的download 根据不同的操作系统进行下载安装,我这里选的是Windows x64 安装好后打开,点击左上角的…
阅读更多...
【知识整理】产研中心岗位评定标准之测试岗位

【知识整理】产研中心岗位评定标准之测试岗位

为贯彻执行集团数字化转型的需要,该知识库将公示集团组织内各产研团队不同角色成员的职务“职级”岗位的评定标准; 一、定级定档目的 通过对公司现有岗位及相应岗位员工的工作能力、工作水平进行客观公正评定,确定各岗位的等级及同等级岗位员工对应的档级,从而为员工以后的晋升…
阅读更多...
Debezium发布历史125

Debezium发布历史125

原文地址&#xff1a; https://debezium.io/blog/2022/07/11/debezium-1-9-5-final-released/ 欢迎关注留言&#xff0c;我是收集整理小能手&#xff0c;工具翻译&#xff0c;仅供参考&#xff0c;笔芯笔芯. Debezium 1.9.5.Final Released July 11, 2022 by Chris Cranford …
阅读更多...
【开源】JAVA+Vue.js实现天然气工程运维系统

【开源】JAVA+Vue.js实现天然气工程运维系统

目录 一、摘要1.1 项目介绍1.2 项目录屏 二、功能模块2.1 系统角色分类2.2 核心功能2.2.1 流程 12.2.2 流程 22.3 各角色功能2.3.1 系统管理员功能2.3.2 用户服务部功能2.3.3 分公司&#xff08;施工单位&#xff09;功能2.3.3.1 技术员角色功能2.3.3.2 材料员角色功能 2.3.4 安…
阅读更多...
九、java 继承

九、java 继承

文章目录 java 继承3.1 根父类Object3.2 方法重写3.3 继承案例&#xff1a;图形类继承体系3.4 继承的细节3.4.1 构造方法3.4.2 重名与静态绑定3.4.3 重载和重写3.4.4 父子类型转换3.4.5 继承访问权限protected3.4.6 可见性重写3.4.7 防止继承final 3.5 继承是把双刃剑3.5.1 继承…
阅读更多...
爬虫-华为云空间备忘录导出到docx-selenium控制浏览器行为-python数据处理

爬虫-华为云空间备忘录导出到docx-selenium控制浏览器行为-python数据处理

背景适用情况介绍 老的荣耀手机属于华为云系统&#xff0c;家里人换了新荣耀手机属于荣耀云系统无法通过云空间将备忘录转移到新手机&#xff0c;不想让他们一个一个搞&#xff0c;于是整了一晚上想办法爬取下来。从网页抓取下来&#xff0c;然后存到docx文档中&#xff08;包…
阅读更多...
福布斯2023年推荐：十佳项目管理软件榜单揭晓

福布斯2023年推荐：十佳项目管理软件榜单揭晓

项目管理软件可以轻松规划项目、分配任务并保持团队井井有条&#xff0c;以便满足截止日期和目标。然而当今市场上有如此多的项目管理系统&#xff0c;选择适合您需求的正确选项可能很困难。为了提供帮助&#xff0c;福布斯小型企业顾问团队分析了数十家领先的提供商&#xff0…
阅读更多...
算法刷题day11

算法刷题day11

目录 引言一、FEB二、区间合并三、管道四、填充 引言 今天是初四&#xff0c;已经休息三天了&#xff0c;今天开始继续学习&#xff0c;然后觉得确实玩久了不太适应&#xff0c;已经有惰性了&#xff0c;不过还好自己喜欢&#xff0c;还是慢慢的坚持了下来&#xff0c;本来今天…
阅读更多...
.NET高级面试指南专题七【SocketWebSocket】

.NET高级面试指南专题七【SocketWebSocket】

Socket&#xff08;套接字&#xff09;是一种在计算机网络中实现通信的一种机制&#xff0c;它提供了一种标准的接口&#xff0c;使不同计算机上的程序能够通过网络进行数据交换。Socket允许在网络中的不同设备之间建立连接&#xff0c;进行双向的数据传输。 Socket通常用于实现…
阅读更多...
Ubuntu如何使用宝塔面板部署开源论坛HadSky并结合内网穿透远程访问？

Ubuntu如何使用宝塔面板部署开源论坛HadSky并结合内网穿透远程访问？

文章目录 前言1. 网站搭建1.1 网页下载和安装1.2 网页测试1.3 cpolar的安装和注册 2. 本地网页发布2.1 Cpolar临时数据隧道2.2 Cpolar稳定隧道&#xff08;云端设置&#xff09;2.3 Cpolar稳定隧道&#xff08;本地设置&#xff09;2.4 公网访问测试 总结 前言 经过多年的基础…
阅读更多...
设计模式-策略模式 Strategy

设计模式-策略模式 Strategy

策略模式 1) 原理和实现1、策略的定义2、策略的创建3、策略的使用 该模式最常见的应用场景是&#xff0c;利用它来避免冗长的 if-else 或 switch 分支判断。不过&#xff0c;它的作用还不止如此。它也可以像模板模式那样&#xff0c;提供框架的扩展点等等。 1) 原理和实现 策…
阅读更多...
Linux基础命令1

Linux基础命令1

Linux基础命令1 持续更新 ‍ 命令参数和作用示例ls列出文件和文件夹当前 ls​指定 ls /home​递归 ls -R​单列 ls -1​逗号分隔 ls -m​显示隐藏 ls -a​显示类型 ls -F​显示权限、所有者等详细信息 ls -l​按字母反序排列 ls -r​按扩展名排序 ls -lX​按日期和时间排序 …
阅读更多...
php数组与字符串函数

php数组与字符串函数

php数组与字符串函数 1. php数组2. 字符串函数 1. php数组 在php中&#xff0c;有三种类型的数组&#xff1a; 数值数组 - 带有数字ID键的数组关联数组 - 带有指定的键的数组&#xff0c;每个键关联一个值多维数组 - 包含一个或多个数组的数组 2. 字符串函数 在PHP中&#xf…
阅读更多...
耳机壳UV树脂制作私模定制耳塞适合什么样的人使用呢？

耳机壳UV树脂制作私模定制耳塞适合什么样的人使用呢？

耳机壳UV树脂制作私模定制耳塞适合以下人群使用&#xff1a; 对音质要求高的人&#xff1a;私模定制耳塞能够完美契合用户的耳朵形状&#xff0c;减少漏音和外部噪音的干扰&#xff0c;提供更好的音质体验。需要长时间佩戴耳机的人&#xff1a;私模定制耳塞能够提高佩戴舒适度…
阅读更多...
＜网络安全＞《30 常用安全标准》

＜网络安全＞《30 常用安全标准》

《常用安全标准》 1 个人信息安全 o《信息安全技术 个人信息安全规范》&#xff08;GB/T35273-2017) o《信息安全技术 个人信息去标识化指南》&#xff08;GB/T37964-2019) 2 工业控制安全 o《信息安全技术 工业控制系统安全检查指南》&#xff08;GB/T 37980-2019) o《信息…
阅读更多...
C++ //练习 6.41 下面的哪个调用是非法的？为什么？哪个调用虽然合法但显然与程序员的初衷不符？为什么？

C++ //练习 6.41 下面的哪个调用是非法的？为什么？哪个调用虽然合法但显然与程序员的初衷不符？为什么？

C Primer&#xff08;第5版&#xff09; 练习 6.41 练习 6.41 下面的哪个调用是非法的&#xff1f;为什么&#xff1f;哪个调用虽然合法但显然与程序员的初衷不符&#xff1f;为什么&#xff1f; char *init(int ht, int wd 80, char bckgrnd ); (a) init(); (b) init(24…
阅读更多...
Git基础命令，分支，标签的使用【快速入门Git】

Git基础命令，分支，标签的使用【快速入门Git】

Git基础命令&#xff0c;分支&#xff0c;标签的使用【快速入门Git】 Git基础常用命令Git工作流程工作区&#xff0c;暂存区和版本库文件状态获取Git仓库 git init | git clone查看文件状态 git status暂存已修改的文件 git add 查看已暂存和未暂存的修改 git diff提交文件更改…
阅读更多...
【JAVA】类加载机制

【JAVA】类加载机制

目录 一、快速梳理JAVA类加载机制 1、JDK8的类加载体系 2、沙箱保护机制 3、Linking链接过程 二、一个用类加载机制加薪的故事 三、通过类加载器引入外部Jar包 四、自定义类加载器实现Class代码混淆 五、自定义类加载器实现热加载 六、打破双亲委派&#xff0c;实现同类…
阅读更多...
最新文章

Copyright @ 2022~2023