该楼层疑似违规已被系统折叠 隐藏此楼查看此楼

muX=[20;10];cvX=[0.25;0.2];

rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX;

aEv=sqrt(6)*sigmaX(1)/pi;uEv=-psi(1)*aEv-muX(1);

sLn=sqrt(log(1+(sigmaX(2)/muX(2))^2));mLn=log(muX(2))-sLn^2/2;

muX1=muX;sigmaX1=sigmaX;

s=rhoX(1,2);t=cvX(2);

r=1.029+0.001*s+0.004*s^2+0.014*t+0.233*t^2-0.197*s*t;

r12=r*rhoX(1,2);rhoX1=[1,r12,r12,1];

x=muX;normX=eps;

while abs(norm(x)-normX)/normX > 1e-6

normX=norm(x);

g=x(2)^2-2*x(1);

gX=[-2;2*x(2)];

cdfX=[1-evcdf(-x(1),uEv,aEv);logncdf(x(2),mLn,sLn)];

pdfX=[evpdf(-x(1),uEv,aEv);lognpdf(x(2),mLn,sLn)];

nc=norminv(cdfX);

sigmaX1=normpdf(nc)./pdfX;

muX1=x-nc.*sigmaX1;

covX=diag(sigmaX1)*rhoX1*diag(sigmaX1);

[a,d]=eig(covX);

muY=a'*muX1;sigmaY=sqrt(diag(d));

y=a'*x;

gY=a'*gX;gs=gY.*sigmaY;

alphaY=-gs/norm(gs);

bbeta=(g+gY'*(muY-y))/norm(gY.*sigmaY)

y=muY+bbeta*sigmaY.*alphaY;

x=a*y;

end

错误使用 *

内部矩阵维度必须一致。

出错 Untitled (line 20)

covX=diag(sigmaX1)*rhoX1*diag(sigmaX1);

把*改成.*还是说错误

有没有大神看看哪的错

matlab lognpdf用法,求助大神相关推荐

  1. vs 2017 建立 php,程序在VS 2017中输出的结果为什么是这样的?求助大神!谢谢!...

    编程论坛 → C语言论坛 程序在VS 2017中输出的结果为什么是这样的?求助大神!谢谢! 菜鸟100学C 发布于 2018-09-29 08:06, 1405 次点击 #include int fu ...

  2. 求助大神:valgrind 出现 address 0x0 is not stack'd or (recently)freee'd怎么解决?linuxclionvalgrind## 标题

    求助大神:valgrind 出现 address 0x0 is not stack'd or (recently)freee'd怎么解决?linux&clion&valgrind## ...

  3. Mysword我的剑圣经破解方法,,,求助大神。

    Mysword我的剑圣经破解方法,,,求助大神. 求大神帮助我给些思路方法破解这款软件,我不知道该从哪里下手,求大家帮助一下,会有一些辛苦费谢谢大家的帮助!!!!2256771988我的qq

  4. matlab 12位 显示不出来,求助大神,为何不同机器运行MATLAB结果不同

    求助:不同机器运行MATLAB结果不同 我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结 ...

  5. 不同cpu matlab,求助大神,为何不同机器运行MATLAB结果不同

    求助:不同机器运行MATLAB结果不同 我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结 ...

  6. evcdf matlab,求助大神

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 muX=[20;10];cvX=[0.25;0.2]; rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX; aEv=sqrt(6)*s ...

  7. matlab cdfx,求助大神

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 muX=[20;10];cvX=[0.25;0.2]; rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX; aEv=sqrt(6)*s ...

  8. matlab模拟出现较大误差是什么原因,求助大神,使用BP神经网络预测数据,为什么误差很大?...

    随机产生了25组数据,利用公式求得结果矩阵,现用20组数据训练BP网络,剩下的5组进行预测,产生的结果误差为什么很大,请大神指教,不甚感激 由于金币只有1.5个,所以没办法发求助帖 %测试BP神经网络 ...

  9. 求助大神:用C#编辑计算器,怎么实现连加、连除功能啊!!!!!

    基础的计算器会做,但是没有连加.连减.连除.连乘等功能就是一个BUG呀!求助哪位大神,怎么实现连加功能啊,谢谢

最新文章

  1. GDAL2.1.1库在Ubuntu14.04下编译时遇到的问题处理方法
  2. (转)MySQL 线程池内幕
  3. Servlet学习笔记(三)之HttpServletRequest
  4. 升级版记事本 Notepad++
  5. FSD键盘钩子框架参考爱写驱动的女装大佬
  6. socket下server端支持多客户端并发访问简单实现
  7. 排序算法之——插入排序
  8. jQuery解决ajax请求的跨域问题
  9. 噩耗!吴孟达去世,终年68岁!
  10. D3D中设备丢失的处理
  11. 收集的关于DLL较全的介绍
  12. 微带贴片天线的分析方法-传输线模型法
  13. Android Accessibility大致解析,通过adb运行纯java代码打开应用的辅助功能
  14. 《深度学习》学习笔记(十二)
  15. mysql错误合集_mysql错误合集
  16. FileNotFoundError: Could not find module ‘C:\Users\Yolanda\anaconda3\envs\python38\lib\site-packages
  17. IDEA右下角git图标上消失,如何重现?
  18. python练习-给数组去重后元素的加序号
  19. Struts2的运行流程
  20. 【Maven】Maven的安装与下载

热门文章

  1. 【虚拟机学习】常用的虚拟化软件
  2. Nginx反向代理静态资源
  3. 十大交易系统源码 (http://bbs.tb18.net/forum.php?mod=viewthreadaction=printabletid=22501)
  4. ADS系列 – 低噪声放大器(LNA)模型下载安装及 LNA仿真设计
  5. 如何在本地进行微信公众号的开发和调试
  6. 服务器负载均衡 链路负载均衡 《CDN技术详解》
  7. python3 zip破解
  8. 数字电路使用555定时器,74163,74194设计跑马灯电路
  9. 分析师分析业务维度,附(个人制作分析思维导图Xmind)
  10. Mysql 分库分表