|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
功能函数Z1=18.46-7.48*X1/X2^3，+ O: |& S# E2 D' r6 v6 @; m9 i, Z0 R
X1～N(10，2)，X2～N(2.5，0.375)，求验算点6 ?) \* h/ K3 t' i: o- \
P*(x1*，x2*)及结构可靠指标r。% p6 t$ ]/ S- y- O' `2 J
编程如下（采用二次多项式算法）& z% |6 ~( X+ i. c- C
Mu=[10 2.5];Sigama=[2 0.235];
/ T. E; [! }/ V* c" l& {( Q( k) RX0=[10 2.5];3 O7 I" l Y% D5 U d
A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];4 `3 }+ {- x5 N. r
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);& B* r! k4 Y8 \ d3 f/ @, u' f
bata=sqrt(fval);%计算可靠度指标值
/ b, m2 Z' ]( kPf=cdf('norm',-bata,0,1);%计算是效率; G# ?5 v7 o/ x( W& Y
function CC=bata2(X)0 n9 k/ a( o: S; h' j |; X
global Mu Sigama' a( V- M5 a- m0 F
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2
9 F r; @5 c4 s& S! r- E$ Qfunction [c,ceq ] =st(X)
2 K6 i) j# ^: ?( O6 Yc=[];$ A& \ o" O, u2 C' J6 w' `
ceq=xym(X);
' G% H1 a2 G$ y6 l$ @$ ]function Z =xym(X)3 P' r; [9 z. S7 e0 T- N
function Z=xym(X)9 J) j. V% ?2 R' A$ G/ i
%Z=18.46-7.48*X(1)/(X(2)^3);
& K9 Z0 _ w T* A; OA1=[1.0000 10.0000 2.5000 100.0000 6.2500;
) s* g6 q6 ^$ _% L. _2 l 1.0000 4.0000 2.5000 16.0000 6.2500;
! F* I4 O2 {# } 1.0000 10.0000 3.6250 100.0000 13.1406;7 {. o) `0 m) [) y/ b) B- F7 ?
1.0000 10.0000 1.3750 100.0000 1.8906];
. r- Q: t ~1 O/ P6 S4 Y" I0 xB1=[13.6728;10.8005;16.5451;16.8897;-10.3136];% ~) [2 e' ^1 s3 f
C=A1\B1;% Z: r1 O6 u. R5 r2 l }
Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);9 y! N6 g z( m( Z' C% j
0 m2 t2 G. ~1 d9 ~0 h' D
运行后出现：??? Error using ==> fmincon at 198 J( s& S! O! M7 H; H" p! T1 u
FMINCON only accepts inputs of data type double.
, B: W2 f5 t* j T4 f2 A8 D) S这是怎么回事？. S/ E% N. b( b- D1 h
|
|
/1 
发表于 2020-2-17 14:38
收藏
淘帖
支持！
反对！