%问题2、按原混凝土设计规范设计的轴心受压构件极限状态方程为：R－SG－SQ = 0，ρ= 2.0；
%求β。
%其中：KR = 1.33，δR = 0.17；
%KSG = 1.06，δSG = 0.07；
%KSQ = 0.70，δSQ = 0.29。
function [bb]=kkdmc
mm=1000;
nn=mm^2;
pp=0.2;%pp=sq/sg荷载效应比值
sg=10;
sq=sg*pp;
r=1.55*(sg+sq);%按一个分项系数设计
u=[1.33 1.06 0.7];%矩阵U第一列代表抗力R,第二列代表SG,第三列代表SQ
p=[0.17 0.07 0.29];%矩阵p的三列分别为R,SG,SQ的标准差
a1=log(lognrnd(u(1)*r,p(1)*r,mm,mm));%对数正态分布随机抽样
a3=normrnd(u(2)*sg,p(2)*sg,mm,mm);%正态分布随机抽样
aaa=pi/(p(3)*sq*sqrt(6));%m为极值I型当量正态化过程中的一个常数
bbb=u(3)*sq-0.5772/aaa;%n为极值I型当量正态化过程中的另一个常数
a=rand(mm);
a4=(aaa*bbb-log(log(1./a)))/aaa;
l=a1-a3-a4;%功能函数
o=find(l<0);%失效的次数的位置
s=length(o);%失效的次数
bb=norminv(1-s/mm^2)%可靠度指标
end
%结果3.8
