|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：
6 ^/ V" z1 y" O5 N
, o( V0 _% k7 Q1 B0 Q3 Y5 P, \function f = adn(x) %主动配电网层目标函数. c( b2 e4 w: w+ Z* [
global A b WT PV j1 PA Psub %全局量8 X, a' a& T/ w) W' X
ADN=0;" k0 o1 d& F5 |& ]
for j=1:24% G2 p& ^8 K, F0 n
ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量0 L" @. }1 \/ f9 k J' S
ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量3 E( O4 k; }& j1 i. [' G
ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量4 a$ K" ^& G) L7 B& E3 a0 U
ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量 / x/ x& \$ n% d
ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 2 `4 G& p; H- g
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
$ g4 I% h+ S$ ?* [2 oend
1 {1 p' ]. g( M3 G1 mf=ADN;
% N! B1 m! L1 M9 ]. w" x+ n0 n: p! g& ]* o9 P3 p: {0 z# w
9 K8 [) s* |& b i* g& c; r
7 K1 R. X/ Q7 h0 @0 }5 }
clc
" _; [. Z! c8 R7 tclear all& R5 U# R) h4 Z! g
global A b WT PV PL j1 PA Psub
- j! K7 u( C; O( ^' QA= [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 + B8 N5 a. `/ K) T# m9 n7 [7 r
b= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh) u+ |: c- B3 z. A* P
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小时实时交易电价 元/MWh
' u5 ~ B$ S) M% y2 Y, u- L7 LWT = [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
+ \" H7 S' }' \. Q8 h9 |: nPV = [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小时输出功率 MW3 W H4 g1 C$ c
WW=WT+PV; %风电、光伏发电出力之和 MW
" b2 j2 \+ D6 XPL2=[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为平衡节点） MW4 i6 q; G" h5 p2 a l0 m% u
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
5 U) `5 N- `3 V" c9 JPL4=[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
6 b/ y2 e3 j" Q4 _, H" U- ?" z+ oPL5=[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+ C1 n. y* J" P5 ?1 r, t* M" _
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
5 B: \6 c: Z2 Q! y8 J5 GPL7=[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
; v# {% T2 T' X( W) u8 o* y% CPL8=[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小时负荷数据 MW0 \( w/ g) o; g+ S: m: g* s5 t
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/ ~9 z. q: b P- E& y* [
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
* h4 b, l% X% Y! g; x: T [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
|. l" L7 C* {- L) {, R/ n/ R5 cPL12=[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 : G9 q; Y" Y+ |$ h* L6 u. E6 {
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
( M% b7 e5 ]$ Z: I" h& y+ E" 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
' q' C& A# u) \1 z! a# QPL15=[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& Y8 b2 X6 N/ h) t& C: u
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. a T* \. w& G! m! f3 H. e9 _
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! ]- V3 d3 H; X5 y# V
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小时负荷数据 MW
* \* Y- U& z) |1 ^ p) r3 @0 {9 P9 {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
9 G7 `1 B# O9 w0 c6 }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小时负荷数据 MW2 U4 ~) U& A% z
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
6 j- L9 \8 ~ n2 k& `0 N$ R) X2 S; e5 q1 UPL22=[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
, T5 H- t- R% s1 p# uPL23=[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. e) R' j( e1 @, z/ P8 `) W3 D
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. A$ p8 ]* U0 f1 U
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小时负荷数据 MW
& W+ W( Z7 F0 v. H8 T- l9 ~4 C# MPL26=[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
3 u9 p& E9 |1 Y6 {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: d. n0 X x/ b9 A" l* x
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
' N1 |3 x9 U+ W! b' D' fPL29=[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 _: }; C3 C% BPL30=[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: y1 \# X$ P5 m# t- P; [3 t& K( I& N
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
1 w; B! Q. j* |( x" }* O5 lPL32=[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! D, y s$ {# B& H
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; s; D5 r% z4 q% b3 F! y ^" z
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
$ Y4 y5 U0 d0 {0 yPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)* I0 s4 i4 W! U6 Z8 l( P0 x& P
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW
$ w. ^2 n" f3 g! uPDG2=zeros(24,1); %初始化分布式电源DG2电量 MW
4 o& F" D' v& a2 OPsub=zeros(24,1); %初始化Psub矩阵+ k S* A5 p/ I( o1 u* r o) T8 `# b
PA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
4 p/ Q( j3 ?& m8 gSA=zeros(24,1); %初始化聚合站电池容量 MWh
# {- w$ e) }4 J8 d: Wx0=zeros(864,1); %变量初始化 , K4 m( ` D. o& f% N0 R. M
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%. ^+ L4 T6 @/ P$ Y! z
Aeq=zeros(793,864); %等式约束系数矩阵初始化& z+ i. G4 |; \; W
beq=zeros(793,1); %等式约束结果矩阵初始化
* q( ^) Z/ m( t0 H) lfor i1=1:793/ M% [1 h4 e, O# u
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%8 n H5 \) K$ S# G
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL! [# \) r$ p& s$ F1 ~$ Y. b4 [6 V
Aeq(i1,i1)=1;
; `3 U. ?7 M5 b6 A beq(i1,1)=PL(i1,1);% M9 `/ {8 ^0 {7 w. c, j/ G4 k
end3 V6 Q2 z6 H; l* g
if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV
" ~5 d" f7 ^* a) l3 I& X9 Q2 _ Aeq(i1,i1)=1;
! y! S I$ _) r$ X. x9 Q* l5 @1 a" s beq(i1,1)=PL(i1,1)-PV(i1-144,1);5 B& J. ]% y5 i5 }
end" V$ W4 P6 D; _/ v o i9 ~
if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
$ @2 N2 M; R( K. L" t, {0 O Aeq(i1,i1)=1;" v- R0 N9 r- j- i: f
beq(i1,1)=PL(i1,1);
- D X! G4 N5 n end
- B+ \1 j6 i; _ if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL
" D0 {% J; P% N. q Aeq(i1,i1)=1;6 `9 b0 {' N: o" A! }2 O' e
Aeq(i1,i1+504)=1;
3 J* H2 E" U, ]' h2 X7 S beq(i1,1)=PL(i1,1);* u9 | |7 b% I' ?$ I/ R" M+ x. }; ~
end* \5 D; @- w1 W; X9 ~' H& _# ~* Y
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL( H! _& M/ ?+ m2 V+ A
Aeq(i1,i1)=1;
+ {7 a A8 F: e c1 x L) \ beq(i1,1)=PL(i1,1);: E) I1 D# }2 N$ F, c' h& n4 R
end4 m" R: c% J a9 `$ V
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL
2 E0 t% r, {* i" P Aeq(i1,i1)=1;1 h1 |% n9 D5 O7 L* N9 t9 u" P
Aeq(i1,i1+384)=1;. W% V. ^) j* N t5 g
beq(i1,1)=PL(i1,1);
) _! i9 g/ n* m4 |: N end
" K2 L% b* o1 _0 w4 W1 x* n if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL% s1 B3 M0 O4 D( u+ h" c n5 @3 {
Aeq(i1,i1)=1;0 j8 ]5 J$ |8 W+ S
beq(i1,1)=PL(i1,1);
7 U3 W! p( v# y: p: T1 o end
0 t, D6 V& v8 D3 ^- K3 c2 n2 j if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT4 T9 Q2 {. ^4 X# t! {; H
Aeq(i1,i1)=1;& y) s7 v$ L, f! ~# j
beq(i1,1)=PL(i1,1)-WT(i1-456,1);
! N& L7 t& D: x end
4 i5 H/ E, f6 V; O# A if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL1 j: x9 f1 T5 M+ U4 ]
Aeq(i1,i1)=1;
' Q0 }$ A3 k" B1 D beq(i1,1)=PL(i1,1);# `/ Y+ A5 i' l$ U
end
! W, |! I9 _3 ], O- M if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL
! {( Q5 p6 x# t+ i( B5 Y3 t5 N$ Z Aeq(i1,i1)=1;
9 x }9 U, Q, L2 a$ E Aeq(i1,i1+144)=1;
) [5 u4 A) T0 y0 \9 }4 b beq(i1,1)=PL(i1,1);, T: r: @ A, M
end0 J$ G2 ^. v$ R+ v2 q
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL
* H" {% Q0 v9 e0 | K& J Aeq(i1,i1)=1;, W7 D# r* E ]% \) p d1 Q5 e% b- {
beq(i1,1)=PL(i1,1);- l( I. O7 s3 @& \( k
end
7 | V! ~" Y- v0 M1 l' Y%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%* k" J; k; D- C; P
if i1>=792&&i1<=792
3 @7 O' v) a; q& {; T Aeq(i1,i1+72)=1; %对应电池容量SA! i8 c$ Z+ p& M. H+ u. Z# E
for i2=1:24
% L0 v3 M! c" _* N if i2<=i1-768" n+ d! t/ t; b, N6 j
Aeq(i1,i2+816)=1;
# [& I+ A! a% j end' q- x" p0 `( f N( l+ B
end
! t; R# F! P$ n beq(i1,1)=5.4; %5.4MWh
) r! E: j" {- p end
# r5 a2 x! K4 ^+ ^5 f%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
' r* M j& s3 l* V if i1>=793&&i1<=793 %对应S24=5.4MWh3 Q* l0 F5 h/ t8 b" v) o
Aeq(i1,i1+71)=1;
, x L1 A! A" u2 A beq(i1,1)=5.4;
# e; H. B0 s9 B0 {1 Y; k end
: ]5 V5 s7 D+ [. [+ `end
6 J( j X+ a" j" W, {1 e%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%' }: z: t. p9 T7 ~. S' X
A=zeros(48,864); %初始化线性不等式约束系数矩阵+ V7 \9 K* F/ N5 r" a: z4 c4 g
b=zeros(48,1); %初始化线性不等式约束结果矩阵 P y' v( n8 q
for i3=1:48 & I) r: L. l) \, a0 v% @9 M
if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
. ^. i( V3 Z+ v$ d4 U A(i3,i3+816)=1; %对应(PA)t的系数
* K- x- J1 j( e4 Q6 q5 { if i3==1
8 v4 v& {1 o5 L3 @" B7 S3 z3 B4 W A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数4 L/ Z( [; Q2 ]3 g" |
b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵: C) Z2 c8 i) @) T/ E [
else( R4 [5 L1 o. F$ e
A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数4 N/ V! r1 s7 Y4 d
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵
, P" O; G) w: M end/ [3 `3 Q' m# ?6 q7 b
end 8 @/ z/ Z* O; I' L( a
if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
) }% Q- R3 c( S, K1 L A(i3,i3+792)=-1; %对应(PA)t的系数# |; s" O- L9 L$ ?( Y* y; x
if i3==25, j9 Q2 ~2 s/ z" R/ t9 v8 h
A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
8 x) G5 V# o9 ` b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵8 x0 M- b7 O% t+ s! ` X. Q
else" s+ ^8 ]8 M0 s( R$ Q
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
- m% W% z, ~9 Z- [3 D) h" s b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵' |, M! j/ b) G0 k* L6 Q+ ]
end
7 [/ P0 j$ E# N3 y+ M" q end+ J* N' ?# c; w+ ~4 q
end
4 v# V2 m$ w" \8 x& r: F, N%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
: l6 E$ Z! x# M0 PIb=zeros(864,1); %初始化下限矩阵8 E# M, o( e; I' }, d# t9 L
Ub=zeros(864,1); %初始化上限矩阵7 L+ x& S* O" \. o. M$ N! W( j
for i4=1:864 %Ib Ub赋值
/ z4 b+ [ l* Y) r4 Z' V if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
, L$ y4 z5 m! b- s Ib(i4,1)=0;! P- { P. l! F+ }
Ub(i4,1)=3.73;: l5 {( X! |0 s U' t `& `
end5 J" q/ {9 Z( |0 _
if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW4 A6 d H y* Y3 l* k. f
Ib(i4,1)=0;
; R! L- ?3 A! ^ Ub(i4,1)=3;
/ w+ J7 M/ V+ e* W' b9 I6 C- R end
: L) Z% k. }6 q. k | if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW
, w" ]& \, d, D% ] Ib(i4,1)=0;' h/ W' h* [ }* s- W1 C9 }8 i
Ub(i4,1)=3;# ~) T0 L" `0 c. |( t3 G
end
1 T$ c, m; r Y' i if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW0 D0 P2 y$ h m; j% Y
Ib(i4,1)=-2;
, c7 b* Y( l$ v3 v" S Ub(i4,1)=2;, } u! R* T2 x8 s( z
end# V3 K! o" K5 v# }- G
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh
0 [" Z' t0 e" a" B+ [/ d# ]5 A. W Ib(i4,1)=1.08;, ` v- h. `. v/ z/ X
Ub(i4,1)=5.4;; g1 K" [# @9 }. n- x
end( u7 c3 g5 e7 |$ H
end, p: _) |' ^" r ?7 ?: t+ ~
%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%* e( M# J5 \ U) f1 H
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);! ^9 U: y8 U& ^! {1 z( v0 d1 V
for i5=1:768
8 `8 d( ~" S, h# _4 o# c Pij(i5,1)=x(i5,1);/ P4 y% C/ U9 T; X" P' Y$ r
end
4 L# `# p& N, Z5 o; {: afor i6=1:24
' x' n" h( L# p; j PDG1(i6,1)=x(i6+768,1);
# z* L) g9 [! A PDG2(i6,1)=x(i6+792,1);
B0 v# d! J( M0 Y3 N5 S/ t9 d PA(i6,1)=x(i6+816,1);+ Y" W; o* Y. H V0 `
SA(i6,1)=x(i5+840,1);8 ?0 t$ K) F: H- d/ ?: T
end2 _: Q' e8 `) j8 L. U% H
for i7=1:240 F0 ^! C0 b6 Z
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);( Y1 m7 E, `! V! k: g
end& V' l; v& B: `. O# Z
disp(['运行成本：',num2str(fval1),'元']);
% w0 F6 h1 t& Qsubplot(2,1,1);; t8 q# U7 ?9 q& [" P9 L
plot(PA,'b-','LineWidth',1.5);
7 E* w ]( e5 [2 e6 C7 y$ B0 |% ]hold on
. P3 I1 t, P9 `0 fgrid on) t! G* S9 q7 v0 E9 u
( z" o$ P. H: _) o
1 \0 @% D4 s8 X: r! r
|
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！