| 查看: 383 | 回复: 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'. 什么原因 我是照着张明编写的结构可靠度分析--方法与程序 中照写的 |
回复此楼» 猜你喜欢
长年满屏的广告，版主太不责任了。基金也等的急
已经有3人回复
-
filecode
已经有12人回复
-
不应该看fileCode
已经有4人回复
-
应该是93bebmhtak前后十一个字符比较关键
已经有33人回复
-
FileCode能看出啥？
已经有19人回复
-
2019年青年基金涵评意见，大家看看几个A，几个B?
已经有5人回复
-
职称评审，求友友推荐见刊最快的期刊
已经有4人回复
-
奇怪，两个人的filecode固定段从头到尾一模一样
已经有10人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有13人回复
-
是这周出结果还是下周出结果？
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求介绍国内外做工程结构可靠度分析做得比较好的高校或学者
已经有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
