|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：
7 d; M0 K( ^1 [, D W# I3 M7 R' \+ U4 Q7 i1 T3 \
function f = adn(x) %主动配电网层目标函数
6 \! `; p! @& Z* u, E9 S/ Nglobal A b WT PV j1 PA Psub %全局量
6 P- a) D: H' @4 l- k2 L& hADN=0;
6 D5 Q+ A2 r) Ufor j=1:24
6 E3 \+ @2 I+ Y1 z( i ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量/ A! x0 X* N& a1 q+ [7 \* e
ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量% v/ Z5 O" e, N
ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量
+ E! o/ K5 p6 J! m+ U2 _ ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量
# H; n# v% F3 E8 i ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 $ g: x5 Z' g9 A+ m; \ c
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
) L0 F; M: M/ l/ u( X+ T8 ?end
+ ]! J7 |* \$ l# k$ lf=ADN;- C" p' k# D" x' E% I( n6 g
* w; H. y9 @- y; j" U& t r8 H9 S1 G
9 T# H/ Z0 r' z, o: V8 [4 j; u! u) \, h+ z
clc
: o* I$ V |$ Oclear all6 l$ _+ n i4 C$ c# h- o7 j
global A b WT PV PL j1 PA Psub/ M1 Q# l0 }9 h! c& N, v
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
# Q6 w$ |: K6 p) C6 Z# z/ l6 \/ }" Db= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh
w4 I2 n6 ?- t$ 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
, I7 e `' {) A5 r7 W9 x4 AWT = [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
" K5 i4 J ]8 p/ N) {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, h" ?. c- s/ C q1 F- \ M
WW=WT+PV; %风电、光伏发电出力之和 MW8 n! u6 l/ b' J; o' O" W
PL2=[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
: e1 ]/ Q8 H1 h) \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
# A/ v: d, s/ E/ I8 l4 ?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
% m; e- m3 Y+ B6 ]2 I4 K/ MPL5=[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
" w& V9 n+ i4 I7 A- CPL6=[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! l d% K% j' Q
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% D8 d' A4 i* f. j+ n- X
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: a& Y# }# q' d) t/ D. ?) h# Y
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
( x$ V4 L$ Q# E, \) S( 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
& h$ {& I# s! q. TPL11=[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$ y. a( j) j) P) W) D h% R6 z
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
' O8 G+ p! b& W9 B$ RPL13=[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
7 l# u, F' j. M/ H! `1 P+ XPL14=[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
% S' L7 D% k8 {3 D& A. q: {* k3 h- TPL15=[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小时负荷数据 MW3 Z* t9 y1 a- W% u4 l4 a
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+ h3 v- v4 Q7 f- v5 I
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
% U0 Y; N- ]: O6 l1 j" }' E7 ZPL18=[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小时负荷数据 MW6 C; I' d- q5 Y
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
: G$ r8 c7 w9 N) K* @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小时负荷数据 MW$ r8 `: w) K9 h, q5 w( P; E
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
4 w& l* B' A& b* zPL22=[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
7 O; I* b; W( ~) gPL23=[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) U' A1 Z" c5 Z* ^& |
PL24=[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
/ U0 T/ v, O1 o1 s# i8 BPL25=[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* Y( w# `( ]* H( j+ P4 C- v
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
' j; d1 q" u+ g8 {! HPL27=[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小时负荷数据 MW2 M1 x5 S% F; Y, m8 n
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小时负荷数据 MW8 i5 X }% T$ c- h
PL29=[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
, K6 c0 l& Z6 l/ v UPL30=[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小时负荷数据 MW7 w; x3 l' J* t/ T
PL31=[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
3 X P9 I% u4 F ]8 d0 ePL32=[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# i9 ?# B! ^, r, U" r2 K
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, O& _7 C+ k' { z9 p/ |
PL=[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
$ W- _ x: {4 _8 f- {1 MPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)( Z n6 X, l3 I' w9 h) l/ j
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW e: P! y) f" O* a- _
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW
/ g; u9 N# {5 c; vPsub=zeros(24,1); %初始化Psub矩阵
# N" Z; q8 Z! k+ \$ N; ~. cPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW$ [( y3 N$ @) @: p& O5 L6 P& F
SA=zeros(24,1); %初始化聚合站电池容量 MWh% Y0 z3 b) J4 T+ a6 I
x0=zeros(864,1); %变量初始化 , h% O/ O( B& w h
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
6 ~4 I# ]: I' s* lAeq=zeros(793,864); %等式约束系数矩阵初始化/ O) B( e% P$ S: W/ `
beq=zeros(793,1); %等式约束结果矩阵初始化, V) p, T E0 }: Z$ b$ s7 u
for i1=1:793
; i6 b% E" B& w. T5 L4 ?3 ?%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%7 {' F; L, _7 s4 F% d5 T1 V
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL2 }9 Q" L- a. G% G: H: l
Aeq(i1,i1)=1;
1 q6 i" B7 F2 g W beq(i1,1)=PL(i1,1);
4 C( }7 Z( b2 g! j4 K- \4 o end
6 c) `1 W0 F1 S* O4 B# E if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV+ s* V: F6 U2 S& M
Aeq(i1,i1)=1;
) j' W& O% b! B6 O8 S3 b beq(i1,1)=PL(i1,1)-PV(i1-144,1);
. a, v0 ?8 J6 O end2 @/ m$ m# H& c" g4 j1 L
if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL! r, `1 V, m- n0 q) |8 q: n& B
Aeq(i1,i1)=1;! F# Q7 J/ ~$ |4 G o" z
beq(i1,1)=PL(i1,1);
+ ?# |: S) i5 A4 E& b+ \' ] end
) e9 Z. F3 F5 Y# e9 e; V if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL; ^2 s& a" D) c% N
Aeq(i1,i1)=1;5 E% P7 h. g; _ f
Aeq(i1,i1+504)=1;; D2 F- D! m4 l- w# o# w
beq(i1,1)=PL(i1,1);
9 y% e$ ?6 n; z9 ` end) M; f# V4 N8 \
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL' w( W, H6 s& o0 a9 @
Aeq(i1,i1)=1;
; `# t6 x0 o% l- h beq(i1,1)=PL(i1,1);7 n( ]! g8 [! x
end6 W! I. G4 F% ^ ?+ ^' s2 {
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL2 R9 ]# Z7 y' P0 ?: ]
Aeq(i1,i1)=1;1 t% U; `3 m# s8 @( M+ q
Aeq(i1,i1+384)=1;
0 }# K9 o, D2 K5 y# m2 a3 x beq(i1,1)=PL(i1,1);$ d4 }; p: X# a( m$ @
end$ w1 M: Q2 `( g* ^; N8 h
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL
- m, W" b# M% V Aeq(i1,i1)=1;" k* b+ Y4 {5 r: W6 y
beq(i1,1)=PL(i1,1);
2 f) d) U7 Y# I2 h6 V. x! \ end
/ @$ ~+ j" M9 I& j7 ~( D8 V if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT. X+ { M+ o) y* P$ [1 \
Aeq(i1,i1)=1;6 l. l# a* N9 n b9 s1 Q( r/ ^, |2 }
beq(i1,1)=PL(i1,1)-WT(i1-456,1);
4 J, w5 I) Z" h' b1 X4 V. x end: G4 z9 L3 U" }4 R% d4 M
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL; v# D. x, N# \3 a
Aeq(i1,i1)=1;. b, n3 f5 j* b! B* Y# _
beq(i1,1)=PL(i1,1);
8 t9 C& P/ P9 r. j; ~ end
h5 h" u1 x" c6 { if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL V \$ T( E2 g- }
Aeq(i1,i1)=1;- U' c% ]3 H- g
Aeq(i1,i1+144)=1;
- Z8 Q) l( F: P" I4 { beq(i1,1)=PL(i1,1);6 K* C" n! n7 s- m1 q, t. L' c
end, s8 w* o7 `, ?" q
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL/ J6 b6 E5 R; ]: A
Aeq(i1,i1)=1;
' o' T$ i3 X" Y. O& I# o2 g beq(i1,1)=PL(i1,1);
1 y$ e' o9 N! X x1 X end
# I& y' d3 F2 d. _& N0 \%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%, @) {% N( `3 _4 z# u8 _1 T+ k5 Y
if i1>=792&&i1<=792- `" b$ b- \# L* a, |+ M
Aeq(i1,i1+72)=1; %对应电池容量SA
5 a( t5 u8 d7 S, W for i2=1:24
W& L5 H& @) h' e7 m if i2<=i1-768
" k4 A4 Q- \1 D; G E( E) e Aeq(i1,i2+816)=1;
. J( U1 B( D" f% D5 Q- A end( j/ G) @5 y& R
end6 s! F; H1 u- R
beq(i1,1)=5.4; %5.4MWh. @$ w7 T& {% D$ h: u: [3 K
end6 n, q" L9 Z) D/ q- G
%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
! K/ l0 E6 k" n if i1>=793&&i1<=793 %对应S24=5.4MWh2 Q# c% e6 F" P
Aeq(i1,i1+71)=1;0 {' `: }: t) U- F
beq(i1,1)=5.4;
+ A& K2 h- y; ^/ e( d end1 D0 k& [, m, f) A0 j/ z& \
end
5 b4 F! J8 K, g5 y- p9 G%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%4 Y1 P9 J& X9 o: E X. Z
A=zeros(48,864); %初始化线性不等式约束系数矩阵
7 t" G, A: H& s, f. t+ pb=zeros(48,1); %初始化线性不等式约束结果矩阵
( ^1 o$ h- L7 c6 H' zfor i3=1:48 5 n' D* m/ b; L' x) K t* h5 R
if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t8 s7 x" @: G6 |
A(i3,i3+816)=1; %对应(PA)t的系数
3 [9 V# D: }# n" q if i3==1
3 u' |2 S$ W4 z$ @( Z. B7 J A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
: U/ r' ]( z& k b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵" m0 p" z) t" J2 ]
else
3 p4 b" G$ e) |7 N A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数9 m) }2 s) O6 \
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵6 r5 l; ]% o2 v& }. s3 J
end% c1 t7 a4 W3 X- D! J& l
end
$ L/ C3 o: F3 Z' \2 t* o5 W/ _ if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1) J. e8 V" s H6 \5 s
A(i3,i3+792)=-1; %对应(PA)t的系数7 T3 g6 c, s: g4 Q7 v% V
if i3==25
' b# [. q# A" |: i5 z6 o! ~4 x A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
, p8 E# o* t4 N7 h b5 A b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
- k9 x" A) i( H. |. _; S else8 X3 f' R9 [/ G% F3 E% e) N4 J0 ~- y
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
* p! U5 X0 j7 C% G% { c0 ]: \9 x" L% X9 Z b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
" Z' l& \- c8 ?/ H4 _ end
( |: b7 t$ m% d! ~) B end
/ g) Y# R, `/ o; Cend3 ]- c7 w& ?8 z* k
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
/ h1 ^5 k8 }% h# m1 z* oIb=zeros(864,1); %初始化下限矩阵5 o) N1 f7 @& ~+ n
Ub=zeros(864,1); %初始化上限矩阵, t/ {4 V9 F& k- o0 p" [- B
for i4=1:864 %Ib Ub赋值+ y. w2 }, o2 o3 I
if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求$ Y8 w d2 e# i: E, U
Ib(i4,1)=0;
) g6 O1 h; U' d7 K2 k Ub(i4,1)=3.73;
) m! C, S/ ? c, L' m3 ]$ B end
/ |8 U, v$ R3 p1 L3 W6 _% o! W if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW% F: y+ k! ^. y$ V7 g
Ib(i4,1)=0;2 e9 Y8 \$ `1 |( W
Ub(i4,1)=3;' L3 `3 P8 }# N! k% K( k6 W1 |
end8 @9 }# }0 z, }1 J; w- |
if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW- g) s: z J8 c
Ib(i4,1)=0;0 M8 k, O9 b7 K) x- ^" H# R2 v
Ub(i4,1)=3;' a r) `( d" z3 |
end$ K( \4 D. d9 b
if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW" ^% C! k) x& w) ?
Ib(i4,1)=-2;
; \7 c: y1 K3 O9 I- c6 c9 O Ub(i4,1)=2;" _; t- l- ]; B6 U* `( y
end
2 g5 k: \: N. d if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh
2 D+ M" b4 C7 `, q/ R Ib(i4,1)=1.08;
8 |0 ]; O7 }9 h0 W6 | Ub(i4,1)=5.4;4 l4 W$ S+ X# x
end! A# v0 H4 z, K
end
* j N4 I) z3 P. n- k%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%. h- u7 J, o, |, R/ c
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);0 i/ U8 b. \) l
for i5=1:768
) A+ H( d" Q& E# j& B9 b: k Pij(i5,1)=x(i5,1);! n: z! y5 N( |: w* z1 I
end
9 ^8 b- n' g- F) \9 A4 lfor i6=1:24
4 E* d+ r! A/ [; t$ v PDG1(i6,1)=x(i6+768,1);
& w$ { M, R! t& e; u PDG2(i6,1)=x(i6+792,1);- X$ x1 D Q! O7 c
PA(i6,1)=x(i6+816,1);
* T3 T" J7 s$ Z K: A! {8 L" ? SA(i6,1)=x(i5+840,1); Z1 f, l7 E6 E! A6 e1 T
end: ]- R+ P+ s. L- R( q% I& [" ]
for i7=1:24( M/ ~7 E8 F5 Q& q& l4 a
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);) r' ~, U5 ~6 ]3 S7 H+ F9 k" t
end
7 j u+ t# J- P0 Vdisp(['运行成本：',num2str(fval1),'元']);
. _; Z& i4 t) D: T! i* }subplot(2,1,1);8 S$ g" M/ r, {3 k
plot(PA,'b-','LineWidth',1.5);: Q3 o, U, a# K* m9 S+ v
hold on3 M8 [3 D$ V8 Z+ @. M3 k
grid on( y: N1 H0 a0 A) x h
# }0 m4 {% e% F/ v8 ~6 w7 T; w& `
; b# W2 a L# j1 ?/ [; n6 M
|
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！