当和时极坐标是如何分布的?


19

令随机点的笛卡尔x,y坐标为st (x,y)U(10,10)×U(10,10)

因此,半径ρ=x2+y2并不是ρpdf所暗示的均匀分布。

尽管如此,我希望θ=arctanyx几乎是均匀的,不包括由于边缘4个残差而导致的假象:

在此处输入图片说明

以下是grafically计算概率密度函数θρ在此处输入图片说明

现在,如果我让分布为st那么似乎是均匀分布的:x,yx,yN(0,202)×N(0,202)θ

在此处输入图片说明

为什么当时不均匀而当时是均匀的吗?θ(x,y)U(10,10)×U(10,10)x,yN(0,202)×N(0,202)

我使用的Matlab代码:

number_of_points = 100000;
rng('shuffle')

a = -10;
b = 10;
r = (b-a).*randn(2,number_of_points);
r = reshape(r, [2,number_of_points]);
I = eye(2);
e1 = I(:,1); e2 = I(:,2);
theta = inf*ones(1,number_of_points);
rho = inf*ones(1,number_of_points);

for i=1:length(r(1,:))
    x = r(:,i);
    [theta(i),rho(i)] = cart2pol(x(1),x(2));        
end

figure
M=3;N=1; bins = 360;
subplot(M,N,1); 
histogram(rad2deg(theta), bins)
title('Polar angle coordinate p.d.f');

subplot(M,N,2); 
histogram(rho, bins);
title('Polar radius coordinate p.d.f');

subplot(M,N,3); 
histogram(r(:));
title('The x-y cooridnates distrbution (p.d.f)');

替换第三行:r = (b-a).*randn(2,number_of_points);with r = (b-a).*randn(2,number_of_points) +a ;将把的分布从正态更改为均匀。(x,y)


5
该问题在每次编辑时看上去都越来越漂亮,并且标题更清晰,更简洁。@ 0x90做得好。
Michael R. Chernick

3
+1。有趣的是,正态分布是唯一导致角度均匀分布(即旋转对称2D分布)的分布,请参阅stats.stackexchange.com/a/255417/28666
变形虫说莫妮卡(Monica)恢复职权

Answers:


13

您指的是从一对独立变量到极坐标表示(半径和角度)的转换,然后查看的边际分布。(X,Y)(R,θ)θ

我将提供一些直观的解释(尽管密度的数学推导本质上是我非正式描述的内容)。

请注意,如果按某个通用比例缩放两个变量X和Y(例如,从U(-1,1)到U(-10,10)或从N(0,1)到N(0,20)同时作用在两个变量上)不会影响角度的分布(它只会影响半径分布的比例)。因此,让我们仅考虑单位案例。

首先考虑统一案例的情况。请注意,分布在单位平方上是均匀的,因此内包含的区域中的概率密度与该区域的面积成比例。具体来说,查看与角度元素在水平附近(近角)和对角线上(在近角)相关的密度:[1,1]2dθθ=0θ=π/4

在此处输入图片说明

显然,当角度接近对角线之一时,与角度()元素相对应的概率元素(即面积)更大。确实考虑在正方形内刻一个圆圈。圆内一个给定的微小角度所覆盖的面积是恒定的,然后,当我们接近对角线时,圆外的部分会增大,即最大。dfθdθ

这完全说明了您在模拟中看到的模式。

实际上,我们可以看到密度必须与从正方形的中心到其边缘的线段的长度成比例;简单的三角法足以从那里导出密度,然后很容易找到使密度积分到1所需的常数。

[编辑:添加了下一点以讨论半径,因为自我的原始答案以来问题已经改变。]

