数值积分与数值微分实验报告.docxVIP

实验项目名称：数值积分与数值微分 实验目的：熟悉数值积分与数值微分步骤，并用计算机实现 实验内容： a. 用三点公式求在处的导数值，并估计误差。的值由下表给出： x1.01.11.2f(x)0.25000.22680.2066 b. 用不同的数值方法计算积分： 取不同的h,分别用复合梯形及复合辛普森求积计算积分，给出误差中关于h的函数，并与积分精确值比较两个公式的精度，是否存在一个最小的h,使得精度不能再被改善？ 用龙贝格求积计算完成问题(1) 用自适应辛普森积分，使其精度达到 实验环境：matlab r2010b 实验过程与分析： a. 一．实验原理： 三点公式有以下三种形式：  = 1 \* ROMAN I： f (x0) =12h(-3y0+4y1-y2)  = 2 \* ROMAN II： f (x1) =12h(y2-y0)  = 3 \* ROMAN III： f (x2) =12h(y0-4y1+3y2) 二．源程序： 1). 主函数(保存为m文件直接点击运行即可得到结果) clear clc X =[1.0 1.1 1.2]; Y=[0.2500 0.2268 0.2066]; df=threepoints(X,Y) 2). 三点公式子函数threepoints.m: function df=threepoints(X,Y) syms x; L=Lagrange(X,Y); L=diff(L,x); df=subs(L,x,X); end 3).求插值函数的子函数Lagrange.m： function L=Lagrange(X,Y,t) if(length(X)~=length(Y)) disp(d do not agree); else k=length(X); syms x; w=1; for j=2:k w=w*(x-X(j))/(X(1)-X(j)); end L=Y(1)*w; for i=2:k-1 w=1; for j=[1:(i-1),(i+1):k] w=w*(x-X(j))/(X(i)-X(j)); end L=L+Y(i)*w; end w=1; for j=1:k-1 w=w*(x-X(j))/(X(k)-X(j)); end L=L+Y(k)*w; end if nargin==3 L=subs(L,x,t); else L=vpa(collect(L,x),6); end end b. 一．实验原理： 复合梯形公式: Tn=h2[fa+2k=1n-1fxk+f(b)] 复合辛普森公式： Sn=h6[fa+4k=0n-1fxk+12+2k=1n-1fxk+f(b)] 然后给出误差公式，并与积分精确值比较两个公式的精度，判断是否存在一 个最小的h，使得精度不能再被改善。 龙贝格公式： Tm(k)=4m4m-1Tm-1(k+1)-14m-1Tm-1k k=1,2,… 计算表如下： khT0(k)T1(k)T2(k)T3(k)T4(k)…0b-aT0(0)1b-a2T0(1)T1(0)2b-a4T0(2)T1(1)T2(0)3b-a8T0(3)T1(2)T2(1)T3(0)4b-a16T0(4)T1(3)T2(2)T3(1)T4(0)……………………可以证明，如果f(x)充分光滑，那么上表每一列的元素及对角线元素均收敛到所求的积分值I，即： limk→∞Tm(k)=I (m固定)， limm→∞Tm(0)=I. 自适应辛普森方法： 设给定精度要求ε0,计算积分 If=abf(x)dx 的近似值，先取步长h=b-a，应用辛普森公式有 If=abfxdx=Sa,b-b-a180(h2)4f4(η), η∈(a,b)  = 1 \* GB3 ① 其中 Sa,b=h6[fa+4fa+b2+f(b)]. 若把区间[a,b]对分，步长h