承受恒載和樓面活荷載的鋼筋混凝土軸心受壓短柱，已知恒載產方法一muX=[4560,765.5,1159.1];%均值cvX=[729.6/4560,222/765.5,81.1/1159.1];%變異系數mLn=log(muX(1)/sqrt(1+cvaEv=pi/(sqrt(6)*sigmaX(2));%求極值I型分布的參數psi=0.5772while%記錄循環次數%可靠度的初始值%可靠度的第二個值%抗力R當量正態化muX1(1)=muX(1)*(1+mLn-log(mt=exp((-aEv)*(muX(2%活載L當量正態化%標準正態概率密度反函數求值%標準正態概率密度求值muX1(2)=muX(2)-a*sigAlphar=-sigmaX1(1)/w;%求方向余弦x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1):%求循環后驗算點的坐標值x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*disp(結果：)fprintf(循環次數m:m=%d\n',m);fprintf(可靠度指標貝塔：bbeta=%1.2f\n',bbeta);fprintf(最后驗算點坐標：muX=[%1.2f;%1.2f;%1.2f]\n’,muX);循環次數m:m=6cvX=[729.6/4560,222/7M-file如下：whileabs(norm(x)-normX)/normX>1e-6%記錄循環次數%向量x二范數的前后兩個值的相對誤差%計算x的二范數%抗力R的當量正態化muX1(1)=muX(1)*(1+mLn-log(mut=exp((-1/aEv)*(muX(2)-uEv));%活載L的當量正態化muX1(2)=muX(2)-a*sigx(1)=muX1(1)+Alphar*bbeta*x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*fprintf(可靠度指標貝塔：bbeta=%1.2f\n',bbeta);fprintf(最后驗算點坐標：muX=[%1.2f;%1.2f;%1.2f]\n',muX);循環次數m:m=7cvX=[729.6/4560,222/7-fmLn=log(muX(1)/sqrt(1+cv%計算標準差以便下面使用%求極值I型分布的參數muX1(1)=muX(1)*(1+mLn-log(mu%對抗力R進行當量正態化%對活載L進行當量正態化x(1)=muX1(1)+Alphar*bbeta*sigmaX1(1);%求循環后驗算點的坐標值x(2)=muX1(2)+Alphal*bbetamuX=x;%將x賦予給muX以便下次循環使用fprintf(可靠度指標貝塔：bbeta=%1.2f\n',bbeta);fprintf(最后驗算點坐標：muX=[%1.2f;%1.2f;%1.2f]\n',muX);方法一.cvX=[0.17,0.29,0.07];%變異系數muX1(2)=muX(2)-a*sig%計算x的二范數%活載L的當量正態化x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbetafprintf(循環次數m:m=%d\n',m);fprintf(最后驗算點坐標：muX=[%1.2f;%1.2f;%1.2f]\n’,muX);循環次數m:m=9方法二cvX=[0.17,0.29,0.07];%變異系數%求極值I型分布的參數%對活載L進行當量正態化muX1(2)=muX(2)-a*sigsigmaX1=[sigmaX1(1)Alphar=-sigmaX1(1)/w;%求方向余弦x(2)=muX1(2)+Alphal*bbetax(3)=muX1(3)+Alphag*bbeta*x(1)=x(2)+x(3);%求循環后驗算點的坐標值x=[x(1);x(2);x(3)];%將計算的驗算點坐標放入數組xmuX=x;