matlab lognpdf用法,求助大神
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
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用法,求助大神相关推荐
- vs 2017 建立 php,程序在VS 2017中输出的结果为什么是这样的?求助大神!谢谢!...
编程论坛 → C语言论坛 程序在VS 2017中输出的结果为什么是这样的?求助大神!谢谢! 菜鸟100学C 发布于 2018-09-29 08:06, 1405 次点击 #include int fu ...
- 求助大神: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## ...
- Mysword我的剑圣经破解方法,,,求助大神。
Mysword我的剑圣经破解方法,,,求助大神. 求大神帮助我给些思路方法破解这款软件,我不知道该从哪里下手,求大家帮助一下,会有一些辛苦费谢谢大家的帮助!!!!2256771988我的qq
- matlab 12位 显示不出来,求助大神,为何不同机器运行MATLAB结果不同
求助:不同机器运行MATLAB结果不同 我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结 ...
- 不同cpu matlab,求助大神,为何不同机器运行MATLAB结果不同
求助:不同机器运行MATLAB结果不同 我调用MATLAB优化工具箱的库函数fmincon,使用相同的初始解(可行解),对同一个问题进行局部搜索(算法为序列二次规划,即SQP),但在不同机器上得到的结 ...
- evcdf matlab,求助大神
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 muX=[20;10];cvX=[0.25;0.2]; rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX; aEv=sqrt(6)*s ...
- matlab cdfx,求助大神
该楼层疑似违规已被系统折叠 隐藏此楼查看此楼 muX=[20;10];cvX=[0.25;0.2]; rhoX=[1,0.5;0.5,1];sigmaX=cvX.*muX; aEv=sqrt(6)*s ...
- matlab模拟出现较大误差是什么原因,求助大神,使用BP神经网络预测数据,为什么误差很大?...
随机产生了25组数据,利用公式求得结果矩阵,现用20组数据训练BP网络,剩下的5组进行预测,产生的结果误差为什么很大,请大神指教,不甚感激 由于金币只有1.5个,所以没办法发求助帖 %测试BP神经网络 ...
- 求助大神:用C#编辑计算器,怎么实现连加、连除功能啊!!!!!
基础的计算器会做,但是没有连加.连减.连除.连乘等功能就是一个BUG呀!求助哪位大神,怎么实现连加功能啊,谢谢
最新文章
- GDAL2.1.1库在Ubuntu14.04下编译时遇到的问题处理方法
- (转)MySQL 线程池内幕
- Servlet学习笔记(三)之HttpServletRequest
- 升级版记事本 Notepad++
- FSD键盘钩子框架参考爱写驱动的女装大佬
- socket下server端支持多客户端并发访问简单实现
- 排序算法之——插入排序
- jQuery解决ajax请求的跨域问题
- 噩耗!吴孟达去世,终年68岁!
- D3D中设备丢失的处理
- 收集的关于DLL较全的介绍
- 微带贴片天线的分析方法-传输线模型法
- Android Accessibility大致解析,通过adb运行纯java代码打开应用的辅助功能
- 《深度学习》学习笔记(十二)
- mysql错误合集_mysql错误合集
- FileNotFoundError: Could not find module ‘C:\Users\Yolanda\anaconda3\envs\python38\lib\site-packages
- IDEA右下角git图标上消失,如何重现?
- python练习-给数组去重后元素的加序号
- Struts2的运行流程
- 【Maven】Maven的安装与下载
热门文章
- 【虚拟机学习】常用的虚拟化软件
- Nginx反向代理静态资源
- 十大交易系统源码 (http://bbs.tb18.net/forum.php?mod=viewthreadaction=printabletid=22501)
- ADS系列 – 低噪声放大器(LNA)模型下载安装及 LNA仿真设计
- 如何在本地进行微信公众号的开发和调试
- 服务器负载均衡 链路负载均衡 《CDN技术详解》
- python3 zip破解
- 数字电路使用555定时器,74163,74194设计跑马灯电路
- 分析师分析业务维度,附(个人制作分析思维导图Xmind)
- Mysql 分库分表