|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：
& I) x/ A) |9 b, T% D' A% I0 ?5 }( d& s) X: x& e
function f = adn(x) %主动配电网层目标函数# @9 u6 D1 ?; r+ M
global A b WT PV j1 PA Psub %全局量" z3 G7 R& B6 `7 u; b9 H
ADN=0;' K: F) y8 d8 w7 }% p8 t
for j=1:24
6 x% W! h! U5 s) O- @( J; n8 N+ j ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量
7 {0 s' E) M+ n& M% G) \ ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量
) k3 x/ { Q( x8 j ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量
7 `! g' Z* h: U( M1 l9 C7 O ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量
) G) h+ c9 p& u1 _9 J* H1 ]* t ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 0 x! Z3 h( u0 y5 t* i* I
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本% \8 G. k# G+ H! F
end
$ R' @$ z# I4 z# Q$ I: nf=ADN;& f) S1 ?8 z; d2 x# D! I+ `, p( ^
# q6 D) H% s3 X4 m/ t9 u) G+ O1 y/ Q7 V6 c
* N* H! S1 K9 t( l, @# _. I' R2 f% C: u# k' N; E
clc. H# M& X f9 s6 G& h
clear all
7 M! ?7 D4 a+ Dglobal A b WT PV PL j1 PA Psub
+ H4 H$ ~: _, m T2 D4 AA= [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 4 [' x! K, ]6 |9 k* ^1 v1 d
b= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh" C2 x8 F+ V8 S3 g5 z: \
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
* ?; e9 i+ K8 J2 Z$ qWT = [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
; X# z1 Q/ ]+ d5 ~$ T* @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小时输出功率 MW4 ^) B$ A* `! ]. {6 c- S. s/ i( d
WW=WT+PV; %风电、光伏发电出力之和 MW6 @ y7 N6 O$ r- I I1 C% _1 v
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 W1 y: J: V: K% [: S
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
% w+ p* `9 B' k9 PPL4=[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 U8 c( j+ Z; l* uPL5=[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
' _1 \6 J3 ]. E9 Q" oPL6=[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
/ k, D6 W J/ G M' E0 F. [! c7 aPL7=[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小时负荷数据 MW5 ] K4 ~& ]! B2 Q% d% ?( U
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
* U- \6 J+ W3 b- g* y @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
/ y& W R' C5 o& lPL10=[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 : n, R- ~ s9 S' N- j
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
2 U* g8 @* D7 K) dPL12=[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 + V# z3 k$ k+ N0 `+ g P$ C
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' q5 `8 g0 a) U S: m f+ s
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小时负荷数据 MW& A$ | L! ] A7 x6 u% A
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# K# A8 \4 m I( o( H, ^
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
- ~! m8 ?! W6 j' yPL17=[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
* D! s0 B* R: q% x/ V$ I. `( xPL18=[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小时负荷数据 MW4 k& f& M7 o8 y; F
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 T+ S% _. y7 g% _3 t3 n6 N
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小时负荷数据 MW9 F6 L& \9 A+ A( k6 Q3 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! \( I3 w4 Z5 R0 ~
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 4 K4 ?$ m6 c' R2 R4 S% s- ]
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
5 j9 C& G; H, C- l0 F' i& hPL24=[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小时负荷数据 MW6 U+ x* _7 U: l/ 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
. v. G0 J& m( B( c5 vPL26=[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+ o1 l7 H/ g8 k4 t4 ?
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
& [- Y% k4 e# V" `- F5 e, ^& B. lPL28=[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 p( V# t1 D& ^
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
* [7 s0 y: B3 c) LPL30=[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+ }& H, o' ]! a. F/ a
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
+ l! ?' D9 H# q. I$ I B, k2 mPL32=[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; f1 {3 C) W, l9 i+ n8 J+ D) d
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
3 P. u2 s, Q5 F5 G# B6 i3 n, tPL=[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
2 g0 y3 o3 E: v# B4 S' ZPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768): V& Z0 q2 |0 @, t
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW0 p2 Q( f! Z+ B5 d- Q% W& u
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW
: { V9 B7 t5 e: [9 [Psub=zeros(24,1); %初始化Psub矩阵
' L3 [* p) @: L1 e5 ^- PPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
2 S% X2 R- p k9 z$ }SA=zeros(24,1); %初始化聚合站电池容量 MWh
+ [2 d' z5 e0 K+ U, ux0=zeros(864,1); %变量初始化 E$ w, ~4 c) q: V+ _0 s& C% N
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
( u0 [9 `! `$ Z) h$ cAeq=zeros(793,864); %等式约束系数矩阵初始化; F4 I! M) s5 C2 Q% O* Y: j8 m
beq=zeros(793,1); %等式约束结果矩阵初始化
2 \/ {5 a. H1 W. |1 O& J Lfor i1=1:7933 `6 P0 M5 l% ^5 Z5 F
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%+ x9 Z, n7 S, I
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL6 \+ F3 I5 Z% i
Aeq(i1,i1)=1;- J( G2 G4 x% `% i. c
beq(i1,1)=PL(i1,1);$ i* [. U% u# c$ D# j2 I* v
end
1 s" r$ f* d1 A' c if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV
, A! R0 {2 B: m! X0 G( X Aeq(i1,i1)=1;7 T% ~- Z% X3 O) E% n7 i
beq(i1,1)=PL(i1,1)-PV(i1-144,1);8 j: U- G& J% b5 C e" v
end
" N1 m$ A" w& C" @' ? if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
/ p" K) E+ z7 o f8 K } Aeq(i1,i1)=1;! S7 Q% p% ~) ?/ {$ W
beq(i1,1)=PL(i1,1);
) ]" y$ j& b2 Q a end! s8 Z3 l& [. s9 Z6 O
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL4 p9 r3 N: w# r
Aeq(i1,i1)=1;5 Z" i+ H% A* K
Aeq(i1,i1+504)=1;
6 E1 Q( p/ {. _ @0 U beq(i1,1)=PL(i1,1);
5 t, L( i& y8 {4 T3 u end; T, L& ~! m! F* J8 {
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL
& I B3 ~. G6 y% y+ [( E$ z9 U Aeq(i1,i1)=1;
9 s) S$ n0 I& A* N6 d x' x beq(i1,1)=PL(i1,1);
; b9 W. {) Q7 m/ P' z u4 |1 ^! I end# q; |& N2 ~: `/ i, l/ _: X2 d
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL
5 ^ F/ g& e" P' y( ? Aeq(i1,i1)=1;
$ Y) ^% I, p! R$ o7 Y+ b/ t Aeq(i1,i1+384)=1;
# c, _. W/ K0 K beq(i1,1)=PL(i1,1);* U0 D( \5 [( } x) h
end- i# M& p1 ]) m6 @9 L
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL2 S. h0 P" {2 X3 r8 J" ]
Aeq(i1,i1)=1;! R6 e% A7 r7 u4 _2 T
beq(i1,1)=PL(i1,1);: F0 O9 K5 T3 _# _6 l
end
- Y4 S; {4 \- L4 {3 n if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT
: C. i# ?/ t/ c Aeq(i1,i1)=1;
4 v& W0 N5 Y" X# b' h, H beq(i1,1)=PL(i1,1)-WT(i1-456,1);0 v" h2 Y _+ \( [
end3 n$ @. V% z, ^' X; u" m
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL
1 z' R r+ e4 s# @3 o Aeq(i1,i1)=1;
( n2 p* b- P9 P U beq(i1,1)=PL(i1,1);5 K; M& Q1 m* K4 P* w" J
end
}' j5 Y! N: |* U9 _# \* k if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL
1 M1 n% V' f% }# o0 ` Aeq(i1,i1)=1;
! u! l: v0 ] `! t Aeq(i1,i1+144)=1;3 g: L- z6 j5 I$ Z
beq(i1,1)=PL(i1,1);1 R: z6 H/ R7 ?7 U% V g3 Z: d
end9 W% p7 n) O% W% T
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL- e2 \0 g- U+ g6 {( i
Aeq(i1,i1)=1;' Q$ U0 D2 u, G) Q
beq(i1,1)=PL(i1,1);
: }8 ^- q ?) N1 t2 D! R end* e* [& N+ w1 Z8 M
%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
( _! q$ |0 F: E0 x if i1>=792&&i1<=792
4 C* g. f+ g7 ?! _5 \ Aeq(i1,i1+72)=1; %对应电池容量SA. g* R! `. A) M8 P5 M; J
for i2=1:24
: v+ [8 \6 `4 L8 X- \ } if i2<=i1-768
/ Z) v/ x& I& F" `' G3 k& h Aeq(i1,i2+816)=1;
: J D' S" Y* S end
4 h9 _& {8 E) `* c6 g+ v end% w, K$ Q( ]$ v2 d' x" _: p' N
beq(i1,1)=5.4; %5.4MWh
2 ^ [1 w7 ^. m. W7 o! b end
( S. y2 I, v8 L% e%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
" `3 d( j# [7 A) [. C if i1>=793&&i1<=793 %对应S24=5.4MWh9 S, V3 H( R+ M9 ]# C5 d
Aeq(i1,i1+71)=1;9 M) m: l& ]. \( A- b9 v. G2 s& x3 f! X
beq(i1,1)=5.4;( U" S9 f3 r. n( k9 k' d/ p5 _
end
2 {& ~, j2 w- |# s! Y7 }: h4 ~+ Fend ; f' `$ d, e' S8 A. M
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%: B+ O7 _& G! y2 y8 @
A=zeros(48,864); %初始化线性不等式约束系数矩阵5 ~5 r. J5 S' a5 y {
b=zeros(48,1); %初始化线性不等式约束结果矩阵
0 c K) n: X* o3 l5 z* E2 x) ^for i3=1:48 9 d0 c& F8 E* ^
if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
9 _# o( q" K7 O A(i3,i3+816)=1; %对应(PA)t的系数7 d/ K* v: b. Y
if i3==1
& Q$ r/ J) c( W+ {* H- L- [" Y P A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数! }# l7 c B7 N, Y( a/ S
b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵4 ?+ I: p* P, a4 u$ s
else8 I& D/ ~" c, k+ t/ |9 W
A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数. r; m/ C0 c7 M, C. o8 C9 ?8 P i8 I
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵
! p7 h* D9 G9 |& m end
/ F* W8 p. [; x+ S! z2 g7 C end * d3 o" d6 I' r- ]
if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
' Y' T/ ?/ L0 j A(i3,i3+792)=-1; %对应(PA)t的系数 v: C: U, u3 v
if i3==25% r) w& C+ o( T' i! a4 \
A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
& N2 n, x4 z! [ z" t b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
. T9 D: \ K& E0 O5 t: v else
/ A6 m' l% S3 j$ m% }* D' i A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
' ~6 ^+ v9 F0 h7 O: y1 y5 ` b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵0 _: C( e4 H& |, \
end6 N9 q# A) s- y
end
3 p) H3 y- u1 T- W9 kend
2 O$ I% i; P6 \: u0 F* H%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%9 M% M. W) s1 j! h3 F
Ib=zeros(864,1); %初始化下限矩阵
* K" x3 m8 h) t L1 C3 r% V- |. |Ub=zeros(864,1); %初始化上限矩阵( _; }' T Y# o
for i4=1:864 %Ib Ub赋值
$ ?) w1 V+ f* b+ V if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
8 c2 p0 b1 n2 a( D2 R8 s. m- z5 [ Ib(i4,1)=0;
7 r# G. C! p& R Ub(i4,1)=3.73;
+ `7 }& `+ f) k8 a! F6 s0 \2 @ end- ?' ?" r( ?$ A; h, S) r- g
if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW, B7 {8 ]+ Y) }4 u0 m
Ib(i4,1)=0;* D7 d& t+ a! V
Ub(i4,1)=3;$ j( L: m; d) _( K- k. Z$ S/ v; H
end$ r: M! ]5 z+ G" h7 l
if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW
6 u8 h: @( }$ k' z+ M0 J Ib(i4,1)=0;
! _7 q K3 { L8 Q2 ?7 u+ e. Z, [ Ub(i4,1)=3;# w* C/ Y/ |/ ~1 ~' @7 K( s) f
end
+ b/ Q K" s& s% S& {2 ~" B* L if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW
6 }' t! t% p5 W/ |5 ?4 N0 A* l: F Ib(i4,1)=-2;+ D* X) h% l9 ^% F- f' ?: Q. X
Ub(i4,1)=2;" ~$ Y Y+ s% r5 r5 w& ^2 \
end
0 U9 q! w8 A" X7 h if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh
4 _ x5 R* }/ C* i Ib(i4,1)=1.08;0 u# w$ t+ l- g }
Ub(i4,1)=5.4;% p. ^6 V8 T3 p' v& u
end
0 W: l) p1 d! z# Mend
) I# g3 C. P8 T0 M+ F) Y0 u' l& m%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
6 b, y' s( E* w. p) I4 ]% \ n& w[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);9 i/ F0 |; i+ O J* _2 g8 Y
for i5=1:7685 u' V- s4 J( E2 |8 F/ a
Pij(i5,1)=x(i5,1);
- q& j- C; }9 ~$ |. [* D+ ?& P! c9 [end
5 _+ E5 ?! p3 I8 Sfor i6=1:24
8 k7 _7 v6 ~0 I+ H8 ?6 n PDG1(i6,1)=x(i6+768,1);3 h: E* e1 |' j3 V) |
PDG2(i6,1)=x(i6+792,1);
% M6 T2 k6 ~+ x5 d6 k3 F* C PA(i6,1)=x(i6+816,1);
( d( g0 ]1 i* M0 K% k SA(i6,1)=x(i5+840,1);
4 D% O! t8 [, Cend
% H( t, ~4 v8 G! i Dfor i7=1:24
( u; v5 ?7 s7 l" P" x& f9 y% X; Y 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);; y" W9 O1 V5 f
end0 j8 _! w1 @ v# G( X! }) D
disp(['运行成本：',num2str(fval1),'元']); - \$ e( k' W1 X3 K1 Q
subplot(2,1,1);
) u. U4 _8 o* S$ ]2 F' B+ N% qplot(PA,'b-','LineWidth',1.5);
5 M" Z1 i# ~' h+ H4 a v6 chold on5 j9 B) ]- I2 r/ ]2 i8 n
grid on
& Q0 w6 |0 N" k2 N% J" n# A9 G, u/ \! u$ e0 E$ J
5 d, }; y/ ?) m# ?; s% d |
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！