y=sinx的概率分布推导

最新推荐文章于 2024-07-29 22:29:57 发布

B417科研笔记

最新推荐文章于 2024-07-29 22:29:57 发布

阅读量4.9k

点赞数 11

CC 4.0 BY-SA版权

分类专栏：通信中的常用数学文章标签：线性代数概率论 matlab

本文链接：https://blog.csdn.net/weixin_39274659/article/details/121305911

通信中的常用数学专栏收录该内容

48 篇文章

订阅专栏

博客探讨了当角度x在[0,2π]上均匀分布时，正弦函数sin(x)的分布并非均匀。通过MATLAB模拟和推导，得出sin(x)在[0,π]和[π,2π]两个区间上的概率分布，并给出精确的概率密度函数公式。最终通过程序验证了推导的正确性。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

近期碰到这样一个小问题，当 $x$ 是在 $2\pi]$ 上均匀分布时，求 $y=sin⁡(x)y=\sin(x)$ 的分布？
然后发现百度出来的答案都是错的…

先上结论，用matlab先看一下概率分布图，代码如下：

for i = 1 : 100000
    a = unifrnd(0, 2 * pi);
    b(i) = sin(a);
end
cdfplot(b)

在这里插入图片描述

结果如图所示，也就是说， $sin⁡(x)\sin(x)$ 并不是如直觉所想，在 $[- 1, 1]$ 上均匀分布！
注意到由于 $sin⁡x=−sin⁡(x+π)\sin x = -\sin(x + \pi)$ ，因此 $sin⁡(x)在[0,π]\sin(x)在[0,\pi]$ 上的分布与 $sin⁡(x)在[π,2π]\sin(x)在[\pi,2\pi]$ 应该刚好对称。
那么我们进行推导, 当 $x$ 在 $\pi]$ 上均匀分布时：
$FY(y)=P(Y≤y)=P(sin⁡X≤y)=P(X≤sin⁡−1y)+P(X≥π−sin⁡−1y)F_Y(y)=P(Y\le y)=P(\sin X\le y)=P(X\le\sin^{-1}y)+P(X\ge\pi-\sin^{-1}y)$
这里我们假定, $arcsin⁡:[−1;1]↦[−π/2;π/2]\arcsin :[-1 ; 1] \mapsto[-\pi / 2 ; \pi / 2]$ 。
接着有：
$P(X\le\sin^{-1}y) =\int_{-\infty}^{\sin ^{-1} y} f_{X}(x) d x=\frac{1}{\pi}\int_{-\infty}^{\sin ^{-1} y} d x=\frac{\sin ^{-1} y}{\pi}$
类似的，
$P(X≥π−sin⁡−1y)=1−P(X≤π−sin⁡−1y)=sin⁡−1yπP(X\ge\pi-\sin^{-1}y)=1 - P(X\le \pi-\sin^{-1}y)=\frac{\sin ^{-1} y}{\pi}$
因此 $FY(y)=2sin⁡−1yπF_Y(y) = \frac{2\sin ^{-1} y}{\pi}$ .
当 $x$ 在 $[π,2π][\pi, 2\pi]$ 上均匀分布时：
$FY(y)=P(Y≤y)=P(sin⁡X≤y)=P(π−sin⁡−1y≤X≤2π+sin⁡−1y)F_Y(y)=P(Y\le y)=P(\sin X\le y)=P(\pi-\sin^{-1}y\le X\le2\pi+\sin^{-1}y)$
接着有：
$P(X\le2\pi+\sin^{-1}y) =\int_{-\infty}^{2\pi + \sin ^{-1} y} f_{X}(x) d x=\frac{1}{\pi}\int_{-\infty}^{2\pi + \sin ^{-1} y} d x=\frac{\pi + \sin ^{-1} y}{\pi}$
类似的,
$P(X≤π−sin⁡−1y)=−sin⁡−1yπP(X\le\pi-\sin^{-1}y)=\frac{-\sin ^{-1} y}{\pi}$ ,
因此 $FY(y)=π+2sin⁡−1yπF_Y(y) = \frac{\pi+2\sin ^{-1} y}{\pi}$ .

综上，
$F_Y(y)= \begin{cases} \frac{\sin ^{-1}y}{\pi} + \frac{1}{2}&y\ge0 \\ \frac{\pi+2\sin ^{-1} y}{2\pi} & y<0\end{cases}.$

用matlab程序验证：

y = -1 : 0.001: 1;
for i = 1 : length(y)
    if y(i) >= 0
        z(i) = asin(y(i))/pi + 0.5;
    else
        z(i) = (pi + 2 * asin(y(i))) / 2 / pi;
    end
end
plot(y,z)

最后两条线完美重合，证明正确！
在这里插入图片描述