nctionrelia%定义主函数
x0= [ 0.6,2.18,3.49];%初始迭代点,这里取均值
%调用优化工具箱求解
options=optimset(′LargeScale′,′off′);
[x,fval]=fmincon(@ obj,x0,@st,options);
fprintf(%可靠度指标为bata=%f n′,fval);
fprintf(%验算点为[%f,%f,%f] n′,x(1),x(2),exp(x
(3)));
functionCC=obj(x) %目标函数子函数
CC=sqrt(x(1)- 0.6)/0.0786) 2+((x(2)-
2.18)/0.0654) 2+… ((x(3)- 3.49)/0.03)
2);%即为β
function[c,ceq]=st(x) %约束条件子函数
c=[];
ceq= 567*x(1)*x(2)- 0.5* exp(x(3)) 2;