找回密码
 注册
关于网站域名变更的通知
查看: 600|回复: 3
打印 上一主题 下一主题

fmincon优化

[复制链接]

该用户从未签到

跳转到指定楼层
1#
发表于 2021-3-2 15:30 | 只看该作者 |只看大图 回帖奖励 |倒序浏览 |阅读模式

EDA365欢迎您登录!

您需要 登录 才可以下载或查看,没有帐号?注册

x
我在用fmincon优化求最小成本,但是出现:User supplied objective function must return a scalar value.求助大神!代码如下:2 s% S+ k  r: q/ u  V. }

, c$ f" ~' M# r9 u6 bfunction f = adn(x)                  %主动配电网层目标函数
# a% T: {5 Y, _% r, M9 Dglobal  A b WT PV j1 PA  Psub          %全局量3 u$ L: _. b; P
ADN=0;% F2 _! N9 X/ B7 T, R
for j=1:243 z1 K( `, B& `" y' S
    ADN=ADN+Psub(j,1)*A(j,1);        %输电网部分的电价*电量2 x3 Z2 P) ]( s6 X
    ADN=ADN+x(j+768,1)*b;            %(输电网+DG1)部分的电价*电量
# j, H" ?" [/ f  M    ADN=ADN+x(j+792,1)*b;            %(输电网+DG1+DG2)部分的电价*电量
) Y( b* U4 D. m' L    ADN=ADN+WT(j,1)*b;               %(输电网+DG1+DG2+风电)部分的电价*电量            
' A* [. |4 b8 r    ADN=ADN+PV(j,1)*b;               %(输电网+DG1+DG2+风电+光伏)部分的电价*电量 & S4 ~& U6 w7 V, \* h" a9 s, W
    ADN=ADN+j1(1,j)*PA(j,1);         %(输电网+DG1+DG2+风电+光伏+聚合站)部分的电价*电量,即配电网运行成本) D  M  o+ K2 L  X5 g- [
end
* A' a' b5 w% Zf=ADN;
) H; L5 ?9 W* }! M
0 N- A9 M" l0 p% L+ w
" e9 o, r/ ]4 g
7 E  E4 p. v4 ~8 Nclc" y+ D. K" p8 I: m$ J
clear all5 ~* b4 t0 c. f( t
global A b WT PV  PL j1 PA Psub
5 d7 e3 J! ]0 |1 N; FA= [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                 
7 _0 q$ t! ?/ V6 j: Zb= 400;                                                                                                                            %可再生分布式电源电价(=上网电价-补贴电价)元/MWh6 m, |" q; ?0 M3 S  J2 q+ g8 }/ 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: Y5 ^# ]7 A9 u# F3 c! g# o
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               + L1 Z2 L- ^3 o$ Z4 v! y: O
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
! q. t' T2 ]6 u- S2 b: z9 GWW=WT+PV;   %风电、光伏发电出力之和 MW) {9 I% B4 N4 r# z  H
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+ u& ?# ]* |, X8 C8 n5 w
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小时负荷数据 MW4 S: M* c' T' h" f9 M; c6 h" r
PL4=[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% r9 {$ C* o9 Z$ Z+ e2 v
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- R7 H8 f) d& q) Q9 `& V. x- k
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小时负荷数据 MW5 \: M0 d9 p4 M7 \0 l
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
  x7 S+ n! S6 sPL8=[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
6 @, u/ ~9 U& V; Q; f! tPL9=[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
+ I) W! _2 r9 z2 \# K4 P; ]5 @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
. \. H7 x& W" zPL11=[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" j5 }9 U6 G6 O4 E- j$ T& s8 ^
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  % K3 N* G# L( e% \6 r6 Q
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
* h$ U, ]/ l" F4 a% z9 S! NPL14=[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% i- P: Q( F+ \; ~7 D( Y: e0 D1 ]1 j1 W
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小时负荷数据 MW9 ^: u# A4 T! {# k3 c! ^
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" W# D5 T" N; z# x
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
4 p8 U4 ]5 y8 I7 o; x% uPL18=[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
! W" u$ d* g# j# c# _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
# j* v) 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
3 o: o# q6 k7 O, X$ Z7 \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% A7 V7 P; y/ G2 f
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
1 Z* `% q+ e  ]9 r9 h" GPL23=[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 w2 Z& Z# z0 C, B1 l# G) D) A1 T1 gPL24=[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: h" {- B  H) 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
) g7 `! O: q* s$ e7 L: s; u( TPL26=[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$ E% A- G+ a# d$ x
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- `  {, X/ ~9 _# Y
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小时负荷数据 MW
: A# ~1 S0 m" B( f' y( XPL29=[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
/ F1 V# n* F) u* T, s+ P+ mPL30=[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- J0 R4 O+ {0 l: u0 N% 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
) Z2 [+ r! A" ]" J1 B# QPL32=[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
5 b, w% B, l/ ?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; s0 A/ J% X  ~% A8 C3 D' \
PL=[PL2L3L4L5;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
% h( h$ ]* `+ K1 b9 j+ zPij=zeros(768,1);         %初始化输电网电量(Psub即Pij)MW(32条支路32*24=768)
, }* r+ Q$ h+ ]( @5 O0 V9 l$ N2 ~+ j  CPDG1=zeros(24,1);         %初始化分布式电源给DG1电量 MW, S# w/ V3 t$ C$ x
PDG2=zeros(24,1);         %初始化分布式电源DG2电量 MW. P/ i' E3 R1 i2 t
Psub=zeros(24,1);         %初始化Psub矩阵
& ~8 ]2 i) N6 X7 r$ BPA=zeros(24,1);           %初始化主动配电网与聚合站交易电量 MW0 J& O% |+ o: g8 I: Q2 S% _2 h+ P
SA=zeros(24,1);           %初始化聚合站电池容量 MWh* |, Q  q% E' Y. ?. I1 e1 M
x0=zeros(864,1);          %变量初始化        . K! n! u; L0 r- G4 F
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
2 {2 H9 ^, S& |" N  C: g( uAeq=zeros(793,864);       %等式约束系数矩阵初始化
/ |/ U1 b' n. b, E! ^2 Nbeq=zeros(793,1);         %等式约束结果矩阵初始化
5 k$ _. ?' j, [3 [( C% ^for i1=1:793
: D7 b; n8 u' A5 F, O) C/ q%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%- }: |& p2 X5 B; L' D8 c1 Z! {
    if i1>=1&&i1<=144          %节点2-节点7这些非特殊节点:Pij=PL) d2 _1 D0 i0 [, w* m% U! c
        Aeq(i1,i1)=1;
0 U/ O- L: R' T        beq(i1,1)=PL(i1,1);
4 H( N/ n  j% _! s2 K    end* e6 ~$ v# w  N' _- U! J  V
    if i1>=145&&i1<=168        %节点8是带PV的节点:Pij+PV=PL即Pij=PL-PV
$ a; ?) h: p& K; J         Aeq(i1,i1)=1;
+ g5 `6 j- I, q, F' Z  h/ Z         beq(i1,1)=PL(i1,1)-PV(i1-144,1);
  V, s( a; P$ e5 F" D4 G: t    end
/ K; g' O0 r/ k: }- G; l    if i1>=169&&i1<=312        %节点9-节点14为非特殊节点:Pij=PL2 B& Z  d4 j  H) R6 D) m4 u
         Aeq(i1,i1)=1;& z; v! e" h6 }) ]
         beq(i1,1)=PL(i1,1);- a) `# l# Y; A: f
    end- K7 P0 P: t4 `& F( c5 @. u
    if i1>=313&&i1<=336        %节点15是带PEV的特殊节点:Pij+PA=PL
# D7 b% V$ B' x2 T& N4 F         Aeq(i1,i1)=1;
0 @$ `8 J- p1 F: M; k: x8 I- m/ h4 U         Aeq(i1,i1+504)=1;  p9 T& z* j. w* z2 s6 L% O% x
         beq(i1,1)=PL(i1,1);
4 a) g3 J: Q7 V* o- Q    end! @" Y. A* B( R* }
    if i1>=337&&i1<=384        %节点16-节点17为非特殊节点:Pij=PL
6 @* g* n+ g2 x) D- j& Z% a        Aeq(i1,i1)=1;# T+ D, f- |+ q- ]: b$ T5 U2 F
        beq(i1,1)=PL(i1,1);& x' J2 q' X: k% n/ c% @4 l! q
    end7 `3 b$ n5 q) }7 q
    if i1>=385&&i1<=408       %节点18是带DG1的特殊节点:Pij+PDG1=PL8 ~6 G3 a' v; C' }
        Aeq(i1,i1)=1;/ x9 N5 |& V; r
        Aeq(i1,i1+384)=1;7 s7 a) m0 [5 r$ G$ O0 j3 F4 b
        beq(i1,1)=PL(i1,1);- j, S, h, c2 P7 g  m  {) p( ~
    end
0 R; z" C2 U6 y! H' Z' C    if i1>=409&&i1<=456       %节点19-节点20是非特殊节点:Pij=PL
, u# v7 n: I. v: h& P! r$ [" x        Aeq(i1,i1)=1;
$ @/ r: G0 l; h8 G+ f        beq(i1,1)=PL(i1,1);
8 {7 b0 g! X( u0 W3 ]    end9 g* G+ G- c" F+ o* k- I  J
    if i1>=457&&i1<=480      %节点21是带WT的特殊节点:Pij=PL-WT8 U' O6 |/ @9 N. i: O) D) \4 H
        Aeq(i1,i1)=1;
: f) I, u  k+ F        beq(i1,1)=PL(i1,1)-WT(i1-456,1);3 i4 E/ G- [6 Q8 J. s/ ]# R5 g
    end; N, Q! D6 A4 u0 I) p' Z) b8 J% g
    if i1>=481&&i1<=648      %节点22-节点28为非特殊节点:Pij=PL3 t) {: t& _1 I
        Aeq(i1,i1)=1;6 y  r/ F9 Y- {8 D
        beq(i1,1)=PL(i1,1);! @7 m% p/ R3 F4 @
    end
0 X+ n1 Z/ @# g$ U! {    if i1>=649&&i1<=672     %节点29为带DG2的特殊节点:Pij+PDG2=PL
' V* o, E$ r, O$ P4 X9 s# a        Aeq(i1,i1)=1;
! |, A! l! A3 g  ?8 s: r        Aeq(i1,i1+144)=1;
  c% N$ g3 A2 k# n        beq(i1,1)=PL(i1,1);7 C/ }3 V* a$ M7 A
    end
1 R+ h( s8 b4 I7 o1 I; y0 {    if i1>=673&&i1<=768      %节点30-节点33为非特殊节点:Pij=PL
" Y& i; c0 i' G; |9 x. N7 t  ^        Aeq(i1,i1)=1;* r9 c: o' x+ u, ^8 H( j
        beq(i1,1)=PL(i1,1);
3 e& s3 A* |) @2 l    end
8 R+ X7 O2 i6 v4 c0 e, F" o%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%. v* J  ]  R4 x, G3 C
   if i1>=792&&i1<=792
( x1 B, [# O; J# x+ q5 v       Aeq(i1,i1+72)=1;      %对应电池容量SA
! V* t, r$ Y  B" z5 A       for i2=1:24
, ]" S" A0 F3 ~* v$ h$ t7 }" b           if i2<=i1-768* ?! @3 R# O. P7 L: L/ H$ a
               Aeq(i1,i2+816)=1;8 Z8 J; j% H: }4 d$ W; Z2 h
           end
) ~# r" @- U( C9 K       end
( ~  M! x3 S/ L8 O: ~9 X       beq(i1,1)=5.4;        %5.4MWh& [! O/ V- A2 Y. }6 G9 b0 ]2 u
   end( j8 g" F0 n9 x3 o. p6 T
%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
+ X$ A% @( a2 c2 V& ?8 j   if i1>=793&&i1<=793       %对应S24=5.4MWh. }! H, M& N0 @: m5 B, N# L
       Aeq(i1,i1+71)=1;8 A) X& r/ i4 A
       beq(i1,1)=5.4;1 [1 \/ f- g( g; l. h1 P
   end7 R* M) R* x% {3 z3 A: e) ^
end & o$ M9 {+ [! @# o- p
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%3 b' y( Q8 I2 K4 M8 j
A=zeros(48,864);         %初始化线性不等式约束系数矩阵
. i: k6 X# |9 F7 u6 y) ?' Kb=zeros(48,1);           %初始化线性不等式约束结果矩阵
  o* `0 ^0 H$ J; c/ Rfor i3=1:48                 8 D) I! D5 [* j
     if i3>0&&i3<=24                                %约束条件VI<0.5,(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t# {/ {2 D& R$ V
         A(i3,i3+816)=1;                           %对应(PA)t的系数
1 Q* g8 A* a) {3 w" k% f         if i3==14 c* q* z& i( W. }2 {
             A(i3,840)=-1.5;                       %1时刻的上一时刻为24时,对应t=1时刻的(PA)t-1的系数; Y1 z; n( K! |
             b(i3,1)=1.5*WW(24,1)-WW(1,1);         %对应t=1时刻的不等式约束结果矩阵# L1 y: r7 ]5 s7 l# I+ _% P
         else
# L! y6 A* F" V& V5 k  f             A(i3,i3+815)=-1.5;                    %对应其他时刻(PA)t-1的系数+ x! `0 |5 b# ?9 R( x$ C; C
             b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1);      %对应其他时刻的不等式约束结果矩阵
/ p& _( p7 n0 a& I" Q, G- }         end& R0 ^* l# U/ X& i# B8 B
     end
* J7 G: B- }- p- p/ l      if i3>24&&i3<=48                              %约束条件VI>-0.5,(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-13 M( w$ o' k. U. L( c
         A(i3,i3+792)=-1;                          %对应(PA)t的系数: ?5 \& i4 x: b7 R" h. f
         if i3==25
+ B5 j- w6 W  I$ _& _5 r& A             A(i3,840)=0.5;                        %1时刻的上一时刻为24时,对应t=1时刻的(PA)t-1的系数
/ Z4 C5 e0 V' _1 Q             b(i3,1)=WW(1,1)-0.5*WW(24,1);         %对应t=1时刻的不等式约束结果矩阵
) _% ^* X7 ?) t/ D         else
7 J" U- o! u; n; Y$ B* S5 Z; N             A(i3,i3+791)=0.5;                     %对应其他时刻(PA)t-1的系数        0 x% k5 Y0 b& T
             b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
' E9 S  y) Q% U* p1 ]         end
0 N, e1 i. J- k. S# }9 p0 |* ~      end
6 V1 c9 R. L9 _; x0 o$ I3 nend3 l- I$ p- ~9 P9 O7 T  e1 U3 u
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
6 [4 c8 j8 h3 }+ ?Ib=zeros(864,1);                    %初始化下限矩阵
2 w  m. \: L# S6 Q; BUb=zeros(864,1);                    %初始化上限矩阵* F% K: z) ~6 r  Z! I/ Z
for i4=1:864                        %Ib Ub赋值. Q) ]6 ]% Y* H5 C$ n5 w
    if i4>0&&i4<=768                %Pij上下限:3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求6 M) y" j/ F" ^; y, C) n( }# t
        Ib(i4,1)=0;7 u" ?$ Z) J5 [
        Ub(i4,1)=3.73;
