一个fmincon优化的问题

dsahdajs

功能函数Z1=18.46-7.48*X1/X2^3，
+ ^% u9 H) X' @8 U4 [- q) lX1～N(10，2)，X2～N(2.5，0.375)，求验算点
P*(x1*，x2*)及结构可靠指标r。
$ ^6 m9 W; B8 l" J" _/ V编程如下（采用二次多项式算法）
Mu=[10 2.5];Sigama=[2 0.235];
( |# {5 P9 M+ g3 p& ?; oX0=[10 2.5];
* l2 \2 T! H4 ]( R( D: c/ TA=[];b=[];Aeq=[];beq=[];lb=[];ub=[];
  E9 U) ~. t. L( {/ T[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);
bata=sqrt(fval);%计算可靠度指标值
Pf=cdf('norm',-bata,0,1);%计算是效率
* I/ |6 l3 l+ `$ E: J! }9 Ufunction CC=bata2(X)
# K2 G! g) K" R' V, }global Mu Sigama
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
; [0 t! V/ _. B* N6 h/ K7 @$ [- efunction [c,ceq ] =st(X)
6 `9 p. L7 X8 `; k3 qc=[];
ceq=xym(X);
. h6 |/ o4 ?7 d0 j& Hfunction Z =xym(X)
function Z=xym(X)
%Z=18.46-7.48*X(1)/(X(2)^3);
A1=[1.0000 10.0000 2.5000 100.0000 6.2500;
/ p; f9 R; K0 f( j  c    1.0000 4.0000 2.5000 16.0000 6.2500;
    1.0000 10.0000 3.6250 100.0000 13.1406;
    1.0000 10.0000 1.3750 100.0000 1.8906];
, v8 c1 X/ R; P  z7 SB1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
C=A1\B1;
3 x$ z6 X. w+ ^; G$ c+ E$ Y2 Y1 }Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);

1 q5 T* P1 k( D& A, ~运行后出现：??? Error using ==> fmincon at 198
FMINCON only accepts inputs of data type double.
这是怎么回事？

yin123

顶一个