试验三 数值积分.doc

《试验三 数值积分.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、佛山科学技术学院实验报告课程名称数值分析实验项目常微分方程问题初值问题数值解法专业班级姓名学号指导教师成绩日期试验三数值积分一、实验目的1、理解如何在计算机上使用数值方法计算定积分的近似值;2、学会复合梯形、复合Simpson和龙贝格求积分公式的编程与应用。3、学会利用Matlab提供的积分函数求二重积分。函数名含义函数名含义trapz梯形求积法dblquad二重积分quadSimpson求积法triplequad三重积分二、实验要求（1）按照题目要求完成实验内容；（2）写出相应的Matlab程序；（3）给出实验结果(可以用表格展示实验

2、结果)；（4）分析和讨论实验结果并提出可能的优化实验。（5）写出实验报告。三、实验步骤1、用不同数值方法计算积分（1）取不同的步长，分别用复合梯形及复合辛普森求积计算积分，给出误差中关于的函数，并与积分精确值比较两公式的精度，是否存在一个最小的，使得精度不能被改善？（2）用龙贝格求积计算完成问题（1）。2、计算二重积分（1）若积分区域；（2）若积分区域求积公式参考代码：1、复化梯形求积公式functionI=T_quad(x,y)%复化梯形求积公式，其中%x---向量，被积函数自变量的等距节点%y---向量，被积函数在节点处的函数值n=

3、length(x);m=length(y);ifn~=merror('ThelengthsofXandYmustbeequal');return;endh=(x(n)-x(1))/(n-1);a=[12*ones(1,n-2)1];I=h/2*sum(a.*y);2、复化Simpson求积公式functionI=S_quad(x,y)%复化Simpson求积公式，其中%x---向量，被积函数自变量的等距节点%y---向量，被积函数在节点处的函数值n=length(x);m=length(y);ifn~=merror('Thelength

4、sofXandYmustbeequal');return;endifrem(n-1,2)~=0%如果n-1不能被2整除，则调用复化梯形公式warning('给出的点数不适用于Simpson公式，调用梯形公式求积分');I=T_quad(x,y);return;endN=(n-1)/2;h=(x(n)-x(1))/N;a=zeros(1,n);%如果等分成N段，且从1开始计算，应该是公式N=(n-1)/2;fork=1:Na(2*k-1)=a(2*k-1)+1;a(2*k)=a(2*k)+4;a(2*k+1)=a(2*k+1)+1;%如果

5、下标k从1开始计算，则每小段区间的下标应该为1、3、5、7、、、，中间节点%的下标分别是2、4、6、、，对应书本P107公式（3.4），它们的系数分别为1，4，1；endI=h/6*sum(a.*y);3、Romberg求积公式functionI=R_quad_iter(fun,a,b,ep)%Romberg求积公式，其中%fun---被积函数，在调用之前需要用inline()函数定义其为内置函数。%a,b---积分区间的端点，要求a

6、in指设定参数的个数ifnargin<4ep=1e-5;end;m=1;h=b-a;I=h/2*(feval(fun,a)+feval(fun,b));T(1,1)=I;while1N=2^(m-1);h=h/2;I=I/2;fori=1:NI=I+h*feval(fun,a+(2*i-1)*h);endT(m+1,1)=I;M=2*N;k=1;whileM>1;T(m+1,k+1)=(4^k*T(m+1,k)-T(m,k))/(4^k-1);M=M/2;k=k+1;endifabs(T(k,k)-T(k-1,k-1))

7、;endm=m+1;endI=T(k,k);实验三参考答案：1.1用复合梯形公式计算积分（用不同的步长h=0.05和0.001,与精确值0.比较）(1)步长h=0.05x=0.0001:0.1255:1;y=sin(x)./x;I=T_quad(x,y)I=0.4063%也可以用eps代替0.0001，表示一个很小的正数（<10-15）x=eps:0.125:1;y=sin(x)./x;>>I=T_quad(x,y)I=0.27011.2用复合辛普森求积公式计算积分。x=0.001:0.0001:1;y=sin(x)./x;I=S_qu

8、ad(x,y)I=0.2740x=0.0001:0.0001:1;y=sin(x)./x;I=S_quad(x,y)Warning:给出的点数不适用于Simpson公式，调用梯形公式求积分>InS_quad