★ 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+y21 /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 ? 上述哪一个值更佳？