|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
功能函数Z1=18.46-7.48*X1/X2^3，
" c9 b! k5 Q2 W/ L# b+ y: YX1～N(10，2)，X2～N(2.5，0.375)，求验算点
' l4 Z; \+ L5 n2 t7 b6 H4 `; i) OP*(x1*，x2*)及结构可靠指标r。% w3 `. T6 L$ J, F0 t
编程如下（采用二次多项式算法）
( l6 H) }8 _, a' [/ M$ \Mu=[10 2.5];Sigama=[2 0.235];
, r2 @4 Y9 h( R9 W& cX0=[10 2.5];
' `. i2 Y2 v. n: ?% x! VA=[];b=[];Aeq=[];beq=[];lb=[];ub=[];% z9 V: [' Q; m+ j$ @3 b
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);: p$ q! ~1 v/ m. H/ L. W
bata=sqrt(fval);%计算可靠度指标值
. ? }( K4 T- zPf=cdf('norm',-bata,0,1);%计算是效率3 N* s- K1 a# R6 A1 h
function CC=bata2(X)
+ r$ [( v2 j/ ~; Sglobal Mu Sigama- ^! J) U7 C. W, a6 p
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^2' y6 O% {; I- A v( G$ {; ?
function [c,ceq ] =st(X)
4 F! O3 q2 t I7 x4 N* p2 pc=[];
0 X! V/ C5 b' g" g5 U$ c: rceq=xym(X);: w: C6 Y% T4 Q
function Z =xym(X): q0 A5 Q( P* F& U6 t1 i
function Z=xym(X)
% E( h+ C( P0 q( ~; |) z%Z=18.46-7.48*X(1)/(X(2)^3);
2 p. S* c. S, v$ W0 ]8 y- k- L4 f0 y3 sA1=[1.0000 10.0000 2.5000 100.0000 6.2500;
8 V. R) H( p8 K 1.0000 4.0000 2.5000 16.0000 6.2500;
. H3 m+ {6 M6 R/ L8 { 1.0000 10.0000 3.6250 100.0000 13.1406;
9 P7 q y' n- e) u 1.0000 10.0000 1.3750 100.0000 1.8906];$ l, ~7 l1 N. L2 F: f6 j
B1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
" O4 \& T6 ]8 c4 F* HC=A1\B1;, o( [3 \' c& s7 K
Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);
' r p+ u1 ~! C, H( d/ y+ i# [% H
运行后出现：??? Error using ==> fmincon at 198
- V( T8 i$ j5 y! l% k) n {, R) QFMINCON only accepts inputs of data type double.* E6 c% N2 E! C& ]" x
这是怎么回事？
# n2 K, Z7 c% v% z |
|
/1 
发表于 2020-2-17 14:38
收藏
淘帖
支持！
反对！