| 查看: 439 | 回复: 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'. 什么原因 我是照着张明编写的结构可靠度分析--方法与程序 中照写的 |
回复此楼» 猜你喜欢
【长期招聘】复旦大学 AI 制药实验室 博士后招聘
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求介绍国内外做工程结构可靠度分析做得比较好的高校或学者
已经有10人回复
liaobudenan
木虫 (初入文坛)
- 应助: 1 (幼儿园)
- 金币: 1821.3
- 红花: 1
- 帖子: 46
- 在线: 119.8小时
- 虫号: 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
50