二维正态随机变量概率密度函数三维图的matlab实现
1.二维正态随机变量
二维正态随机变量是最常见的一种二维随机变量分布。其联合概率密度函数为:
p
(
x
,
y
)
=
1
2
π
σ
X
σ
Y
1
−
r
2
⋅
e
x
p
{
−
1
2
(
1
−
r
2
)
[
(
x
−
m
X
2
)
σ
X
2
−
2
r
(
x
−
m
X
)
(
y
−
m
Y
)
σ
X
σ
Y
+
(
y
−
m
Y
2
)
σ
Y
2
]
}
p(x,y)=\frac{1}{2\pi \sigma _X\sigma _Y\sqrt{1-r^2}}\cdot exp\{ -\frac{1}{2(1-r^2)}[\frac{(x-m_X^2)}{\sigma_X ^2}-\frac{2r(x-m_X)(y-m_Y)}{\sigma_X\sigma_Y}+\frac{(y-m_Y^2)}{\sigma_Y^2}]\}
p(x,y)=2πσXσY1−r21⋅exp{−2(1−r2)1[σX2(x−mX2)−σXσY2r(x−mX)(y−mY)+σY2(y−mY2)]}
变量 | 含义 |
---|---|
σ X \sigma_X σX | 随机变量X的方差 |
σ Y \sigma_Y σY | 随机变量Y的方差 |
m X m_X mX | 随机变量X的方差 |
m Y m_Y mY | 随机变量Y的方差 |
r | 随机变量X、Y相关系数 |
2.Mtalab画联合概率密度三维图
取 σ X = σ Y = 1 , m X = m Y = 5 , r = 0 \sigma_X=\sigma_Y=1,m_X=m_Y=5,r=0 σX=σY=1,mX=mY=5,r=0,画联合概率密度的三维曲面如下:
- 三维视图
- X-Z视图:
- Y-Z视图:
- 任意视图(体验视觉冲击力)
3.matlab代码
clc
close all
clearvars
Dx=1;%方差
Dy=1;%方差
mx=5;
my=5;
r=0;
x=0:0.05:10;
y=0:0.05:10;
[X,Y]=meshgrid(x,y);
p2=(1/(2*pi*Dx*Dy*sqrt(1-r^2)))*exp((-1/(2*(1-r^2)))*((X-mx).^2/Dx^2)-(2*r*(X-mx).*(Y-my)/(Dx*Dy)+(Y-my).^2/Dy^2));
mesh(X,Y,p2)
title('随机变量X、Y的联合概率密度')
xlabel('X')
ylabel('Y')
zlabel('联合概率密度')