| 查看: 286 | 回复: 2 | ||
xmszkx铁虫 (著名写手)
|
[求助]
结构可靠度分析中极值Ⅰ型分布总是出现错误提示已有1人参与
|
|
clear;clc; muX = [22;10;2];sigmaX = [2;0.9;0.6]; 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 =[logncdf(x(1),mLn,sLn);1-evcdf(-x(3),uEv,aEv)]; pdfX =[lognpdf(x(1),mLn,sLn);evpdf(-x(3),uEv,aEv)]; nc = norminv(cdfX); sigmaX1(1:2:3) =normpdf(nc)./pdfX; muX1(1:2:3) = [x(1:2:3)-nc.*sigmaX1(1:2:3)]; gs = gX.*sigmaX1; alphaX = -gs./norm(gs); bbeta = (g+gX'*(muX1-x))/norm(gs); x =muX1+bbeta*sigmaX1.*alphaX; end 程序运算总是出现 ??? Undefined function or variable 'evcdf'. 什么原因 我是照着张明编写的结构可靠度分析--方法与程序 中照写的 |
回复此楼» 猜你喜欢
三无产品还有机会吗
已经有6人回复
-
垃圾破二本职称评审标准
已经有7人回复
-
投稿返修后收到这样的回复，还有希望吗
已经有7人回复
-
博士申请都是内定的吗？
已经有14人回复
-
谈谈两天一夜的“延安行”
已经有13人回复
-
氨基封端PDMS和HDI反应快速固化
已经有11人回复
-
之前让一硕士生水了7个发明专利，现在这7个获批发明专利的维护费可从哪儿支出哈？
已经有11人回复
-
论文投稿求助
已经有4人回复
-
Applied Surface Science 这个期刊。有哪位虫友投过的能把word模板发给我参考一下嘛
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求介绍国内外做工程结构可靠度分析做得比较好的高校或学者
已经有10人回复
liaobudenan
木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 1821.3
- 红花: 1
- 帖子: 46
- 在线: 119.6小时
- 虫号: 932613
- 注册: 2009-12-24
- 专业: 结构工程
2楼2014-10-21 21:53:08
xmszkx
铁虫 (著名写手)
- 应助: 5 (幼儿园)
- 金币: 3373.1
- 散金: 126
- 红花: 8
- 帖子: 1868
- 在线: 315.7小时
- 虫号: 1857713
- 注册: 2012-06-13
- 性别: GG
- 专业: 结构工程
3楼2014-10-22 19:20:53