电力系统机组组合（Matlab代码实现）

 💥💥💥💞💞💞欢迎来到本博客❤️❤️❤️💥💥💥

🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者，博主专门做了一个专栏目录，整个专栏只放了一篇文章，足见我对其重视程度：博主专栏目录。做到极度细致，方便大家进行学习！亲民！！！还有我开了一个专栏给女朋友的，很浪漫的喔，代码学累的时候去瞧一瞧，看一看：女朋友的浪漫邂逅。有问题可以私密博主，博主看到会在第一时间回复。
📝目前更新：🌟🌟🌟电力系统相关知识，期刊论文，算法，机器学习和人工智能学习。
🚀支持：🎁🎁🎁如果觉得博主的文章还不错或者您用得到的话，可以关注一下博主，如果三连收藏支持就更好啦！这就是给予我最大的支持！

👨‍🎓博主课外兴趣：中西方哲学，送予读者：

👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。在我这个专栏记录我有空时的一些哲学思考和科研笔记：科研和哲思。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“真理”上的尘埃吧。

     或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎

📋📋📋本文目录如下：⛳️⛳️⛳️

目录

1 概述

2 数学模型

2.1 目标函数

2.2 约束条件 

3 Matlab代码实现

4 运行结果

---

1 概述

机组组合问题是指在满足系统安全约束和电能质量要求的前提下，利用合理的机制控制机组的开/关状态，调整各时段的出力，使某一调度周期（通常指一天）内的总运行成本最小化。机组组合是解决电力系统最优潮流和经济调度的前提。

机组组合问题是NP-hard型问题，具有高维、非凸、离散、非线性等特点，到目前为止，还没有好的解决办法。目前，机组组合的求解算法主要是两类。一类是传统算法：包括优化级表法、混合整数规划、动态规划法、拉格朗日算法等。

2 数学模型

2.1 目标函数

2.2 约束条件 

本文目标函数和约束比较简单，才考虑三台机组，主要是方便入门学习。

3 Matlab代码实现

%% 机组组合
clc
lower_lim =[ 150 100 50 ];  %机组下限(MW)
upper_lim=[600  400  200];  %机组上限
%% ====三台机组机组参数a,b,c======
a=[ 510.0  310.0 78.0];  
b=[7.20 7.85 7.97];  
c=[0.00142  0.00194 0.00482]; 

fc=[1.1  1.0   1.2];%三台机组燃料成本

load=550;  %总负荷

%% ====机组出力====
c1= [ 0     400 150
      550   0   0
      500   0   50
      295   255 0
      267   233 50];
  
f1=zeros(5,1);

for i1=1:length(c1(:,1))
    for j1=1:length(a)
        if(c1(i1,j1))==0  %如果机组出力为0
            f1(i1,j1)=0;
        else
            f1(i1,j1)=(a(j1)+b(j1)*c1(i1,j1)+c(j1)*c1(i1,j1)^2)*fc(j1);  %机组运行成本
        end
    end
end

%% 总成本
t_cost=sum(f1')
[a1,b1]=min(t_cost);
p_comb=c1((b1),:)  %机组出力

4 运行结果