| 查看: 3140 | 回复: 1 | ||
| 【悬赏金币】回答本帖问题，作者雨天小木子将赠送您 5 个金币 | ||
[求助]
无法执行赋值，因为左侧和右侧元素数目不同已有1人参与
|
||
|
clear all clc muX=[300;550];cvX=[0.12;0.10]; sigmaX=cvX.*muX; sLn=sqrt(log(1+(sigmaX(1)/muX(1))^2)); mLn=log(muX(1))-sLn^2/2; aEv=sqrt(6)*sigmaX (2)/pi; uEv=-psi(1)*aEv-muX(2); muX1=muX;sigmaX1=sigmaX; x=muX;normX=eps; while abs(norm(x)-normX)/normX>1e-6 ??? normX=norm(x); ??? g=3*x(1)-x(2); ??? gX=[3;-1]; ?? ?cdfX=[logncdf(x(1),mLn,sLn);1-evcdf(-x(2),uEv,aEv)]; ??? pdfX=[lognpdf(x(1),mLn,sLn); evcdf(-x(2),uEv,aEv)]; ??? nc=norminv(cdfX); ??? sigmaX1(1:2)=normpdf(nc)/pdfX; ??? muX1(1:2)=[x(1:2)-nc.*sigmaX1(1:2)]; ??? gs=gX.*sigmaX1; ??? alphaX=-gs/norm(gs); ??? bbeta=(g+gX'*(muX1-x))/norm(gs) ??? x=muX1+bbeta*sigmaX1.*alphaX end 发自小木虫Android客户端 |
回复此楼» 猜你喜欢
刚博士毕业，入职普通一本高校
已经有21人回复
-
博士白读了
已经有50人回复
-
with editor 两个月了，什么原因？
已经有16人回复
-
普通院校药学硕士，做合成的，感觉找不到工作
已经有17人回复
-
如何分离原料与产品
已经有13人回复
-
LWT投
已经有3人回复
-
真的好想读博！
已经有22人回复
-
基金开始函评了吗？
已经有16人回复
-
没学上
已经有12人回复
-
研二光催化6月底四篇二区什么水平
已经有15人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16178.8
- 红花: 52
- 帖子: 2570
- 在线: 572.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
2楼2020-06-07 13:54:10
5