定积分的运用.docxVIP

实验四 定积分及应用 实验的目的 1、掌握利用 Matlab 进行积分运算； 2、掌握积分在计算面积、体积等问题中的应用； 3、掌握各种积分指令的区别与特点。 实验的基本理论与方法 1、定积分定义： 函数 f ( x) 在区间 [ a,b] 上的定积分定义为： 设函数 f (x) 在 [ a, b] 上有界， 在区间 [ a,b] 上任取 n 1 个分点： a x0 x1 x2 xn 1 xn b ，把 [ a, b] 分成 n 个小区间 i [ xi 1 , xi ] ， i 1,2, , n 。这些分点构成对区间 [ a,b] 的一个分割，用 T 表示。小区 间 i 的长度为 xi xi xi 1。记 T max xi ，称为分割 T 的模。在区间 i [ xi 1 , xi ] 1 i n 上 取 点 i ( xi 1 i xi ) ， 做 函 数 值 f ( i ) 与 小 区 间 长 度 xi 的 乘 积 f ( i ) xi (i 1,2, n) ，并作和 S n xi 。当 T 0 时，和 S 总趋于确定的极限， f ( i ) i 1 f (x) 在区间 [ a, b] 上的定积分，记作 b 这时这个极限为函数 f ( x)dx 。即 a b n i ) xi 。 f (x) dx lim 0 i 1 f ( a T 2、定积分的应用 ①计算平面图形的面积： 由连续曲线 y f ( x)( f ( x) 0) ，直线 x a, x b(a b) 及 S b x 轴所围成的曲边梯形面积为 f ( x) dx； a ②计算旋转体的体积： 由连续曲线 y f ( x) ，直线 x a, x b(a b) 及 x 轴所围成的 曲边梯形绕 x 轴旋转一周所成立体的体积为 V b [ f (x)] 2 dx ； a ③计算平面曲线的弧长：设曲线弧由直线坐标方程 y f ( x)(a x b) 给出，其中 f ( x) 在 [ a,b] 上具有一阶连续导数，则曲线弧长 l b 1 y 2 dx ； a x x(t) t ) 给出，其中 x(t), y(t ) 在 [ , ] 上具有连续 设曲线弧由参数方程 y , ( y(t ) 导数，则曲线弧长 l x t 2 y t 2 dt ； ( ) ( ) 设曲线弧由极坐标方程 r r ( )( ) 给出，其中 r ( ) 在 [ , ] 上具有连续导 数，则曲线弧长 l r ( ) 2 r ( )2 d 。 ④定积分的近似计算： 用分点 a x0 x1 xn b 将区间 [ a, b] 分成 n 个长度相等的小区间，每个小区 间长度为 xi b na xi a (b a) ni , xi 1 xi b na yi f ( xi ) 。 b f ( x)dx b na ( y0 y1 yn 1 ) 矩形法公式： a b b na ( y1 f ( x)dx y2 yn ) a b 梯形法公式： f ( x)dx b n a 21 ( y0 yn ) y1 y2 yn 1 a 抛物线法公式： b b a f (x) dx ( y0 yn ) 2( y2 y4 yn 2 ) 4( y1 y3 yn 1 ) 3n a 实验使用的函数与命令 1、符号积分命令 int int(fun) ：求函数 fun 的不定积分； int(fun,var) ：求函数 fun 关于变量 var 的不定积分； int(fun, var, a ， b,)：求函数 fun 的在 [a,b] 间的定积分或广义积分； 2、数值计算定积分 quad， quad1，trapz ① quad(fun,a,b,tol) ： fun 为被积的函数名， a, b 微积分上下限， tol 为精度，若缺省，其 缺省值为 1.0e-6； ② quad(fun,a,b,tol,trace) ：参数 fun,a,b,tol 用法与上面相同， 而输入第五个非零参数 trace， 是对积分过程通过被积函数上的图像进行跟踪。 对于 quad 使用自适应步长 Simpson 法，而 quadl 的调用格式与 quad 一致，但它使用 Lobbato 算法，其精度和速度要比 quad 的精度高。 ③ Z = trapz(X,Y) ：用于进行梯形积分，精度低，适用于数值函数和光滑性不好的函数。 其中 X 表示积分区间的离散化变量， Y 与 X 同维的向量，表示被积函数， Z 返回积分的近 似值。 quad， quad1， trapz 都不能用于广义积分，此外由于数值积分的特点，对一些假奇异 积分也不能直接求解。 3、 Inline 实验指导 例 1 用符号