实验报告7—数值积分.docVIP

标题：积分方程的数值积分计算 1.实验描述： 数值积分最突出的优点是它可以计算无法解析求解的积分问题。根据节点的选择方法可将数值积分分为常见的：组合梯形公式法、组合辛普生公式法、龙贝格积分法、自适应积分法、高斯—勒让德积分法。本实验利用5种方法计算同一积分，通过误差分析比较各种方法的优缺点。 2.实验内容： 计算，并进行误差分析。具体内容如下： 1.用组合梯形公式计算。 2.用组合辛普生公式计算。 3.用龙贝格积分计算，本次实验中采用4阶公式计算。 4.用自适应积分方法计算，本次实验中起始容差：。 5.用5点高斯—勒让德积分计算。 通过误差分析比较各种方法的优缺点。 3.实验原理及分析： 数值积分的目的是：通过在有限采样点上计算在区间上的定积分。设，若有：，其中形如：，则称为面积公式，为截断误差，为面积节点，为权。根据节点的选择方法可将积分方法分为：组合梯形公式法、组合辛普生公式法、龙贝格积分法、自适应积分法、高斯—勒让德积分法。下面着重介绍5种方法的原理： ①组合梯形公式法及误差分析： 设等距节点，将区间划分为宽度为的个子区间，个子区间的组合梯形积分公式有3种等价表示方法： ②组合辛普生公式法误差分析： 设等距节点，将区间分为个宽度为的子区间,个子区间的组合辛普生积分公式也有3种等价表示方法： ③龙贝格积分法及误差分析： 龙贝格积分法是利用理查森外推法来提高精度的，下面给出一般公式： 其中 ,为梯形公式；,为辛普生公式；，为布尔公式。 ④自适应积分法及误差分析： 自适应积分法节点非等间距，它能在函数值变化较大的区间采用较小步长进行计算，其基础是辛普生公式。计算的关键是结果能否通过精度测试。 ⑤高斯—勒让德积分法及误差分析： 高斯—勒让德积分法是一种非常高效的计算方法，其理论依据是函数基于勒让德多项式的表示，一般点高斯—勒让德公式对次数小于的多项式是精确的，公式为： 其中定义在区间上 对于非规范情形可根据高斯—勒让德变换进行计算：设定义在上，做变换，则积分公式为： 4.实验结果及结论分析： ①组合梯形公式法及误差分析： 计算结果：T =0.1697，ET = -0.0301 ②组合辛普生公式法误差分析： 计算结果：S =0.1401，ES = -0.0596 ③龙贝格积分法及误差分析： 计算结果： r =0.1997，Er =-2.3430e-005 ④自适应积分法及误差分析： 计算结果：Y=0.1783,EY=-0.0214 ⑤高斯—勒让德积分法及误差分析： 计算结果：G =0.1970，EG =-0.0027 附件1（代码）： ①组合梯形公式法及误差分析： ②组合辛普生公式法误差分析： ③龙贝格积分法及误差分析： ④自适应积分法及误差分析： ⑤高斯—勒让德积分法及误差分析： ?数值计算方法?实验报告 X=linspace(0,3,11); Y=sin(4*X).*exp(-2*X); h=0.3; T=0.5*0.3*(Y(1)+Y(11))+0.3*sum(Y(2:10)) syms x y=sin(4*x)*exp(-2*x); T0=int(y,0,3); ET=T-T0; ET=double(ET) X=linspace(0,3,11); Y=sin(4*X).*exp(-2*X); h=0.3; A=[]; for k=1:4 A=[A Y(2*k)]; end B=[]; for i=1:5 B=[B Y(2*i+1)]; end S=1/3*0.3*(Y(1)+Y(11))+2/3*0.3