实验5蒙特卡罗法求PI排序算法解答.doc

浙江大学城市学院实验报告 课程名称 多核与并行程序设计 实验项目名称 实验五 蒙特卡罗法求PI及排序算法 学生姓名 专业班级 学号 实验成绩 指导老师（签名 ） 日期 【实验环境】 硬件平台：联想4核,4GZ内存 编译器：Microsoft Visual Studio C++ 6.0 操作系统：Windows 2003 server sp2 测试数据集合：由随机数函数产生的数据集合 【实验1】 一、问题描述 蒙特卡洛算法可理解为通过大量实验，模拟实际行为，来收集统计数据。本例中，算法随机产生一系列点，模拟这些点落在如下图所示的正方形区域内的情况。其几何解释如下 图1 如图1所示，正方形边长为1，左下顶点与原点重合，两边分别与x，y轴重合。曲线为1/4圆弧，圆心位于原点，与正方形左下定点重合，半径为1。正方形面积S1=1，圆弧内面积S2=。算法模拟大量点随机落在此正方形区域内，落在圆弧内的点的数量（n2）与点的总数（n1）的比例与面积成正比关系。即 （1） 由此可得 （2） 因此，只要计算出落在圆弧内的点的数量在点总数中所占的比例，就能求出的值。 由图1可知，所有点均落在正方形范围内，因此点的x坐标满足。又，当点落在圆弧范围内，则点的二维坐标关系满足。检验每一个点是否满足此关系即可判定改点是否落在圆弧内。 二、串行算法描述 本项目中使用了标准C语言库中的产生随机数函数。该函数原型为： int rand( void ); 此函数产生随机数列，每次调用时均返回0到RAND_MAX之间的一个整数。 void srand( unsigned int seed ); 此函数为rand（）函数所生成的伪随机数序列设置起始点，使之产生不同的伪随机数。 算法： 产生2n个随机数据，范围[0，1]，对每个数据点计算其坐标是否满足，统计满足此关系的点的数量count，则 三、并行算法 3.1 并行算法描述 算法步骤： 1、确定需要产生的点的个数n，参与运行的处理器数m； 2、对每一个处理器，生成两个随机数x，y，范围[0，1]； 3、判断两个随机数x，y是否满足； 4、若满足，则变量COUNTi++； 5、重复步骤2-4，直至每个处理器均生成n/m个随机点； 6、收集COUNTi的值，并累加至变量COUNT中，此即为随机点落在圆弧内的数量； 7、通过（2）式计算的值。 3.2 并行算法在Windows下的一个例子 #includestdio.h #includewindows.h #include time.h //#include process.h #include iostream #include fstream #includestdlib.h using namespace std; HANDLE evFinish; long cs=0; //总循环次数 long count=0; //主线程有效次数 long count_thread=0; //thread线程有效次数 time_t start, finish; //定义开始结束时间 //thread线程计算量为总数的一半 DWORD WINAPI thread(LPVOID param) //这个函数在PDF里 WaitForSingleObject(evFinish,INFINITE);//两线程同步 count+=count_thread; finish=time(NULL); //记录结束时间 printf(并行情况：\n\n); printf(用时=%f 秒\n,difftime(finish,start)); //计算时间差 printf(总共的循环次数=%d次\n,cs); printf( 线程有效次数=%d次\n,count); printf(pi= %f \n,4*(double)count/(double)cs); printf(串行行情况：\n); count=0; start=time(NULL); //记录开始时间 for(i=0;ics;i++) { x=(long double)ra