王向东数学实验课本(可编辑)2-4.doc

实验四　积分 【实验目的】 了解积分的基本概念。 掌握积分在计算面积、体积等问题中的应用。 了解变步长积分和广义积分 学习掌握MATLAB软件有关的命令。 【实验内容】 用符号积分命令int计算积分 计算数值积分 计算数值积分 【实验准备】 1．积分的有关理论 定积分：积分是微分的无限和，函数在区间上的积分定义为 其中从几何意义上说，上非负函数，记分值是曲线与直线及轴所围的曲边梯形的面积。有界连续（或几何处处连续）函数的积分总是存在的。 微积分基本定理（Newton-Leibniz公式）：在上连续，且，则有 这个公式表明导数与积分是一对互逆运算，它也提供了求积分的解析方法：为了求的定积分，需要找到一个函数，使的导数正好是，我们称是的原函数或不定积分。不定积分的求法有学多数学技巧，常用的有换元积分和分部积分法。从理论上讲，可积函数的原函数总是存在的，但很多被积函数的原函数不能用初等函数表示，也就是说这些积分不能用解析方法求解，需用数值积分法解决。 在应用问题中，常常是利用微分进行分析，而问题最终归结为微分的和（即积分）。一些更复杂的问题是含微分的方程，不能直接积分求解。 多元函数的积分称为多重积分。二重积分的定义为 当非负时，积分值表示曲顶柱体的体积。二重积分的计算主要是转换为两次单积分来解决，无论是解析方法还是数值方法，如何实现这种转换，是解决问题的关键。 平面曲线的长度 空间曲线的长度 曲面的面积 2．积分的数值方法 梯形法：将划分为若干小区间则 在每一小区间上近似为一直线，用弦代替，有 从而 称为梯形公式。通常将区间等分，， 可以证明，当时由上述公式给出的梯形法是收敛的。 变步长积分法：以上介绍的梯形法中，划分是给定的。在实际应用中，等分数往往是难以确定的。下面介绍变步长梯形法，其思路是对于给定的精度，从开始，等分数逐次加倍，直至为此，一般说来，这样的做法会浪费计算量，幸运的是，与有如下的递推关系 其中。 重积分：重积分的数值计算可通过若干次单积分的组合实现，如对于二重积分 先化为二次计分 利用梯形法，先将区间等分，利用梯形积分公式可得 再将区间等分，利用梯形积分公式可得 3．积分的MATLAB命令 MATLAB中主要用int进行符号积分，用trapz,dblquad,quad,quad8等进行数值积分。 int(s) 符号表达式s的不定积分 int(s,x) 符号表达式s关于变量x的不定积分 int(s,a,b) 符号表达式s的定积分，a,b分别为积分的上、下限 int(s,x,a,b) 符号表达式s关于变量x的定积分，a,b分别为积分的上、下限 trapz(x,y) 梯形积分法，x时表示积分区间的离散化向量，y是与x同维数的向量，表示被积函数，z返回积分值。 quad8(‘fun’,a,b,tol) 变步长数值积分，fun表示被积函数的M函数名，a,b分别为积分上、下限，tol为精度，缺省至为1e-3. fblquad(‘fun’,a,b,c,d) 矩形区域二重数值积分，fun表示被积函数的M函数名，a,b分别为x的上、下限，c,d分别为y的上、下限. 可以用help int, help trapz, help quad等查阅有关这些命令的详细信息 【实验方法与步骤】 练习1 用符号积分命令int计算积分，MATLAB代码为： clear; syms x; int(x^2*sin(x)) 结果为 ans =-x^2*cos(x)+2*cos(x)+2*x*sin(x) 如果用微分命令diff验证积分正确性，MATLAB代码为： clear; syms x; diff(-x^2*cos(x)+2*cos(x)+2*x*sin(x)) 结果为 ans =x^2*sin(x) 练习2 计算数值积分.先用梯形积分法命令trapz计算积分，MATLAB代码为： clear; x=-2:0.1:2; y=x.^4; %积分步长为0.1 trapz(x,y) 结果为 ans = 12.8533 实际上，积分的精确值为。如果取积分步长为0.01, MATLAB代码为： clear; x=-2:0.01:2; y=x.^4; %积分步长为0.01 trapz(x,y) 结果为 ans =12.8005 可用不同的步长进行计算，考虑步长和精度之间的关系。一般说来,trapz是最基本的数值积分方法,精度低,适用于数值函数和光滑性不好的函数. 如果用quad8命令计算积分，得先编写被积函数的M函数 %M函数fun1.m function y=fun1(x) y=x.^4; 保存后,在命令窗口用MATLAB代码: clear; quad8(fun1,-2,2) vpa(quad8(fun1,-2,2),10) %以10位有效