|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
功能函数Z1=18.46-7.48*X1/X2^3，
$ O1 ^" e- g% n4 x$ L7 xX1～N(10，2)，X2～N(2.5，0.375)，求验算点. B% m8 O6 R) O! m$ r% i' T& o
P*(x1*，x2*)及结构可靠指标r。6 z' b7 n t' Y& m# e. q8 e! \
编程如下（采用二次多项式算法）
1 b5 v% Q, n* j2 e. r7 AMu=[10 2.5];Sigama=[2 0.235];" N. t. q) {; S
X0=[10 2.5];. Z, i, ~- e5 F9 Y' b
A=[];b=[];Aeq=[];beq=[];lb=[];ub=[];7 a/ a+ M, r( o
[X,fval,exitflag,output]=fmincon(@bata2,X0,A,Aeq,beq,lb,ub,@st);, F+ p* a2 `6 w/ Q: R
bata=sqrt(fval);%计算可靠度指标值
' Q5 m2 `8 F9 d/ EPf=cdf('norm',-bata,0,1);%计算是效率
% w& |4 v- o/ F! Cfunction CC=bata2(X). `+ E$ } {4 T& B+ G# c7 q& A6 Y/ m
global Mu Sigama$ [1 M4 D5 S+ x
CC=((X(1)-Mu(1))/Sigama(1))^2+((X(2)-Mu(2))/Sigama(2))^21 u# A& `! n- v. `$ {& G
function [c,ceq ] =st(X)
6 A, j, p% Z3 Y0 o# u5 U6 Cc=[];
5 R) X! g3 f! O9 h0 @5 l; Mceq=xym(X);
$ o, m$ Z; F& ~2 R1 tfunction Z =xym(X)5 [8 ^& n0 l; v" F" p: }
function Z=xym(X)
/ K7 R0 ~ j* X1 [- @$ D6 G%Z=18.46-7.48*X(1)/(X(2)^3);
! s2 b+ W' V- {A1=[1.0000 10.0000 2.5000 100.0000 6.2500;, E5 j& l6 y1 u7 e6 S& @
1.0000 4.0000 2.5000 16.0000 6.2500;$ R' Q4 c9 i) v4 j$ ^' X
1.0000 10.0000 3.6250 100.0000 13.1406; M6 u9 N3 w* z& W9 ~
1.0000 10.0000 1.3750 100.0000 1.8906];
$ W w, {0 z& O$ z+ rB1=[13.6728;10.8005;16.5451;16.8897;-10.3136];
+ ?1 e. ^; C$ J6 S0 E# EC=A1\B1;0 @" z0 g5 Z6 W' k0 y
Z=C(1)+C(2)*X(1)+C(3)*X(2)+C(4)*X(1)+C(5)*X(2);. z, \4 v3 [% x/ r# h, V2 }
8 U9 V* B( S. o& A$ n
运行后出现：??? Error using ==> fmincon at 1987 r' J; e+ }" l/ y+ p- m2 M( {- `! ?
FMINCON only accepts inputs of data type double.
/ G' N3 W3 ^/ D. f这是怎么回事？& b8 F" g5 A6 J( E5 v8 a8 _& e
|
|
/1 
发表于 2020-2-17 14:38
收藏
淘帖
支持！
反对！