注意,如果我们有超过单位圆(即,我们在以前的正方形内切的一个),则半径的密度,这将是与半径成比例的均匀分布(考虑宽度的小的环形元件的面积在半径即在和 -具有与成正比的面积。然后,当我们经过圆外时,具有较大半径的新环形区域仅从正方形中的零件获得密度贡献,因此密度在到之间减小(最初相当快,然后更慢)。(再次,相当简单的几何概念足以获得密度的功能形式,如果需要的话。)drrrr+drr12


相反,如果关节分布关于原点旋转对称,则概率元素在某个角度不依赖于角度(这实质上是重言式!)。两个独立标准高斯的双变量分布关于原点旋转对称:

在此处输入图片说明

(此图片的代码基于此处的 Elan Cohen的代码此处有一个不错的选择 ,而此处两者之间也有区别

因此,该体积中含有一些角度是每个相同,所以与相关联角度的密度上均匀。dθθ[0,2π)

[通常用于对实线上的法线密度进行积分的极点技巧可用于确定半径平方的密度为负指数,并且从此处可以通过简单的变换自变量轻松确定半径的密度。分布函数]


4
分布中的四个峰值确实是由于正方形四个角。请注意,任何球对称分布都会导致上的均匀分布,从以为中心的球和圆上的均匀性开始。- 10 10 2 θ 0 0 θ(10,10)2θ(0,0)
西安

2
+1。有趣的是,正态分布是唯一导致旋转对称2D分布的分布,请参阅stats.stackexchange.com/a/255417/28666。这让我感到惊讶。
变形虫说莫妮卡(

3
@amoeba是的,这是唯一的圆形对称分布,是独立边距的乘积。
Glen_b-恢复莫妮卡(Monica)

2
我认为这真是太神奇了。考虑在答案中提及它!
变形虫说莫妮卡(Monica)恢复职权

6

我将回答有关导致统一分配正常情况的问题 众所周知,如果Y独立且正态分布,则恒定概率密度的轮廓在x - y平面上是一个圆。半径R = XYxy具有瑞利分布。为了对此进行很好的讨论,维基百科文章标题为“瑞利分布”。R=X2+Y2

现在让我们使用极坐标看随机变量Y。XY

Y = r sin θ 。注意, X 2 + Y 2 = r 2。如果 θ是均匀的上0 2 π - [R具有瑞利分布 X ÿ将独立法线每个具有 0均值和方差常见。反之亦然。相反的证明是我认为OP希望作为问题第二部分的答案。X=rcos(θ)Y=rsin(θ)X2+Y2=r2θ(0,2π)rXY0

这是证明的草图。不失一般性,我们可以假设,是分布式Ñ 0 1 ÿ分布Ñ 0 1 和相互独立的。XN(0,1)YN(0,1)

然后将联合密度。使用对极坐标的变换获得g r θ 。由于x = r sin θ y = r cos θ 。所以,[Rf(x,y)=(1/2π)exp[([x2+y2])/2]g(r,θ)x=rsin(θ)y=rcos(θ)θ=反正切x/y。计算变换的雅可比行列式,并将其适当代入fxy。作为结果[R θ将是- [REXP[-[R2/2π][R00θr=x2+y2θ=arctan(x/y)f(x,y)g(r,θ)rexp[(r2)/(2π)]r0。由此可见, [R和θ是独立的 - [R具有瑞利分布和θ具有恒定密度 1 /2 π 0θ2πrr1/(2π)


这就是说,如果您在距中心(本例中为原点)的固定径向距离处查看双变量密度的高度,则该圆上所有点的值将相同。
Michael R. Chernick


@ 0x90是,您的链接显示了一种查看此现象的方法,即查看密度指数的二次形式。因此,通常对于双变量正态设置,指数等于常数定义了常数密度的轮廓,并且该方程式是椭圆之一。在特殊情况下,当协方差矩阵是可缩放的单位矩阵时,椭圆简化为一个圆。
Michael R. Chernick

2
我觉得有看到均匀性更简单的方法:独立正常的均值0,很容易证明他们比柯西0 1 。由于柯西的CDF只是缩放和转换arctan,因此通过概率积分变换arctan X / Y 只是移位和缩放的标准均匀随机变量。X,Y0Cauchy(0,1)arctan arctan(X/Y)
弗朗西斯

1
@Francis大多数情况下,我很感谢您对我所有方程式的详尽编辑。我还想说,您在上面的评论绝对显示了一种解决theta均匀性问题的富有想象力的方法。我敢肯定,有些人会同意这样做会更容易。
Michael R. Chernick

6

为了完成格伦和迈克尔给出的还算不错的答案,我只是计算的密度时的分布X Y ^ 是均匀的广场上[ - 1 1 ] × [ - 1 1 ]。该均匀密度为1θ(X,Y)[1,1]×[1,1]在这个正方形上为 4,在其他地方为0,也就是说,在正方形的给定区域中对一个点进行采样的概率为1140该区域的面积。14

这个问题的关注区域是此图上的红色部分: 带阴影部分的正方形

θθ+dθθθ+dθθ

θ[π4,π4]π2

1cosθ

1cos(θ+dθ)=1cosθ+sinθcos2θdθ.

abα12absinα

12(1cosθ)(1cosθ+sinθcos2θdθ)sindθ=dθ2cos2θ
dθsindθ=dθ

θ

18cos2θ
θ[π4,π4]π2

验证:

x <- runif(1e6, -1, 1)
y <- runif(1e6, -1, 1)
hist( atan2(y,x), freq=FALSE, breaks=100)
theta <- seq(-pi, pi, length=500)
lines(theta, 0.125/cos((theta + pi/4)%%(pi/2) - pi/4)**2, col="red" )

直方图+密度

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.