第七讲-定积分的近似计算.pptVIP

例：用int函数计算定积分：解：symsx;f=1/(1+x^2);int(f,x,0,1)f=sym(1/(1+x^2));int(f,x,0,1)int(1/(1+x^2),x,0,1)或int(1/(1+x^2),0,1)或或int举例第29页，共39页，星期日，2025年，2月5日double(a)将a转化为双精度型，若a是字符，则取对应的ASCII码a=3;double(a)double(a)例：ans=3ans=97其它相关函数第30页，共39页，星期日，2025年，2月5日x=1:0.001:2;y=exp(x.^(-2));trapz(x,y)梯形法：抛物线法：quad(exp(x.^(-2)),1,2,10e-10)符号积分法：symsxint(exp(x^(-2)),x,1,2)例1：用Matlab函数近似计算积分数值实验第31页，共39页，星期日，2025年，2月5日第1页，共39页，星期日，2025年，2月5日问题背景和实验目的定积分的近似计算定积分计算的基本公式是牛顿－莱布尼兹公式。但当被积函数的原函数不知道时，如何计算？这时就需要利用近似计算。特别是在许多实际应用中，被积函数甚至没有解析表达式，而是一条实验记录曲线，或一组离散的采样值，此时只能用近似方法计算定积分。本实验主要研究定积分的三种近似计算算法：矩形法、梯形法和抛物线法。同时介绍Matlab计算定积分的相关函数。第2页，共39页，星期日，2025年，2月5日矩形法梯形法抛物线法数值积分的常见算法主要内容Matlab求积分函数数值积分函数：trapz、quad、dblquad、triplequad符号积分函数：int第3页，共39页，星期日，2025年，2月5日矩形法定积分的定义：定积分的近似计算第4页，共39页，星期日，2025年，2月5日矩形法d充分小定积分的近似：通常我们取左点法右点法中点法点可以任意选取，常见的取法有：

左端点，右端点和中点。第5页，共39页，星期日，2025年，2月5日步长节点右点法：中点法：左点法：左点法、右点法和中点法第6页，共39页，星期日，2025年，2月5日解：矩形法举例==h=1/100=0.01,xi=i*h,a=0,b=1,n=100例：用不同的矩形法计算下面的定积分(取n=100)，

并比较这三种方法的相对误差。左点法：右点法：中点法：(i=0,1,2,...,100)第7页，共39页，星期日，2025年，2月5日理论值：左点法相对误差：误差分析矩形法举例右点法相对误差：中点法相对误差：不同的方法有不同的计算精度有没有更好的近似计算定积分的方法?第8页，共39页，星期日，2025年，2月5日定积分几何意义第9页，共39页，星期日，2025年，2月5日曲边小梯形的面积可以由直边小梯形的面积来近似整个曲边梯形的面积：梯形法第10页，共39页，星期日，2025年，2月5日如果我们n等分区间[a,b]，即令：则==梯形公式梯形法梯形公式与中点公式有什么区别?第11页，共39页，星期日，2025年，2月5日解：==例：用梯形法计算下面定积分(取n=100)，

并计算相对误差梯形法举例a=0,b=1,n=100,f(x)=1/(1+x2)==h=1/100=0.01,xi=i*h,yi=f(xi)相对误差：第12页，共39页，星期日，2025年，2月5日2n等分区间[a,b]，得该直线用抛物线代替，计算精度是否会更好？计算每个节点上的函数值：抛物线法在区间[x0,x2]上，用过以下三点的抛物线来近似原函数f(x)。第13页，共39页，星期日，2025年，2月5日设过以上三点的抛物线方程为：则在区间[x0,x2]上，有y=?x2+?x+?=p1(x)抛物线法第14页，共39页，星期日，2025年，2月5日同理可得：相加即得：抛物线法第15页，共39页，星期日，2025年，