计算方法数值积分演示文稿.pptVIP

其中 (k=0,1,2,…,4N)， 3.复合Cotes公式 当前第30页\共有67页\编于星期五\14点 4、复合Simpson公式算法 (1) 输入a,b,N (2) (3) 当 i=1,2, …,N时 做循环 ① x=x+h ② s=s+4f(x) ③ x=x+h ④ s=s+2f(x) (4) 当前第31页\共有67页\编于星期五\14点 例 5.2：利用数据表 xk 0 1/8 1/4 3/8 1/2 5/8 3/4 7/8 1 f (xk) 4 3.93846 3.76470 3.50685 3.20000 2.87640 2.46000 2.26549 2 计算积分 这个问题有明显的答案 取n = 8用复合梯形公式 当前第32页\共有67页\编于星期五\14点 取n=4,用辛普森公式 当前第33页\共有67页\编于星期五\14点 二、复合求积公式的余项 梯形公式的余项为 对于复合梯形公式则有 若 在[a,b]上连续，则存在 ，使 1、复合梯形公式的余项 当前第34页\共有67页\编于星期五\14点 所以 由 在[a,b]上连续可知， 在[a,b]上有界，于是存在常数M2，使 ， 1、复合梯形公式的余项 故 当前第35页\共有67页\编于星期五\14点 2、复合Simpson公式的余项 同理 由 在[a,b]上连续可知， 在[a,b]上有界，于是存在常数M4，使 故 当前第36页\共有67页\编于星期五\14点 3、复合Cotes公式的余项 由 在[a,b]上连续可知， 在[a,b]上有界，于是存在常数M6，使 同理 故 当前第37页\共有67页\编于星期五\14点 当 时， ，于是从这些余项公式可以看出， 当时，复合求积公式TN ,SN , CN都收敛于定积分值I，而且收敛速度一个比一个快。 二、复合求积公式的余项 当前第38页\共有67页\编于星期五\14点 例5.3 用复合梯形公式、复合Simpson公式、复合Cotes公式在取相同节点的情况下，计算定积分 的近似值。设把区间8等分。 解： 把区间[0,1] 8等分， ，共有9个节点， 节点表示为 (k=0,1,2,…,8)。 当前第39页\共有67页\编于星期五\14点 (1) 用复合梯形公式计算，相当于取 当前第40页\共有67页\编于星期五\14点 (2) 用复合Simpson 公式计算，相当于取N=4，把区间[0,1]N等分，然后在每个子区间上使用Simpson公式， 当前第41页\共有67页\编于星期五\14点 (3) 用复合Cotes 公式计算，相当于取N=2，把区间[0,1]N等分，然后在每个子区间上使用Cotes公式， 的准确值为0.9460831…… 当前第42页\共有67页\编于星期五\14点 回顾：复合求积公式的余项 1. 复合梯形公式的余项 2. 复合辛普森公式的余项 3. 复合柯特斯公式的余项 当前第43页\共有67页\编于星期五\14点 一、变步长梯形公式 1.当把区间[a,b] 等分时，步长 复合梯形公式为 2.当把区间[a,b] 等分时,步长 复合梯形公式为 §5.3 变步长求积公式 当前第44页\共有67页\编于星期五\14点 改写上式得： 复合梯形公式的递推公式 当前第45页\共有67页\编于星期五\14点 二、变步长梯形公式算法 1. 输入a,b,精度eps; 2. h=b-a 3. 做循环 ④ ① T1=T ② S=0 ③ 对 x=a+h/2(初值), b(终值) h (步长)做循环 s=s+f(x) ⑤ h=h/2 4. 则返回 3 5. 输出T 当前第46页\共有67页\编于星期五\14点 1. 将[a,b]N等分后复合梯形公式的余项，h= 2. 将[a,b]2N等分后复合梯形公式的余项，h= 设 在[a,b]上变化不大，即有 于是 整理得 当前第47页\共有67页\编于星期五\14点 同