数值课程设计详解.docVIP

数值分析课程设计报告 用复梯形公式、复Simpson公式和Cotes公式用Romberg积分法计算积分上机完成 （1）用复化梯形公式、复化Simpson公式Cotes公式计算积分，自己设置不同精度要求，对结果进行比较分析。 （2）用Romberg积分法计算积分，自己设置不同精度要求，对结果进行比较分析。 （3）记，在上面的计算中只取4位有效数字或7位有效数字，计算结果又有什么不同。 （5）上面计算精度可达8-20位有效数字吗？若可以，请说明实现过程，并举例。划分为等份，步长，选取等距节点构造出的插值型求积公式 In = (1.1) 称为牛顿-柯特斯公式，式中称为柯特斯系数，引进变换，则有 (1.2) 由于是多项式的积分，柯特斯系数的计算不会遇到实质性的困难，则当时， , 这时的求积公式就是我们所熟悉的梯形公式。 当时，按（1.2）式，这时的柯特斯系数为 ==, =- , = , 相应的求积公式时下列的辛普森公式： ， （1.3） 而n=4时的牛顿-柯特斯公式则特别称为柯特斯公式，而形式是 (1.4) 这里, () 2.复合梯形公式 将区间划分为n等份，分点，,,在每个子区间 ()上采用梯形公式（1.1），则得 == (2.1) 记 = (2.2) 称为复合梯形公式，其余项可由（1.10）式得 ,, 由于，且 ， 所以使 于是复合梯形公式的余项为 （2.3） 3. 复合辛普森求积公式 将区间 分为n等份，在每个子区间上采用辛普森公式，若,则得 == = （3.1） 记 = = +2] （3.2） 称为复合辛普森求积公式，其余项由（2.5）得 于是当时，与复合梯形公式相似有 （3.3） 由（3.6）式看出，误差阶为h4，收敛性是显然的，实际上，只要则可得到收敛性，即 此外，由于Sn中的求积系数均为正数，故知复合辛普森公式计算稳定。 ：算法程序 1、对积分 ，,使其精度达到 解： 先建立三种复化公式的函数文件,它们分别为复化梯形公式trap.m、复化公式为simpson.m、公式为cotes.m,三个函数的源程序如下: 复化梯形公式trap.m function T=trap(f,a,b,n) %trap.m %复化梯形公式求积分值 %f为积分函数 %[a,b]为积分区间 %n是等分区间份数 h=(b-a)/n;%步长 T=0; for k=1:(n-1) x0=a+h*k; T=T+limit(f,x0); end T=h*(limit(f,a)+limit(f,b))/2+h*T; T=double(T); (2) 复化公式simpson.m: function S=simpson(f,a,b,n) %simpson.m %Simpson公式求积分值 %f为积分函数 %[a,b]为积分区间 %n是等分区间份数 h=(b-a)/(2*n);%步长 s1=0; s2=0; for k=1:n x0=a+h*(2*k-1); s1=s1+limit(f,x0); end for k=1:(n-1) x0=a+h*2*k; s2=s2+limit(f,x0); end S=h*(limit(f,a)+limit(f,b)+4*s1+2*s2)/3; S=double(S); (3) 复化公式cotes.m: function C=cote(f,a,b,n) %cote.m %复化cotes公式求积分值 %f为积分函数 %[a,b]为积分区间 %n是等分区间份数 h=(b-a)/n;%步长 C=0; for i=1:(n-1) x0=a+i*h; C=C+14*limit(f,x0); end for k=0:(n-1) x0=a+h*k; s=32*limit(f,x0+h*1/4)+12*limit(f,x0+h*1/2)+32*limit(f,x0+h*3/4); C=C+s; end C=C+7*(limit(f,a)+limit(f,b)); C=C*h/90; C=double(C); 输入命令其结果显示如下： 计算积分 积分区间分别分为1,10,20等份 b.计算积分 2、用积分法计算,，精度. 解： 编写积分法的函数M文件,源程序如下(romberg.m): 将上述源程序另存为romberg.m后