第24课-计算圆周率π的近似值课件.ppt

1、授课人：授课人：杨鹏杨鹏第第24课课 计算圆周率计算圆周率的近似值的近似值授课人：授课人：杨鹏杨鹏1.用迭代法 求的近似值。2.用蒙特卡洛法求的值。3.设计一个算法计算1900年1月1日之后的任何一天是星期几？用VB程序实现这一算法。授课人：授课人：杨鹏杨鹏复习引入复习引入 解析算法基本思想：用解析的方法找出表示问题的前提条件与所求结果之间关系的数学表达式，并通过表达式的计算来实现问题求解。在设计解析算法时，必须注意保证计算过程描述的正确性。授课人：授课人：杨鹏杨鹏问题提出问题提出 圆周率圆周率=3.1415926有无穷多位，历史上 有很多科学家计算过这个数，古希腊学者阿基米德、我国杰出的数学

2、家和天文学家祖冲之等，都提出许多求解的方法。公元前3世纪，古希腊学者阿基米德用圆内接或外切正多边形来近似代替圆，得到近似圆周率。到了公元5世纪，我国杰出的数学家和天文学家祖冲之在刘徽“割圆术”的基础上，使用更加精密的方法，精确到小数点后七位，求出结果在3.14159263.1415927间，这个的圆周率，在当时是非常了不起的成就。直到16世纪，才由中亚人计算出精确度为小数16位。随着计算机的出现，的计算取得了新的突破。利用计算机技术将值计算到小数点后20多亿位，甚至到上千亿位。授课人：授课人：杨鹏杨鹏古今中外的数学家们计算的近似值主要有三类方法：1.利用“正多边形逼近正多边形逼近”的方法。2.

3、使用迭代法迭代法。可利用下列式子计算。3.使用蒙特卡洛法蒙特卡洛法（Monte-Carlo method）。蒙特卡洛法是使用随机模拟实验结果进行统计来求得的近似值的方法。授课人：授课人：杨鹏杨鹏实践活动实践活动用级数公式，计算圆周率的近似值。授课人：授课人：杨鹏杨鹏1.算法分析按照级数公式，由若干项求和得到，第1项为，第2项为，第3项为，第4项为分析可得第n项为an与前一项an-1的关系是。2S11a312a52313a ,7352314a1211nnaann授课人：授课人：杨鹏杨鹏2.算法流程图：变量n:输入计算的项数变量pi:计算的结果变量an:第n项的值授课人：授课人：杨鹏杨鹏3.程序界

4、面设计：授课人：授课人：杨鹏杨鹏控件属性值的设置参考如下表：对象名属性名属性值说明Form1Caption计算计算的近似值的近似值说明程序的功能Text2Text空串空串输出计算的结果（的近似值）BackColor蓝色蓝色Text1Text空串空串输入各个电阻的阻值（每次输入一个）BackColor黄色黄色Command1Caption计算计算说明命令按钮的作用Label1Caption输入项数输入项数n：说明文本框Text1的作用Label2Caption计算计算的近似值的近似值说明文本框Text2的作用Image1Picturepi.jpg显示图片pi.jpg（级数公式）授课人：授课人：杨

5、鹏杨鹏4.分析算法流程，编写代码PrivateSubCommand1_Click()DimPiAsDouble,anAsDoubleDimm,nAsIntegern=Val(Text1.Text)Pi=1an=1Form=2Tonan=an*(m-1)/(2*m-1)计算第n项的值Pi=Pi+an进行累加NextText2.Text=Str(Pi*2)EndSub授课人：授课人：杨鹏杨鹏5.运行调试程序;授课人：授课人：杨鹏杨鹏课堂练习1.使用蒙特卡洛法蒙特卡洛法求的近似值。利用求单位正方形与内接圆面积的比例关系来求得的近似值。单位圆的1/4面积是一个扇形，它是边长为1单位正方形的一部分。如果

6、能求出扇形面积s1在正方形面积s中占的比例k=s1/s，它的值也等于/4，从而就计算得到的值。怎样求出扇形面积在正方形面积中占的比例k呢？蒙特卡洛法是在正方形中随机投入很多点，使所投的点落在正方形中每一个位置的机会相等。有些点将落在扇形内，而另一些点将会落在扇形外，落在扇形内的点数m与所投点的总数n之间比m/n即为k的近似值。蒙特卡洛(MonteCarlo)方法，或称计算机随机模拟方法，是一种基于“随机数”的计算方法。这一方法源于美国在第二次世界大战时研制原子弹的“曼哈顿计划”。该计划的主持人之一、数学家冯诺伊曼用驰名世界的赌城摩纳哥的MonteCarlo来命名这种方法，为它蒙上了一层神秘色彩

7、。授课人：授课人：杨鹏杨鹏使用随机函数Rnd随机产生两个小数x、y,可构成一个坐标点(x,y），判断是否落在扇形内的充要条件是x2+y2=1，总计投放n次，根据条件判断统计得到落在扇形内的次数m，再根据两者的比例m/n，求得=4*m/n的近似值。授课人：授课人：杨鹏杨鹏1）参考程序代码：PrivateSubCommand1_Click()DimPiAsDouble,xAsDouble,yAsDoubleDimmAsLong,nAsLongRandomizeTimer随机数初始化n=Val(Text1.Text)读入投放次数nIfn=0ThenMsgBox请输入投放次数nExitSubEndIf

8、m=0ForI=1Tonx=Rnd()y=Rnd()Ifx2+y2=1 And day=1 And day=1 Or day=1 Or day=1 And day=1 And day=1 And day=1 And day=28)Then ok=1Then ok=1 If ok=0 Then If ok=0 Then Text4.Text=Text4.Text=输入日期有错误输入日期有错误 授课人：授课人：杨鹏杨鹏ElseElse If month=1 Or month=2 Then If month=1 Or month=2 Then year=year year=year-1-1 mont

9、h=month month=month+12+12 End If End If c=year 100 c=year 100 y=year Mod 100 y=year Mod 100 m=month m=month d=day d=day w=Int(c/4)-2 w=Int(c/4)-2*c+y+Int(y c+y+Int(y/4)+Int(13 /4)+Int(13*(m+1)(m+1)/5)+d 1 /5)+d 1 蔡勒公式蔡勒公式 Text4.Text=Str(c)+Str(y)+Str(m)+Str(d)+Str(wText4.Text=Str(c)+Str(y)+Str(m)+St

10、r(d)+Str(w)w=(w+700)Mod 7+1 w=(w+700)Mod 7+1 求除以求除以7 7的余数的余数,w,w加上加上700700保保证余数一定是整数证余数一定是整数 Text4.Text=WeekdayName(wText4.Text=WeekdayName(w)转换为星期形式转换为星期形式 End IfEnd IfEnd SubEnd Sub授课人：授课人：杨鹏杨鹏PrivateSubText1_Click()Text1.Text=Text2.Text=Text3.Text=“利用WeekdayName函数可将值转换成星期形式：值1234567返回值星期日星期日星期一星期一星期二星期二星期三星期三星期四星期四星期五星期五星期六星期六授课人：授课人：杨鹏杨鹏程序运行结果参考如下：星期计算问题测试数据表测 试 数 据 运 行 结 果第一组19051905年月年月1 1日日星期日星期日第二组19881988年年8 8月月8 8日日星期一星期一第三组20052005年年5 5月日月日星期日星期日第四组20052005年年1 1月月5454日日输入有误输入有误第五组20052005年年1313月月4 4日日输入有误输入有误第六组20052005年年2 2月月2929日日输入有误输入有误