一个fmincon优化的问题

dsahdajs

功能函数Z1=18.46-7.48*X1/X2^3，
X1～N(10，2)，X2～N(2.5，0.375)，求验算点
P*(x1*，x2*)及结构可靠指标r。
% S# r6 i, o9 G# M1 }6 c. T编程如下（采用二次多项式算法）
Mu=[10 2.5];Sigama=[2 0.235];
6 p2 t# J" v( [X0=[10 2.5];
A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];
! P! n6 n, E  I7 j6 V  a! s[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);
( w9 M2 u  U: a0 }6 z. T, f& l; [bata=sqrt(fval);%计算可靠度指标值
Pf=cdf('norm',-bata,0,1);%计算是效率
function CC=bata2(X)
global Mu Sigama
# l- w! j; d! u' v6 SCC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
2 c' x/ l8 A& K" \' X8 [function [c,ceq ] =st(X)
( n. x# [  W9 P% L: Ic=[];
) v$ P5 q0 d0 X' vceq=xym(X);
' Y( S, F1 {2 N; {, m+ Ofunction Z =xym(X)
5 F4 S' a8 @# f  Z1 Afunction Z=xym(X)
%Z=18.46-7.48*X(1)/(X(2)^3);
( c0 ?% Z/ O  n, lA1=[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];
B1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
9 y( n! a2 Y. L% A$ sC=A1\B1;
Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);

运行后出现：??? Error using ==> fmincon at 198
; f$ L! X/ m9 D! P0 ^0 U( ZFMINCON only accepts inputs of data type double.
这是怎么回事？

yin123

顶一个