实习手册 二 的计算.ppt

课程主页 /ligangliu/Courses/MathSoftware_2006/default.htm 实习手册 实验二: ? 的计算 想一想：怎样算 ? ？ 圆的面积 单位圆的面积等于 ? 计算第一象限内的单位圆的面积，方法为：把它分成n个窄的曲边梯形，计算S大，S小，其中n可以为1000, 10000, ... 无穷级数法 arctg x =x-x3/3+x5/5-x7/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(x,y) P在圆内?x2+y2?1 ?/4的近似值为落在圆内 的点数与总投点数的 比值 实验三: 最佳分数近似值 问题： 祖冲之的圆周率为 A. 3.141592… B. 355/113 ? 上述哪一个值更佳？ 寻找最佳 最佳近似值：误差小，分母小。 准确值 z= ? , 分数近似值 u/v， 误差 d=|z-u/v|。 选法1：递增选v，使d递减。 选法2：递增选v，使 vd 递减。 选法3：以v2d小为佳，选十佳。 * * 数值积分法 梯形公式 辛普森公式 A S C B O 某次随机投点的结果 * * *