|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：
3 V, p; |7 p; f! Y/ c. O
6 Q- T* S7 Z( l1 gfunction f = adn(x) %主动配电网层目标函数
; x" p1 n2 b' @global A b WT PV j1 PA Psub %全局量# _% ?# @% E$ a
ADN=0;! P0 {6 k6 Z0 \2 M, J
for j=1:24
# V# M$ L! {" `5 e! F8 _* p ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量
% @+ T, ?3 ^3 k" S% S ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量. @- d+ S) N: a% O
ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量( t4 E- R: D! W9 K
ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量 # t5 i$ f; ]3 g& u @
ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 $ W# L9 u) I9 r) J4 h$ F' y
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本- u1 v) O$ H! C
end+ z7 b; o7 |1 Z& x* U( y' f
f=ADN;
! f9 n& v5 s+ d4 d* T% Q% p' J5 H' A1 b+ }4 _9 J
" M: `/ I! N0 }. m, p1 V3 h
5 ~; x" p3 w3 T- h& d2 Eclc
4 o: h# T7 u8 s: h: Sclear all
# o# D3 ?9 j* L% R) Hglobal A b WT PV PL j1 PA Psub# x+ o% ?9 c7 l3 o7 v/ n9 Z$ _/ _3 y# ]
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
9 z3 m3 w* |% F( c, y9 D$ H4 fb= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh
. P. }2 s Q' u" Sj1=[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; d. I; m2 s$ J2 X: c
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 ) \4 |9 J% Y5 {+ J
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小时输出功率 MW5 }( R# [* [9 Y
WW=WT+PV; %风电、光伏发电出力之和 MW
! Z7 D6 z+ I+ p O& [; ^. jPL2=[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% C1 M0 z. L6 r j8 W4 {
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' ~6 t6 u# ~& g* 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
2 Y9 P7 k4 L: NPL5=[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
+ L) E3 @* u# ^0 HPL6=[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小时负荷数据 MW3 v4 y1 C& i% L# Y# k' s
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
6 Z6 P& u) k% K- E& r) ^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% h: F/ ^. a( d3 F8 ]( S
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
' |$ [- `7 B9 s( _9 u, 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+ N& Y/ ~. V- F. d5 a( ?
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小时负荷数据 MW5 |+ c8 }% L" a9 P
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
* f& j6 W$ \4 H+ s' v/ APL13=[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小时负荷数据 MW9 e0 k* T& o* s) V6 t5 V Y
PL14=[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小时负荷数据 MW1 P1 _9 ]! U5 k
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% \& s: c8 _. J9 h4 a/ G4 e' t5 f
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
! i1 ` O) o# x; x6 NPL17=[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 Q, {1 f5 [$ p2 iPL18=[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$ I8 u( q2 r. N9 r9 p- p# C0 D. I; W
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
) v) L2 Z; ^4 Z: T5 UPL20=[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
1 }7 a! e5 a: ?5 T- O: R+ o; jPL21=[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, e2 ~4 _8 Q3 D/ @3 F! l. x
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 # Q. I7 j. v4 L5 O- `( z3 b- o
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. z' p( }2 j1 R4 S8 G+ U) s
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
8 P! s4 t* r% UPL25=[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
* B2 V9 X( o' A, c' h2 M- JPL26=[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
. H( [/ |3 H/ @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
4 u. K( {. ~: D6 m: I2 PPL28=[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; S/ g( ]. j) J3 g
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小时负荷数据 MW9 v3 L& X, D3 m2 X
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
. [& w9 e* F4 }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: d. Q' ]( t- C2 a
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
; e$ g% E2 A! N. a* S$ Q+ c# h8 ?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
' x+ T: Q$ g& j; R8 f/ }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
8 T; }6 K/ W N+ g0 L' G( V* tPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)* J- k- X q: |1 X
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW
8 Y9 f( N- r, B+ bPDG2=zeros(24,1); %初始化分布式电源DG2电量 MW/ Q) ?8 a1 }8 P. J3 i
Psub=zeros(24,1); %初始化Psub矩阵
1 I+ l+ ~+ O, VPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
8 b S) d- o+ p* l+ cSA=zeros(24,1); %初始化聚合站电池容量 MWh, h+ x* p4 {+ g# ~. u
x0=zeros(864,1); %变量初始化 1 w p% F4 Q% A& U
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
/ h: ]6 i3 U7 H W4 CAeq=zeros(793,864); %等式约束系数矩阵初始化
I8 {& R& c$ H0 G) Z2 ]beq=zeros(793,1); %等式约束结果矩阵初始化
7 U; ^" c' {5 k2 W, t1 H, t* Afor i1=1:7939 y0 V$ U3 e$ }) X B5 L, g+ [8 E" N
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%
) _& r- ^3 Z$ ]9 e7 a+ `9 d' i. V if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL8 X4 D' n$ T+ z8 ], H8 S$ x
Aeq(i1,i1)=1; P5 a, l, u) o3 A u3 M0 b
beq(i1,1)=PL(i1,1);: h% Y- t3 j' L( ?( B; |4 {* D
end
6 |, u9 K1 X) }$ p0 H if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV- d M+ f1 [3 G7 W
Aeq(i1,i1)=1;. l. b. ]( ]9 ~; |! n
beq(i1,1)=PL(i1,1)-PV(i1-144,1);
% a6 h, j6 H7 X end
2 p+ v% ]6 h1 `; S if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
$ s l9 o2 J4 y' z# y Aeq(i1,i1)=1;
. F2 F8 ~0 U; f6 f: m beq(i1,1)=PL(i1,1);( S/ J- d- s8 G5 r Z6 D) M
end; x: {) ?( @9 z) c3 Q
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL1 b8 X8 ~5 O! W
Aeq(i1,i1)=1;
. ?% Q. L T2 _1 K Aeq(i1,i1+504)=1;, H, Y9 V+ t$ G1 M6 G
beq(i1,1)=PL(i1,1);
" n- D7 S. X# h end
9 O; F, D7 Q6 e! g6 ` if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL- y3 f$ ?! s/ }. \3 F/ k
Aeq(i1,i1)=1;7 Q. y) F( P7 m7 x* ]0 _
beq(i1,1)=PL(i1,1);. ~" K. ^7 `7 h" c
end2 L* j {4 K; e1 M' i4 c v2 G$ }
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL
: V: @% V0 |% _8 Y: e: U" N Aeq(i1,i1)=1;
4 t @8 A+ u4 q8 {* x7 N( j* u6 j Aeq(i1,i1+384)=1;7 W5 m+ ?4 }/ h
beq(i1,1)=PL(i1,1);
m" ?6 Z k! Q1 w9 l9 Y# S& G/ y) B! p end/ R. N7 `( ?6 ]* w6 O( r4 t; ?" }. l) U
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL6 _! b4 e1 ~- x# ?9 j
Aeq(i1,i1)=1;
2 O$ X" Y" f2 U" {+ T5 x beq(i1,1)=PL(i1,1);
0 ^' L6 s4 _% p# h. W# B end+ F3 t; V) a& g: b
if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT
0 a" n; h$ m6 K" z. x Aeq(i1,i1)=1;
6 K G: Y/ x6 ?/ s5 b+ p beq(i1,1)=PL(i1,1)-WT(i1-456,1);
. j4 L5 P+ y, j end
5 q4 }# M J: m if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL
; Z$ d' l% `8 Z" x Aeq(i1,i1)=1;7 o) v, L5 n, ?/ W# v ~# c& j8 z
beq(i1,1)=PL(i1,1);
$ K: f' A. {+ o4 _* d end
; [1 w2 L" D3 m l8 G if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL B# |* ?% K/ A/ ]* u2 I
Aeq(i1,i1)=1;
2 \4 K( h- l/ l9 ?3 F; z Aeq(i1,i1+144)=1;
- E3 Q5 `; @/ N7 r beq(i1,1)=PL(i1,1);
/ j9 g; ^' @ [4 \5 p end
: a/ z3 M8 i: P1 q% z4 P if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL
. g0 N$ |0 P2 m Aeq(i1,i1)=1;4 G! ]$ a; V0 t) X+ j0 Z& L9 _
beq(i1,1)=PL(i1,1);0 R) w9 o0 y" L3 x% ~: L' J- f' c
end
4 f P# I5 x7 G* o2 M%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
~. b1 n. j: L I3 { if i1>=792&&i1<=792
6 X- N7 L3 W" B6 Y: E c Aeq(i1,i1+72)=1; %对应电池容量SA |2 ~ ?6 i6 N& J# [6 R
for i2=1:24
' n6 v! _& C) z% p i if i2<=i1-768$ u4 o9 _( q5 W* f
Aeq(i1,i2+816)=1;% X. K& S0 s. ?/ J; M
end
4 x; n1 U5 [, T+ V: A: R" S, \ end
2 Y* U1 U7 v+ f4 Z r beq(i1,1)=5.4; %5.4MWh1 I2 z7 O5 j% l4 q x. H
end
, d5 H+ Y3 w0 v) ~4 ]0 f% ~%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%& _8 m9 ^$ f' R8 o6 N6 k' q
if i1>=793&&i1<=793 %对应S24=5.4MWh
% |" m# M$ s" X. i0 t" B7 F4 u2 L( k Aeq(i1,i1+71)=1;0 }; x" g9 y2 {/ Q
beq(i1,1)=5.4; B0 r" ^; d# H0 I, O
end) `+ h, S c3 Y, w' f4 q
end
7 d* t9 V8 g J" h$ u C; h%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
, d$ ^( s& H7 Y: p$ TA=zeros(48,864); %初始化线性不等式约束系数矩阵, Z* F. h1 U2 T
b=zeros(48,1); %初始化线性不等式约束结果矩阵
* _5 Z+ w- X* B: W m8 Tfor i3=1:48 . l! b4 N3 {" I9 d+ F% k+ }
if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
: r2 f b( ~/ A2 M! O A(i3,i3+816)=1; %对应(PA)t的系数. h( u+ O [6 W# I+ W0 G
if i3==1- I0 a' S* K" a; g0 M* {8 k
A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
, M& w$ w% G0 M& E# P6 X; I b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵9 ^+ y9 d$ _# b2 ]! [" ?' ?$ l
else% N4 `8 c1 Z) }5 m
A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数
) u1 M0 F1 I7 H/ Q3 B- e b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵
6 N& R9 f- q# b3 W7 ^ end. _8 K- V$ y# v+ r' Q) r, y4 J
end % ]- ^4 w8 S& \0 Q' n% u" H) B
if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-19 F8 ?* ?5 I" A, U5 r& {+ G
A(i3,i3+792)=-1; %对应(PA)t的系数
3 d6 G7 v* r* X( a if i3==25
7 ^" z( U; L# u" F/ q' Z A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
" `' X0 H# H1 e) U( y b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵1 u& x; q% r# G" `, b+ Z6 U& N
else* {9 D( d/ j5 E O9 h6 p r
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
3 k" b+ z, i8 i' j1 o b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵% V' o; M0 r1 ^7 w+ @
end. R5 J( Q2 u1 j5 M* Y
end3 |4 o& y3 m/ {( j: A1 _! R
end1 W* P5 @0 G w. _
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
7 I. c2 C5 o$ J0 v( W2 f7 _& NIb=zeros(864,1); %初始化下限矩阵
& L1 L/ q' c- g/ h; R4 @* IUb=zeros(864,1); %初始化上限矩阵. B8 {) q4 m4 g9 G( h; U: o% G+ w
for i4=1:864 %Ib Ub赋值: j8 g, h/ J+ ?& I7 D% {
if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
8 e. I0 i6 C2 |8 O Ib(i4,1)=0;
' B* e* h1 m* S3 O Ub(i4,1)=3.73;6 }9 R1 l5 S7 R' ^* X" `
end, D& [- \! I3 ]. s) f, G" ? D- i
if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW
/ D$ Z( u3 Z+ R! A6 a5 p' A9 b4 ^* k Ib(i4,1)=0;6 J9 p% Q, m. R: n" Z5 m
Ub(i4,1)=3;
1 g/ o8 }$ j' c1 r* T0 ? end
; q8 k5 V4 R1 T3 R& U6 ^( A0 \6 q if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW
. v( Q( s- ]5 G* a. W Ib(i4,1)=0;# k& q6 ]$ N7 m$ q- q A0 I
Ub(i4,1)=3;& }" c; ^8 S0 M
end
+ m4 i' ~$ i1 ]1 Y# a3 Q if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW
3 ~* g- h2 ^" a' F2 P& f Ib(i4,1)=-2;8 \) h5 K8 u7 i( N1 [3 a$ \; h8 D0 O
Ub(i4,1)=2;# G A) I8 B, W. l% Z& l
end6 v1 \* I Z1 u
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh7 N2 [! B8 Y/ k: K
Ib(i4,1)=1.08;
4 C8 |8 U4 I& o" ? Ub(i4,1)=5.4;+ E% G- _# M/ Y4 V3 Q5 ?0 R# ^
end
9 F R# `) Q( p" W3 l8 p1 g' r* Kend; I3 q; J2 r, u Y! L0 e2 [
%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%4 B# _# D% t, \8 l z
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
: y) u( O7 u1 `) Zfor i5=1:768
/ b) V* J, y# ?" \0 @9 i3 E$ Y5 V Pij(i5,1)=x(i5,1);
' w# I6 W& o a% Xend4 U- A L# {; G0 k" }
for i6=1:24% q- i8 }* w- J7 Z
PDG1(i6,1)=x(i6+768,1);" |: ~+ a$ Z- Z' ^( b( @# ^
PDG2(i6,1)=x(i6+792,1);: [. N# V6 T1 y- y- A' _
PA(i6,1)=x(i6+816,1);* z4 t/ m0 C s5 f% P8 h
SA(i6,1)=x(i5+840,1);
. N( O$ G/ ~$ z$ S# W2 B! aend
; H: _) F* n+ a- k& J5 p3 Tfor i7=1:24
# d0 q4 X, U& i. C/ s6 w0 f" 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);7 f! v% N, u7 R- b, I: I: F
end, p, F# p3 v) e! o7 a
disp(['运行成本：',num2str(fval1),'元']); ( C5 p6 l- w1 {4 v9 d, C
subplot(2,1,1);
9 G v {7 j- y2 X; j/ D3 z/ H" Pplot(PA,'b-','LineWidth',1.5);! R% _1 Y' R( x) [# i6 M
hold on6 ?8 M1 k% e6 h) R! K5 P3 n2 \
grid on! J. y! U- D/ U. I6 V, e3 [
, J2 k& `3 a: N0 Y% ~+ G; S
/ v3 a6 q! K1 [% D( x- p$ B A
|
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！