一个fmincon优化的问题

dsahdajs

功能函数Z1=18.46-7.48*X1/X2^3，
X1～N(10，2)，X2～N(2.5，0.375)，求验算点
P*(x1*，x2*)及结构可靠指标r。
编程如下（采用二次多项式算法）
Mu=[10 2.5];Sigama=[2 0.235];
X0=[10 2.5];
A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);
bata=sqrt(fval);%计算可靠度指标值
' S0 Y" x: |; m; p2 u& TPf=cdf('norm',-bata,0,1);%计算是效率
) k: ]2 @# V- Z. V% t3 Jfunction CC=bata2(X)
global Mu Sigama
* b  k. a/ K+ j1 RCC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
" j7 U8 F- t  ~: n. J8 u  dfunction [c,ceq ] =st(X)
$ ]5 o1 h; D- A) C) Yc=[];
$ D3 w- z9 G; w! rceq=xym(X);
# s- l! o; C4 vfunction Z =xym(X)
2 m2 x' K5 G$ T" x) E7 e' yfunction Z=xym(X)
%Z=18.46-7.48*X(1)/(X(2)^3);
- \$ ~4 Z. l' _A1=[1.0000 10.0000 2.5000 100.0000 6.2500;
    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];
) G" L% F' M$ xB1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
- q9 d2 h$ f1 H1 v1 Z" U$ dC=A1\B1;
' h! `/ I+ e/ pZ=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);
$ N" ?6 L3 \6 a% U, }
* f% F" q" k" J" E6 t+ L6 q运行后出现：??? Error using ==> fmincon at 198
FMINCON only accepts inputs of data type double.
这是怎么回事？

yin123

顶一个