24小时热门版块排行榜    

查看: 984  |  回复: 6

cheey

银虫 (小有名气)

[求助] 求助:可靠度程序

极限状态方程为G=1.67*0.01^(-0.2)*5^(0.4)*600^(0.5)*d-5.36*0.01^(-0.4)*5^(0.6)e-f-600*g*j/h/5

参数        分布情况
f        对数正态分布(100,60)
g        正态分布(1,0.5)
h        正态分布(1.0,0.6)
d        正态分布(1,0.1)
e        正态分布(1,0.2)
j        正态分布(1,0.3)
这是我写的用蒙特法求解可靠度的程序,算出来的结果不太理想,可能程序那地方写错了,请大家帮我看看,谢谢
mm=1000;
nn=mm^2;
muX=[1;1;100;1;1;1];
sigmaX=[0.1;0.2;60;0.5;0.6;0.3];
d=normrnd(muX(1),sigmaX(1),mm,mm);
e=normrnd(muX(2),sigmaX(2),mm,mm);
f=log(lognrnd(muX(3),sigmaX(3),mm,mm));
g=normrnd(muX(4),sigmaX(4),mm,mm);
h=normrnd(muX(5),sigmaX(5),mm,mm);
j=normrnd(muX(6),sigmaX(6),mm,mm);
g=1.6*0.01.^(-0.2).*5.^0.4.*600.^0.5.*d-5.3*0.01.^(-0.4).*5.^0.6.*e-f-g.*600.*j./h./5;   o=find(g<0);
   s=length(o);
   bb=norminv(1-s/nn)
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cheey

银虫 (小有名气)

木有人回复啊
2楼2011-11-16 13:32:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiumo

新虫 (初入文坛)

写得很复杂,看不懂
3楼2011-11-23 17:31:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

geozb

木虫 (正式写手)

【答案】应助回帖

没有指出自己要计算的目标以及计算的远离,别人也没法给你指导的啊
要幸福的度过每一天,爱身边的人
4楼2011-11-26 09:12:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

jiumo

新虫 (初入文坛)

用matlab写的吗  感觉挺好的  算出来怎么不理想了
5楼2011-12-25 09:17:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lihai209

禁虫 (初入文坛)

本帖内容被屏蔽

6楼2013-12-15 12:53:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zxqqd

新虫 (小有名气)

你好,知道问题在哪里了吗,我最近也在计算一个可靠度,程序和你差不多,计算出来是负无穷大,不知道为什么
7楼2014-02-21 00:32:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cheey 的主题更新
信息提示
请填处理意见