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

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);

把*改成.*还是说错误

有没有大神看看哪的错

evcdf matlab,求助大神相关推荐

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

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

  2. bg算法 matlab,求助大神,有关BG算法

    程序如下 function [ FLAGS ] = BGSA( x,y,P0,L0 ) % [ FLAGS ] = BGSA( x,y,P0,L0 ) %  x: 序列的x坐标(仅用于绘图,如果不使用 ...

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

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

  4. edge robert matlab,哪位熟悉matlab的大神路过瞄一眼哈

    cxfx(believe truth believe me) UID240430 帖子100 精华积分1755 蛋蛋币1755 枚 威望0 BT积分0 阅读权限60 性别男 在线时间125 小时 注册 ...

  5. 求助大神: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## ...

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

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

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

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

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

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

  9. matlab cdfx,求助大神

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

最新文章

  1. Android — 使用 SharedPreferences 本地保存 key-value 数据
  2. linux 生成o 文件,急,linux为什么生成不了.o驱动程序文件
  3. java学习(129):hashmap的方法
  4. java 最小生成树_图的最小生成树(java实现)
  5. 腾讯开源国内首个H.266/VVC视频播放器
  6. JDK 环境变量配置
  7. 【渝粤教育】国家开放大学2018年秋季 7389-22T劳动与社会保障法 参考试题
  8. 如何拯救「大泥球」架构?人人车Android客户端架构演进实录
  9. html 双 直线 公文,公文格式标准页面设置
  10. SVN下载安装汉化教程
  11. 良心安利游戏音效素材网站
  12. 【Python游戏】Python实现一个推箱子小游戏 | 附带源码
  13. java-多线程编程
  14. 【计算机组成原理】定点乘法运算之补码一位乘法(Booth算法)(对初学者的步骤详解)
  15. 极虎病毒结束度假,感染量再攀高峰
  16. CSS3炫酷的发光文字 可自定义文字色彩
  17. python处理点云数据_python将指定点云文件(asc)转换为PCD格式
  18. IT、TT、TN系统,你真的了解吗?
  19. 汇编语言中xor指令_这个汇编代码有什么作用? (TEST,XOR,JNZ)
  20. HALCON图像格式

热门文章

  1. 网页游戏架设_这10年来手机游戏的迭代,也是一部硬件发展史丨触乐
  2. java 克隆_Java实现对象克隆的方法
  3. mysql sql variant_SQL 数据类型
  4. python unix时间戳格式化输出_python正常时间和unix时间戳相互转换的方法
  5. 怎样设置html的元素属性,Jquery怎么设置元素的属性和样式?
  6. MySQL面试三连杀:如何实现可重复读、又为什么会出现幻读、是否解决了幻读问题?...
  7. 皮一皮:所以说...哪里的程序猿最能吃辣?
  8. Spring Boot 2.x基础教程:多文件的上传
  9. 面试:知道 CopyOnWriteArrayList 吗?
  10. 如果你不习惯新版的 Github 的 UI 界面,可以试试这款插件