一个fmincon优化的问题

dsahdajs

功能函数Z1=18.46-7.48*X1/X2^3，
X1～N(10，2)，X2～N(2.5，0.375)，求验算点
P*(x1*，x2*)及结构可靠指标r。
; }& M/ i  O" _! Q% ~! _& g编程如下（采用二次多项式算法）
Mu=[10 2.5];Sigama=[2 0.235];
/ o, M9 [& V3 a1 |+ g, b! J) _X0=[10 2.5];
A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];
! ^+ }; X8 M9 v3 f/ D* u3 i[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);
bata=sqrt(fval);%计算可靠度指标值
7 I4 `- b  `  o0 i1 @- w( gPf=cdf('norm',-bata,0,1);%计算是效率
! L! v6 t) j& J! d0 dfunction CC=bata2(X)
global Mu Sigama
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
function [c,ceq ] =st(X)
c=[];
% o8 J$ c& o: ^+ l4 Kceq=xym(X);
function Z =xym(X)
function Z=xym(X)
* ^" i  V: Y& h# `, R" A) ~%Z=18.46-7.48*X(1)/(X(2)^3);
% p  v% Q1 z( y7 cA1=[1.0000 10.0000 2.5000 100.0000 6.2500;
    1.0000 4.0000 2.5000 16.0000 6.2500;
2 D/ z0 F/ d5 N) s9 y# g    1.0000 10.0000 3.6250 100.0000 13.1406;
    1.0000 10.0000 1.3750 100.0000 1.8906];
B1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
C=A1\B1;
% E! C5 l9 m% g$ ~, h" O8 WZ=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);
" z0 |4 x& m7 X' n( J: v
* z$ l, O0 w$ @4 n9 E+ S3 {. W% {, d运行后出现：??? Error using ==> fmincon at 198
FMINCON only accepts inputs of data type double.
这是怎么回事？
+ U: _( [; ]1 B$ n, O

yin123

顶一个