|
|
楼主 |
发表于 2007-4-30 12:39
|
显示全部楼层
%蒙特卡罗直接抽样
function mcs_1
syms i x s mm k
mm=1000;
nn=mm^2;
s=0;
i=1;
u=[450 162 135];%矩阵U第一列代表抗力R,第二列代表恒载G,第三列代表活载L
p=[90.4 12.96 21.6];%矩阵p的三列分别为R,G,L的标准差
m=pi/(p(1,3)*sqrt(6));%m为极值I型当量正态化过程中的一个常数
n=u(1,1)-0.5772/m;%n为极值I型当量正态化过程中的另一个常数
uu=log(u(1,1)/sqrt(1+(p(1,1)/u(1,1))^2));
pp=sqrt(log(1+(p(1,1)/u(1,1))^2));
a=lognrnd(u(1,1),p(1,1),mm,mm);%对数正态分布随机抽样
b=normrnd(u(1,2),p(1,2),mm,mm);%正态分布随机抽样
c=norminv(rand(mm))*pp+uu*ones(mm);%极值I型分布随机抽样
for i=1:nn
if a(i)-b(i)-c(i)<0
s=s+1;
end
i=i+1;
end
bb=norminv(1-s/nn)%可靠度指标
disp(s)
我用的直接抽样，为什么一个失效点也找不到，应该是有解的，请大家看看那里应该修改啊？ |
|
发表于 2007-4-17 12:02
提升卡
变色卡
显身卡
发表于 2007-10-26 23:31
