圆周率的几种C语言编程代码

1、-. z.计算的方法一、蒙特卡罗法Monte Carlo蒙特卡洛方法是一种以概率统计理论为指导的一类非常重要的数值计算方法。蒙特卡罗是摩纳哥公国的一座城市赌城。Buffon投针就是蒙特卡洛方法的思想，但是Buffon是蒙特卡洛方法的起源。蒙特卡洛方法同样可以通过随机通过产生随机数的方式来解决计算问题。比方估算和计算不规则图形的面积。图1 蒙特卡罗方法模拟这种方法是一种利用计算机随机数的功能基于随机数的算法，通过计算落在单位圆内的点与落在正方形内的点的比值求PI。图2 蒙特卡罗方法由于图形的对称性，我们靠考虑该图的四分之一局部。假定一点能够均匀地扔到一个正方形中，计算落入其中的点个数。通过计数其

2、中落入内切圆的点的个数；如果一共投入N个点，其中有M个落入圆中，则只要点均匀，假定圆周的半径为R，则：式1该方法得到的要得到的精度与投入点的个数有关，一般个数较大时精度比拟高。/Monte Carlo Methoddoublerand_pi(intn)intnumInCircle=0;double*,y;doublepi;for(inti=0;in;i+)*= rand()*1.0/RAND_MA*;y=rand()*1.0/RAND_MA*;if(*+y*y1)numInCircle+;pi=(4.0*numInCircle)/n;returnpi;图3 蒙特卡罗例如代码写出来你的代码实际测

3、试结果10 Marks表1 蒙特卡罗方法实际测试结果次数10102103104105106107108109计算值运行时间ms误差二、数学公式级数由数学公式：式2当n取1000时就很接近真实值了n=1000时，pi=3.59946。/Numerical Series/Sigma(1/n2) = pi/6doublemath_Pi(intn)intnumInCircle=0;doublesum=0;doublepi;for(inti=1;i=n;i+)sum+=1.0/(double)i*i);pi=sqrt(sum*6);returnpi;图4 级数方法例如代码写出来你的代码实际测试结果10

4、Marks表2 级数方法实际测试结果次数10102103104105106107108109计算值运行时间ms误差三、划分网格计算图5 蒙特卡罗例如代码将图片分为n*n个小方形，统计落在圆内的个数占所有方形的比例。/Grid Countdoublegrid_Pi(intn)inti;doublesum=0;for(i=0;in;i+)sum+=(int)sqrt(n*(double)n-i*(double)i);return(4.0*sum)/n/n;图6 网格划分例如代码写出来你的代码实际测试结果10 Marks表3 网格划分方法实际测试结果次数10102103104105106107108

5、109计算值运行时间ms误差Appendi* 时间统计方法一：clock统计在待测试程序的开场局部和完毕局部分别参加clock统计，两个之间的差值就是该局部代码运行时间。需要注意的是Windows下和Linu*下头文件不一样。clock_tstart;clock_tend;/get current clock countstart=clock();/Insert your code here/get current clock countend=clock();printf(The time duration is %d ms.n,end-start);方法二：系统API准确统计时间Windo

6、ws 和Linu*下均有准确计时方法。在Windows下可以用系统提供的API函数QueryPerformanceFrequency 和QueryPerformanceCounter来进展高精度的计时，现在的机器根本上都提供这种高精度的计时啦，所以不用担忧。利用该函数可以准确的计时到微妙级别。QueryPerformanceFrequency() 可以得到CPU的时钟频率。QueryPerformanceCounter() 可以通过两次的差值来得到CPU的时钟周期差值。这样就可以计算出两次差值之间所花费的时间，当然这里计算出来的时间是秒，一般都换算为毫秒或者微秒来表示。Windows#ifde

7、f_WIN32LARGE_INTEGERnum;longlongstart,end,freq,time1,time2;#endif/开场计时#ifdef_WIN32QueryPerformanceFrequency(&num);freq=num.QuadPart;QueryPerformanceCounter(&num);start=num.QuadPart;#endif/Insert your code here/插入待测时间代码/ Sleep(1000);/完毕计时时间#ifdef_WIN32QueryPerformanceCounter(&num);end=num.QuadPart;ti

8、me1=(end-start)*1000/freq;printf(nDuration time=%d msn,(end-start)*1000/freq);#endifgettimeofday 准确到微秒级别计时函数int gettimeofday(structtimeval *tv,struct timezone *tz);说明：gettimeofday将时间保存在构造tv之中，strut timeval gives the number of seconds and microsecondssince the Epoch，(00:00:00 UTC, January 1, 1970), m

9、easured inmicroseconds.tz(即timezone，很少用到)一般我们使用null来代替。保存时间的构造体strut timeval long tv_sec;long tv_usec;注：需要头文件#include Linu*#ifdef linu*struct timeval ta, tb;double mseca, msecb;#endif/开场计时#ifdef linu*gettimeofday(&ta, NULL);mseca = ta.tv_sec*1000000.0+ ta.tv_usec;#endif/插入运行代码/Insert your code here/ Sleep(1000);