《计算方法 第五章 数值积分.pptVIP

程序5-1: 利用复化Simpson公式求被积函数f(x)在给定区间上的积分值。 文件: FSimpson.m function S = FSimpson(f,a,b,N) h = (b-a)/N; fa = feval(f,a); fb = feval(f,b); S = fa+fb; x = a; for i=1:N x=x+h/2; fx=feval(f,x); S=S+4*fx; x=x+h/2; fx=feval(f,x); S=S+2*fx; end S = h*S/6; 实例: 利用复化Simpson公式计算积分S = ∫ x/(4+x^2) dx 输入: f=@fq; a=0;b=1;N=256; S=FSimpson(f,a,b,N); 上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 积分第二中值定理 梯形(trapezia)公式具有 次代数精度。 故 均差性质 1 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 设在区间[a, b]上函数f (x)连续，而函数?(x)可积且不变号，则在开区间(a, b)内至少存在一点ξ，使 积分第二中值定理： 均差性质： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 2. Simpson公式及其余项 Cotes系数为 求积公式为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 上式称为Simpson求积公式，也称三点公式或抛物线公式 记为 Simpson公式的余项为 Simpson公式具有 次代数精度。 3 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 3. Cotes公式及其余项 Cotes系数为 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 求积公式为 上式称为Cotes求积公式，也称五点公式 记为 Cotes公式的余项为 Cotes公式具有5次代数精度。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 第3节 复化求积公式 高次插值有Runge 现象，故采用分段低次插值 ? 分段低次合成的 Newton-Cotes 复化求积公式。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. ? 复化梯形公式： 在每个 上用梯形公式： = Tn 由介值定理知： 使 即有： 余项： Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. ? 复化 Simpson 公式： 4 4 4 4 4 = Sn 注：为方便编程，可采用另一记法：令 n’ = 2