y=sinx的概率分布推导

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

B417科研笔记

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

阅读量4.7k

点赞数 9

分类专栏：通信中的常用数学

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

版权

概率分布均匀分布正弦函数 MATLAB 概率密度函数

关键词由CSDN通过智能技术生成

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

48 篇文章

订阅专栏

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

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

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

在这里插入图片描述

结果如图所示，也就是说， $\sin(x)$ 并不是如直觉所想，在 $[- 1, 1]$ 上均匀分布！
注意到由于 $\sin x = -\sin(x + \pi)$ ，因此 $\sin(x)在[0,\pi]$ 上的分布与 $\sin(x)在[\pi,2\pi]$ 应该刚好对称。
那么我们进行推导, 当 $x$ 在 $\pi]$ 上均匀分布时：
$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] \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\ge\pi-\sin^{-1}y)=1 - P(X\le \pi-\sin^{-1}y)=\frac{\sin ^{-1} y}{\pi}$
因此 $F_Y(y) = \frac{2\sin ^{-1} y}{\pi}$ .
当 $x$ 在 $[\pi, 2\pi]$ 上均匀分布时：
$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\le\pi-\sin^{-1}y)=\frac{-\sin ^{-1} y}{\pi}$ ,
因此 $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)

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