蒙特卡洛模拟——蒲丰投针问题

蒲丰投针简述：

        蒲丰投针法（Buffon's Needle）背后有一个有趣的故事，与法国数学家蒲丰（Georges-Louis Leclerc, Comte de Buffon）的名字相关。这个故事源于18世纪，是数学领域的一个经典问题，涉及到概率和几何。

        蒲丰于1707年出生在法国，他是一位多才多艺的自然科学家，数学家和哲学家。在他的一本著作《自然历史》（Histoire naturelle）中，他提出了一个有趣的问题，探讨了投针在平行线上的分布。这个问题后来就以他的名字命名为蒲丰投针法。

        蒲丰的问题可以简化为：如果你有一根长度为d的针，在一组间距为t的平行线上随机投掷针，那么针与线相交的概率是多少？这个问题看似简单，但涉及了几何学和概率学的复杂关系。

        蒲丰用这个问题来探讨圆周率 π，他发现通过统计投针与线相交的次数和总投针次数，可以估计出圆周率 π 的值。这个方法成为了一种早期的蒙特卡洛方法示例，因为它依赖于随机抽样。

        蒲丰投针法在数学和概率教育中常被用来介绍概率统计方法，以及如何使用蒙特卡洛模拟来估计数学常数。它代表了一种将数学与实际问题相结合的方法，展示了数学可以用于解决各种实际世界的难题。

蒲丰投针的原理：

1. 几何关系：假设有两组平行线，一组距离为 d （针的长度），另一组距离为 t（两相邻平行线之间的间距，也可以称为“格子宽度”）。我们将针随机投放在这些线之间。

2. 相交情况：当针与平行线相交时，可以分成两种情况：

   • 情况一：针的中心点位于两条相邻线之间，但是有一条线在针的两个端点之间。
   • 情况二：针横跨了两条相邻线，中心点距离其中一条线的距离小于针的一半。

3. 概率计算：如果将针随机投放，并且假设针的长度不大于格子宽度，那么情况一和情况二出现的概率分别为：

   • 情况一概率：P(针中心位于两条线之间) = 1
   • 情况二概率：P(针横跨两条线) = 针的长度 / 格子宽度

4. 计算圆周率：现在，我们可以使用概率的思想来计算圆周率。假设 N 次投针后，有 M 次针与平行线相交（即情况二出现的次数）。那么根据上面的概率计算，我们有：

   • P(情况二出现) = M / N = 针的长度 / 格子宽度
   • 由于 P(情况二出现) = 针横跨两条线的概率，所以可以得到 π 的估计：π ≈ (2 * 针的长度) / 格子宽度

        通过进行足够多次的投针实验，记录情况二出现的次数 M 和总投针次数 N，就可以使用上述公式估计圆周率 π。

蒲丰投针的MATLAB实现：

% 参数设置
numNeedles = 1000000;  % 增加投针次数
needleLength = 0.8;    % 调整针的长度
gridWidth = 2;         % 两相邻平行线之间的间距

% 初始化计数器
countCrossings = 0;

% 模拟投针过程
for i = 1:numNeedles
    % 随机生成针的中心点位置和角度
    center = rand * gridWidth;
    angle = rand * pi;
    
    % 判断针是否与平行线相交
    if center < (needleLength/2) * sin(angle) || (gridWidth - center) < (needleLength/2) * sin(angle)
        countCrossings = countCrossings + 1;
    end
end

% 计算估计的圆周率
estimatedPi = (2 * needleLength * numNeedles) / (gridWidth * countCrossings);

% 显示估计的圆周率和实际圆周率
disp(['估计的圆周率: ', num2str(estimatedPi)]);
disp(['实际的圆周率: ', num2str(pi)]);