复化求积公式求数值积分.docVIP

应用软件课程设计 ——复化求积公式求数值积分 数学0801班 *** 一、实验目的 程序1：利用复化梯形公式以及复化辛普森公式求解定积分的数值解。 程序2：分析剖分区间个数对复化梯形公式精度的影响。 程序3：比较MATLAB的quad、quadl命令与上述两种方法的精度；比较在相同剖分区间下两种求数值积分方法的精度。 分析与探讨两种方法精度不同的原因。 说明：原题目给出的积分精确值I=4.006994稍过粗糙，所以通过计算机求解得到更为精确的解析解。详见测试结果。 二、算法说明 自定义函数有： 积分函数（在本题中为）：hanshu.m（附录1） 复化梯形公式求定积分数值解tixing.m（附录2） 复化辛普森公式求定积分数值解xps.m（附录3） 积分原函数解析解jxj.m（附录4） 程序1：p1.m（附录5） 1、n=2000;%确定积分区间分割份数 2、tx=tixing(0,2,n);%用复化梯形公式求解 3、xps=xps(0,2,n);%用复化辛普森公式求解 4、显示结果 disp([积分区间分割分数为：,num2str(n)]) disp(复化梯形公式的求解结果：),disp(tx) disp(复化辛普森公式的求解结果：),disp(xps) 程序2-1：p21.m（附录6） 1、jxj=jxj(2)-jxj(0);%求出解析解。 2、tx=zeros(5,1);%给数值解向量赋值。 d=zeros(5,1);%给误差向量赋值。 3、i=1,2,3,4,5 4、n=10^i;%定义剖分份数。 5、tx(i)=tixing(0,2,n);%将剖分份数n代入，求出该n下的数值解。 d(i)=tx(i)-jxj;%求出误差。 6、结束循环 7、输出结果 程序2-2：p22.m（附录7） 1、jxj=jxj(2)-jxj(0);%求出解析解。 2、n=0;l=2000;%迭代次数。 dd=5;%迭代步长。 d=zeros(1,l); 3、i=1,2,…,l; 4、n=n+dd;%积分区间分割份数。 d(i)=tixing(0,2,n)-jxj;%求出误差。 5、结束循环 6、figure，x=dd:dd:l*dd;%构造横坐标向量。 7、plot(x,d,-m,linewidth,3),xlabel(积分区间等分份数),ylabel(误差大小),title(误差曲线),grid on 8、%axis([200 10000 0 5*10^-7])%用以控制横纵坐标范围，便于考察不同范围内的误差曲线情况，可以置于命令窗口中运行。 程序3：p3.m（附录8） 1、a=0;%给定定积分下界。 b=2;%给定定积分上界。 l=5;%给定数值积分个数。 c=1:l;%给定10的次方数。 n=10.^c;%给定积分区间的等分份数。 jxj=jxj(b)-jxj(a);%求该积分的解析解。 tx=zeros(l,1);xp=zeros(l,1);wucha=zeros(l,2);%分别定义两种求数值解的方法不同迭代次数下的解矩阵以及误差矩阵。 2、qwc=abs(quad(hanshu,0,2)-jxj);%求quad方法的误差。 qlwc=abs(quadl(hanshu,0,2)-jxj);%求quadl方法的误差。 3、i=c%进行循环 4、tx(i)=tixing(a,b,n(i));%以复化梯形公式求数值解。 xp(i)=xps(a,b,n(i));%以复化辛普森公式求数值解。 wucha(i,1)=abs(tx(i)-jxj);%求复化梯形公式的误差绝对值。 wucha(i,2)=abs(xp(i)-jxj);%求复化辛普森公式的误差绝对值。 5、结束循环 6、wc=num2str([c wucha]); 7、显示运行结果 disp(wc) disp(说明：误差矩阵第一列代表迭代次数10^n。) disp( 第二列代表复化梯形公式求得数值解的误差绝对值。) disp( 第三列代表复化辛普森公式求得数值解的误差绝对值。) disp([ quad方法求得的解析解误差绝对值为：,num2str(qwc)]) disp([ quadl方法求得的解析解误差绝对值为：,num2str(qlwc)]) 三、测试结果 计算的解析解。利用MATLAB软件对原函数求不定积分，得到原函数符号解。在MATLAB中format long的格式下，显示I的值为：4.006994223254704（可在MATLAB中直接通过命令“jxj(2-jxj(0)”驱动jxj.m显示该值。），精度高于题目所给。为保证计算中保持该数值的最大精度，在求解全过程中，凡是需要此解析解进行运算的程序，均通过直接调用函数