结构可靠度分析matlab编程命令流如下： function bbeta=qqjj(muX,sigmaX) [muX,sigmaX]=input;
sLn=sqrt(log(1+(sigmaX(1)/muX(1))^2));mLn=log(muX(1))-sLn^2/2; aEv=sqrt(6)*sigmaX(3)/pi;uEv=-psi(1)*aEv-muX(3); muX1=muX;sigmaX1=sigmaX; x=muX;normX=eps;
while abs(norm(x)-normX)/normX>1e-6
normX=norm(x); g=x(1)-x(2)-x(3); gX=[1;-1;-1];
cdfX=1-evcdf(-x(3),uEv,aEv); pdfX=evpdf(-x(3),uEv,aEv); nc=norminv(cdfX);
sigmaX1(3)=normpdf(nc)/pdfX; muX1(3)=x(3)-nc*sigmaX1(3); muX1(1)=x(1)*(1+mLn-log(x(1))); sigmaX1(1)=sLn*x(1);
gs=gX.*sigmaX1;alphaX=-gs/norm(gs); a=gX'*(muX1-x); bbeta=(g+a)/norm(gs)
x=muX1+bbeta*sigmaX1.*alphaX
end
注：需要参照可靠度的相关理论进行理解与分析哦～～