数值分析-第6章数值积分与数值微分.pptVIP

第6章 数值积分 与数值微分 内容概要 求积公式简介 牛顿-柯特斯公式 复化求积公式 龙贝格求积公式 高斯型求积公式 数值微分 6．1 求积公式简介 6．1 ．1 求积公式 6．1．2 求积公式的余项与代数精度 6．1．2 求积公式的余项与代数精度 6．1．3 插值型求积公式 6．1．3 插值型求积公式 插值型求积公式的代数精度 6．1．4 内插求积公式 6．2 牛顿-柯特斯公式 6．2．1 梯形公式 梯形公式的余项与代数精度 6．2．2 抛物形公式 抛物形公式的几何意义 抛物形公式的余项与代数精度 抛物形公式的MATLAB实现 程序6-1 function Z=srule(f,a0,b0) %f是被积函数。 %a0和b0是被积区间的上下限。 %输出的值是所求的积分值。 h=(b0-a0)/2; C=zeros(1,3); C=[feval(f,a0),feval(f,(a0+b0)/2),feval(f,b0)]; S=h*(C(1)+4*C(2)+C(3))/3; Z=S; 6．2．3 牛顿-柯特斯公式 6．2．3 牛顿-柯特斯公式 牛顿-柯特斯公式的余项与代数精度 牛顿-柯特斯公式 特例 牛顿-柯特斯公式 的几何意义 牛顿-柯特斯公式 的几何意义 6．3 复化求积公式 复化求积公式的基本思想：将积分区间划分成若干个小区间，在各小区间上采用低次的求积公式（梯形公式或抛物形公式），然后再利用积分的区间可加性，把各区间上的积分加起来，得到新的求积公式。 6．3．1 复化梯形公式 复化梯形公式的误差估计 复化梯形公式的MATLAB实现 程序6-2 function s=trapr1(f, a, b, n) %f是被积函数。a, b分别为积分的上下限。 %n是子区间的个数。s是梯形总面积。 h=(b-a)/n; s=0; for k=1:(n-1) x=a+h*k; s=s+feval(f, x); end s=h*(feval(f, a)+feval(f, b))/2+h*s; 6．3．2 复化抛物形公式 复化抛物形公式的误差估计 复化抛物形公式的MATLAB实现 程序6-3 function s=simpr1(f, a, b, n) %f是被积函数。a，b分别为积分的上下限。 %n是子区间的个数。s是梯形总面积。 h=(b-a)/(2*n); s1=0; s2=0; for k=1:n x=a+h*(2*k-1); s1=s1+feval(f, x); end for k=1:(n-1); x=a+h*2*k; s2=s2+feval(f, x); end s=h*(feval(f, a)+feval(f, b)+4*s1+2*s2)/3; 6．4 龙贝格求积公式 6．4．1 龙贝格求积公式 6．4．1 龙贝格求积公式 6．4．1 龙贝格求积公式 数表 程序6-4（第二部分） while((errtol)(Jn))|(J4) J=J+1; h=h/2; s=0; for p=1 : M x=a+h*(2*p-1); s=s+feval(f, x); end R(J+1,1)=R(J, 1)/2+h*s; M=2*M; for K=1:J R(J+1,K+1)=R(J+1,K)+(R(J+1,K)-R(J,K))/(4^K-1); end err=abs(R(J,J)-R(J+1,K+1))； end quad=R(J+1,J+1)； 6．5 高斯型求积公式 2. 最高代数精度的求法 3. 高斯型求积公式 高斯型求积公式的截断误差 6．5．2 几种常用的高斯型求积公式 高斯-勒让德求积公式的截断误差 2．高斯-勒让德求积公式的MATLAB实现 程序6-5 function quad=gauss(f,a,b,X,A) %f是被积函数。 %a和b分别为积分的上下限。 %X是节点的横坐标构成的向量。 %A是系数构成的向量。 %输出所求积分值。 N=length(X); T=zeros(1,N); T=(a+b)/2+((b-a)/2)*X; quad=((b-a)/2)*sum(A.*feval(f,T)); 3. 高斯-拉盖尔（Gauss-Laguerre）求积公式 3. 高斯-拉盖尔（Gauss-Laguerre）求积公式