|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：7 i6 E& O' ?. x* }, O, I& S7 i
. \; f% Q n: j* b
function f = adn(x) %主动配电网层目标函数, e7 V F8 B0 k: J9 G
global A b WT PV j1 PA Psub %全局量 k# G& b. a0 z! t0 ~9 x
ADN=0;
3 N3 t, L0 M& {% j3 L I% j& bfor j=1:247 F9 E$ S1 Q% m9 ]* `
ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量' k# _' ]6 C2 g7 N* S
ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量" C6 p( X: C" i4 P8 X
ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量2 d3 Z/ ]# n2 z0 n
ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量
! k- u3 ~$ j. X! N/ z+ S ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 ! R1 |2 c; z' S
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
" ^9 b% O' e' Z2 U" \! V, m* mend
& |" P9 D7 z; F' yf=ADN;
; d) x6 r8 C8 c( P* ?5 ]1 V# N5 a8 d8 q- d- b' D. c
4 k5 a. ?0 l( R! `% a/ E6 X$ x H) r% ]
5 g! H8 Y, J( H8 J/ y8 \1 bclc
4 I6 ^7 i3 b1 V0 M" vclear all
+ x$ k8 j6 F- ]" dglobal A b WT PV PL j1 PA Psub
$ D* O9 m# o& \8 T5 V. wA= [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
: N" w. S2 ]% y: u& Vb= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh2 A* R* e5 L1 m" x. [, X
j1=[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小时实时交易电价 元/MWh7 N/ w; Z* W/ g0 \/ X) O& E
WT = [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 " ?' W% ]% F& S; j, |% W
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
& w; h8 b1 G5 Z5 N6 ~! xWW=WT+PV; %风电、光伏发电出力之和 MW* G! c% f3 h3 @0 f# p" h w4 O* s* ]$ r4 D
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
( K, a+ R& ]5 T* Z9 Q( L4 ZPL3=[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" \' N4 d& l3 _( y" h
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
0 @' m; }# Z" Q3 `; e' j7 K0 ZPL5=[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小时负荷数据 MW7 T! x& l1 r% ^6 i' C: ]: k$ k2 {9 o6 Y
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小时负荷数据 MW8 J5 K) ?9 F7 b* h1 o
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
" ]/ a1 p4 n! e/ APL8=[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
' O# V0 Y9 |. r3 gPL9=[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
- ]) H3 ]1 u% w/ v! W1 x. H6 ^PL10=[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
$ e8 I! z1 a" X I6 v* e; W% jPL11=[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" l+ n# Z5 h; C. k( G% a; `- @ 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
- a0 M6 y$ U, r0 H8 u' d* hPL13=[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
1 h! z7 _& S2 s( _# a2 A' uPL14=[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
: }& A6 M* W1 q8 APL15=[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! N& s* _# q( l( b3 r
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
& I$ A' V& x0 P2 b* x3 U- GPL17=[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
# `+ K, E/ @4 E' VPL18=[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小时负荷数据 MW9 G7 e2 G1 p9 L. k2 h: D
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- d+ K @5 p% {) }
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小时负荷数据 MW3 h) X# A7 s3 d" b" ]- Q( V
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. l3 Z9 r% ]: P+ ~# P4 z
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
. S1 {: ?7 \$ t9 QPL23=[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 {- ~0 K2 X" }9 z) A- u4 u4 y0 j
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% M" `9 b% X) w4 Z6 w2 l7 A
PL25=[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小时负荷数据 MW0 p1 S! y2 v4 n6 `8 M% h
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; E( Z' c0 ^, v; O+ M
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小时负荷数据 MW
' X4 \ x% d+ B* O2 b+ GPL28=[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 h3 T/ \) ~2 U# U B& w
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
! m/ |6 K# U5 H' hPL30=[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
$ G" O1 h; C7 D7 \/ ZPL31=[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
0 r. G( u# X) e8 M5 y$ X/ G5 d4 Z1 T; N0 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小时负荷数据 MW2 {6 x/ @1 m H6 p2 K/ s
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
2 R5 W: ~* y V# n8 MPL=[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小时各节点负荷数据 MW5 P- \: z4 J/ H3 J4 a/ R2 z& k
Pij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)6 [1 }/ X; k& F8 c5 @
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW: z y+ n( h5 j! q y" \" g0 n
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW+ R3 b$ o' {0 i# ~& f
Psub=zeros(24,1); %初始化Psub矩阵
1 o5 Q" r8 M$ F( N" WPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
2 q! o; ` N0 A% u+ ^! c! ASA=zeros(24,1); %初始化聚合站电池容量 MWh
# E9 s4 F) Q; Q' Dx0=zeros(864,1); %变量初始化 7 H7 g. U* R% v' [0 X/ s. @/ o* p
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%9 n: T8 d" h5 K8 V: K4 ~; p$ F
Aeq=zeros(793,864); %等式约束系数矩阵初始化
8 Q& b+ }( p1 m$ v& r9 qbeq=zeros(793,1); %等式约束结果矩阵初始化$ G: R d G @0 p
for i1=1:793
' `: O% g( j, V5 C4 `) Q%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%/ \2 C O, f4 Y8 X# R {; w5 G% H
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL
6 h E% ?; k7 Z Aeq(i1,i1)=1;
, D8 s- w: T7 }2 u1 |/ { beq(i1,1)=PL(i1,1);6 L' c6 n& q+ z0 g+ T
end
. g4 |5 N' o8 G. v6 p. ^6 Y if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV
% `" q$ w, ~% N$ e# Z {3 \3 d- W Aeq(i1,i1)=1;. r! k& @7 B; l7 c/ |/ Q
beq(i1,1)=PL(i1,1)-PV(i1-144,1);
% [9 {- E+ a$ z# }! k end1 B D. R5 r1 o
if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
7 W e2 K1 ^. R' p, [- i0 v Aeq(i1,i1)=1;
5 y9 M& w! A3 U5 a6 c* j1 c* G beq(i1,1)=PL(i1,1);9 _& x: Q+ X& Z; B* _
end
" H0 V' B+ y& P& ? if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL- Y, A8 J6 v7 }+ U
Aeq(i1,i1)=1;0 ?9 b) s3 {4 h" I3 O
Aeq(i1,i1+504)=1;" p3 |! Y$ |, f4 n# _0 z5 r( ^) J
beq(i1,1)=PL(i1,1);' \6 h% [* g' \+ `, R/ E
end
- s, F: _! e1 G4 z x# L if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL" m6 ^6 F) J% j- P! S O
Aeq(i1,i1)=1;' g6 ~; V& Q! j
beq(i1,1)=PL(i1,1);
# ~/ Y1 y2 e9 T* t end
* k% e4 O" m1 D j& n if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL, ~6 I8 \0 {" r$ B( s7 m
Aeq(i1,i1)=1;
2 m7 x% r$ V* U" e Aeq(i1,i1+384)=1;
. l, ]# e5 _4 M) X: Z6 ^ beq(i1,1)=PL(i1,1);
! {8 W$ F* `! d! K end' C! i. x) ~) h8 o
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL0 f9 J; B! X' x( C# \
Aeq(i1,i1)=1;
h4 r# D5 L: A" W( Q& T beq(i1,1)=PL(i1,1);
1 t" U3 }( \4 ?, @, ^+ T: k end
8 X7 n8 n# m. R; P9 E& `$ i( d if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT
, j1 w' \- x: ~2 E3 ] Aeq(i1,i1)=1;
& ^* `- m: R$ J beq(i1,1)=PL(i1,1)-WT(i1-456,1);
4 L( s" ]' d0 e1 }+ }+ | end3 N6 k, N3 Q% {# O
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL- Q; f* ~0 m9 H+ M# O& `$ [
Aeq(i1,i1)=1;1 [( g+ k# y5 A8 p. |
beq(i1,1)=PL(i1,1);
7 H, B5 S4 z' s# _5 v D end
" o. D' W& l4 E if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL8 i! e$ k9 _/ h
Aeq(i1,i1)=1;: w4 z$ ~7 B4 S) ^9 `5 N9 c5 b
Aeq(i1,i1+144)=1;1 @5 }9 }( V0 b' Q# `& L; Q
beq(i1,1)=PL(i1,1);
% h4 M+ [: M# p% O% ^% u% ` end
# r( r( G7 S0 q! m: k. | if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL" C! v8 `0 ~( F- q& h; q, x
Aeq(i1,i1)=1;
$ g+ W9 C/ P' B, U beq(i1,1)=PL(i1,1);
% d: V8 X# u$ E& T0 y end
& m2 P9 i }& i9 T0 y%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%' f$ ~% h8 _5 g2 O5 t F$ x5 `
if i1>=792&&i1<=792
+ O# x( L* V. r& [* r$ F; R Aeq(i1,i1+72)=1; %对应电池容量SA# ~' z4 Z, x1 [% ]3 m$ h% p6 v( J
for i2=1:24 P _9 w6 [/ u k
if i2<=i1-768 `$ o) C1 R9 g5 a: N' p; O) P Z
Aeq(i1,i2+816)=1;
' F( u% V- G' j3 B end
3 v( |9 d3 _9 P- F$ Q end
. ~; T |+ `/ z2 N1 K5 p2 Q beq(i1,1)=5.4; %5.4MWh! H; R: J/ ?* h! ^; E
end+ b+ ?- K2 y! q# Y# } [- x+ y
%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%: d! s* W* U! g! B9 W: G
if i1>=793&&i1<=793 %对应S24=5.4MWh0 k. q* k) P6 c# w5 b- b5 m" A
Aeq(i1,i1+71)=1;& k$ H5 U$ Z9 ?* ^3 C
beq(i1,1)=5.4;% g' p" |4 j& K
end7 |5 X* _9 `. Y. M) O
end 5 ^% @4 G+ s. G: B) l4 ?
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
' Q3 L% i1 ?- zA=zeros(48,864); %初始化线性不等式约束系数矩阵; b! j/ _) s7 r; J) h6 y% o
b=zeros(48,1); %初始化线性不等式约束结果矩阵- _4 |8 X# u7 n0 [" ~4 i( |0 o5 M
for i3=1:48 U( O2 U; O3 C. Q3 f/ }: K! T
if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t" T$ c6 Q3 O; I( B9 G- U
A(i3,i3+816)=1; %对应(PA)t的系数
9 c/ T1 Q3 N9 B if i3==16 |* @. m2 [; D, B% m1 ?
A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数6 I2 e' s4 _9 v# Q* f
b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵
5 k& {' Q1 Q0 X else
, R4 G0 ^4 J Y/ x' X, p6 C7 I A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数
) T! x7 u2 L9 @/ L6 U b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵6 f K L+ _& g- c8 X W. v0 k
end
( a- j+ B3 Y: U, y end
: y' l; @! Q- A& D$ |0 ~; |" S if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1( [. v2 v* K0 ^; _
A(i3,i3+792)=-1; %对应(PA)t的系数
! g8 j: R: x2 u' y5 M/ T, F if i3==25
. `5 z( L6 T% x$ U$ {+ | A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数+ \: e4 ~ @$ M- k0 V: P& Q9 A8 c% M
b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
3 @8 F- B, \$ E. c& U" F* U else7 _' w8 u T: e. x0 g
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数 ; Z. M! D9 v1 X p$ [/ E2 j5 ]
b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵" x$ |6 k( U/ C0 @9 W
end( F, V" u& x4 k
end1 Z+ x4 D- `3 j
end' s5 N s9 ^6 r. K1 j$ n' q1 H. C
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%* |) Q8 N8 p# f3 r6 ?) X
Ib=zeros(864,1); %初始化下限矩阵
7 m9 f7 T5 \+ x. f9 wUb=zeros(864,1); %初始化上限矩阵
1 ~# n& Y0 i! n/ A+ w- k1 u8 e2 lfor i4=1:864 %Ib Ub赋值
/ z* _2 y( H1 s" F3 B, E if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
. b& T4 s, v! J+ J. C Ib(i4,1)=0;8 C9 a2 W- F8 W9 w' I) u
Ub(i4,1)=3.73;+ i9 Y4 p. h( `7 ]+ R) T; q
end
0 v( i9 D9 [1 ~9 H l1 a, ? if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW
1 E( b; j2 u/ O( f4 {! }3 Q/ m! a Ib(i4,1)=0;
Z5 M/ }& F0 ^& z Ub(i4,1)=3;3 I) ^# }* |3 J* I
end
: @5 k+ O, ]& f$ L1 S if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW
3 n0 {' |5 O! F Ib(i4,1)=0;
9 ?. A* Q) P% q! H2 \& A/ s8 z Ub(i4,1)=3;0 d& |; n" Y' L0 J
end
/ p* C) W U. M4 R3 J, Z if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW
& u1 M/ T3 x1 ] Ib(i4,1)=-2;- K `- P) `# s, U0 e
Ub(i4,1)=2;( x0 f* V3 K" t5 h: s$ W% X! k
end
v2 e" r2 ?8 h if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh
! z& D) A5 l: l+ x: o) B Ib(i4,1)=1.08;% F- h" U8 N' ]3 Z: { K
Ub(i4,1)=5.4;5 I% d0 O7 @, Y# ~; p! |+ M n
end/ K4 H8 |6 E |) `6 y; E
end
* S0 h! `1 L, {7 k7 l0 T4 [" j%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%/ [ I8 k# `" `3 M: f' N2 s
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
+ i# x# l' p) o5 X9 b+ u. z7 qfor i5=1:768, O/ ^4 N+ u0 m) y0 [7 R
Pij(i5,1)=x(i5,1);! T1 R9 ?' Y7 [& M n3 Q0 r
end$ W0 l# \# l7 s
for i6=1:24 M! i( @( W' w, B. ]+ H7 M h
PDG1(i6,1)=x(i6+768,1);
$ r9 p z, z2 w9 T0 Q PDG2(i6,1)=x(i6+792,1);* |3 q; H/ k: d$ V9 \7 u, q
PA(i6,1)=x(i6+816,1);
, k* d: e1 b. T8 H A SA(i6,1)=x(i5+840,1);
& ~8 b; o) }! ?& nend
2 x8 y7 m/ u1 n8 Afor i7=1:24
# o7 B- t1 n7 Z1 W 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);; s7 p9 \/ Z# f8 S6 w
end) B( L" ^! O) J, l4 q
disp(['运行成本：',num2str(fval1),'元']);
% ]1 ?4 z4 P. z; T& X7 |/ F: X( esubplot(2,1,1);3 V' o- ?- n% }6 R# _0 g" J; g
plot(PA,'b-','LineWidth',1.5);
& x& n! T8 S4 p! {) Hhold on s7 ~: J& B5 }# O
grid on9 z- P! L# s6 @8 @6 ?% N
' T4 W) U/ ^& y, u
( C; ^( l, B. _3 z+ l4 n |
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！