fmincon优化

Uifhjvv

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

' b) S2 U" ?1 Dfunction f = adn(x)                  %主动配电网层目标函数
- O, y/ X; q2 ~# Y- mglobal  A b WT PV j1 PA  Psub          %全局量
* x8 [5 p7 S: C+ K; Y1 dADN=0;
for j=1:24
- y  m* \4 z& m6 f3 ?: U- E    ADN=ADN+Psub(j,1)*A(j,1);        %输电网部分的电价*电量
% k% J' j& O4 J: a$ W    ADN=ADN+x(j+768,1)*b;            %（输电网+DG1）部分的电价*电量
  ?$ A; x2 k9 W1 N( s  R- C2 V. X    ADN=ADN+x(j+792,1)*b;            %（输电网+DG1+DG2）部分的电价*电量
: C, b7 f: B$ Q: p    ADN=ADN+WT(j,1)*b;               %（输电网+DG1+DG2+风电）部分的电价*电量            
    ADN=ADN+PV(j,1)*b;               %（输电网+DG1+DG2+风电+光伏）部分的电价*电量
- I3 M$ `. [. k9 i# s, R' n9 S    ADN=ADN+j1(1,j)*PA(j,1);         %（输电网+DG1+DG2+风电+光伏+聚合站）部分的电价*电量，即配电网运行成本
2 A. h0 m/ k) T; Y) `# f- iend
f=ADN;
; ^$ _0 V: ~+ p: }: p
) B( P, U1 }- g( |

clc
: w# f' \5 X- P! B+ A$ T1 Cclear all
global A b WT PV  PL j1 PA Psub
) {  r7 S8 [8 M. JA= [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                 
2 s2 @6 D5 V( B7 m. X6 [b= 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
, W5 T- ^; m# ^; l; B* ^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               
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
WW=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
: O9 L$ U' I* N9 ~! C3 q. d+ q  oPL4=[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
+ l" y' d3 i( c2 e5 R" APL5=[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
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
5 P9 M4 o: N6 Q, TPL7=[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
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
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
* ?0 B+ [: k, {7 s( ^/ VPL10=[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
6 }. v- t3 Z  A7 jPL11=[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
+ B1 n6 j$ X9 \( Y: j% N8 {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  
# q. V8 D8 X8 ]! R, `7 x. }! [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
# M2 Q  t- O# c: |% F2 X+ w. RPL14=[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
9 U0 k( p7 r4 E& y3 nPL15=[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
& e. C* J2 y  E" [3 ?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
9 L5 F) E' `5 UPL17=[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
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
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
; X6 [' }+ ~0 U* z0 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
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
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
  e! S' y: j$ L' TPL25=[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
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
4 G( p9 \8 z  q+ A( x/ uPL27=[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
. v( B3 Y$ F# a/ b& n* XPL28=[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
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
" v3 @( h8 [5 w+ V& ]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
# J5 _1 W) M. \! T7 p  Y9 vPL31=[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
6 v" H+ z2 F& j5 n5 WPL33=[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
2 C( R; [& i" V; l8 _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
, c  |) c7 ]8 G  k8 sPij=zeros(768,1);         %初始化输电网电量(Psub即Pij)MW（32条支路32*24=768)
, h; m2 x9 i7 w7 d) kPDG1=zeros(24,1);         %初始化分布式电源给DG1电量 MW
" g- T: \, g: L5 ?, i+ EPDG2=zeros(24,1);         %初始化分布式电源DG2电量 MW
- n% k( c- z9 ^' _Psub=zeros(24,1);         %初始化Psub矩阵
4 N8 g, r0 X1 lPA=zeros(24,1);           %初始化主动配电网与聚合站交易电量 MW
4 u. P  m; r2 |  W/ j0 W( Z$ cSA=zeros(24,1);           %初始化聚合站电池容量 MWh
- c) q- q  L. ox0=zeros(864,1);          %变量初始化        
%%%%%%%%%%%%%%%等式约束%%%%%%%%%%%%%%%%%%%%
5 U! i: g! c0 N' x' Z2 HAeq=zeros(793,864);       %等式约束系数矩阵初始化
! a; ]9 \/ U; ?( }beq=zeros(793,1);         %等式约束结果矩阵初始化
for i1=1:793
%%%%%%%%%%%%有功功率平衡等式约束%%%%%%%%%%%%%%%%%%%%%
+ A0 Y4 Z, {/ S# ^+ T    if i1>=1&&i1<=144          %节点2-节点7这些非特殊节点：Pij=PL
3 P, c- Y6 S/ N, Y& `( X* F        Aeq(i1,i1)=1;
+ ^0 R9 B- V2 V0 h        beq(i1,1)=PL(i1,1);
2 D( x( t, E/ y, x3 Q) |3 q, u. P( A5 \    end
* E. ^3 Z2 i5 R& z- v+ Q' j' q    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);
    end
) A; k* `  ^; `- z    if i1>=169&&i1<=312        %节点9-节点14为非特殊节点：Pij=PL
+ ^' D, L1 p* Y         Aeq(i1,i1)=1;
: g; k  k; C4 N0 ^3 B         beq(i1,1)=PL(i1,1);
+ T; v( ~- ^: @1 E  b    end
    if i1>=313&&i1<=336        %节点15是带PEV的特殊节点：Pij+PA=PL
% Y; k9 n* X, R3 N( p         Aeq(i1,i1)=1;
8 K( i% b8 @6 {+ O* ~1 H& L         Aeq(i1,i1+504)=1;
         beq(i1,1)=PL(i1,1);
7 u8 ^" m: U7 t( a/ ^, G( ^2 ^) B    end
    if i1>=337&&i1<=384        %节点16-节点17为非特殊节点：Pij=PL
& Q7 w8 j3 N7 f) O2 Q0 W! E9 }        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
    end
, m. f) t5 a  F0 r    if i1>=385&&i1<=408       %节点18是带DG1的特殊节点：Pij+PDG1=PL
( G5 X: K  C; D% S4 f. ~+ }        Aeq(i1,i1)=1;
7 g1 b; N4 ]& Z; b        Aeq(i1,i1+384)=1;
& X5 I; ^  ?1 Q$ Z; I1 r; _        beq(i1,1)=PL(i1,1);
    end
, {' C0 G' S6 u6 T7 ~    if i1>=409&&i1<=456       %节点19-节点20是非特殊节点：Pij=PL
) x- B6 T; i5 O, N5 R% ?) {4 R        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
9 w8 s$ `7 D3 R" W    end
3 p$ p* _& {; F8 N" I% Q/ {    if i1>=457&&i1<=480      %节点21是带WT的特殊节点：Pij=PL-WT
        Aeq(i1,i1)=1;
8 C. c* o" H7 P5 U( ]$ R/ u        beq(i1,1)=PL(i1,1)-WT(i1-456,1);
% F! L3 X4 W/ s7 j0 Q4 Z+ q  I    end
    if i1>=481&&i1<=648      %节点22-节点28为非特殊节点：Pij=PL
        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
% a6 F! `7 ?0 Y: f; x    end
$ e0 {& j1 D  F& j) b  X( ?    if i1>=649&&i1<=672     %节点29为带DG2的特殊节点：Pij+PDG2=PL
        Aeq(i1,i1)=1;
        Aeq(i1,i1+144)=1;
        beq(i1,1)=PL(i1,1);
. I: E3 C/ w5 [+ E: }/ t' e- p    end
    if i1>=673&&i1<=768      %节点30-节点33为非特殊节点：Pij=PL
        Aeq(i1,i1)=1;
        beq(i1,1)=PL(i1,1);
    end
' ?% `6 W6 a# s/ T%%%%%%%%%%%%%%%%%充放电电池容量等式约束%%%%%%%%%%%%%%%%%%%%
! D) t3 {7 L; m8 I) ^. y3 N& w   if i1>=792&&i1<=792
       Aeq(i1,i1+72)=1;      %对应电池容量SA
       for i2=1:24
) g- y1 {$ h$ G8 L, J. Q* ?  c           if i2<=i1-768
+ |% h- I+ ?$ k$ r* |               Aeq(i1,i2+816)=1;
           end
! {# Q$ a3 T$ I       end
       beq(i1,1)=5.4;        %5.4MWh
( l2 w$ B( L( ]  k& F# ?   end
  C1 u) I- R! ?. _9 }2 g& ^%%%%%%%%%%%%%%充放电容量平衡约束%%%%%%%%%%%%%%%%%%%%%%%%
9 j  S" w2 x7 C+ ]   if i1>=793&&i1<=793       %对应S24=5.4MWh
       Aeq(i1,i1+71)=1;
0 Z+ X' u8 Y/ d5 J       beq(i1,1)=5.4;
   end
/ W8 ^0 H" ?* W8 V0 uend
%%%%%%%%%%不等式约束%%%%%%%%%%%%%%%%%%
) K$ y0 y. p9 C! \6 `+ @, p3 |) OA=zeros(48,864);         %初始化线性不等式约束系数矩阵
0 i: r* T5 G9 U9 ^. A( Ab=zeros(48,1);           %初始化线性不等式约束结果矩阵
for i3=1:48                 
) _! _- [8 T  N" [     if i3>0&&i3<=24                                %约束条件VI<0.5，(PA)t-1.5(PA)t-1<1.5(WW)t-1-(WW)t
         A(i3,i3+816)=1;                           %对应(PA)t的系数
! X% h* p" f# L2 L9 _6 J7 N         if i3==1
2 q$ L% |- L- d3 [! w( G( A             A(i3,840)=-1.5;                       %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
* k# W- ^* ^: g: l9 Q; O9 ]/ E             b(i3,1)=1.5*WW(24,1)-WW(1,1);         %对应t=1时刻的不等式约束结果矩阵
         else
! L; G3 c4 v, ?$ P$ Z0 I             A(i3,i3+815)=-1.5;                    %对应其他时刻(PA)t-1的系数
( v$ |, l% }6 K; R! M9 M9 E  t             b(i3,1)=1.5*WW(i3-1,1)-WW(i3,1);      %对应其他时刻的不等式约束结果矩阵
         end
. ~6 ]8 |  M' j     end
      if i3>24&&i3<=48                              %约束条件VI>-0.5，(-PA)t+0.5(PA)t-1<(WW)t-0.5(WW)t-1
         A(i3,i3+792)=-1;                          %对应(PA)t的系数
( D7 F  W1 k6 @# H8 C6 F         if i3==25
! ^+ u& v5 e5 u             A(i3,840)=0.5;                        %1时刻的上一时刻为24时，对应t=1时刻的(PA)t-1的系数
             b(i3,1)=WW(1,1)-0.5*WW(24,1);         %对应t=1时刻的不等式约束结果矩阵
9 E3 e0 U& a& Q9 v; M         else
             A(i3,i3+791)=0.5;                     %对应其他时刻(PA)t-1的系数        
1 m9 w& c& C9 X: `: {             b(i3,1)=WW(i3-24,1)-0.5*WW(i3-25,1); %对应其他时刻的不等式约束结果矩阵
6 G* y9 g! w" b) E* y1 D         end
9 Z- D; ~9 L1 O9 J( u2 z3 x      end
end
%%%%%%%%%%%%%%%%%%上下限约束%%%%%%%%%%%%%%%%%%%%
Ib=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为配电网中日最大有功负荷需求
        Ib(i4,1)=0;
        Ub(i4,1)=3.73;
$ |3 K; Y8 ~2 X+ o  \2 g    end
/ m/ d# k' ]$ a& U9 J% z# h    if i4>=769&&i4<=792              %DG1功率上下限：3>=PDG1>=0 MW
5 p6 [7 H$ R7 N, y2 m5 Q        Ib(i4,1)=0;
  [; t- {) h# u$ ]" F        Ub(i4,1)=3;
    end
  N5 N& m9 ~7 F1 `5 u    if i4>=793&&i4<=816              %DG2功率上下限：3>=PDG2>=0 MW
        Ib(i4,1)=0;
        Ub(i4,1)=3;
    end
4 h& X0 Y7 k, h8 Q! p6 E% L: I    if i4>=817&&i4<=840              %聚合站充电功率上下限：2>=PA>=-2  MW
        Ib(i4,1)=-2;
6 H/ A5 ^  R2 @2 D  I7 [1 A        Ub(i4,1)=2;
    end
    if i4>=841&&i4<=864              %电池容量上下限:5.4>=SA>=1.08 MWh
" L9 r% ~% J7 [/ Z& E        Ib(i4,1)=1.08;
- c2 M, k1 m2 S. _- P        Ub(i4,1)=5.4;
    end
: o% `' k. K+ ~& r* j! x7 r- yend
1 J6 h! P7 u+ N, v( z%%%%%%%%%%%%%%%%成本优化%%%%%%%%%%%%%%%%%%%%%%%
1 C  r5 I) I% t& e[x,fval1] = fmincon(@adn,x0,A,b,Aeq,beq,Ib,Ub);
for i5=1:768
  b0 K  F0 `/ K    Pij(i5,1)=x(i5,1);
end
; N5 y$ I+ |/ R( Efor i6=1:24
5 w! L/ H  J1 g3 ?" l$ A$ g4 s- ]    PDG1(i6,1)=x(i6+768,1);
9 ]- h5 s  l; t! f5 z9 I' H    PDG2(i6,1)=x(i6+792,1);
5 |; k6 t% v# s9 X* l0 m/ _    PA(i6,1)=x(i6+816,1);
+ @- l3 r, ^, c9 j( g    SA(i6,1)=x(i5+840,1);
end
! X5 ?, \# ^5 Z* Q( c% N) ^, \9 Kfor i7=1:24
* w) q8 K, G/ j/ G# J  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);
0 ~6 e( B0 I; |! e- a8 a6 Mend
disp(['运行成本：',num2str(fval1),'元']);  
subplot(2,1,1);
: @  |) d( U; k- Y* fplot(PA,'b-','LineWidth',1.5);
hold on
grid on
% U) c  V+ {/ ]* N. k

1 I: H* v- w, S1 F

shuddkk

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

勇往直前11

楼主解决了吗

llbnmo

来学习一下