注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

Mr.Right

不顾一切的去想,于是我们有了梦想。脚踏实地的去做,于是梦想成了现实。

 
 
 

日志

 
 
关于我

人生一年又一年,只要每年都有所积累,有所成长,都有那么一次自己认为满意的花开时刻就好。即使一时不顺,也要敞开胸怀。生命的荣枯并不是简单的重复,一时的得失不是成败的尺度。花开不是荣耀,而是一个美丽的结束,花谢也不是耻辱,而是一个低调的开始。

网易考拉推荐

2维高斯随机变量的pdf -- matlab  

2012-08-27 11:38:33|  分类: 绘图 |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
调用:gauss_2d_pdf([1;1], 0.9)


function [C,p]=gauss_2d_pdf(sigmas,rho)
%
% USAGE: p=gauss_2d_pdf(sigmas,rho);
%
%  Inputs: sigmas = 1x2 row vector of the variances of X and Y
%          rho = correlation coefficient (between -1 and 1) of X and Y
%
% Outputs:    C = covariance matrix
%             p = 2-D pdf
%             Surface Plot of PDF
% Example: gauss_2d_pdf([1;1], 0.9)

if abs(rho)>1
    error('rho must be between -1 and 1')
end

MAX=max(sigmas);
x=3*MAX*(-100:100)/100;
y=x;

cov_xy=rho*sqrt(sigmas(1)*sigmas(2));

C=[sigmas(1) cov_xy;cov_xy sigmas(2)];


Const=1/(2*pi*sqrt(det(C)));

[X,Y]=meshgrid(x,y);

C_inv=inv(C);
quad_form = C_inv(1,1)*(X.^2) + C_inv(2,2)*(Y.^2)  + 2*C_inv(1,2)*(X.*Y);


p=Const*exp(-0.5*quad_form);

subplot(2,1,1)
surfl(x,y,p); shading interp;
subplot(2,1,2)
contour(x,y,p)
axis square



  评论这张
 
阅读(843)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2016