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

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 cdfx,求助大神相关推荐

  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 lognpdf用法,求助大神

    该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 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. 施一公:西湖大学评审终身副教授不看文章数,而是看你能不能讲好一个故事...
  2. linux chown 命令 更改设置文件对应的用户和组
  3. python编程员工资-Python 程序员的工资能有多高?
  4. Eclipse,myeclipse开发中常用技巧总结
  5. matlab中max和find函数详解
  6. [蓝桥杯][算法提高]能量项链(区间dp)
  7. html插入图片和文字,HTML第三课文字图片插入
  8. git 回退上一个版本
  9. python程序加密_Python简单的加密程序:如何循环Z回到A
  10. SSH config 文件的作用
  11. DBeaver连接GBase数据库
  12. SLT2021: CONVERSATIONAL END-TO-END TTS FOR VOICE AGENTS
  13. AG9311/AG9310 Type-C转HDMI设计方案|替代AG9310/AG9311芯片|GSV2201可完全替代兼容AG9310/AG9311
  14. 前端好还是后端好,看看7年前端和后端怎么说
  15. ActiveMQ双机热备时的连接跳动问题解决
  16. 【LeetCode】935. Knight Dialer 解题报告(Python)
  17. JCJC错别字检测系统-年底收官
  18. KVM环境初搭建【工作笔记】
  19. 如何解决css动画抖动的问题
  20. 快速给视频添加背景和水印的步骤

热门文章

  1. 北理工嵩天Python学习笔记
  2. python:计划持有基金n年,求n年的每年复利_利率
  3. 怎么让抖音快速涨粉之抖音数据分析
  4. python格式和JSON格式转换
  5. 罗技 Logitech flow 连接
  6. Unity --- 导航网格 与 导航的使用
  7. 【Visual C++】游戏开发笔记之九 游戏地图制作(一)平面地图贴图
  8. Python3 juypterlab pyecharts 饼图制作
  9. 计算机面试工作计划,信息技术部面试自我介绍
  10. DDLMS-DFE算法