|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：" s) E) T, d8 @# [+ C
' b) S2 U" ?1 Dfunction f = adn(x) %主动配电网层目标函数
- O, y/ X; q2 ~# Y- mglobal A b WT PV j1 PA Psub %全局量
* x8 [5 p7 S: C+ K; Y1 dADN=0;/ ?8 k& W& I; P/ `; ~9 a
for j=1:24
- y m* \4 z& m6 f3 ?: U- E ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量
% k% J' j& O4 J: a$ W ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量
?$ A; x2 k9 W1 N( s R- C2 V. X ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量
: C, b7 f: B$ Q: p ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量 ! R: }* f6 l9 t2 v( K
ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量
- I3 M$ `. [. k9 i# s, R' n9 S ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
2 A. h0 m/ k) T; Y) `# f- iend: L2 z u+ H4 @7 y7 A7 Y) q: Z
f=ADN;
; ^$ _0 V: ~+ p: }: p
) B( P, U1 }- g( |, o6 \# V5 h5 h O- z
r- u$ i0 _9 _# J+ G
clc
: w# f' \5 X- P! B+ A$ T1 Cclear all6 Q9 b, X$ b' `& [( f; O1 ~7 \$ ]
global A b WT PV PL j1 PA Psub
) { r7 S8 [8 M. JA= [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
2 s2 @6 D5 V( B7 m. X6 [b= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh) H. t! [5 Y% r( ~9 t
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
, W5 T- ^; m# ^; l; B* ^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 ! d4 z0 s+ x8 ?1 C. T7 s
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小时输出功率 MW0 H; M9 q2 L: J' q0 n1 X. |$ ]' I
WW=WT+PV; %风电、光伏发电出力之和 MW" _6 Z( _ H4 x0 ]3 h6 V" P/ b/ U5 q& [
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* T2 E8 Y& {- ^
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
: O9 L$ U' I* N9 ~! C3 q. d+ q oPL4=[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
+ l" y' d3 i( c2 e5 R" APL5=[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( t& q; p6 c) 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 P9 M4 o: N6 Q, TPL7=[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% b4 u3 d" o* @5 ~, \( z
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 V$ z0 j1 w* R, R3 a3 d
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
* ?0 B+ [: k, {7 s( ^/ VPL10=[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
6 }. v- t3 Z A7 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
+ B1 n6 j$ X9 \( Y: j% N8 {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
# q. V8 D8 X8 ]! R, `7 x. }! [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
# M2 Q t- O# c: |% F2 X+ w. RPL14=[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
9 U0 k( p7 r4 E& y3 nPL15=[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
& e. C* J2 y E" [3 ?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
9 L5 F) E' `5 UPL17=[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$ a- d- A- R/ X% a/ K
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) q. K# V1 p6 D3 ?7 r
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小时负荷数据 MW0 P4 J! @3 H# S- 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小时负荷数据 MW
; X6 [' }+ ~0 U* z0 xPL21=[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% E+ G8 d1 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 ) B3 I& z. u" x2 B
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- A8 Q! F0 @* w, F5 B/ y6 @, X
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
e! S' y: j$ L' TPL25=[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小时负荷数据 MW5 r5 h, d( _3 v( \8 K5 G
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
4 G( p9 \8 z q+ A( x/ uPL27=[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
. v( B3 Y$ F# a/ b& n* XPL28=[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: \0 D2 ]9 n8 o- Q# M
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
" v3 @( h8 [5 w+ V& ]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
# J5 _1 W) M. \! T7 p Y9 vPL31=[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. n8 H, K4 G$ d6 S
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
6 v" H+ z2 F& j5 n5 WPL33=[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 C( R; [& i" V; l8 _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
, c |) c7 ]8 G k8 sPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)
, h; m2 x9 i7 w7 d) kPDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW
" g- T: \, g: L5 ?, i+ EPDG2=zeros(24,1); %初始化分布式电源DG2电量 MW
- n% k( c- z9 ^' _Psub=zeros(24,1); %初始化Psub矩阵
4 N8 g, r0 X1 lPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
4 u. P m; r2 | W/ j0 W( Z$ cSA=zeros(24,1); %初始化聚合站电池容量 MWh
- c) q- q L. ox0=zeros(864,1); %变量初始化 ; @ p, C0 U0 o$ F
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
5 U! i: g! c0 N' x' Z2 HAeq=zeros(793,864); %等式约束系数矩阵初始化
! a; ]9 \/ U; ?( }beq=zeros(793,1); %等式约束结果矩阵初始化5 [1 |% s* W7 t" {0 D
for i1=1:793$ J* `0 c y! s& W% {- ^) g; d
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%
+ A0 Y4 Z, {/ S# ^+ T if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL
3 P, c- Y6 S/ N, Y& `( X* F Aeq(i1,i1)=1;
+ ^0 R9 B- V2 V0 h beq(i1,1)=PL(i1,1);
2 D( x( t, E/ y, x3 Q) |3 q, u. P( A5 \ end
* E. ^3 Z2 i5 R& z- v+ Q' j' q if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV) S3 L. R7 \, h
Aeq(i1,i1)=1;' \* i3 h/ \" J- [2 C
beq(i1,1)=PL(i1,1)-PV(i1-144,1);& Y1 f, b0 r8 V1 u
end
) A; k* ` ^; `- z if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
+ ^' D, L1 p* Y Aeq(i1,i1)=1;
: g; k k; C4 N0 ^3 B beq(i1,1)=PL(i1,1);
+ T; v( ~- ^: @1 E b end, z6 F: d* |) u) h3 g* a
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL
% Y; k9 n* X, R3 N( p Aeq(i1,i1)=1;
8 K( i% b8 @6 {+ O* ~1 H& L Aeq(i1,i1+504)=1;0 V' q, {7 u r) _$ q1 [
beq(i1,1)=PL(i1,1);
7 u8 ^" m: U7 t( a/ ^, G( ^2 ^) B end7 n/ @" f l# J, \6 o; N9 _
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL
& Q7 w8 j3 N7 f) O2 Q0 W! E9 } Aeq(i1,i1)=1;6 a% n8 J+ M+ v( X4 H
beq(i1,1)=PL(i1,1);5 p. R( d- ~$ m& _+ Z
end
, m. f) t5 a F0 r if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL
( G5 X: K C; D% S4 f. ~+ } Aeq(i1,i1)=1;
7 g1 b; N4 ]& Z; b Aeq(i1,i1+384)=1;
& X5 I; ^ ?1 Q$ Z; I1 r; _ beq(i1,1)=PL(i1,1);- w% ~. ~9 w W1 J* W+ y; Y
end
, {' C0 G' S6 u6 T7 ~ if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL
) x- B6 T; i5 O, N5 R% ?) {4 R Aeq(i1,i1)=1;. v* r- @$ D0 x( h( H2 P3 [- S
beq(i1,1)=PL(i1,1);
9 w8 s$ `7 D3 R" W end
3 p$ p* _& {; F8 N" I% Q/ { if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT) P2 N8 a$ W! y# J0 ^" U* C
Aeq(i1,i1)=1;
8 C. c* o" H7 P5 U( ]$ R/ u beq(i1,1)=PL(i1,1)-WT(i1-456,1);
% F! L3 X4 W/ s7 j0 Q4 Z+ q I end3 X7 m; d1 I( @! Q
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL H" i8 _% j7 H# S
Aeq(i1,i1)=1; t+ X" F$ b5 M
beq(i1,1)=PL(i1,1);
% a6 F! `7 ?0 Y: f; x end
$ e0 {& j1 D F& j) b X( ? if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL4 y0 d" q: ]& ^0 @+ f, \( ]
Aeq(i1,i1)=1;; H- `5 e9 T: \$ j" U
Aeq(i1,i1+144)=1;3 ?+ F: B/ e3 k+ C
beq(i1,1)=PL(i1,1);
. I: E3 C/ w5 [+ E: }/ t' e- p end5 r4 _' n5 Z5 E' ?
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL6 r; ^# {3 X% O; T5 s! C
Aeq(i1,i1)=1;% T/ X) Q6 R/ D9 W; Q t! P
beq(i1,1)=PL(i1,1);9 N: A4 Q9 u# y
end
' ?% `6 W6 a# s/ T%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
! D) t3 {7 L; m8 I) ^. y3 N& w if i1>=792&&i1<=792( C& R' U' ]) A) ]5 c2 K8 G
Aeq(i1,i1+72)=1; %对应电池容量SA8 y6 z: ?% s4 Z! x3 F
for i2=1:24
) g- y1 {$ h$ G8 L, J. Q* ? c if i2<=i1-768
+ |% h- I+ ?$ k$ r* | Aeq(i1,i2+816)=1;+ Z6 I g5 X3 t7 N3 {. J2 c
end
! {# Q$ a3 T$ I end" U% W, |9 f4 Z/ v. D
beq(i1,1)=5.4; %5.4MWh
( l2 w$ B( L( ] k& F# ? end
C1 u) I- R! ?. _9 }2 g& ^%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
9 j S" w2 x7 C+ ] if i1>=793&&i1<=793 %对应S24=5.4MWh. v( j/ ]) u2 I7 d) K
Aeq(i1,i1+71)=1;
0 Z+ X' u8 Y/ d5 J beq(i1,1)=5.4;( O; C" e7 v9 m( W! g3 [8 w7 F
end
/ W8 ^0 H" ?* W8 V0 uend + I: X9 O4 G Q
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
) K$ y0 y. p9 C! \6 `+ @, p3 |) OA=zeros(48,864); %初始化线性不等式约束系数矩阵
0 i: r* T5 G9 U9 ^. A( Ab=zeros(48,1); %初始化线性不等式约束结果矩阵. ^) W- K- ~1 c, O- e7 o* v
for i3=1:48
) _! _- [8 T N" [ if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t, Y7 M2 D( q+ `: i
A(i3,i3+816)=1; %对应(PA)t的系数
! X% h* p" f# L2 L9 _6 J7 N if i3==1
2 q$ L% |- L- d3 [! w( G( A A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
* k# W- ^* ^: g: l9 Q; O9 ]/ E b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵$ [. p: O7 d" q' f7 ]" M5 W/ `
else
! L; G3 c4 v, ?$ P$ Z0 I A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数
( v$ |, l% }6 K; R! M9 M9 E t b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵( h; m4 y5 b( `
end
. ~6 ]8 | M' j end 3 k) A1 c7 h4 w+ _
if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-12 ]3 p( @" u3 }1 d# I
A(i3,i3+792)=-1; %对应(PA)t的系数
( D7 F W1 k6 @# H8 C6 F if i3==25
! ^+ u& v5 e5 u A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数$ i7 C1 f9 }- V* K" X
b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
9 E3 e0 U& a& Q9 v; M else9 a. q5 H* L, X" |
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
1 m9 w& c& C9 X: `: { b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
6 G* y9 g! w" b) E* y1 D end
9 Z- D; ~9 L1 O9 J( u2 z3 x end/ {! k ~% e+ G. W" c) k
end5 `! F, a$ G9 E t( F L, s
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%# m, b! o5 s+ ], c
Ib=zeros(864,1); %初始化下限矩阵 Z5 `2 X; t0 K! \; C& M9 D
Ub=zeros(864,1); %初始化上限矩阵1 l8 q& F- {1 Z) D+ n2 h! I) f# R
for i4=1:864 %Ib Ub赋值& b: }& P/ K% N- ~3 U
if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求: f- V/ H" W' D& h' o; ~" t& Y
Ib(i4,1)=0;0 \/ K- `3 f, |, |. g
Ub(i4,1)=3.73;
$ |3 K; Y8 ~2 X+ o \2 g end
/ m/ d# k' ]$ a& U9 J% z# h if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW
5 p6 [7 H$ R7 N, y2 m5 Q Ib(i4,1)=0;
[; t- {) h# u$ ]" F Ub(i4,1)=3;* U- Q- A6 P; c4 O" n3 r
end
N5 N& m9 ~7 F1 `5 u if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW u2 O: ]: h6 @! n" x6 ]9 u
Ib(i4,1)=0;$ o& E- k! b' {$ a0 S, x5 N) \
Ub(i4,1)=3;% O1 N$ U1 c1 J# p o& n$ B o
end
4 h& X0 Y7 k, h8 Q! p6 E% L: I if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW) y3 ~4 ?! z/ v
Ib(i4,1)=-2;
6 H/ A5 ^ R2 @2 D I7 [1 A Ub(i4,1)=2;4 h: p- s8 M+ Q T b
end' w" j/ U6 T' v" W, M- A \+ U
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh
" L9 r% ~% J7 [/ Z& E Ib(i4,1)=1.08;
- c2 M, k1 m2 S. _- P Ub(i4,1)=5.4;0 B; Q) R) p! d2 y
end
: o% `' k. K+ ~& r* j! x7 r- yend
1 J6 h! P7 u+ N, v( z%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
1 C r5 I) I% t& e[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);$ u; a4 s. ?/ n1 K9 Z
for i5=1:768
b0 K F0 `/ K Pij(i5,1)=x(i5,1);& _; f+ W' v& R' X) a
end
; N5 y$ I+ |/ R( Efor i6=1:24
5 w! L/ H J1 g3 ?" l$ A$ g4 s- ] PDG1(i6,1)=x(i6+768,1);
9 ]- h5 s l; t! f5 z9 I' H PDG2(i6,1)=x(i6+792,1);
5 |; k6 t% v# s9 X* l0 m/ _ PA(i6,1)=x(i6+816,1);
+ @- l3 r, ^, c9 j( g SA(i6,1)=x(i5+840,1);2 C& l1 T( P8 U6 N) R
end
! X5 ?, \# ^5 Z* Q( c% N) ^, \9 Kfor i7=1:24
* w) q8 K, G/ j/ G# J 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);
0 ~6 e( B0 I; |! e- a8 a6 Mend$ L/ B8 a3 G, @3 b: M9 a. j
disp(['运行成本：',num2str(fval1),'元']); 1 I3 z/ i( |' |& _: U2 ~
subplot(2,1,1);
: @ |) d( U; k- Y* fplot(PA,'b-','LineWidth',1.5); W- b! Q' O0 D# w% B/ @. T& F8 K
hold on5 W; n5 |6 Y/ z5 g5 T7 e/ A
grid on
% U) c V+ {/ ]* N. k* M$ d9 b# r$ ^/ s
1 I: H* v- w, S1 F |
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！