[数学]计算方法 第五章 数值积分.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); 上式称为梯形求积公式,也称两点公式，记为 梯形公式的余项为 积分第二中值定理 梯形(trapezia)公式具有 次代数精度。 故 均差性质 1 设在区间[a, b]上函数f (x)连续，而函数?(x)可积且不变号，则在开区间(a, b)内至少存在一点ξ，使 积分第二中值定理： 均差性质： 2. Simpson公式及其余项 Cotes系数为 求积公式为 上式称为Simpson求积