结构可靠度分析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
注:需要参照可靠度的相关理论进行理解与分析哦~~
百度搜索“77cn”或“免费范文网”即可找到本站免费阅读全部范文。收藏本站方便下次阅读,免费范文网,提供经典小说综合文库结构可靠度编程在线全文阅读。
相关推荐: