贝朗特悖论的解决.doc

《贝朗特悖论的解决.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、理学院SchoolofScience课程设计报告学生姓名：李凡学生学号：所在班级：07数学1所在专业：数学与应用数学指导教师：樊嵘实习场所：青岛理工大学实习时间：第六学期课程设计成绩总评学习态度报告质量使用SAS统计模拟方法解决Bertrand’sparadoxBertand’sparadox是法国数学家Bertrand于1889提出的一个概率悖论：在圆内任作一弦，其长度超过圆内内接正三角形边长的概率是多少？他在提出问题之后，给出了三种不同的解法，得到了三个不同的结果，是为悖论。第一种解法如下：

2、由于弦交圆于两点。我们先固定弦的一个端点。以此端点作一个等边三角形(如图)。显然，只有穿过此三角形内的弦才符合要求。而符合条件的弦的另一端正好占整个圆弧的1/3。并且，不论固定的那个端点在圆上的哪个位置，情况都是一样的。所以结果为1/3。第二种解法如下：由于弦长只和圆心到它的距离有关。所以固定圆内一条半径。当且仅当圆心到它的距离小于1/2才满足条件。并且，不论固定的是哪条半径，情况都是一样的。所以结果为1/2。第三种解法如下；弦被其中点唯一确定（除了圆心）。当且仅当其中点在半径为1/2的圆内时才

3、满足条件。此小圆面积为大圆的1/4。所以结果为1/4。三种看似都有道理的解法却得到了不同的答案，所以被称为悖论。在以前对这问题的分析中，倾向于认为得到三种结果的原因是因为采用了不同的等可能性假定。解法一假定端点在圆上均匀分布。解法二假定半径在圆内均匀分布以及弦的中点在半径上均匀分布。解法三假定弦的中点在圆内均匀分布。先不论他们的假设是否合理，从这个问题的提法来看，问题考察的是圆内的随机弦问题。我们应该从弦的本质定义出发，即圆上任意两点的连线为弦。从这个思路，我们可以使用SAS进行统计模拟，确定问

4、题的答案。具体思路如下：1.先进行1000次试验，每次试验进行1000次模拟，每次模拟从圆上随机取两点，计算距离，记录的次数。如此得到1000个数据，数据集为cs，其中的变量只有一个x。对此数据进行分析，得到其方差与均值，可以求出概率。2.为了得到弦长的分布，我们进行1000次模拟，每次模拟从圆上随机取两点，计算距离并记录。如此得到数据集为strx，其中的变量有三个，分别记录两点的角度参数x，y与两点之间距离d。3.从圆进行推广，得到椭圆内随机弦长的分布，思路同上。4.从得到的结果进行理论分析。

5、数据的得到与数据集的建立：使用matlab编程可以得到模拟需要的数据，在SAS中建立各数据集的程序如下：cs数据集：Datalf.cs;Inputx@@;/*x为1000中满足条件的弦的个数*/Cards;331329346337333338355348319328333341318315349307320327337371353341325348329323332319316341341348340312353330335317318315341334318358346351362353290

6、332368335359346……；strx数据集：Datalf.strx；Inputxyz；/*x，y为随机角度，z为弦长*/Cards；5.11915.69130.564420.797885.73891.24373.97320.612861.9881.74993.43621.49356.01626.06260.……；；strx1数据集：Datalf.strx1；/*椭圆的数据*/Inputxyz；/*x，y为随机角度，z为弦长，长轴为2，短轴为1*/Cards；4.68015.60631.1

7、0261.52430.814280.839511.41412.19920.825351.80385.82763.04310.322413.72382.4873……；对数据的分析与结果解读：对于cs数据集中的数据，我们根据林德贝格-勒维中心极限定理，记xn为第n次试验中，满足弦长平方大于3的弦的个数，则不管xn的分布如何，只要n充分大，就可以用正态分布去逼近。于是我们先对数据进行正态性检验，使用Solutions-Analysis-GuidedDataAnalysis，对数据进行分析，得到下面的结

8、果：图1图2从图2中可以看到数据的均值为333.89，标准偏差为14.7。其中Q1，与Q3分别为四分之一和四分之三分位。P：normal=0.25025为正态性检验的概率值。图1为数据直方图与正态曲线，图3为正态概率检验图，从两个图可以看出来，数据是服从正态分布的。且可以估计其期望为334次，于是可以得到结论，圆内随机弦长度大于圆内内接正三角形边长长度的概率为334/1000=1/3。对于strx数据集中的数据，我们的目的是得到弦长的分布，即绘制其密度函数曲线和分布函数的曲线。首先是对弦长数据的