|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
我在用fmincon优化求最小成本,但是出现:User supplied objective function must return a scalar value.求助大神!代码如下:2 s% S+ k r: q/ u V. }
, c$ f" ~' M# r9 u6 bfunction f = adn(x) %主动配电网层目标函数
# a% T: {5 Y, _% r, M9 Dglobal A b WT PV j1 PA Psub %全局量3 u$ L: _. b; P
ADN=0;% F2 _! N9 X/ B7 T, R
for j=1:243 z1 K( `, B& `" y' S
ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量2 x3 Z2 P) ]( s6 X
ADN=ADN+x(j+768,1)*b; %(输电网+DG1)部分的电价*电量
# j, H" ?" [/ f M ADN=ADN+x(j+792,1)*b; %(输电网+DG1+DG2)部分的电价*电量
) Y( b* U4 D. m' L ADN=ADN+WT(j,1)*b; %(输电网+DG1+DG2+风电)部分的电价*电量
' A* [. |4 b8 r ADN=ADN+PV(j,1)*b; %(输电网+DG1+DG2+风电+光伏)部分的电价*电量 & S4 ~& U6 w7 V, \* h" a9 s, W
ADN=ADN+j1(1,j)*PA(j,1); %(输电网+DG1+DG2+风电+光伏+聚合站)部分的电价*电量,即配电网运行成本) D M o+ K2 L X5 g- [
end
* A' a' b5 w% Zf=ADN;
) H; L5 ?9 W* }! M
0 N- A9 M" l0 p% L+ w
" e9 o, r/ ]4 g
7 E E4 p. v4 ~8 Nclc" y+ D. K" p8 I: m$ J
clear all5 ~* b4 t0 c. f( t
global A b WT PV PL j1 PA Psub
5 d7 e3 J! ]0 |1 N; FA= [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
7 _0 q$ t! ?/ V6 j: Zb= 400; %可再生分布式电源电价(=上网电价-补贴电价)元/MWh6 m, |" q; ?0 M3 S J2 q+ g8 }/ I
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: Y5 ^# ]7 A9 u# F3 c! g# o
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 + L1 Z2 L- ^3 o$ Z4 v! y: O
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
! q. t' T2 ]6 u- S2 b: z9 GWW=WT+PV; %风电、光伏发电出力之和 MW) {9 I% B4 N4 r# z H
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+ u& ?# ]* |, X8 C8 n5 w
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小时负荷数据 MW4 S: M* c' T' h" f9 M; c6 h" r
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% r9 {$ C* o9 Z$ Z+ e2 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- R7 H8 f) d& q) Q9 `& V. x- k
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小时负荷数据 MW5 \: M0 d9 p4 M7 \0 l
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
x7 S+ n! S6 sPL8=[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
6 @, u/ ~9 U& V; Q; f! tPL9=[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
+ I) W! _2 r9 z2 \# K4 P; ]5 @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
. \. H7 x& W" zPL11=[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" j5 }9 U6 G6 O4 E- j$ T& s8 ^
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 % K3 N* G# L( e% \6 r6 Q
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
* h$ U, ]/ l" F4 a% z9 S! NPL14=[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% i- P: Q( F+ \; ~7 D( Y: e0 D1 ]1 j1 W
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小时负荷数据 MW9 ^: u# A4 T! {# k3 c! ^
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" W# D5 T" N; z# x
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
4 p8 U4 ]5 y8 I7 o; x% uPL18=[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
! W" u$ d* g# j# c# _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
# j* v) W* |* [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
3 o: o# q6 k7 O, X$ Z7 \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% A7 V7 P; y/ G2 f
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
1 Z* `% q+ e ]9 r9 h" 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
5 w2 Z& Z# z0 C, B1 l# G) D) A1 T1 gPL24=[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: h" {- B H) W
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
) g7 `! O: q* s$ e7 L: s; u( TPL26=[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% A- G+ a# d$ x
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- ` {, X/ ~9 _# Y
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
: A# ~1 S0 m" B( f' y( XPL29=[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
/ F1 V# n* F) u* T, s+ P+ mPL30=[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- J0 R4 O+ {0 l: u0 N% Y
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
) Z2 [+ r! A" ]" J1 B# QPL32=[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
5 b, w% B, l/ ?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; s0 A/ J% X ~% A8 C3 D' \
PL=[PL2 L3 L4 L5;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
% h( h$ ]* `+ K1 b9 j+ zPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW(32条支路32*24=768)
, }* r+ Q$ h+ ]( @5 O0 V9 l$ N2 ~+ j CPDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW, S# w/ V3 t$ C$ x
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW. P/ i' E3 R1 i2 t
Psub=zeros(24,1); %初始化Psub矩阵
& ~8 ]2 i) N6 X7 r$ BPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW0 J& O% |+ o: g8 I: Q2 S% _2 h+ P
SA=zeros(24,1); %初始化聚合站电池容量 MWh* |, Q q% E' Y. ?. I1 e1 M
x0=zeros(864,1); %变量初始化 . K! n! u; L0 r- G4 F
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
2 {2 H9 ^, S& |" N C: g( uAeq=zeros(793,864); %等式约束系数矩阵初始化
/ |/ U1 b' n. b, E! ^2 Nbeq=zeros(793,1); %等式约束结果矩阵初始化
5 k$ _. ?' j, [3 [( C% ^for i1=1:793
: D7 b; n8 u' A5 F, O) C/ q%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%- }: |& p2 X5 B; L' D8 c1 Z! {
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点:Pij=PL) d2 _1 D0 i0 [, w* m% U! c
Aeq(i1,i1)=1;
0 U/ O- L: R' T beq(i1,1)=PL(i1,1);
4 H( N/ n j% _! s2 K end* e6 ~$ v# w N' _- U! J V
if i1>=145&&i1<=168 %节点8是带PV的节点:Pij+PV=PL即Pij=PL-PV
$ a; ?) h: p& K; J Aeq(i1,i1)=1;
+ g5 `6 j- I, q, F' Z h/ Z beq(i1,1)=PL(i1,1)-PV(i1-144,1);
V, s( a; P$ e5 F" D4 G: t end
/ K; g' O0 r/ k: }- G; l if i1>=169&&i1<=312 %节点9-节点14为非特殊节点:Pij=PL2 B& Z d4 j H) R6 D) m4 u
Aeq(i1,i1)=1;& z; v! e" h6 }) ]
beq(i1,1)=PL(i1,1);- a) `# l# Y; A: f
end- K7 P0 P: t4 `& F( c5 @. u
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点:Pij+PA=PL
# D7 b% V$ B' x2 T& N4 F Aeq(i1,i1)=1;
0 @$ `8 J- p1 F: M; k: x8 I- m/ h4 U Aeq(i1,i1+504)=1; p9 T& z* j. w* z2 s6 L% O% x
beq(i1,1)=PL(i1,1);
4 a) g3 J: Q7 V* o- Q end! @" Y. A* B( R* }
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点:Pij=PL
6 @* g* n+ g2 x) D- j& Z% a Aeq(i1,i1)=1;# T+ D, f- |+ q- ]: b$ T5 U2 F
beq(i1,1)=PL(i1,1);& x' J2 q' X: k% n/ c% @4 l! q
end7 `3 b$ n5 q) }7 q
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点:Pij+PDG1=PL8 ~6 G3 a' v; C' }
Aeq(i1,i1)=1;/ x9 N5 |& V; r
Aeq(i1,i1+384)=1;7 s7 a) m0 [5 r$ G$ O0 j3 F4 b
beq(i1,1)=PL(i1,1);- j, S, h, c2 P7 g m {) p( ~
end
0 R; z" C2 U6 y! H' Z' C if i1>=409&&i1<=456 %节点19-节点20是非特殊节点:Pij=PL
, u# v7 n: I. v: h& P! r$ [" x Aeq(i1,i1)=1;
$ @/ r: G0 l; h8 G+ f beq(i1,1)=PL(i1,1);
8 {7 b0 g! X( u0 W3 ] end9 g* G+ G- c" F+ o* k- I J
if i1>=457&&i1<=480 %节点21是带WT的特殊节点:Pij=PL-WT8 U' O6 |/ @9 N. i: O) D) \4 H
Aeq(i1,i1)=1;
: f) I, u k+ F beq(i1,1)=PL(i1,1)-WT(i1-456,1);3 i4 E/ G- [6 Q8 J. s/ ]# R5 g
end; N, Q! D6 A4 u0 I) p' Z) b8 J% g
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点:Pij=PL3 t) {: t& _1 I
Aeq(i1,i1)=1;6 y r/ F9 Y- {8 D
beq(i1,1)=PL(i1,1);! @7 m% p/ R3 F4 @
end
0 X+ n1 Z/ @# g$ U! { if i1>=649&&i1<=672 %节点29为带DG2的特殊节点:Pij+PDG2=PL
' V* o, E$ r, O$ P4 X9 s# a Aeq(i1,i1)=1;
! |, A! l! A3 g ?8 s: r Aeq(i1,i1+144)=1;
c% N$ g3 A2 k# n beq(i1,1)=PL(i1,1);7 C/ }3 V* a$ M7 A
end
1 R+ h( s8 b4 I7 o1 I; y0 { if i1>=673&&i1<=768 %节点30-节点33为非特殊节点:Pij=PL
" Y& i; c0 i' G; |9 x. N7 t ^ Aeq(i1,i1)=1;* r9 c: o' x+ u, ^8 H( j
beq(i1,1)=PL(i1,1);
3 e& s3 A* |) @2 l end
8 R+ X7 O2 i6 v4 c0 e, F" o%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%. v* J ] R4 x, G3 C
if i1>=792&&i1<=792
( x1 B, [# O; J# x+ q5 v Aeq(i1,i1+72)=1; %对应电池容量SA
! V* t, r$ Y B" z5 A for i2=1:24
, ]" S" A0 F3 ~* v$ h$ t7 }" b if i2<=i1-768* ?! @3 R# O. P7 L: L/ H$ a
Aeq(i1,i2+816)=1;8 Z8 J; j% H: }4 d$ W; Z2 h
end
) ~# r" @- U( C9 K end
( ~ M! x3 S/ L8 O: ~9 X beq(i1,1)=5.4; %5.4MWh& [! O/ V- A2 Y. }6 G9 b0 ]2 u
end( j8 g" F0 n9 x3 o. p6 T
%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
+ X$ A% @( a2 c2 V& ?8 j if i1>=793&&i1<=793 %对应S24=5.4MWh. }! H, M& N0 @: m5 B, N# L
Aeq(i1,i1+71)=1;8 A) X& r/ i4 A
beq(i1,1)=5.4;1 [1 \/ f- g( g; l. h1 P
end7 R* M) R* x% {3 z3 A: e) ^
end & o$ M9 {+ [! @# o- p
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%3 b' y( Q8 I2 K4 M8 j
A=zeros(48,864); %初始化线性不等式约束系数矩阵
. i: k6 X# |9 F7 u6 y) ?' Kb=zeros(48,1); %初始化线性不等式约束结果矩阵
o* `0 ^0 H$ J; c/ Rfor i3=1:48 8 D) I! D5 [* j
if i3>0&&i3<=24 %约束条件VI<0.5,(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t# {/ {2 D& R$ V
A(i3,i3+816)=1; %对应(PA)t的系数
1 Q* g8 A* a) {3 w" k% f if i3==14 c* q* z& i( W. }2 {
A(i3,840)=-1.5; %1时刻的上一时刻为24时,对应t=1时刻的(PA)t-1的系数; Y1 z; n( K! |
b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵# L1 y: r7 ]5 s7 l# I+ _% P
else
# L! y6 A* F" V& V5 k f A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数+ x! `0 |5 b# ?9 R( x$ C; C
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵
/ p& _( p7 n0 a& I" Q, G- } end& R0 ^* l# U/ X& i# B8 B
end
* J7 G: B- }- p- p/ l if i3>24&&i3<=48 %约束条件VI>-0.5,(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-13 M( w$ o' k. U. L( c
A(i3,i3+792)=-1; %对应(PA)t的系数: ?5 \& i4 x: b7 R" h. f
if i3==25
+ B5 j- w6 W I$ _& _5 r& A A(i3,840)=0.5; %1时刻的上一时刻为24时,对应t=1时刻的(PA)t-1的系数
/ Z4 C5 e0 V' _1 Q b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
) _% ^* X7 ?) t/ D else
7 J" U- o! u; n; Y$ B* S5 Z; N A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数 0 x% k5 Y0 b& T
b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
' E9 S y) Q% U* p1 ] end
0 N, e1 i. J- k. S# }9 p0 |* ~ end
6 V1 c9 R. L9 _; x0 o$ I3 nend3 l- I$ p- ~9 P9 O7 T e1 U3 u
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
6 [4 c8 j8 h3 }+ ?Ib=zeros(864,1); %初始化下限矩阵
2 w m. \: L# S6 Q; BUb=zeros(864,1); %初始化上限矩阵* F% K: z) ~6 r Z! I/ Z
for i4=1:864 %Ib Ub赋值. Q) ]6 ]% Y* H5 C$ n5 w
if i4>0&&i4<=768 %Pij上下限:3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求6 M) y" j/ F" ^; y, C) n( }# t
Ib(i4,1)=0;7 u" ?$ Z) J5 [
Ub(i4,1)=3.73;
% O( f. H; ?% u1 Z6 g- C$ w end7 h- M- A: x3 |' c
if i4>=769&&i4<=792 %DG1功率上下限:3>=PDG1>=0 MW9 D8 t: \3 ? q, p
Ib(i4,1)=0;
0 @" R/ y4 G1 v z Ub(i4,1)=3;8 [3 E% \% R, H9 W1 K
end
: y s- l/ N7 ?! m) W if i4>=793&&i4<=816 %DG2功率上下限:3>=PDG2>=0 MW1 ]$ k4 l' Y M: {* R% E
Ib(i4,1)=0;
( O! Q) x+ Z' p& s% P3 z) {2 e Ub(i4,1)=3;3 ?$ L$ S6 } p6 x6 d
end
+ O, h) I/ l6 e. ~ }1 t if i4>=817&&i4<=840 %聚合站充电功率上下限:2>=PA>=-2 MW
. |: ?- L% k/ Y0 Y Ib(i4,1)=-2;4 b( T: G# z0 L) {6 M q4 e
Ub(i4,1)=2;
0 \- c! h, |3 w end& O; v2 e/ s0 K* @7 U
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh* ?2 q% _. N9 \3 E X: m9 j+ [
Ib(i4,1)=1.08;- E$ C( Z a! h8 i
Ub(i4,1)=5.4;
1 b2 D4 R+ a# r# h8 q8 l. \% o# E end- J" k% H, h, a' O; t2 A+ s
end
; y: }' | M8 n: ?7 I1 [%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%% R! y3 U% r: n
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
0 Y, f& a9 J2 q. z2 L' n ufor i5=1:768; u8 t4 M/ N. |% b5 A5 q
Pij(i5,1)=x(i5,1);
& @+ W' L2 I% @. W6 a2 Qend
% j- [6 T+ J8 Z9 Q/ Ofor i6=1:24+ j9 x% R+ B5 L6 C2 w7 j9 J7 \
PDG1(i6,1)=x(i6+768,1);
4 F6 t) H0 s( W" B* [ PDG2(i6,1)=x(i6+792,1);8 J) \) b _9 k |9 K
PA(i6,1)=x(i6+816,1);
9 h O( s: v, o+ L- ]4 m SA(i6,1)=x(i5+840,1);2 [0 E$ ^! u- x6 M8 G: X# [
end" \0 A( I' M. _2 }
for i7=1:244 ~1 n; J5 _4 s9 o- Y1 c
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);: u* i! q) x/ ?1 B, F3 s1 E
end
1 l/ g/ D: r' t$ a8 T6 k9 Edisp(['运行成本:',num2str(fval1),'元']);
" K$ t' v2 c T, ~1 u7 lsubplot(2,1,1);* L: K# j+ q) e6 N! b+ k
plot(PA,'b-','LineWidth',1.5);
' T# Z7 Q. T O. P( u% `hold on4 X* L0 O+ s/ m. G
grid on- ?3 o7 u- |; x" k" g& {) Z
1 e% j# M0 [2 I5 K# x/ ^- t
2 ~4 a* _- H9 ^ |
|