在查找程序问题时,发现第一时步赋值没有问题,但是对后面时步进行约束时出现问题
问题函数应该在01决策变量约束上面。
问题1:01决策变量v(i,m)>=v(i,m-1)
可能性:m-1时步 m时步
0 0
0 1
1 1
问题2:对时间的约束,implies表述是否有问题
第一时步的时间都无法表示 T(1,1)=85
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