|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：6 J3 \4 A5 Z' u$ ^
: j' w8 | E* o& Ufunction f = adn(x) %主动配电网层目标函数1 o4 W$ C8 H2 i
global A b WT PV j1 PA Psub %全局量
) M( _+ ]! M. Z3 `1 U5 DADN=0;
6 u5 v: m& b( e7 G; Pfor j=1:24
1 o% a1 k6 R' g$ P8 C ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量) e5 E" d' d9 F
ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量
) a9 D( {- ^/ z0 V$ M ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量
# c9 X. @, r9 c& w. e n1 o ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量
' h# e k0 K, Y ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 / F) M5 x1 ~0 E+ C& T' q: M
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本& a# u+ |$ H6 m% d9 b- ]& Z* C
end. P9 t! r+ m: T! L# n
f=ADN;
- N3 u+ V4 [! d3 w( V# t! H- i6 Q
( j; E# v5 r4 s# F6 F+ Y6 v' A q& w6 j" }8 X; s
clc2 i2 I1 @/ H; m5 R4 y$ c
clear all
+ e& w) r/ j2 Cglobal A b WT PV PL j1 PA Psub/ \- h1 W7 a1 S9 [3 C1 {: i# F
A= [680; 675; 650; 645; 640; 635 ;710; 820 ;870; 920; 975; 960; 910; 860; 795 ;830; 925; 1240; 1210; 930 ;860 ;800 ;740; 685]; %输电网24小时趸售电价 元/MWh ' X' M5 G# ?- O* W5 c* C/ r
b= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh
5 K$ t' G% d/ N5 H2 c I c1 r# C0 Xj1=[500 475 450 443 437 469 563 669 700 744 781 763 713 650 600 625 863 1063 1000 738 663 600 563 475]; %聚合站24小时实时交易电价 元/MWh
' a0 [8 Y; _- _( MWT = [0.92 ;0.84; 1.04; 1.10 ;1.04 ;1.17; 1.17; 1.20; 1.12; 0.97 ;1.00; 0.90; 1.13 ;1.27 ;0.84 ;1.12; 1.04 ;1.16; 0.94 ;1.40; 1.40; 1.40; 1.36; 1.40];%风力发电24小时输出功率 MW / K; ?; ~1 W4 ?$ { d F
PV = [0; 0;0; 0; 0; 0.08; 0.27; 0.43; 0.70; 0.84; 0.96; 0.99; 0.99; 0.95; 0.85; 0.58; 0.37; 0.13; 0; 0; 0; 0; 0; 0]; %光伏发电24小时输出功率 MW
# z( ]4 r3 r( PWW=WT+PV; %风电、光伏发电出力之和 MW
& U2 x2 H* d/ oPL2=[0.066;0.066;0.067;0.071;0.078;0.088;0.09;0.09;0.075;0.09;0.084;0.085;0.084;0.083;0.085;0.09;0.096;0.097;0.101;0.102;0.093;0.086;0.078;0.072]; %节点2的24小时负荷数据（节点1为平衡节点） MW: [' F. I- q! Z% a
PL3=[0.057;0.056;0.057;0.061;0.068;0.078;0.081;0.08;0.065;0.08;0.074;0.075;0.074;0.073;0.075;0.081;0.086;0.087;0.091;0.092;0.083;0.076;0.068;0.062]; %节点3的24小时负荷数据 MW* H" n- S( }* W8 J9 h: s
PL4=[0.087;0.086;0.087;0.091;0.098;0.108;0.111;0.11;0.095;0.111;0.104;0.105;0.104;0.103;0.105;0.111;0.116;0.117;0.121;0.122;0.113;0.106;0.098;0.092]; %节点4的24小时负荷数据 MW! g' G! K4 v/ @
PL5=[0.027;0.026;0.027;0.031;0.038;0.048;0.051;0.05;0.035;0.051;0.044;0.045;0.044; 0.043;0.045;0.051;0.056;0.057;0.061;0.062;0.053;0.046;0.038;0.032];%节点5的24小时负荷数据 MW; Z$ n) Q! S6 T1 }$ n' k, X
PL6=[0.027;0.026;0.027;0.031;0.037;0.048;0.051;0.05;0.035;0.051;0.044;0.045;0.044; 0.043;0.045;0.051;0.056;0.057;0.061;0.062;0.053;0.046;0.038;0.032];%节点6的24小时负荷数据 MW- }3 l6 F9 w- [! J9 I) K2 k; A: |
PL7=[0.167;0.166;0.167;0.171;0.177;0.188;0.191;0.19;0.175;0.191;0.184;0.184;0.184;0.183;0.185;0.191;0.196;0.207;0.201;0.202;0.193;0.186;0.178;0.172]; %节点7的24小时负荷数据 MW
/ s1 g; H7 `' v4 [PL8=[0.167;0.166;0.167;0.171;0.177;0.188;0.191;0.189;0.175;0.191;0.184;0.184;0.184;0.183;0.185;0.191;0.196;0.207;0.201;0.202;0.193;0.186;0.178;0.172];%节点8的24小时负荷数据 MW$ y* `& ~: w' R$ H
PL9=[0.027;0.026;0.027;0.031;0.037;0.048;0.051;0.049;0.035;0.051;0.044;0.044;0.044;0.043;0.045;0.05;0.056;0.067;0.061;0.062;0.053;0.046;0.038;0.032]; %节点9的24小时负荷数据 MW
6 a3 p( q' ?& w- A6 S6 HPL10=[0.027;0.026;0.027;0.031;0.037;0.048;0.051;0.049;0.035;0.051;0.044;0.044;0.044;0.043;0.045;0.05;0.056;0.067;0.061;0.062;0.053;0.046;0.038;0.032];%节点10的24小时负荷数据 MW , G. V- v G: ^
PL11=[0.012;0.011;0.012;0.016;0.022;0.033;0.036;0.034;0.02;0.036;0.029;0.029;0.029;0.028;0.03;0.035;0.041;0.052;0.046;0.047;0.038;0.031;0.023;0.017]; %节点11的24小时负荷数据 MW* R% `3 _# K' p' S4 y6 h
PL12=[0.027;0.026;0.027;0.031;0.037;0.048;0.051;0.049;0.035;0.051;0.044;0.044;0.045;0.043;0.045;0.05;0.056;0.067;0.061;0.062;0.053;0.046;0.038;0.032];%节点12的24小时负荷数据 MW / k$ @5 Z' ?5 T
PL13=[0.027;0.026;0.027;0.031;0.037;0.048;0.051;0.049;0.035;0.051;0.044;0.045;0.045;0.043;0.045;0.05;0.056;0.067;0.061;0.062;0.053;0.046;0.038;0.032];%节点13的24小时负荷数据 MW
2 U: Y7 @! x8 D- RPL14=[0.087;0.086;0.087;0.091;0.097;0.108;0.111;0.109;0.095;0.111;0.104;0.105;0.105;0.103;0.105;0.11;0.116;0.127;0.121;0.122;0.114;0.106;0.099;0.092];%节点14的24小时负荷数据 MW
& x8 G" O0 v- Y4 @PL15=[0.027;0.026;0.027;0.03;0.037;0.049;0.051;0.049;0.035;0.051;0.043;0.045;0.045;0.043;0.045;0.05;0.056;0.067;0.061;0.062;0.054;0.046;0.039;0.032]; %节点15的24小时负荷数据 MW, D0 a. d" K% h. e E
PL16=[0.028;0.027;0.028;0.03;0.037;0.049;0.051;0.049;0.035;0.051;0.043;0.045;0.045;0.043;0.045;0.05;0.056;0.067;0.061;0.061;0.054;0.046;0.039;0.031]; %节点16的24小时负荷数据 MW/ [2 r ?1 j* H
PL17=[0.028;0.027;0.028;0.03;0.037;0.049;0.051;0.049;0.036;0.051;0.043;0.045;0.044;0.044;0.045;0.05;0.056;0.067;0.061;0.061;0.054;0.046;0.039;0.031]; %节点17的24小时负荷数据 MW& G* y# W. m1 Q
PL18=[0.058;0.057;0.058;0.06;0.067;0.079;0.081;0.079;0.066;0.071;0.073;0.075;0.074;0.074;0.075;0.08;0.086;0.097;0.091;0.091;0.084;0.076;0.069;0.061]; %节点18的24小时负荷数据 MW5 f/ B/ f7 I, q' h3 `8 N
PL19=[0.057;0.056;0.057;0.059;0.066;0.078;0.08;0.078;0.065;0.07;0.072;0.074;0.073;0.072;0.074;0.079;0.085;0.095;0.09;0.09;0.083;0.075;0.068;0.06]; %节点19的24小时负荷数据 MW& q: o! g/ o+ x$ v" Z
PL20=[0.057;0.056;0.057;0.059;0.067;0.078;0.08;0.088;0.065;0.07;0.072;0.074;0.073;0.072;0.074;0.08;0.085;0.09;0.09;0.09;0.083;0.075;0.068;0.06]; %节点20的24小时负荷数据 MW0 S1 K `; U; R, F) `
PL21=[0.056;0.056;0.057;0.059;0.067;0.078;0.08;0.088;0.065;0.07;0.072;0.074;0.073;0.072;0.074;0.08;0.085;0.085;0.09;0.09;0.083;0.075;0.068;0.06]; %节点21的24小时负荷数据 MW, m9 c6 W6 ]2 k5 w
PL22=[0.056;0.056;0.057;0.059;0.067;0.078;0.08;0.088;0.065;0.07;0.072;0.074;0.073;0.072;0.074;0.08;0.085;0.085;0.09;0.09;0.083;0.075;0.068;0.06]; %节点22的24小时负荷数据 MW * c5 U) U+ M2 @1 P7 ^
PL23=[0.056;0.056;0.058;0.059;0.067;0.078;0.08;0.088;0.065;0.07;0.082;0.074;0.073;0.072;0.074;0.08;0.086;0.085;0.09;0.09;0.083;0.075;0.068;0.06]; %节点23的24小时负荷数据 MW
# O& [( w0 a. wPL24=[0.386;0.386;0.388;0.389;0.397;0.408;0.41;0.41;0.395;0.4;0.402;0.404;0.403;0.402;0.404;0.41;0.416;0.415;0.42;0.42;0.413;0.405;0.398;0.39]; %节点24的24小时负荷数据 MW
% u$ X$ }* E( [ aPL25=[0.386;0.386;0.388;0.389;0.397;0.408;0.41;0.41;0.395;0.4;0.402;0.404;0.403;0.402;0.404;0.41;0.416;0.416;0.42;0.42;0.413;0.405;0.398;0.39]; %节点25的24小时负荷数据 MW) F: ]; W* D+ h6 i
PL26=[0.026;0.026;0.028;0.029;0.037;0.048;0.05;0.05;0.035;0.04;0.042;0.044;0.043;0.042;0.044;0.05;0.066;0.056;0.06;0.05;0.053;0.045;0.038;0.03]; %节点26的24小时负荷数据 MW( P2 n0 ]# q# S2 ]9 }
PL27=[0.026;0.026;0.028;0.029;0.037;0.048;0.05;0.05;0.035;0.04;0.042;0.044;0.043;0.042;0.044;0.05;0.066;0.056;0.06;0.05;0.053;0.045;0.038;0.03]; %节点27的24小时负荷数据 MW5 B0 |) k. ~4 T- P
PL28=[0.026;0.026;0.028;0.029;0.037;0.048;0.05;0.048;0.035;0.04;0.042;0.044;0.043;0.042;0.044;0.05;0.066;0.056;0.06;0.05;0.053;0.045;0.038;0.029]; %节点28的24小时负荷数据 MW
4 g# M2 v b% D$ O1 w7 |- bPL29=[0.086;0.086;0.088;0.089;0.097;0.108;0.11;0.108;0.095;0.1;0.102;0.104;0.103;0.093;0.104;0.11;0.126;0.116;0.119;0.11;0.113;0.106;0.098;0.089]; %节点29的24小时负荷数据 MW
2 s, \ [/ s0 C: ?* D1 @PL30=[0.166;0.166;0.168;0.169;0.177;0.188;0.21;0.188;0.175;0.18;0.182;0.184;0.173;0.173;0.184;0.19;0.206;0.196;0.199;0.19;0.193;0.186;0.178;0.169]; %节点30的24小时负荷数据 MW
% g0 t# R( {% K5 o4 hPL31=[0.116;0.116;0.118;0.119;0.127;0.138;0.15;0.138;0.125;0.13;0.142;0.134;0.123;0.123;0.134;0.14;0.156;0.146;0.149;0.14;0.143;0.136;0.128;0.119]; %节点31的24小时负荷数据 MW) w( I3 P$ J: d* ]3 F* e' K Z
PL32=[0.176;0.176;0.178;0.179;0.187;0.198;0.21;0.198;0.185;0.19;0.192;0.194;0.193;0.193;0.194;0.2;0.216;0.206;0.209;0.2;0.203;0.196;0.188;0.179]; %节点32的24小时负荷数据 MW! U& E. K) d4 x% ?: k6 w3 ^
PL33=[0.026;0.026;0.028;0.029;0.037;0.048;0.06;0.048;0.035;0.04;0.042;0.044;0.043;0.043;0.044;0.05;0.056;0.059;0.059;0.05;0.053;0.046;0.038;0.03]; %节点33的24小时负荷数据 MW
+ b- x* D" {$ e2 VPL=[PL2 L3L4L5;PL6;PL7;PL8;PL9;PL10;PL11;PL12;PL13;PL14;PL15;PL16;PL17;PL18;PL19;PL20;PL21;PL22;PL23;PL24;PL25;PL26;PL27;PL28;PL29;PL30;PL31;PL32;PL33]; %IEEE33节点系统24小时各节点负荷数据 MW
/ W3 T& i$ f' U: T$ f# zPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)
" q8 J* e9 I' Y" _+ l, D' LPDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW& W; L$ i0 B: @
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW) Y7 ?2 Z8 }( n! S
Psub=zeros(24,1); %初始化Psub矩阵+ P( G( Z* P" `5 i. h9 @
PA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
' L; r' m& d0 k0 V8 s# ?SA=zeros(24,1); %初始化聚合站电池容量 MWh
5 m+ J+ w' a. F! O1 _$ {x0=zeros(864,1); %变量初始化
6 U' D o$ i2 Y& L%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
, d, W0 j# \. H/ T) ]+ b& o/ M% B6 }Aeq=zeros(793,864); %等式约束系数矩阵初始化
" V5 I* B9 Q$ h0 abeq=zeros(793,1); %等式约束结果矩阵初始化6 n2 d9 }5 c! E
for i1=1:793
* M7 }. w2 J! V5 }%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%+ w% r' F* x7 p
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL. Z1 q% O7 y1 x8 X
Aeq(i1,i1)=1;" J( ~/ k; o9 _) q' n+ ^' U
beq(i1,1)=PL(i1,1);
4 P) N2 J7 [. O0 l) n end4 ?$ t6 L/ n: E. s, Y; N
if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV
" y: l1 F' k# X: k; I$ G# P Aeq(i1,i1)=1;) [. S, p6 P. c d4 m! V& l6 j
beq(i1,1)=PL(i1,1)-PV(i1-144,1);3 Q x% s2 ^% x! q2 q
end
6 W; [/ \+ P `0 _9 I if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL3 p. G4 B" X6 d, ?& d
Aeq(i1,i1)=1;. Q. h$ s" _5 U* |, y
beq(i1,1)=PL(i1,1);- j$ P! A1 G3 n' H5 L) t9 D
end6 c( a) ?" h1 C! w: T
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL
4 R6 {$ W) o: Y- P Aeq(i1,i1)=1;. O4 D; K0 C/ r" B, ^: f& r) O- M
Aeq(i1,i1+504)=1;! X0 t Q1 o; ?! }/ `- _
beq(i1,1)=PL(i1,1);) q" m' T) O4 j* R& h0 q$ p3 N) Q8 k3 K
end, q3 @3 E; O. M
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL. o6 @( E" T+ E. M' z w
Aeq(i1,i1)=1;8 A* u) }' n) O7 r; {) W7 F
beq(i1,1)=PL(i1,1);
& ~: p8 y7 y: ~, Z end2 b( ~+ r! D8 a
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL
8 r! a7 u* Y6 _ Aeq(i1,i1)=1;
2 u7 ^0 @! m9 l( V! U) p Aeq(i1,i1+384)=1;
4 u/ ^ K# ]# n4 T+ e8 }/ H q beq(i1,1)=PL(i1,1);
' V/ k+ t! S) Z7 q end% Y' V" L( w! A0 ~) F! |- i
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL
% e# s l5 r8 C Aeq(i1,i1)=1;/ M' Q$ ^' Q7 T3 K3 M
beq(i1,1)=PL(i1,1);% G7 H2 H) [! ], O9 R7 ^
end
" [7 j6 i/ Q. r if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT
+ E( f8 H! v' [0 q5 @2 I* E b Aeq(i1,i1)=1;( R8 H2 ~) P% ?
beq(i1,1)=PL(i1,1)-WT(i1-456,1);
/ m5 `& D: B3 w1 @ end% }( |; {3 q' l4 A& L: E: w
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL
: m4 I1 i5 U# r5 }) s$ Z! e/ h; f Aeq(i1,i1)=1;
% T+ A* ]0 v; j3 Z4 J, ~" k j beq(i1,1)=PL(i1,1);
% k$ E* P* k N/ Z; K# B: ^ end3 l4 n$ z# V0 I+ X
if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL
, V6 a, Z( K/ x Aeq(i1,i1)=1;
6 O F) ^% W8 D) h7 }" q Aeq(i1,i1+144)=1;
$ S: \! o/ H( V9 n8 ? beq(i1,1)=PL(i1,1);2 D1 ?2 T9 L' F* @
end) h2 o, @) p& n- S5 m4 @4 E
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL
n6 R" m& ^9 N ~0 i# D% Y Aeq(i1,i1)=1;
$ Q+ a& E5 u' r$ \/ a beq(i1,1)=PL(i1,1);
4 r8 m6 v( r/ V end! Y [; W5 z' F+ h1 H0 b: N
%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
2 A/ I8 ]* ~8 d. W if i1>=792&&i1<=7922 \* i1 q5 S! z& v
Aeq(i1,i1+72)=1; %对应电池容量SA# H; N! p6 y% X# }$ ?
for i2=1:24) T; \; I# u8 [; K0 A/ e2 L
if i2<=i1-7688 [; k6 k V+ J; [( Z+ A# d$ B0 c
Aeq(i1,i2+816)=1;2 I/ y, w8 T! G, K0 Q
end
7 d* @' c6 }6 c$ @2 G end
P: e, U' f5 F9 ~5 {* N beq(i1,1)=5.4; %5.4MWh) I' ^( Y. _! [# _8 Z: `/ g% u
end& ^8 p4 o; M1 R4 C; G/ X
%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%. W( G o7 w. O6 p; e
if i1>=793&&i1<=793 %对应S24=5.4MWh. c! m/ B" Q' v8 p5 y" ~9 F" B) g
Aeq(i1,i1+71)=1;
4 K- p1 D& N3 q) ]/ X3 k beq(i1,1)=5.4;
) [" W/ @; P9 B1 G9 V end
2 o6 }! S; ]2 \) p/ Mend
- D1 }; D0 m C; @%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
& Z) U9 ]% w: Z( [: NA=zeros(48,864); %初始化线性不等式约束系数矩阵
/ u2 M6 X8 {/ @: p" j2 x0 a2 hb=zeros(48,1); %初始化线性不等式约束结果矩阵
3 g9 G5 `3 r" R F! w0 Afor i3=1:48
) ` m2 S2 Y( n2 v if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t, c* ~8 X$ o/ |0 F. }# L) S
A(i3,i3+816)=1; %对应(PA)t的系数
! m4 ^$ m/ k9 H1 v( u if i3==1$ P2 N* A7 J ]8 ^. q) |3 I
A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数: j J4 A" j& H6 q/ P
b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵
$ ]9 S; C9 c8 z4 y, W/ o2 q6 F: k else+ a$ _& m) u0 b/ E: t$ m
A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数6 R! }; L& T) ^: ?3 @' c9 O& E
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵
' f9 c3 g% E+ f, G end3 S; L: f/ _" y4 ~' C4 F
end " r% T' n, @0 I6 T/ G& e
if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
% O+ i1 G1 ~& B A(i3,i3+792)=-1; %对应(PA)t的系数. O @, K/ B, H0 r- r+ W, i
if i3==25& t! ~% ^8 q9 P( g2 Z1 K
A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数, G5 ?$ i/ O7 R- O _: q
b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵; D, F# ?* ?+ t. C. v/ o
else
5 V7 W- m% f. `0 t0 }: x0 s A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
9 U: ~5 q$ l. G b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵, t3 C1 @3 @% J1 F+ O
end
& A; N2 C( y4 l4 _3 k) Z9 s end% F( F' n, k' C/ q1 x t/ r
end
9 z W7 k: i' ^$ g%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%3 Q; |6 f; Q7 j5 {- {8 D
Ib=zeros(864,1); %初始化下限矩阵
* H: B# V6 g( A$ d% _/ t2 aUb=zeros(864,1); %初始化上限矩阵7 A# S* k0 M0 u9 f0 ?( N
for i4=1:864 %Ib Ub赋值8 q4 ^" V7 U# l7 C
if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
/ v/ J0 a( n8 N6 `8 F' s, T Ib(i4,1)=0;& j' o) x+ {7 u. G! }; }
Ub(i4,1)=3.73;4 m4 b5 q3 e8 y# _- i$ l& P
end) }& U9 z. P& t6 t& V( M
if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW
4 }! w" r# | }+ F/ p1 ]2 S7 j Ib(i4,1)=0;
, c, R @2 H' ~6 @/ Y# I2 k Ub(i4,1)=3;
3 x/ |' @3 n* B: H" I- R+ ] end9 f. k7 F8 o' [) @0 a
if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW
' ^; M" e2 }$ ~0 | Ib(i4,1)=0;
4 o* x5 I' r- j8 T3 u3 W Ub(i4,1)=3;
% d: S: S+ H/ {1 n end, ^( B6 r% z+ {9 ]' ^+ [" t
if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW4 k- I; Y: f5 Y5 Z$ P1 s [
Ib(i4,1)=-2;4 o8 M. E6 U5 c# b4 |+ |% ?' o
Ub(i4,1)=2;& a: ^! G8 d& i5 }
end* Y5 M) {2 x3 J4 W7 j8 [4 S
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh. c. o0 U; N2 J% f& [! w
Ib(i4,1)=1.08;
2 s" q. [" Y, k' E. [ Ub(i4,1)=5.4;3 i6 o( u* |3 _' x3 g0 r8 U$ G* ?
end
( O/ E0 r7 y& g% O6 }( Oend
( Z6 x$ ?% w/ Z) h# w4 B%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%- P$ H1 L" P+ L n& J, n# |# G
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
8 C- e# Q5 \6 Z2 lfor i5=1:768
4 U: C9 c4 O0 M( j0 ^! y% Q Pij(i5,1)=x(i5,1);+ L7 m p! U6 h( \
end' L/ O3 s2 p0 @% r: @
for i6=1:24
' I2 }, u! A. c& ]9 I6 E( L PDG1(i6,1)=x(i6+768,1);
# ~ b9 X! r9 I: C" d9 l. C1 y PDG2(i6,1)=x(i6+792,1);
7 `+ \: M/ P- q5 j PA(i6,1)=x(i6+816,1);
: B' ` S+ D" ]+ [- A9 P; \ SA(i6,1)=x(i5+840,1);) o8 x" J- h# Y8 V5 ?
end1 i; W m% z3 r. x; y
for i7=1:24
7 r, s3 c6 c$ S; q6 q8 q Psub(i7,1)=x(i7,1)+x(i7+24,1)+x(i7+48,1)+x(i7+72,1)+x(i7+96,1)+x(i7+120,1)+x(i7+144,1)+x(i7+168,1)+x(i7+192,1)+x(i7+216,1)+x(i7+240,1)+ x(i7+264,1)+x(i7+288,1)+x(i7+312,1)+x(i7+336,1)+x(i7+360,1)+x(i7+384,1)+x(i7+408,1)+x(i7+432,1)+x(i7+456,1)+x(i7+480,1)+x(i7+504,1)+x(i7+528,1)+x(i7+552,1)+x(i7+576,1)+x(i7+600,1)+x(i7+624,1)+x(i7+648,1)+x(i7+672,1)+x(i7+696,1)+x(i7+720,1)+x(i7+744,1);7 G4 w% b. a8 }+ V* M" H
end5 B" w3 q% L0 l( h
disp(['运行成本：',num2str(fval1),'元']);
3 Z3 z' \9 H% x3 ~subplot(2,1,1);
9 c) [8 a. {: B+ Z. O Tplot(PA,'b-','LineWidth',1.5);
8 t, X$ B& k7 C# [hold on
3 L& e: g" h$ r' g) N, {7 P9 {: b, Igrid on
' c2 N3 U6 g+ `2 N2 L5 w/ Y4 h) f E2 i
2 W1 d) A# y1 n) M0 x8 |2 b' E! J
|
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！