一个fmincon优化的问题

dsahdajs

功能函数Z1=18.46-7.48*X1/X2^3，
6 K; Q8 c: b$ S- W2 a: I' T# hX1～N(10，2)，X2～N(2.5，0.375)，求验算点
. A0 K# T7 _  s$ h1 OP*(x1*，x2*)及结构可靠指标r。
编程如下（采用二次多项式算法）
Mu=[10 2.5];Sigama=[2 0.235];
X0=[10 2.5];
A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];
; A; V5 M* @( B. ~" K' C+ a[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);
; W) l/ s/ b1 R1 \bata=sqrt(fval);%计算可靠度指标值
Pf=cdf('norm',-bata,0,1);%计算是效率
$ Q! u7 C1 N* j+ P  kfunction CC=bata2(X)
global Mu Sigama
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
6 b7 q2 p' O8 K9 G) ~5 d* x+ B4 I7 qfunction [c,ceq ] =st(X)
8 u! X+ Q1 W. ~# @8 f$ ?c=[];
ceq=xym(X);
function Z =xym(X)
function Z=xym(X)
# b, N8 r2 @* K4 ^4 {%Z=18.46-7.48*X(1)/(X(2)^3);
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];
0 \# b) a- {6 w, y/ o% ?) a1 jB1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
/ ?/ Z- c1 l( h- ^2 o( x# eC=A1\B1;
Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);
/ n- H. t1 t4 }2 S9 S
运行后出现：??? Error using ==> fmincon at 198
4 U" \0 h- s; v5 F/ x' YFMINCON only accepts inputs of data type double.
- E$ L* V9 A/ c1 c. }0 ~) C/ g这是怎么回事？
4 j6 N$ r& h7 y  u

yin123

顶一个