fmincon优化

Uifhjvv

我在用fmincon优化求最小成本，但是出现：User supplied objective function must return a scalar value.求助大神！代码如下：

$ r" R/ W6 X5 g7 G5 J/ ~1 }: p. C; Nfunction f = adn(x)                  %主动配电网层目标函数
; e6 H, J9 l+ s/ o$ u$ Kglobal  A b WT PV j1 PA  Psub          %全局量
2 y4 [! u/ \& w& OADN=0;
for j=1:24
    ADN=ADN+Psub(j,1)*A(j,1);        %输电网部分的电价*电量
    ADN=ADN+x(j+768,1)*b;            %（输电网+DG1）部分的电价*电量
% o7 T) l' x& G2 j    ADN=ADN+x(j+792,1)*b;            %（输电网+DG1+DG2）部分的电价*电量
    ADN=ADN+WT(j,1)*b;               %（输电网+DG1+DG2+风电）部分的电价*电量            
    ADN=ADN+PV(j,1)*b;               %（输电网+DG1+DG2+风电+光伏）部分的电价*电量
0 Y! O  Y& v, ?' P    ADN=ADN+j1(1,j)*PA(j,1);         %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
end
! ~/ S9 u9 m8 T; Z0 Bf=ADN;



clc
+ z, I( i/ V( iclear all
& @, i4 C" h5 kglobal A b WT PV  PL j1 PA Psub
1 s, \5 j0 Q) ^7 a/ m; h+ qA= [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                 
) `: |5 c5 y1 d, M7 c( L8 G) m$ Eb= 400;                                                                                                                            %可再生分布式电源电价（=上网电价-补贴电价）元/MWh
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
1 d! c5 R% y" ^9 H* L. vWT = [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               
  S: z0 {) Z% G% oPV = [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
, L! s1 {1 v7 F0 `4 L3 u+ qWW=WT+PV;   %风电、光伏发电出力之和 MW
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
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
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
) C6 {, y+ x$ R  n0 Y2 }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
) E+ E( e: X( D- v- C( TPL6=[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
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
, M! I) l( b9 F1 k! h% HPL8=[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
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
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
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
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  
. x/ Z/ j1 f8 m$ u6 a5 ]2 \# V6 }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
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
% V8 N& }  r7 F' H, d( ZPL15=[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
' G( m4 Z7 F  `  `4 i8 M# IPL16=[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
+ e' z. L1 E4 D, @& W- L  ~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
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
9 S2 G6 x7 @# P0 d& EPL19=[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 Z* R3 ?2 T# s1 D+ D- OPL20=[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
- e' O( U' d4 E5 x. W% w  HPL21=[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
. Q" x+ w' `" I# h" p8 l' ^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
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
' O* @+ T4 j: E4 d+ BPL24=[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
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
1 g; Q% p* B( b* e" b& ^* J/ `5 HPL26=[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
4 ]1 T! {  m6 J+ J" ~( a" y. zPL27=[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
, D# s# J5 F& M+ S: W4 [5 rPL28=[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
) I* c6 c: }: G& {9 R4 \7 j3 |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
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
! C, g1 F/ g# NPL31=[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
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
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
( d, J9 d: _1 z: A" h4 ~* nPL=[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
Pij=zeros(768,1);         %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)
" \! B9 W2 F) q  y9 oPDG1=zeros(24,1);         %初始化分布式电源给DG1电量 MW
, I: c9 I* r8 l, FPDG2=zeros(24,1);         %初始化分布式电源DG2电量 MW
Psub=zeros(24,1);         %初始化Psub矩阵
PA=zeros(24,1);           %初始化主动配电网与聚合站交易电量 MW
5 s( [7 H; i! M1 zSA=zeros(24,1);           %初始化聚合站电池容量 MWh
  E" u0 u2 ^5 J- dx0=zeros(864,1);          %变量初始化        
/ I* S+ _- ~) Q%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
Aeq=zeros(793,864);       %等式约束系数矩阵初始化
beq=zeros(793,1);         %等式约束结果矩阵初始化
0 z% S9 f5 q& Q8 J2 e; m/ M4 s% Sfor i1=1:793
! O8 f" M! ~* m! Q1 u# a4 [- j%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%
: x1 H: i- J8 D7 H% e% X4 `5 {    if i1>=1&&i1<=144          %节点2-节点7这些非特殊节点：Pij=PL
        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
# o( h. `) s+ h    end
- N- G1 P! n: @! g4 P    if i1>=145&&i1<=168        %节点8是带PV的节点：Pij+PV=PL即Pij=PL-PV
         Aeq(i1,i1)=1;
         beq(i1,1)=PL(i1,1)-PV(i1-144,1);
! P" v: [" J. h. S6 a9 g    end
    if i1>=169&&i1<=312        %节点9-节点14为非特殊节点：Pij=PL
         Aeq(i1,i1)=1;
+ c* R+ O9 c6 q5 m1 |         beq(i1,1)=PL(i1,1);
+ v5 ~" b4 N) m% D    end
; U6 p+ n/ e& h9 y$ q9 }0 t    if i1>=313&&i1<=336        %节点15是带PEV的特殊节点：Pij+PA=PL
         Aeq(i1,i1)=1;
3 `" d8 A$ S. I9 Q! Z0 B7 b% B         Aeq(i1,i1+504)=1;
         beq(i1,1)=PL(i1,1);
- v# F  j' [  }% Y9 M0 G    end
    if i1>=337&&i1<=384        %节点16-节点17为非特殊节点：Pij=PL
        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
2 J8 z! Y' M# N( v' k    end
    if i1>=385&&i1<=408       %节点18是带DG1的特殊节点：Pij+PDG1=PL
/ `8 }: l% K) G& \* `1 z( n        Aeq(i1,i1)=1;
        Aeq(i1,i1+384)=1;
        beq(i1,1)=PL(i1,1);
/ |4 B2 ?( Z3 ?9 j  e2 f+ c% O9 n9 a- L" B    end
    if i1>=409&&i1<=456       %节点19-节点20是非特殊节点：Pij=PL
        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
    end
    if i1>=457&&i1<=480      %节点21是带WT的特殊节点：Pij=PL-WT
& ]: h2 O% b8 Z# B! D% l# t% V6 }        Aeq(i1,i1)=1;
, d8 K+ W/ q+ {/ l; _$ Q2 q        beq(i1,1)=PL(i1,1)-WT(i1-456,1);
0 K; d7 G: Z4 V, b* F$ y    end
% c9 w- h' c8 @! V4 K( @3 F4 z    if i1>=481&&i1<=648      %节点22-节点28为非特殊节点：Pij=PL
( u& ?! }7 N" e! j$ `% C        Aeq(i1,i1)=1;
+ q" J5 G' W2 C' P        beq(i1,1)=PL(i1,1);
" l; g; x8 }+ k4 ^9 {    end
- N; U1 }8 F' G/ u2 G1 @% V1 w7 E    if i1>=649&&i1<=672     %节点29为带DG2的特殊节点：Pij+PDG2=PL
8 b4 L( c! \# I* u- o2 S( D        Aeq(i1,i1)=1;
        Aeq(i1,i1+144)=1;
5 \! c+ l/ l5 l/ M- B; y        beq(i1,1)=PL(i1,1);
+ u, M* C/ H# ~6 u5 [8 P    end
    if i1>=673&&i1<=768      %节点30-节点33为非特殊节点：Pij=PL
        Aeq(i1,i1)=1;
  d4 c2 e* V3 ]4 C9 u% z& M0 `        beq(i1,1)=PL(i1,1);
# l$ r% T1 C- n8 O# O, \  _    end
. x% Z" Y% i( B$ W4 @( W+ ?/ G%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
   if i1>=792&&i1<=792
$ W4 q/ q. S) M       Aeq(i1,i1+72)=1;      %对应电池容量SA
       for i2=1:24
3 R- m9 J3 y; p8 r           if i2<=i1-768
               Aeq(i1,i2+816)=1;
           end
' @0 v3 W1 T) H) k5 y' w# x# b       end
       beq(i1,1)=5.4;        %5.4MWh
   end
  K5 x& ^1 j0 e: d# O1 \9 V%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
   if i1>=793&&i1<=793       %对应S24=5.4MWh
- u% M' u) w- Y" P       Aeq(i1,i1+71)=1;
       beq(i1,1)=5.4;
   end
5 h$ O* @! a: S0 M* T! C. ~9 Qend
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
A=zeros(48,864);         %初始化线性不等式约束系数矩阵
9 }8 `3 `. n3 H2 G! qb=zeros(48,1);           %初始化线性不等式约束结果矩阵
for i3=1:48                 
$ I0 c/ e: x+ J: G( n$ {- D: {9 D     if i3>0&&i3<=24                                %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
# H' E5 i# M( O. z         A(i3,i3+816)=1;                           %对应(PA)t的系数
         if i3==1
; q  C% T( \# m" C7 E7 A. [             A(i3,840)=-1.5;                       %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
) m$ Y$ e6 ~" B) Q# a/ n8 T5 }             b(i3,1)=1.5*WW(24,1)-WW(1,1);         %对应t=1时刻的不等式约束结果矩阵
: t; E5 o4 w9 L( [; h         else
* v; w8 j- i8 ~. n# C( s             A(i3,i3+815)=-1.5;                    %对应其他时刻(PA)t-1的系数
             b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1);      %对应其他时刻的不等式约束结果矩阵
, \" b+ p, F7 }& Z/ b         end
& V+ T- x* W+ A5 [     end
8 U$ F. W6 ^6 N* ^7 q+ y      if i3>24&&i3<=48                              %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
. @6 {3 `' e; S* o% c1 G2 o         A(i3,i3+792)=-1;                          %对应(PA)t的系数
         if i3==25
6 \: B# v6 Y9 |7 W( R2 h             A(i3,840)=0.5;                        %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
( c4 h+ l( S& l' O/ I             b(i3,1)=WW(1,1)-0.5*WW(24,1);         %对应t=1时刻的不等式约束结果矩阵
         else
             A(i3,i3+791)=0.5;                     %对应其他时刻(PA)t-1的系数        
" K" c% c* E3 {+ e             b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
         end
3 T" l7 ?* n+ a) Q      end
end
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
( x' R: Q0 ^: l% G$ cIb=zeros(864,1);                    %初始化下限矩阵
Ub=zeros(864,1);                    %初始化上限矩阵
for i4=1:864                        %Ib Ub赋值
    if i4>0&&i4<=768                %Pij上下限：3.73>=Pij>=0 MW,其中3.73MW为配电网中日最大有功负荷需求
, O# X' r8 |0 ]  _3 g+ m        Ib(i4,1)=0;
. R8 j2 l) q5 a        Ub(i4,1)=3.73;
    end
3 e  z0 D0 l1 A2 m, f$ g+ Q5 r    if i4>=769&&i4<=792              %DG1功率上下限：3>=PDG1>=0 MW
        Ib(i4,1)=0;
/ i) w* B' J( W3 z        Ub(i4,1)=3;
    end
    if i4>=793&&i4<=816              %DG2功率上下限：3>=PDG2>=0 MW
) U2 {1 I, G7 G        Ib(i4,1)=0;
        Ub(i4,1)=3;
    end
2 A' @- l/ g: G) ?6 N; O! i    if i4>=817&&i4<=840              %聚合站充电功率上下限：2>=PA>=-2  MW
        Ib(i4,1)=-2;
3 x% |' e5 c6 }$ X        Ub(i4,1)=2;
    end
% A) \  U% j8 d* [( R    if i4>=841&&i4<=864              %电池容量上下限:5.4>=SA>=1.08 MWh
" z& X! N3 i8 m4 e! ]$ {; B        Ib(i4,1)=1.08;
6 o1 i- v. N' B+ q  }        Ub(i4,1)=5.4;
9 D6 Y- J/ p8 f    end
) l( P8 `, B! j2 Tend
%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
1 E; r4 ]7 V5 `[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
" t$ g7 O2 Q' `' k" |for i5=1:768
* }$ Z/ R- P  u; a" H    Pij(i5,1)=x(i5,1);
& f' u- ~. W, U5 D& bend
for i6=1:24
9 f3 M6 w/ ~8 e6 v& n$ T( @6 v    PDG1(i6,1)=x(i6+768,1);
% b& g( b& i2 u" W    PDG2(i6,1)=x(i6+792,1);
    PA(i6,1)=x(i6+816,1);
    SA(i6,1)=x(i5+840,1);
" p7 c/ y  }. I- t% p8 U0 Tend
3 M8 j$ }& p$ z7 I) P& Qfor i7=1:24
  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);
end
disp(['运行成本：',num2str(fval1),'元']);  
subplot(2,1,1);
6 J0 {# b) i% c) Q: x2 ?, kplot(PA,'b-','LineWidth',1.5);
1 m) ?/ H9 h$ ]& P% _+ p" r! Xhold on
/ l" E6 s* h5 X  Ygrid on
. ?- Q- X# x5 l6 n- U

4 y' M. C+ }9 e3 z" m1 N( t4 L7 [

shuddkk

应该是fmincon函数里面的变量参数代入错误了，你再检查一下，可能是里面某个需要输入数值的变量你输成了向量

勇往直前11

楼主解决了吗

llbnmo

来学习一下