MATLAB模拟布丰投针实验

标签(空格分隔)： 算法

Buffon's Needle

桌面上有距离为a的若干平行线，将长度为L的针随机丢在桌面上，则这根针与平行线相交的概率是多少？假定L < a.

思路：从针据横线的距离与夹角得出。

解决：

假设针的中点到最近横线的距离为y，则\(y\in[0,\frac{a}{2}]\);

因为投针是随机的，所以y服从均匀分布：

\[ f(y) =

\begin{cases}

\frac{2}{a}, & \text{$0 \leq y \leq \frac{a}{2}$} \\

0, & \text{others}

\end{cases}

\]

假定横线向右为正向，记投针与横线正向的角为\(\theta\)，则\(\theta \in[0, \pi]\),为均匀分布。

\[f(\theta) =

\begin{cases}

\frac{1}{\pi}, & \text{$0 \leq \theta \leq \pi$} \\

0, & \text{others}

\end{cases}

\]

投针与横线有交点，即\(y \leq \frac{L}{2}sin\theta\)

布丰投针估算\(\pi\) -- 蒙特卡罗模拟

针与横线有交点的概率：

\(P(x) = \int_{0}^{\pi}\int_{0}^{\frac{L}{2}sin\theta}f(y,\theta)dyd\theta = \int_{0}^{\pi}\int_{0}^{\frac{L}{2}sin\theta}f(y)f(\theta)dyd\theta \\

= \int_{0}^{\pi}\int_{0}^{\frac{L}{2}sin\theta}\frac{2}{a} * \frac{1}{\pi}dyd\theta = \frac{2L}{a\pi}\)

如果做n次投针实验，其中有k次针与横线相交，则针与横线相交的频率为：\(\frac{k}{n}\),根据大数定理，频率也就为概率。

$ \frac{2L}{a\pi} \approx \frac{k}{n}$ 即， \(\pi \approx \frac{2Ln}{ak}\)

MATLAB模拟实验

用布丰投针实验近似计算\(pi\)的值：

代码如下：

l = 0.6; %针的长度

a = 1; %平行线的宽度

n = 1000000; %做n次投针试验

k = 0; %记录针与平行线相交的次数

y = unifrnd(0, a/2, 1, n); %在[0, a/2]内服从均匀分布随机产生n个数

theta = unifrnd(0, pi, 1, n); %在[0, pi]内服从均匀分布随机产生n个数

for i=1:n

if y(i) < (l/2)*sin(theta(i))

k = k + 1;

end

end

f = k / n;

Pi = (2*l*n)/(a*k);

结果如图所示：

如此进行多次实验，进行估计。

如图为进行100次重复投针实验，每次投针1000000次，结果如图所示：

多普勒失真信号采样Matlab模拟分析

多普勒失真信号采样Matlab模拟分析 方案 水声通信指的是使用声信号在水中数据传输. 相对而言.电磁信号在水中吸收严重衰减过快,光信号受水中悬浮颗粒的影响,也无法完毕远距离传输. 这两种信号的传播距 ...

模拟生产搭建Standby RAC实验环境&lpar;11&period;2&period;0&period;4 DG&rpar;

模拟生产搭建Standby RAC实验环境(11.2.0.4 DG) 环境:RHEL 6.5 + Oracle 11.2.0.4 GI.DB 1.需求背景介绍 2.准备工作 3.主库配置 4.备库配置 ...

csdn vip文章：使用matlab模拟镜头失真