matlab+yalmip+cplex/gurobi风电场站内风机启动优化程序结果错误以及0-1变量为0.5的问题

最新推荐文章于 2023-07-14 14:44:59 发布

三年过后去爬杆

最新推荐文章于 2023-07-14 14:44:59 发布

阅读量1.4k

点赞数

文章标签： matlab

本文链接：https://blog.csdn.net/u013317992/article/details/126197633

版权

在风电场站内风机启动优化问题中，使用MATLAB结合YALMIP和CPLEX/GUROBI求解器进行0-1整数规划。目标是最大化恢复时间内的发电量，同时满足稳态频率和惯量响应约束。问题涉及27台风机，分为5组，每组在不同时步只能启动风机或辅机。然而，程序运行出现警告，指出连续松弛可能非凸，可能导致全局最优解无法找到。最终，优化问题被标记为不可行，得到的0-1变量结果为0.5。

摘要由CSDN通过智能技术生成

目标函数：给定恢复时间内发电量最大

约束条件：稳态频率约束+惯量响应约束

风电场内有27台风机，共分为5组。

每个操作时步有两种可能：辅机启动、风机启动；

每组辅机或风机一起启动，且某组辅机启动后下一时步风机启动。

第一时步预估后，设置启动功率。

warning('off');
clear
clc
yalmip;

%新能源场站内节点上的辅机及新能源机组信息
%第1列：节点号  第2列：辅机功率  第3列：额定容量  第4列：功频静特性系数（机组）  第5列：功频静特性系数（辅机）  第6列：惯性时间常数
Ebus=[
1,18000,2000000,25,1,10;
2,18000,2000000,25,1,10;
3,18000,2000000,25,1,10;
4,18000,2000000,25,1,10;
5,18000,2000000,25,1,10;
6,18000,2000000,25,1,10;
7,18000,2000000,25,1,10;
8,18000,2000000,25,1,10;
9,18000,2000000,25,1,10;
10,18000,2000000,25,1,10;
11,18000,2000000,25,1,10;
12,18000,2000000,25,1,10;
13,18000,2000000,25,1,10;
14,18000,2000000,25,1,10;
15,18000,2000000,25,1,10;
16,18000,2000000,25,1,10;
17,18000,2000000,25,1,10;
18,18000,2000000,25,1,10;
19,18000,2000000,25,1,10;
20,18000,2000000,25,1,10;
21,18000,2000000,25,1,10;
22,18000,2000000,25,1,10;
23,18000,2000000,25,1,10;
24,18000,2000000,25,1,10;
25,18000,2000000,25,1,10;
26,18000,2000000,25,1,10;
27,18000,2000000,25,1,10;
];
ND=length(Ebus);%有27个节点

%新能源场站内线路信息
%第1列：线路号  第2列：首节点  第3列：末节点  第4列：rij  第5列：xij
Ebranch=[1,0,1;
2,1,2;
3,2,3;
4,3,4;
5,4,5;
6,0,6;
7,6,7;
8,7,8;
9,8,9;
10,9,10;
11,0,11;
12,11,12;
13,12,13;
14,13,14;
15,14,15;
16,0,16;
17,16,17;
18,17,18;
19,18,19;
20,19,20;
21,20,21;
22,0,22;
23,22,23;
24,23,24;
25,24,25;
26,25,26;
27,26,27;
];
NC=length(find(0==Ebranch(:,2)));%有5组

%% 额定常数
NT=2*NC;%最大总时步为2倍分组数
fmax=0.2/50;%最大频率偏差；0.2/50P.U.
R=0.5/50;%频率变化量最大值  0.5HZ/s
SG=2500000;%外部带电域或储能的额定容量  2.5MW
KG=50;%外部带电域或储能的功频静特性系数  50p.u.
HG=5;%外部带电域或储能的惯性时间常数  5s
PL=Ebus(:,2);
SW=Ebus(:,3);
KW=Ebus(:,4);
KL=Ebus(:,5);
HW=Ebus(:,6);

%% 决策变量
K&