利用复化梯形公式、复化simpson公式计算积分.pdfVIP

要 实 1、利用复化梯形公式、复化 simpson 公式计算积分 求 验 2、比较计算误差与实际误差 目 的 或 取 n=2,3,…,10 分别利用复化梯形公式、复化 simpson 公式计算积分 1 2 ，并与真值进行比较， 并画出计算误差与实际误差之间的曲线。 I x dx 0 利用复化梯形公式的程序代码如下： function f=fx(x) f=x.^2; ％首先建立被积函数，以便于计算真实值。 a=0; ％积分下线 实 b=1; ％积分上线 验 T=[]; ％用来装不同 n 值所计算出的结果 原 理 for n=2:10; （ h=(b-a)/n; ％步长 算 x=zeros(1,n+1); ％给节点定初值 法 for i=1:n+1 流 程 x(i)=a+(i-1)*h; ％给节点赋值 图 end 或 y=x.^2; ％给相应节点处的函数值赋值 者 含 t=0; 注 for i=1:n 释 t=t+h/2*(y(i)+y(i+1)); ％利用复化梯形公式求值 的 end 源 代 T=[T,t]; ％把不同 n 值所计算出的结果装入 T 中 码 end ） R=ones(1,9)*(-(b-a)/12*h.^ 2*2); ％积分余项（计算误差） true=quad(@fx,0,1); ％积分的真实值 A=T-true; ％计算的值与真实值之差（实际误差） x=linspace(0,1,9); plot(x,A,r,x,R,*) ％将计算误差与实际误差用图像画出来 注：由于被积函数是 x.^2 ，它的二阶倒数为 2，所以积分余项为： (-(b-a)/12*h.^ 2*2) — 利用 复化 simpson 公式的程序代码如下： 同样首先建立被积函数的函数文件： function f=fx1(x) f=x.^4; 实 验 原 理 （ a=0; ％积分下线 算 b=1; ％积分上线 法 T=[]; ％用来装不同 n值所计算出的结果 流 for n=2:10 程 h=(b-a)/(2*n); ％步长 图 x=zeros(1,2*n+1); ％给节点定初值 或 for i=1:2*n+1 者 x(i)=a+(i-1)*h; ％给节点赋值 含 end 注 y=x.^4; ％给相应节点处的函数值赋值 释 t=0; 的 for i=1:n 源 t=t+h/3*(y(2*i-1)+4*y(2*i)+y(2*i+1)); ％利用复化 simpson公式求值