一个fmincon优化的问题

dsahdajs

功能函数Z1=18.46-7.48*X1/X2^3，
3 G( K+ n; b3 }. @0 I6 B; jX1～N(10，2)，X2～N(2.5，0.375)，求验算点
8 {0 S/ l' _, }P*(x1*，x2*)及结构可靠指标r。
: v/ o8 {& @4 ]3 y# u0 }编程如下（采用二次多项式算法）
) \0 l3 H; E$ e1 u2 dMu=[10 2.5];Sigama=[2 0.235];
+ B1 q! ?8 \8 g. uX0=[10 2.5];
  X+ j$ ], `+ K2 a- g& {A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);
/ ^& N5 R4 K" n2 i9 c7 }4 C  I* \# Mbata=sqrt(fval);%计算可靠度指标值
, v& D  X* O. QPf=cdf('norm',-bata,0,1);%计算是效率
3 Y; A- H6 R) I' ~% K6 x: ~function CC=bata2(X)
7 a% w" y! O9 b0 ]: Bglobal Mu Sigama
* I  f8 E3 V) ^9 B8 @# oCC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
function [c,ceq ] =st(X)
c=[];
. B: k4 G/ e) J/ S! {- @3 x2 bceq=xym(X);
; W) D' ~7 Q% `function Z =xym(X)
) W2 U0 f$ }1 h, H3 e( t- P: ifunction Z=xym(X)
% P' `4 V8 ?! k1 S%Z=18.46-7.48*X(1)/(X(2)^3);
- t# b$ V5 j7 T; iA1=[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;
. U7 E" s+ R; Y3 M: p5 c% j; \$ G    1.0000 10.0000 1.3750 100.0000 1.8906];
B1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
! }. {9 E' W9 }- K1 p( a( [# _8 E# LC=A1\B1;
) [. B; ~6 m" q3 B# r9 {Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);

4 [; H3 E. @- c, j% ~8 C运行后出现：??? Error using ==> fmincon at 198
' N. k4 @% D" c9 g" oFMINCON only accepts inputs of data type double.
这是怎么回事？

yin123

顶一个