绘制二元正态分布的密度函数图和等高线图MATLAB

  题目:绘制二元正态分布的密度函数图和等高线图,其中,均值u=(0;0),协方差矩阵为[4,3;3,9]

  分析:采用MATLAB中的mesh函数(密度图)和contour函数(等高线)

  二元正态分布密度函数公式:



  代码如下:

%u=[4;3];
%v=[1,-1;-1,2];
u=[0;0];%均值
v=[4,3;3,9];%协方差阵
x=-7:0.05:7;
y=-7:0.05:7;
%x=-2:0.05:10;
%y=-3:0.05:9;
[X,Y]=meshgrid(x,y);
s2x=v(1,1)%x的方差
s2y=v(2,2)
sx=sqrt(s2x)%标准差
sy=sqrt(s2y)
Cov=v(1,2)
r=Cov/(sx*sy)
a=1/(2*pi*sx*sy*sqrt(1-r^2));
b1=-1/(2*(1-r^2));
b2=((X-u(1))./sx).^2;
b3=((Y-u(2))./sy).^2;
b4=2*r.*(X-u(1)).*(Y-u(2))./(sx*sy)
Z=a*exp(b1*(b2+b3-b4));%也就是f(x1,x2)的表达式
%mesh(x,y,Z);
%figure
%grid on
mesh(X,Y,Z),title('密度函数图')
figure
contour(X,Y,Z),title('等高线图')

结果如下:







  • 12
    点赞
  • 70
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值