|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：7 J. P& e( ~) q) R& C
) N2 I( E8 c) {2 d# i- e9 u
function f = adn(x) %主动配电网层目标函数
1 w" W% r& _4 H0 H) w9 w, N. Aglobal A b WT PV j1 PA Psub %全局量) @. U, Z+ [3 _3 Y
ADN=0;: Z) u* K4 @3 u5 N# c
for j=1:24
! b( s9 n8 e* N! D* H% K# @ ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量* v8 K8 h! I; [) P; O Q
ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量$ {$ n8 z/ f# ^% ^. h. d4 o
ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量
1 D8 r) E/ B; V: n- E! W! b ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量
' |0 i' F1 Q, t8 N: x( V: g: h ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 $ y( S0 w: ?# O: H- G
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
, ^, q) T2 D- Zend& @& H/ b" C6 A4 i- s, ^
f=ADN;
* K1 _4 A, e3 Q7 M+ A' }+ T
) |- K, l; y4 O. j: A, G/ u+ m% [" [5 K# K9 X$ |
+ x! z# w3 u# s1 K9 j
clc' @" t6 b) n& J5 `
clear all2 _8 u* v% T/ g: A
global A b WT PV PL j1 PA Psub0 {' |/ u5 e$ x# n+ ~6 U/ P. \
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
0 v ?' N F. l# Q- s! X5 bb= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh
1 d% {+ s6 Q+ r. h, lj1=[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; y- i' }, x6 Q* }
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
9 V% Y+ y9 Z' A7 z* j0 t" mPV = [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 k# n0 _1 y5 L3 W: o* d1 P, e
WW=WT+PV; %风电、光伏发电出力之和 MW
! q+ r; ]% S) g# {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
' g2 C3 E$ l& g) [+ xPL3=[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
+ s) t3 p0 Q9 A+ u3 l3 w, mPL4=[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
/ O9 Y: A2 F9 g9 gPL5=[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- D6 X1 a0 O" J. B* 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
: n1 @6 L, `; O, xPL7=[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
1 p# R/ n/ e% e- ]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) o# F' b% P% J# Z/ Z
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
+ ?5 w" r6 G3 Q6 Y! Z. 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 8 l. Z* O% w Z" D, e8 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小时负荷数据 MW
- y R# n( n) S5 r( \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 / l& m. z- F2 _7 H
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
( X6 _" w- g3 JPL14=[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
8 h+ G7 o1 J( v S3 M5 XPL15=[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) p0 P' ], Z1 k- Q
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
( [- k- x* Y7 C6 B0 IPL17=[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小时负荷数据 MW6 `+ v8 ^( G, n1 W& _. C
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
5 }% x/ `1 d; A6 QPL19=[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" w, |( X9 K* ~+ c+ H) k. X2 I
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
1 e8 d3 p3 V8 |0 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 Y. X$ s* j a8 R$ ^! o! \
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 l+ t3 t: M1 }* F3 @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( ~/ u. D2 q) l5 D
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
* w3 S# q- [" f& MPL25=[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
" ^* C3 f( D7 B& f4 r' `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% `2 s2 i: j, E
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
A8 E; T6 c' WPL28=[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小时负荷数据 MW5 D/ q7 J7 F( f
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小时负荷数据 MW5 x! L( ]; ^# a+ x! y7 o" 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小时负荷数据 MW3 Y1 {" t! R' R& 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. _. |) F6 _, P+ q2 U0 E. o' F
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
$ ~, I. X, `9 T/ z Y$ tPL33=[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
# s/ @& C( [8 f% a9 e* ?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
6 x6 ~4 o3 @7 F7 I, O: d7 k6 QPij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)2 C0 p6 _; p& v7 y* D
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW7 i6 w6 y. P* k
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW7 q* A u1 `4 ^; @- h9 Z& Q
Psub=zeros(24,1); %初始化Psub矩阵% [% k/ {) u; c
PA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW7 @1 F3 E9 U: I
SA=zeros(24,1); %初始化聚合站电池容量 MWh2 W/ X% g$ j$ D0 q9 o
x0=zeros(864,1); %变量初始化 % g4 F- n+ U# P6 S6 J
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
; }) @2 q$ i- XAeq=zeros(793,864); %等式约束系数矩阵初始化) W7 A. H! ]) G* M' M3 I- `
beq=zeros(793,1); %等式约束结果矩阵初始化
' q& d( ]0 J$ c b! L# P ofor i1=1:793* J& N8 d. _4 Q% @
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%" P/ F. G3 r* @9 O1 M' c: N
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL
. z$ `2 @2 V p( X& v! J/ O8 d Aeq(i1,i1)=1;
/ i# k' t8 l* l) [: s- G beq(i1,1)=PL(i1,1);& N# k6 Y7 u. C( L0 w2 `5 |. L: U
end# M! K) Q: G$ R: \* S0 L$ v- [
if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV
1 j1 ?$ s2 ]+ j' | Aeq(i1,i1)=1;
+ y) E; k( g& l$ v beq(i1,1)=PL(i1,1)-PV(i1-144,1);) O" Q! i$ C( d( Z
end: I! [7 @# J8 ?$ p5 f
if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
1 B6 H( { F5 R7 K- }' |9 t% a Aeq(i1,i1)=1;$ [/ k( C( b: I; v# t$ r' u l- Y
beq(i1,1)=PL(i1,1);) P7 t# R6 K0 \0 t
end# k: o! [& L! J
if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL
# z, p: @1 T1 r0 m# q Aeq(i1,i1)=1;) s9 n3 m, K9 J( X6 X: [3 f
Aeq(i1,i1+504)=1;
* w: C3 d1 y8 N# Q1 I0 ~ beq(i1,1)=PL(i1,1);9 U. r* V7 _9 |) ~
end+ |; y# M* V+ u/ L# P4 O
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL9 u2 g! x9 i5 q
Aeq(i1,i1)=1;/ Q! J, Z' F7 c Z" n' `
beq(i1,1)=PL(i1,1);) U) N) @- G1 r$ ^- i/ p3 D& x
end" c. i F$ I0 G; j- F2 S) T
if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL
, g; l! ^4 c( S! l5 }# r9 }, V Aeq(i1,i1)=1;
# |4 p% F! I# W, B5 z7 S/ x Aeq(i1,i1+384)=1;! k, B: }# v$ P$ J$ F
beq(i1,1)=PL(i1,1);; t+ g" ?3 Y/ P" v9 L
end2 V- W8 r0 o& I7 T/ X- `
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL
5 Z' z9 D' @7 y) G Aeq(i1,i1)=1;
/ U1 `* o/ g+ y6 l6 e4 u beq(i1,1)=PL(i1,1);
- k( L$ V4 `; q. c: a8 S3 Z4 B end) V" X8 G6 s8 P& D8 F+ b; n
if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT
: W v+ G' h3 z& Y2 E0 @, n Aeq(i1,i1)=1;
' X% _' S0 A% X7 z2 F beq(i1,1)=PL(i1,1)-WT(i1-456,1);7 h8 _2 ^7 i1 g3 ~
end, _: J7 h: g* J2 V" M
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL
/ g# t9 H) x4 T; e, W2 k Aeq(i1,i1)=1;
: O: E2 ~: K# h& M8 g( z beq(i1,1)=PL(i1,1);
' Q& r, T0 |. h: `6 ? end" E4 D3 H2 D: Q- f
if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL
7 h( F0 i- R; T; T0 a6 J- A Aeq(i1,i1)=1;' c" Q/ ^) {3 O* d5 o: [
Aeq(i1,i1+144)=1;
& Z* m" M! D: Y, N# J beq(i1,1)=PL(i1,1);
% C4 ~8 Q4 P# }! P1 n/ Y end0 X/ s. C; ^9 z, { g
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL8 Y8 |9 Z5 `& y! [6 H6 L& f
Aeq(i1,i1)=1;
+ |1 u# v1 ]# }5 Q% V beq(i1,1)=PL(i1,1);- s, L1 ^% Q1 I0 ~( O
end) n2 p, ]3 S: D& l5 r0 X$ ^
%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%9 V( p; A6 y, G% Q4 G+ ^" y( F
if i1>=792&&i1<=792
+ ^7 Z. J4 s: f7 f1 [% a/ } Aeq(i1,i1+72)=1; %对应电池容量SA
( `' H+ L5 r* i7 R' M for i2=1:24
- \2 m* f9 i, m% r8 s) q) | if i2<=i1-768
- h% ~4 T, Z9 R5 W Aeq(i1,i2+816)=1;; v$ U2 R9 n' H4 V' @4 Y
end
- z) r& T! X. j9 z4 g end7 |5 N4 X& k% k& i& o6 U
beq(i1,1)=5.4; %5.4MWh/ A H" Z2 c2 E- z
end2 W9 s$ E2 p( m+ X& \- U) K
%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
! C& k9 p, J0 V( [9 U8 E if i1>=793&&i1<=793 %对应S24=5.4MWh" b7 v7 H4 `( n# q
Aeq(i1,i1+71)=1;& k' |4 A: P7 ?1 g i H4 P
beq(i1,1)=5.4;% U6 L2 u6 h+ r/ K: x
end
5 n8 ^- u$ X d7 C2 Y0 Xend
. D9 d# S& L, q) ^$ Q%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
7 a. `" W3 j6 E @6 [A=zeros(48,864); %初始化线性不等式约束系数矩阵
" h) {1 C! L, `, K+ a/ F8 T1 q4 rb=zeros(48,1); %初始化线性不等式约束结果矩阵
. E" O2 A U1 C j' O) C# mfor i3=1:48
# f, {5 u9 r' E% j' \) v: G if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t I- n' J7 ], H1 @$ K
A(i3,i3+816)=1; %对应(PA)t的系数. S/ ]) `4 L; i7 h% v v
if i3==1
6 F+ @, w1 q/ x5 N$ }$ u) p A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
/ u7 {, L) w# T, a$ n$ o. G b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵
/ k6 E. _7 K# v/ [7 k" D/ R1 G8 T7 \ else
% d" B# w5 `' _# }3 S& T; ?) n A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数2 l! `! n m% R9 z4 c
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵6 e, J: a5 J* {
end
* J7 L. s c2 v+ f$ C5 [$ N end
' v0 q3 f( U) m# u: ?- w$ ^ if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
, ]- t. @/ D( b' y3 [5 m A(i3,i3+792)=-1; %对应(PA)t的系数
q! ]6 F9 N% K% w( D9 W if i3==254 O5 s/ B& N- ]
A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数, ]# U V: q" K3 R: b7 l
b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵
) Y% P( e2 u: R4 h0 ] else/ w7 Q% \" ?$ d
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
v. T+ n2 {% d# u' ~ b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵# r; X2 f0 H k
end
: F2 T" @. L8 x end
% k* }, {+ q: j) G, ^end# Q1 d, {& R! u5 }: M: `
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%, L& q! w% x( x
Ib=zeros(864,1); %初始化下限矩阵/ J- q3 ?2 V" F/ @6 k* ^& v
Ub=zeros(864,1); %初始化上限矩阵' b. I: g0 y h! K
for i4=1:864 %Ib Ub赋值
( g! a) r6 @! u9 I" D; w if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求# {3 j* ~, o% }1 d
Ib(i4,1)=0;; g5 x: r `0 G' D3 y
Ub(i4,1)=3.73;, }9 t0 t# x( c' B
end
5 j: M# b0 }4 `# x, w3 ~ if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW/ `& ]; X8 o/ Z
Ib(i4,1)=0;! M, C8 t s( T: s0 x' p: P
Ub(i4,1)=3;% O, y4 _# [) k0 g4 u
end4 X: J; ]; `. R5 ~# H) y
if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW
4 c+ Y0 O3 }* n( X9 g# Z5 Q) n Ib(i4,1)=0;
/ a) ?3 s8 |3 w0 t& f b Ub(i4,1)=3;, ?2 r8 B9 G/ k6 C6 l& v
end, _( e) J/ i" e9 Z
if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW
, U9 C' \( ?( H1 w3 X Ib(i4,1)=-2;3 d+ g- D! `' J, U8 J; \
Ub(i4,1)=2;; N" t. j; u! {$ e" q8 @) r" Z
end
( P" T$ J4 P9 ^: Q2 J4 V if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh; p% o; \( M: p5 q
Ib(i4,1)=1.08;/ }+ o- N8 ~$ E+ |: f( d7 q4 b0 A
Ub(i4,1)=5.4;% k" A4 v8 k$ ^. H) R! p
end
* X* L' s+ B" Dend
. |( L+ v' m9 J7 J% P# Y%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
5 j0 X0 }& d4 j$ ~) u4 W[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);: N& @' k1 n z4 H% S* U( _6 h
for i5=1:768. D& q+ q6 V. [* h: ?
Pij(i5,1)=x(i5,1);
! x6 s2 R$ i, f7 ^9 \end
$ p; J& y; U% y8 g) ^1 Nfor i6=1:24% M1 ?- U8 B% O# A
PDG1(i6,1)=x(i6+768,1);
0 Y1 c2 R a1 ^2 @4 W$ F @ PDG2(i6,1)=x(i6+792,1);2 A/ J; Q G/ A/ @) z
PA(i6,1)=x(i6+816,1);
% }4 U4 t2 v: ?# ^" M SA(i6,1)=x(i5+840,1);" Z+ D: d; l8 h7 ^4 a# F) z( }
end/ _! J ?5 d0 a ?. }. d
for i7=1:24
6 R t9 k* m3 k- r2 O5 h6 V 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);
' v% Q4 P+ C. o9 ?. Q7 I1 R+ J5 ^end
2 c8 u: c# V* M9 I8 b. Udisp(['运行成本：',num2str(fval1),'元']); 3 }$ S8 p; g. V/ `# @7 F
subplot(2,1,1);
/ {7 g/ D6 N1 e* Iplot(PA,'b-','LineWidth',1.5);
5 \6 }! z/ k9 Z9 V6 S: P& \! n0 b' uhold on
' ?# ^0 U* J( ?: Q1 }* [grid on
) q' E* W1 z( w( {9 h5 m: ?/ T
4 O% k) ~. x6 |5 v; S
, @" [. t+ S+ N% s: e/ S) I
|
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！