% O( f. H; ?% u1 Z6 g- C$ w    end7 h- M- A: x3 |' c
    if i4>=769&&i4<=792              %DG1功率上下限:3>=PDG1>=0 MW9 D8 t: \3 ?  q, p
        Ib(i4,1)=0;
0 @" R/ y4 G1 v  z        Ub(i4,1)=3;8 [3 E% \% R, H9 W1 K
    end
: y  s- l/ N7 ?! m) W    if i4>=793&&i4<=816              %DG2功率上下限:3>=PDG2>=0 MW1 ]$ k4 l' Y  M: {* R% E
        Ib(i4,1)=0;
( O! Q) x+ Z' p& s% P3 z) {2 e        Ub(i4,1)=3;3 ?$ L$ S6 }  p6 x6 d
    end
+ O, h) I/ l6 e. ~  }1 t    if i4>=817&&i4<=840              %聚合站充电功率上下限:2>=PA>=-2  MW
. |: ?- L% k/ Y0 Y        Ib(i4,1)=-2;4 b( T: G# z0 L) {6 M  q4 e
        Ub(i4,1)=2;
0 \- c! h, |3 w    end& O; v2 e/ s0 K* @7 U
    if i4>=841&&i4<=864              %电池容量上下限:5.4>=SA>=1.08 MWh* ?2 q% _. N9 \3 E  X: m9 j+ [
        Ib(i4,1)=1.08;- E$ C( Z  a! h8 i
        Ub(i4,1)=5.4;
1 b2 D4 R+ a# r# h8 q8 l. \% o# E    end- J" k% H, h, a' O; t2 A+ s
end
; y: }' |  M8 n: ?7 I1 [%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%  R! y3 U% r: n
[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
0 Y, f& a9 J2 q. z2 L' n  ufor i5=1:768; u8 t4 M/ N. |% b5 A5 q
    Pij(i5,1)=x(i5,1);
& @+ W' L2 I% @. W6 a2 Qend
% j- [6 T+ J8 Z9 Q/ Ofor i6=1:24+ j9 x% R+ B5 L6 C2 w7 j9 J7 \
    PDG1(i6,1)=x(i6+768,1);
4 F6 t) H0 s( W" B* [    PDG2(i6,1)=x(i6+792,1);8 J) \) b  _9 k  |9 K
    PA(i6,1)=x(i6+816,1);
9 h  O( s: v, o+ L- ]4 m    SA(i6,1)=x(i5+840,1);2 [0 E$ ^! u- x6 M8 G: X# [
end" \0 A( I' M. _2 }
for i7=1:244 ~1 n; J5 _4 s9 o- Y1 c
  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);: u* i! q) x/ ?1 B, F3 s1 E
end
1 l/ g/ D: r' t$ a8 T6 k9 Edisp(['运行成本:',num2str(fval1),'元']);  
" K$ t' v2 c  T, ~1 u7 lsubplot(2,1,1);* L: K# j+ q) e6 N! b+ k
plot(PA,'b-','LineWidth',1.5);
' T# Z7 Q. T  O. P( u% `hold on4 X* L0 O+ s/ m. G
grid on- ?3 o7 u- |; x" k" g& {) Z
1 e% j# M0 [2 I5 K# x/ ^- t

2 ~4 a* _- H9 ^

该用户从未签到

2#
发表于 2021-3-2 16:05 | 只看该作者
应该是fmincon函数里面的变量参数代入错误了,你再检查一下,可能是里面某个需要输入数值的变量你输成了向量

该用户从未签到

3#
发表于 2021-3-3 10:58 | 只看该作者
楼主解决了吗0 g5 h3 s2 Z, y+ N

该用户从未签到

4#
发表于 2021-3-3 11:02 | 只看该作者
来学习一下
您需要登录后才可以回帖 登录 | 注册

本版积分规则

关闭

推荐内容上一条 /1 下一条

EDA365公众号

关于我们|手机版|EDA365电子论坛网 ( 粤ICP备18020198号-1 )

GMT+8, 2025-10-6 09:20 , Processed in 0.171875 second(s), 26 queries , Gzip On.

深圳市墨知创新科技有限公司

地址:深圳市南山区科技生态园2栋A座805 电话:19926409050

快速回复 返回顶部 返回列表