中国科学技术大学：《数学实验》课程教学资源（实验讲稿PPT）实验二：π的计算 实验三：最佳分数近似值

★ NERSIN i oF ScIENCE TECHNOLOGY CHINA 数学实验 李尚志教授 中国科学技术大学 数学系

数学实验 李尚志 教授 中国科学技术大学 数学系

实验二:π的计算 想一想:怎样算π? 当一回祖冲之!

实验二:  的计算 想一想：怎样算  ？ 当一回祖冲之！

员的面积 单位圆的面积等于兀 计算第一象限内的单位圆的面积, 方法为:把它分成n个窄的曲边梯形, 计算S大S小,其中n可以为1000 10000

圆的面积 单位圆的面积等于  计算第一象限内的单位圆的面积， 方法为：把它分成n个窄的曲边梯形， 计算S大，S小，其中n可以为1000, 10000,

数值积分法 2 ∫1/(1+x2)4x=4 20212/19

2021/2/19 数值积分法 4 1 (1 ) 4 1 1 0 2 1 0 2   + = − =   x dx x dx

梯形公式 b f(x)dx af(a)+f(b)/2+ ∑f(x),x=a+ 20212/19

2021/2/19 梯形公式 i n b a x a n b a f x f x dx f a f b i n i i b a − = + −  + +   − = ( )] , ( ) [( ( ) ( ))/ 2 1 1

辛普森公式 f(dx elf(a)+ f(b)+ 4∑f(x105)+2∑f(x) b i=0 6n 20212/19

2021/2/19 辛普森公式 n b a f x f x f x dx f a f b n i i n i i b a 6 4 ( ) 2 ( )] ( ) [ ( ) ( ) 1 1 1 0 0.5 − +  + +    − = − = +

无穷级数法 arct X=X×3/3+×5/5-×7/7+×9/9 ■■■ π/4=arcg1=1-1/3+1/5-1/7+1/9 收敛太慢! 应当比1小很多,级数收敛才快。 /4= arct 1/2+ arct 1/3 /4=4 arct 1 5-arctg 1/239

无穷级数法 arctg x =x-x 3 /3+x5 /5-x 7 /7+x9 /9-...  /4 = arctg 1 =1-1/3+1/5-1/7+1/9-… 收敛太慢！ |x| 应当比 1 小很多，级数收敛才快。 /4 = arctg 1/2 + arctg 1/3 /4 = 4 arctg 1/5 - arctg 1/239

蒙特卡洛方法 求出S=/4占S正方形ACBo的比例 随机投点P(Xy) P在圆内台×2+y2≤1 丌/4近似值为落在圆内B 的点数与总投点数的 比值 O A

蒙特卡洛方法 求出S=/4占S正方形ACBO的比例。 随机投点P(x,y) P在圆内x 2+y21 /4的近似值为落在圆内 的点数与总投点数的 比值 A S C B O

某次随机投点的结果 E上0000000000008,E 20212/19

2021/2/19 某次随机投点的结果 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 3.14800000000000013` 3.14800000000000013` 3.14800000000000013`

实验三:最佳分数近似值 问题: 祖冲之的圆周率为 A.3.141592..B.355/113? 上述哪一个值更佳?

实验三: 最佳分数近似值 问题： 祖冲之的圆周率为 A. 3.141592… B. 355/113 ? 上述哪一个值更佳？