|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：
- O! U; ?1 ^4 F# w3 B0 g
' X5 g9 {+ ^3 }1 g) p/ ofunction f = adn(x) %主动配电网层目标函数' C+ a0 q, |2 [: s3 s- e+ }4 m& B g; n
global A b WT PV j1 PA Psub %全局量
' U/ A. S5 u( s$ Q5 ~6 mADN=0;
( L% k! Y o# i+ l( nfor j=1:24
* e- z7 m7 T: k1 e ADN=ADN+Psub(j,1)*A(j,1); %输电网部分的电价*电量1 n+ d; X+ q, O; E9 o8 K
ADN=ADN+x(j+768,1)*b; %（输电网+DG1）部分的电价*电量
9 J5 K+ C; u) z' @/ h ADN=ADN+x(j+792,1)*b; %（输电网+DG1+DG2）部分的电价*电量/ f2 |/ k+ s$ e0 K3 W
ADN=ADN+WT(j,1)*b; %（输电网+DG1+DG2+风电）部分的电价*电量 Y' x' Q- d3 K( a$ @, w, p
ADN=ADN+PV(j,1)*b; %（输电网+DG1+DG2+风电+光伏）部分的电价*电量 % u; Y; b* m8 T0 R% ]2 m2 e0 B
ADN=ADN+j1(1,j)*PA(j,1); %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
: o3 Z$ |* V; M; s; O1 w$ Wend
: ?! y+ G% k Wf=ADN;1 s7 ^# |' G( S! G
( u0 f$ j, C3 K/ @! b* M h
; ~2 A( d: a+ M1 L; C3 v' j6 c% J$ v U5 O5 m
clc1 B9 D' C6 f* \; e5 \
clear all3 q: e& a7 g+ d" X* y
global A b WT PV PL j1 PA Psub
8 ~& }, a8 h2 ]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 $ M9 R3 Q" @# E$ E1 y5 n/ |3 R# j6 Q
b= 400; %可再生分布式电源电价（=上网电价-补贴电价）元/MWh3 l J7 w( s$ y7 a& q- 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* ~+ S+ J- h1 {+ @. {# i/ M
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
- n. l; t* ~, E+ 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小时输出功率 MW
6 c |$ M0 j, c1 @0 w6 WWW=WT+PV; %风电、光伏发电出力之和 MW9 [& T, [- W: W, q: R
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
. D9 p5 a2 i1 d v9 Z* fPL3=[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
9 @& x' u3 _; N BPL4=[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
3 W4 ~) n5 \, ]3 {: ~8 E4 WPL5=[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- ]5 i2 _/ \. R& g# `
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( ~. K; Y. o0 Z, t: g
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
, K/ h8 e6 r- Q9 nPL8=[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小时负荷数据 MW1 r6 I+ v; h. w: X# Y$ W" q
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
( k' R5 E5 d# K0 j: [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
! T' `% N" H1 }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
. ~3 E3 {8 C w0 R5 A; q' ePL12=[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
3 I! u% ~. S: C+ M* ~0 mPL13=[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' Y* @) \' [0 O, T
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) n$ w( u. L- H [8 _
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
8 c. Y/ B) i6 i) X1 A' _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( d2 @- | F1 l: b3 l6 C
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
. V% w7 N4 L8 {2 J3 o4 _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. Q7 x; H4 j- A! `
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 h% r' G+ S" {+ o/ p4 w6 iPL20=[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
, f) C/ j* Y! F0 r5 QPL21=[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
9 {. |& R% B+ D ]0 e C/ z4 jPL22=[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
R, K6 X/ x C! l# v8 l, pPL23=[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小时负荷数据 MW9 l p: K+ {' c
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
" @1 _5 t% k! t& b% u/ k' H, k7 hPL25=[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
$ c: O* Q6 X& x7 [+ D! ]+ n# W( d1 KPL26=[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# X! O8 `5 B$ m' o
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% [$ T; Y- o z: D* ^
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小时负荷数据 MW4 M) `3 ^5 h% y( T, X1 @
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" A8 u/ k) A. l9 F8 p7 L6 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小时负荷数据 MW
) \1 G* C, x* @% _8 V- X; ]; {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
+ G; U3 _% |0 Q% j" q- jPL32=[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* g/ q( H8 x1 f- x0 q
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 @8 [, H6 k0 x+ n: h/ t F
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& Z& l' A# `5 R' D" R: Y8 p
Pij=zeros(768,1); %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)9 u1 Z1 D6 ~, o( N2 o
PDG1=zeros(24,1); %初始化分布式电源给DG1电量 MW: C( F j7 A5 X! ?9 V
PDG2=zeros(24,1); %初始化分布式电源DG2电量 MW5 D) C) R& ~/ g& {! [
Psub=zeros(24,1); %初始化Psub矩阵
( v; K$ @: \+ K) c! ^+ E+ E DPA=zeros(24,1); %初始化主动配电网与聚合站交易电量 MW
, e8 Z: f' Q8 w& JSA=zeros(24,1); %初始化聚合站电池容量 MWh0 e( d# @, H) K+ t: P
x0=zeros(864,1); %变量初始化
% W8 i: D4 Q, O2 `& I3 L%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
; f3 N4 o2 Q/ B- C3 BAeq=zeros(793,864); %等式约束系数矩阵初始化! q' Q8 Z! h B. w8 \, j
beq=zeros(793,1); %等式约束结果矩阵初始化# l5 {5 Q. z5 M0 W: m9 P
for i1=1:7936 u8 Q6 Y+ Z4 P5 k) J
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%( S0 V% ^7 H8 X- g5 |/ S
if i1>=1&&i1<=144 %节点2-节点7这些非特殊节点：Pij=PL! b' Q( Z# P. }6 F* r
Aeq(i1,i1)=1;' r6 o( ]' A9 `1 q6 Y/ `6 y) g
beq(i1,1)=PL(i1,1);: J: q4 \* Y$ `$ P% \) m
end
7 n$ X- M9 z) c. L+ ~( V$ ^8 L" ~/ l; M# } if i1>=145&&i1<=168 %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV- _6 @/ \# z1 T# g6 M
Aeq(i1,i1)=1;+ q6 |# y8 J/ r" a' D: P
beq(i1,1)=PL(i1,1)-PV(i1-144,1);) U3 p( u, W# ~8 |/ o# \( ?, V
end& v& f% k) ?3 t! K- a, C6 t
if i1>=169&&i1<=312 %节点9-节点14为非特殊节点：Pij=PL
5 f. K1 G- d. o8 G& p Aeq(i1,i1)=1;
: g% Z5 U |# N5 D8 P2 v/ z3 J beq(i1,1)=PL(i1,1);$ C0 j* e. v! f
end
& \% s7 |) r }) f/ ?0 A% R if i1>=313&&i1<=336 %节点15是带PEV的特殊节点：Pij+PA=PL0 M7 a' y/ b( F* Z
Aeq(i1,i1)=1;9 ^+ u7 q1 D0 q0 j9 ~8 o3 J% Y
Aeq(i1,i1+504)=1;
7 J0 P! b% g7 u1 V7 @ beq(i1,1)=PL(i1,1);
; f+ V% T( H4 t; @$ ^+ ^" q end o ?0 N) X3 W t
if i1>=337&&i1<=384 %节点16-节点17为非特殊节点：Pij=PL8 d8 D7 \/ V1 i( S0 D7 Y6 l( `! l" c" u
Aeq(i1,i1)=1;
2 w; E. X% x3 }5 J beq(i1,1)=PL(i1,1);
- Z1 k" k2 v1 l/ h end
) J4 U2 Y5 q9 C! s- e* L if i1>=385&&i1<=408 %节点18是带DG1的特殊节点：Pij+PDG1=PL0 G+ q$ E( t3 N/ ]/ a8 y
Aeq(i1,i1)=1;
3 G, [3 B. s5 A Aeq(i1,i1+384)=1;
& c/ E$ j+ S* n' Z* ~: o y4 ` beq(i1,1)=PL(i1,1);) l" R8 M4 a _; @
end) Y/ Z7 @( }+ V; p& I( I
if i1>=409&&i1<=456 %节点19-节点20是非特殊节点：Pij=PL
& d2 l2 t0 N. @: `" x% S. K) S Aeq(i1,i1)=1;
" g3 [. E; D! ^+ \2 z beq(i1,1)=PL(i1,1);6 Q7 T; S$ ~) G; u/ ~/ c
end
/ e* V0 M/ r: f if i1>=457&&i1<=480 %节点21是带WT的特殊节点：Pij=PL-WT7 B, o5 K, H& q& s# d
Aeq(i1,i1)=1;
; v0 z, k) }. ?" V4 y beq(i1,1)=PL(i1,1)-WT(i1-456,1);
' |! v& \/ i' X end; q9 _2 K# [9 q4 F; H6 {0 W& s$ m# r
if i1>=481&&i1<=648 %节点22-节点28为非特殊节点：Pij=PL
5 B9 M) g0 @" M Aeq(i1,i1)=1;
" E" k/ h& S. |. @7 r/ C3 r, t beq(i1,1)=PL(i1,1);8 V) ^# ?) d& X& G) M- R1 T9 D
end
0 n9 G' ]! j. \% |+ p0 g( X& h if i1>=649&&i1<=672 %节点29为带DG2的特殊节点：Pij+PDG2=PL# O* e, Y9 b6 [* ~6 |3 Z5 B
Aeq(i1,i1)=1;
8 o3 h4 r- i: w7 f* M Aeq(i1,i1+144)=1;4 V$ @1 T/ ~6 f6 N
beq(i1,1)=PL(i1,1);
+ ]) L, i; U7 d end0 `; T* `2 z2 F. z
if i1>=673&&i1<=768 %节点30-节点33为非特殊节点：Pij=PL
^0 b6 K5 q" _% q/ h# o Aeq(i1,i1)=1;
0 ?+ K+ B# b1 H8 Y) [$ z beq(i1,1)=PL(i1,1);
& ?& Q! p( G, ?" } end2 E9 v2 q! |9 W$ B- ]
%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
5 B+ J4 ^2 i7 C( }! @# I if i1>=792&&i1<=792
; X7 l/ e# `/ g$ o c1 G Aeq(i1,i1+72)=1; %对应电池容量SA9 m: s6 x+ ] N2 k: P
for i2=1:24
. L. Y! {$ R6 ~5 X) I- C- w& i; e if i2<=i1-768
( T. V+ S/ C" V Aeq(i1,i2+816)=1;+ O1 f ]0 Q6 h$ d% d$ J$ M: r1 v
end X3 _+ ^: }8 P9 |
end
" T! M% L4 q' H9 r5 ^* h* G beq(i1,1)=5.4; %5.4MWh
: H$ l) p* ]( F" l; k) n4 h end
. ?4 V5 q0 }' a+ o" K+ r%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%5 I. y! Z' ~& g. C* A3 h
if i1>=793&&i1<=793 %对应S24=5.4MWh6 ?3 l. C# U0 I
Aeq(i1,i1+71)=1;- S' @) [- d7 e. L; I* B- I
beq(i1,1)=5.4;
! g4 [' q; z+ M, C, ^ end
& n- z' ?" L( q* Z! n2 [ [& }end ( E# U' F: |9 }9 ^7 [3 _! A# K
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
1 [0 T/ b m6 KA=zeros(48,864); %初始化线性不等式约束系数矩阵
9 O+ o/ J' t* m2 L& I/ ^b=zeros(48,1); %初始化线性不等式约束结果矩阵7 z% l$ @6 D: ^# Q% Q: n$ c8 O' ~
for i3=1:48 ' j" H4 U8 T1 [4 ?2 y
if i3>0&&i3<=24 %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
# `: b% G4 F; x9 W% q9 u' k A(i3,i3+816)=1; %对应(PA)t的系数
% `; \; P$ N' ^2 P0 q* G6 l* x" L if i3==1
# k+ j. b6 e. [' o. H5 X' v; U& S A(i3,840)=-1.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
$ p- m- F8 Z. p b(i3,1)=1.5*WW(24,1)-WW(1,1); %对应t=1时刻的不等式约束结果矩阵; F' H& c4 g# c9 V4 g1 N* z3 \0 v e
else& H. K' g/ r% U
A(i3,i3+815)=-1.5; %对应其他时刻(PA)t-1的系数7 q; A* j( J/ Q8 Y) Z6 S
b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1); %对应其他时刻的不等式约束结果矩阵$ a4 R2 h! H1 [/ d2 q9 ]1 H
end' W$ z; W1 V: p" X8 [) r- z7 o
end ! u& i" W+ B& ~, i* `' X
if i3>24&&i3<=48 %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1/ Y( K2 V5 Y7 J* H& c+ j
A(i3,i3+792)=-1; %对应(PA)t的系数7 n8 _/ K5 [3 e# \
if i3==25- q V- X7 m1 S. K% W2 y% F& X
A(i3,840)=0.5; %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数* ]# Q3 k* J' n1 ~) q) Y; _4 |7 Z: q
b(i3,1)=WW(1,1)-0.5*WW(24,1); %对应t=1时刻的不等式约束结果矩阵+ Q: j3 Q, [8 M# Y _
else! \( _8 J& M# J& g, I
A(i3,i3+791)=0.5; %对应其他时刻(PA)t-1的系数
8 {, I6 q# {( s S I b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵3 g6 c a) w) D6 W4 q5 M
end
' L; A$ P7 ?* c" q X/ D! Q {: | end9 j4 g9 M/ E- m1 `& w. n
end
0 k3 z) ?& e' E. ~) L9 N* ^%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
{% ]- L% W1 y0 c5 I( a TIb=zeros(864,1); %初始化下限矩阵6 G7 e( s* C+ D5 H
Ub=zeros(864,1); %初始化上限矩阵9 \1 `. `. T& [, W( X
for i4=1:864 %Ib Ub赋值
* Q. k7 N ?: q9 R8 i. e/ B7 [8 y1 { if i4>0&&i4<=768 %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求8 i) @4 P. ~$ V# R
Ib(i4,1)=0;( ?) H! ^% w' q
Ub(i4,1)=3.73;" u. a* u1 L5 q
end
( V# v' k8 H5 O V if i4>=769&&i4<=792 %DG1功率上下限：3>=PDG1>=0 MW
2 V6 b' ], `6 l& T, ?9 [; ^ Ib(i4,1)=0;
# J) D8 d, B5 n7 O3 x: S Ub(i4,1)=3;
7 q4 A6 f- y+ o% V. }( ` end& n) L y* }6 n
if i4>=793&&i4<=816 %DG2功率上下限：3>=PDG2>=0 MW' D% B) ~, K% P1 w& Q) x0 |; b
Ib(i4,1)=0;6 k: K9 ^, c4 w+ K
Ub(i4,1)=3;% v1 C9 A1 E4 B8 b9 R( B6 o/ V
end1 h3 G i% @8 C9 X6 ]
if i4>=817&&i4<=840 %聚合站充电功率上下限：2>=PA>=-2 MW( r _7 C/ u, x4 i; m' n& }' J
Ib(i4,1)=-2;% F5 t9 s }1 t1 F
Ub(i4,1)=2;6 r- ]6 w, B4 U; _5 U
end3 A, ^! J3 M, g4 g4 G. \5 }( J M
if i4>=841&&i4<=864 %电池容量上下限:5.4>=SA>=1.08 MWh
: ?: ?: K+ H F Ib(i4,1)=1.08;
4 z- B8 N% k5 Y Ub(i4,1)=5.4;) v% u6 w& g1 E' F | |0 v
end" U# g: k- k' _7 r$ K" Q1 f% [
end' \$ M' n: P# V6 X8 H$ b- X P
%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
7 D8 g f3 r* B+ [+ U+ ]4 ^[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
$ t/ x$ C5 z! z/ Q% O4 |for i5=1:768" K( u) P/ m3 s4 H5 g# @* v
Pij(i5,1)=x(i5,1);$ B2 F! ?2 E$ s
end
; p. c: p' w9 D' ]for i6=1:241 ~1 n7 L, W5 r/ m/ u7 m
PDG1(i6,1)=x(i6+768,1);1 Q( W7 C) ^5 }+ U, ~4 ] R( v
PDG2(i6,1)=x(i6+792,1);
w. A% L1 M$ y* m( g PA(i6,1)=x(i6+816,1);# w7 t) z0 s; R9 o" j( l8 L
SA(i6,1)=x(i5+840,1); I* k+ o8 z+ B% M5 s* Y
end6 z7 X$ n& o5 M$ [$ p: b( Y
for i7=1:248 X+ U, Z" y u9 H" e0 N9 O! r
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! n7 V: H# n& U
end
" ?. F. d3 e( Y. V* Rdisp(['运行成本：',num2str(fval1),'元']);
e3 o8 |: y# K0 ^' F% G: ?subplot(2,1,1);
0 v7 Q& ^, c0 X! Lplot(PA,'b-','LineWidth',1.5);
1 h% Y7 s+ l; G& B* Khold on" t& b- a8 O" Q6 @
grid on
0 L# B$ F9 A) Y" @/ p6 ]
$ j4 k% ~/ D3 J1 K, W! E
2 h3 ]1 w8 q0 v2 [1 ?4 Y
|
|
/1 
发表于 2021-3-2 15:30
收藏
淘帖
支持！
反对！