|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
我在用fmincon优化求最小成本,但是出现:User supplied objective function must return a scalar value.求助大神!代码如下:
+ _. V2 X) C5 R2 z) Y( f" j d
0 Z" Q9 `. c: Y* F+ J/ Tfunction f = adn(x) %主动配电网层目标函数
2 r8 b9 J F5 U; S4 }global A b WT PV j1 PA Psub %全局量
, M8 n$ J3 y8 ZADN=0;
# n( z6 t: `7 `6 }3 j* `for j=1:24
. X9 I5 o, K- x5 j/ T) n ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量" _5 q6 _; A2 Q& x# R' ]; D
ADN=ADN+x(j+768,1)*b; %(输电网+DG1)部分的电价*电量
) A% D( T6 y3 o ADN=ADN+x(j+792,1)*b; %(输电网+DG1+DG2)部分的电价*电量+ x* @: y1 n* |' r7 Q/ `) P
ADN=ADN+WT(j,1)*b; %(输电网+DG1+DG2+风电)部分的电价*电量
5 q) X5 I$ \ k1 r9 e! y6 Z0 m ADN=ADN+PV(j,1)*b; %(输电网+DG1+DG2+风电+光伏)部分的电价*电量 3 c$ J; c M' C, h: i4 H* Y
ADN=ADN+j1(1,j)*PA(j,1); %(输电网+DG1+DG2+风电+光伏+聚合站)部分的电价*电量,即配电网运行成本! [4 A8 B1 a8 }
end" j# S# M' ?! R; G
f=ADN;
7 T. I% L- }8 Z% O0 S( d+ y, @& C/ N0 }# T& s5 o. e6 K7 W
$ S* F/ n( H: S5 b* F
1 d3 L) d# T' V2 m7 }clc$ {; U7 ?1 }. f% o' F! {7 U
clear all' c3 \2 |( @) K' w2 x0 I- G
global A b WT PV PL j1 PA Psub6 R9 y- L' @3 y! I# ]4 F
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 ! d9 N% o% ^' |8 c! H3 {3 {
b= 400; %可再生分布式电源电价(=上网电价-补贴电价)元/MWh
+ D' _) [- q& q9 W/ Ej1=[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! x) A v* D. M9 [
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
' {5 c0 X6 ^. J" K8 cPV = [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, p1 z" j% ^% [ T* q j
WW=WT+PV; %风电、光伏发电出力之和 MW
/ H! ~; A0 @) P2 A0 z$ U. rPL2=[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
4 X2 x; w" l- [ }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
! o2 x2 e" ~- w& qPL4=[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+ D; v6 X7 i" M3 w: n
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 ]) ?0 Z3 P: j6 H0 Y9 z
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
o5 t' g, v1 d& X2 sPL7=[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 N$ a* u& X0 E, R/ `* ~! d
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+ Q: D3 y/ H+ F0 ?1 t! h
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! V/ @/ s" y4 L6 q3 X! t/ D' ?: v
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 & M5 Y% j$ D, {" C7 n8 k/ B6 Y
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* N( ~( }0 e. O5 c2 H. S4 `
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
7 a' }: o6 {" x" HPL13=[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" [# D, z. c7 J/ N! N, O: k7 j3 X
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
, r% V! ~$ i5 a W9 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
- g2 b% ]; h& ]; Y9 T T. N; J( Y$ MPL16=[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
6 v5 [8 y x! o* b$ @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
* a% `9 F1 j ` JPL18=[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. p8 s1 L# `2 [! o$ S, W% [' Z
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 p8 c, R0 k! [4 o9 b$ d- v0 ?, 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& y9 F' H0 L5 C3 ]# l7 d# f
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小时负荷数据 MW4 @9 N+ i: R: p
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 ; v# A) u3 b$ Q# c" ^
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+ ^# [- L, w& y5 G# P3 I! p
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+ j% x9 ~1 J H, @; M
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( I: t6 Y: l) p' ~7 l& y1 ~, x! X
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
7 b5 k. ?* C5 ^8 m3 IPL27=[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
0 m3 e, F, q4 V ePL28=[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
8 b1 ^% h8 N& ?: UPL29=[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% y: J7 _- _5 ?6 |- T
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
% m) @* z3 x7 R" T; t0 h( q0 OPL31=[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% B8 K: U+ m# _" C' J
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小时负荷数据 MW9 s; _8 E; A7 U' n
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小时负荷数据 MW4 T9 R1 o1 o! J: |/ Y4 O7 E
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* T# R" E9 G5 H" D
Pij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW(32条支路32*24=768)
$ h4 T' }# B. T. k3 C) Q# f N6 {PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW
; A( r) `* p5 q2 m& U3 l7 f1 D! QPDG2=zeros(24,1); %初始化分布式电源DG2电量 MW
5 J7 `% b' w A c* \Psub=zeros(24,1); %初始化Psub矩阵( @( m. K1 R3 g6 d+ w8 x" \
PA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
, S. O7 I* J' x3 SSA=zeros(24,1); %初始化聚合站电池容量 MWh
+ ]) e% } B- g% s( q) u, j5 ox0=zeros(864,1); %变量初始化
- w4 o5 A5 f- l' m7 m5 x, O Q%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
. {% n' c: M% n0 |, }Aeq=zeros(793,864); %等式约束系数矩阵初始化# \7 d2 Y6 r- k+ f9 w7 {
beq=zeros(793,1); %等式约束结果矩阵初始化
; a( T5 a6 ?% P8 B8 J& ~7 `$ {for i1=1:793
5 i8 Q" R5 ~8 i1 U! R( f; {4 s%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%9 d' g4 O0 g4 ~: r7 h
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点:Pij=PL
9 ]) N7 i+ w# T# D) J6 Q3 ] Aeq(i1,i1)=1;& F, ?/ M. U5 t) ^8 K* u: V* x
beq(i1,1)=PL(i1,1);' @3 F! u6 s3 v
end( b7 q6 Z* I; h. A2 i, J1 c6 r/ p
if i1>=145&&i1<=168 %节点8是带PV的节点:Pij+PV=PL即Pij=PL-PV y# n, {' W O+ M+ w/ J n+ x! P
Aeq(i1,i1)=1;; I6 {, l( [( `9 `" z/ m3 c
beq(i1,1)=PL(i1,1)-PV(i1-144,1);- w( {' r2 h+ O/ u. i+ I$ A/ \7 s
end
7 ]- b8 P+ y+ v3 {$ e' { if i1>=169&&i1<=312 %节点9-节点14为非特殊节点:Pij=PL0 \6 A, |6 ~! Z: b# Z6 t
Aeq(i1,i1)=1;
# [0 W5 |" v/ o$ w, p8 \ beq(i1,1)=PL(i1,1);
3 b- `* P9 V: \/ G- i5 [9 Q. C' { end- y' n) z, K$ H4 n5 A
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点:Pij+PA=PL
* i2 e! w h) A+ G6 u0 u: l Aeq(i1,i1)=1;
, F2 P" q2 H7 c; F+ y0 n Aeq(i1,i1+504)=1;. n' G! l F& S3 u) r+ T0 w; k M
beq(i1,1)=PL(i1,1);/ ?% `" z, D* z+ }* F
end( t: G- o6 q" P
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点:Pij=PL
) h U2 ^" g0 M C# Z; ^* g Aeq(i1,i1)=1;
/ u$ a1 z$ c) U& a$ ~: R beq(i1,1)=PL(i1,1);5 \2 f( b5 V" _9 E+ d
end
; r' | j' j% t if i1>=385&&i1<=408 %节点18是带DG1的特殊节点:Pij+PDG1=PL" n- h- _- b" D
Aeq(i1,i1)=1;
/ W) W5 A5 d/ y( k$ |3 `; H Aeq(i1,i1+384)=1;
R2 {7 N: M% K+ `2 c' Y beq(i1,1)=PL(i1,1);/ e: a7 c8 p/ m4 G: I4 A4 y8 m( k
end
4 L6 Q0 d5 V! D) i& @ if i1>=409&&i1<=456 %节点19-节点20是非特殊节点:Pij=PL
. o! X; S; f( z/ k3 C, v Aeq(i1,i1)=1;, `2 h8 J+ A2 a- M
beq(i1,1)=PL(i1,1);
1 K) R8 D6 l7 X% C) F6 Y end# T, D, ^: {3 X* F$ c8 v5 V
if i1>=457&&i1<=480 %节点21是带WT的特殊节点:Pij=PL-WT7 K$ G' C6 |* f) D4 N" H$ c
Aeq(i1,i1)=1;( P: m' i8 U2 o* Y8 ?0 g
beq(i1,1)=PL(i1,1)-WT(i1-456,1);7 u; I' B k6 P. H2 M, \% L3 i
end
4 r. p( p( C* h% k if i1>=481&&i1<=648 %节点22-节点28为非特殊节点:Pij=PL1 L9 k6 v: u7 f! I6 Y
Aeq(i1,i1)=1;
- f% u @- _" ]9 H* J. q0 t9 e beq(i1,1)=PL(i1,1);+ c4 Y( F4 S* g( A
end* G- ?8 s3 k6 M) G5 d- ^7 Z
if i1>=649&&i1<=672 %节点29为带DG2的特殊节点:Pij+PDG2=PL8 B1 \& {6 X/ s9 G
Aeq(i1,i1)=1;' }% } u2 R6 v6 J. a. p( H
Aeq(i1,i1+144)=1;
" m0 D& V7 t) N1 d4 r6 {2 ] beq(i1,1)=PL(i1,1);; [% A. }6 t% s: Q! [' ^
end
- J Q- z$ x$ H if i1>=673&&i1<=768 %节点30-节点33为非特殊节点:Pij=PL
* n% `# b0 S# \1 A2 m Aeq(i1,i1)=1;% x Z0 G5 n; u; V8 _
beq(i1,1)=PL(i1,1);
n4 F0 w. ]' ~) F, u end0 B. g+ W% T) _3 [5 d0 C/ U8 l
%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
4 U* L$ Y4 i5 h- o8 @ if i1>=792&&i1<=792
% |) q7 U( q, C# x7 k+ n Aeq(i1,i1+72)=1; %对应电池容量SA& [9 j$ ?5 Z4 }6 _3 p
for i2=1:24+ J0 S1 a9 [0 z, u/ j
if i2<=i1-768, ]" ~& f- e, B
Aeq(i1,i2+816)=1;8 }. f3 P5 r2 @# e" M3 O, s
end
# }5 q* ?" K! n0 B/ B end6 d6 k! o! f: B5 U% W
beq(i1,1)=5.4; %5.4MWh
/ {. e3 q5 d7 @. g end
( B3 L3 U) F6 d9 E%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%% J( o/ i- g& w9 o+ J! u
if i1>=793&&i1<=793 %对应S24=5.4MWh
; Z- e1 E7 a% t6 h7 S Aeq(i1,i1+71)=1;9 ?. I& o0 d& O" k8 {0 \
beq(i1,1)=5.4;
# }/ f9 M- I( w: ]8 F end8 |7 o, a5 F' D9 {) k, | M6 C2 j
end
1 G6 y" P9 M6 b. V: r9 \% |+ d: w1 x%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
* [+ P; \' y) ]. N" RA=zeros(48,864); %初始化线性不等式约束系数矩阵& M* ~6 V9 n+ M9 q8 `: X" w
b=zeros(48,1); %初始化线性不等式约束结果矩阵* ?2 J, U5 {% D+ s, F
for i3=1:48
* M2 D. ^: s2 ~% Y if i3>0&&i3<=24 %约束条件VI<0.5,(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
4 {" T0 O& V5 y8 e A(i3,i3+816)=1; %对应(PA)t的系数
, ~7 ?$ C& b' A if i3==1
9 C' W' X7 Z1 U) o/ M A(i3,840)=-1.5; %1时刻的上一时刻为24时,对应t=1时刻的(PA)t-1的系数
7 q, K Q4 b, }0 N b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵1 l7 t* [/ w2 J. ~$ I- ` P6 n
else
5 d1 @ d& `0 s$ @ A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数
( a" J0 a, e/ C) _* h1 J b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵
# J3 n9 p6 F" `8 c. W0 z end
. T( o# F/ n* r# P8 Y: \ end
8 Q; ?$ C: L& \, c3 u$ E if i3>24&&i3<=48 %约束条件VI>-0.5,(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
8 x0 e1 L$ o1 s; o; [ A(i3,i3+792)=-1; %对应(PA)t的系数! y! J* R) z; m+ F
if i3==25
: E' b M0 ?! A7 v. D A(i3,840)=0.5; %1时刻的上一时刻为24时,对应t=1时刻的(PA)t-1的系数6 C0 H# M" M1 i) p) f3 {6 ]
b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
. I6 a& l9 \% {$ l5 A else+ s2 u; K! h# x; p a" ]* g
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
; B9 g. |3 I( e+ _0 R b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
! ?- w. o* _' Q) `6 R end, P$ q `- B h- {
end, d" \- G: `6 f% Z" `5 M
end
7 R& J/ i7 G! I%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%% q2 x( k/ S/ g+ w2 p4 q# O; j
Ib=zeros(864,1); %初始化下限矩阵# ~% G7 T }& \- b, U
Ub=zeros(864,1); %初始化上限矩阵
0 E+ ]& Q g/ Tfor i4=1:864 %Ib Ub赋值8 `8 e, D2 u$ U0 z' f0 z+ }. X
if i4>0&&i4<=768 %Pij上下限:3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
5 S6 H: e4 f( M$ h/ |- o) @5 o Ib(i4,1)=0;* O/ }) }8 B" p; a/ y% r
Ub(i4,1)=3.73;% Z: X+ F/ Q& I1 k- T+ A4 ]
end
* W7 B" `. {9 v, @8 s' Q' @ if i4>=769&&i4<=792 %DG1功率上下限:3>=PDG1>=0 MW" Z; m0 p7 V6 F8 T t$ t! y
Ib(i4,1)=0;
; Y' f/ B; n3 {4 ? Ub(i4,1)=3;1 Y4 a" H& }/ Y# R( |# V" P
end a- {" }! X/ p( M. @. u4 E
if i4>=793&&i4<=816 %DG2功率上下限:3>=PDG2>=0 MW7 |, Y V) J+ V+ T0 [" b
Ib(i4,1)=0;
( _- Y1 F* x& G$ S/ t! A6 { Ub(i4,1)=3;
/ [! e% f0 D( n5 E( O1 L2 x( [ end
; K# ?* n q( F$ J: K1 X7 g" u if i4>=817&&i4<=840 %聚合站充电功率上下限:2>=PA>=-2 MW
. B( n* l) }9 K, h4 c2 Z Ib(i4,1)=-2;
. R$ z5 S3 j! p) f9 M9 }$ T) J' q) R Ub(i4,1)=2;
+ j5 C$ I- p) A* K end% q6 @/ }9 M& S1 A
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh! }6 [1 h/ ]) Q) [
Ib(i4,1)=1.08;
$ R7 m. C# W4 n* w' W% v% D# ] Ub(i4,1)=5.4;
: z+ A4 M( S! f& f% i r1 \: b$ y end
4 t8 Y' z( B0 Y; d3 u; eend
0 a# a2 k; T$ _" c; B%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
# x6 u; m% Z0 h4 E' y' d4 D[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);6 c* a- Y9 ^8 w7 X7 _ }
for i5=1:768
* K7 e( `& A! g. J Pij(i5,1)=x(i5,1);
- Y" O, ~: J/ o8 U1 O. {end
, f e R! B' i' C2 Qfor i6=1:24
! P) s+ p9 d9 f" ~; @ PDG1(i6,1)=x(i6+768,1);
8 V! f9 q% b% v S3 T; W) m PDG2(i6,1)=x(i6+792,1);
6 |& `+ y' b' _7 R( [ PA(i6,1)=x(i6+816,1);
5 o* D. d2 N7 t3 C# Y3 g: Q3 F SA(i6,1)=x(i5+840,1);' \* L" U' ?" S) j
end
4 \, p; }9 U5 Z( X* w- B/ Q$ ]8 M) Xfor i7=1:24) W& o, z: b; J1 }1 ^4 p7 H
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);4 R2 L. F( K/ K+ w
end
o0 i* V7 s2 Mdisp(['运行成本:',num2str(fval1),'元']); $ g( @3 d0 t! h; ]% R
subplot(2,1,1);
; E5 ?) W7 Y2 N3 r5 _/ ? ]plot(PA,'b-','LineWidth',1.5);
/ m5 o; l A( H3 C) a. ]3 e/ j2 khold on
) G; i' r' C0 f, |9 Y, O) Zgrid on3 p- s+ _* u1 d* r+ c0 m
# _' ?. L5 W+ {; ^, p% k% K# u
! K5 B' r4 [( q. L" P( j) \2 \3 y |
|