圆周率蒙特卡洛程序.docVIP

实验名称： 圆周率的计算 姓名 李志梅 学号实验地点 T5-207 实验类型 综合 设计 实验要求 选修 学时量 6 所用知识 数学建模 数学软件 运筹学 第一部分 割圆术 割圆方式 面积 面积 面积 周长 周长 分割程度 12 24 1572864 1536 1572864 所用工具 matlab matlab matlab matlab matlab 圆周率值 3.000000000000000 3.105828541230248 3.141592653694816 3.141590463228051 3.141592653587705 第二部分 数值积分 所用函数 积分方法 矩形 复化辛普森 复化辛普森 复化梯形 复化辛普森 分割程度 10000 1876 10 914 5 所用工具 Matlab Matlab Matlab Matlab Matlab 圆周率值 3.141592651923118 3.141590653696272 3.141592652423170 3.141592453646339 3.141592502458706 第三部分 蒙特卡洛方法 所用函数 循环次数 100 1000 10000 100000 满足的次数 74 769 7846 78419 7853665 所用工具 matlab matlab matlab matlab matlab 圆周率值 2.960000000000000 3.076000000000000 3.138400000000000 3.136760000000000 3.141466000000000 第四部分 级数方法 所用级数 计算项数 10 1000 10000 10 1000 精确位数 20 20 20 20 20 所用工具 Matlab Matlab Matlab Matlab Matlab 圆周率值 3.0418396189294032439 3.140592653839794135 3.1414926535900344895 3.1415926535897932385 3.1415926535897932385 第五部分 蒙特卡洛方法的应用 计算积分 积分结果 0.77364454279011131791 0.77364454279011131791 0.77364454279011131791 0.77364454279011131791 循环次数 10000 100000 1000000 100000000 满足的次数 7768 77331 773836 所用工具 matlab matlab matlab matlab 积分的估计值 0.776800000000000 0.773310000000000 0.773836000000000 0.773614250000000 第六部分 心得体会 通过这次数学实验体会到了其中的一些乐趣，也重新学习了一些与matlab有关的东西，同时也体会到了数值分析的重要性，利用数值计算，可以很好地计算一些积分的近似值，可以把复杂的问题简单化，在写某些算法的过程中遇到了好多小问题，但是在网络和同学的帮助下很好地解决了那些问题，最后都得到了令人满意的结果，所以也很喜欢这样的课程。 割圆术的算法： 面积： function s=f(n) p=1; format long; for k=1:n p=sqrt((1-sqrt(1-p^2/4))^2+p^2/4); end b=(1+1+p)/2; a=6*2^n s=a*sqrt(b*(b-1)^2*(b-p)) 周长： function c=f(n) p=1; format long; for k=1:n p=sqrt((1-sqrt(1-p^2/4))^2+p^2/4); end a=6*2^n c=p*a/2 矩形 function s=juxing(n) y=@(x)1/(1+x^2); x=0:(1/n):1; s=0; y=1./(1+x.^2); digits(22) for k=1:length(x)-1 s=s+((y(k)+y(k+1))/2)/n; end format long 4*s vpa(s,20) 的复化辛普森 function t=fx1 f=@(x)sqrt(1-x.^2); a=0; b=1; n=1;A=1; while(A0.9*10^(-7)) h=