第11章 积分和微分方程组.pdfVIP

下载 第11章 积分和微分方程组 在有效的M AT L A B命令帮助下，可以求解出定积分和普通微分方程的数字解并绘制出其 图形。 11.1 积分 在M AT L A B 中能求解如下形式的定积分并给出数字解： 有许多方法都可以能够解决积分问题 (又叫做求面积 ) 。如果要用M AT L A B监控整个计算 过程，可以使用 q u a d命令。同样能计算出被积函数 g 的值，并且让M AT L A B使用梯形规则和 t r a p z命令计算出积分。当只有离散的数据点和被积函数的数学表达式为未知时，这种方法 是非常有效的。 命令集1 0 7 定积分计算 t r a p z ( x , y ) 计算出函数x 的积分并将结果返回到y 。向量x和y 有相同的长度， (xi, yi)代表曲线上的一点。曲线上点的距离不一定相等，x值也 不一定有序。然而，负值间距和子区间被认为是负值积分。 t r a p z ( y ) 计算方法同上，但x 值间隔为 1。 t r a p z ( x , A ) 将A 中每列的值带入x 的函数算出其积分，并返回一组包含 积分结果的向量。A 的列向量必须和向量x 的长度相同。 Z = t r a p z ( x , A , d i m ) 在矩阵A 中d i m指定的维内进行数据积分。如果给定向量 x ， 则x 的长度必须与size(A, dim)相同。 c u m t r a p z ( A , d i m ) 返回大小和A相同的数组，包含的是将矩阵 A进行梯形积分 的累积值。如果d i m 已给定，则在d i m维内进行计算。 q u a d ( f c n , a , b ) 返回在区间[a, b]上g的积分近似值。字符串f c n包含一个与g相 对应的M AT L A B 函数名，也就是预定义函数或者是 M文件。 这个函数接收一个向量参数，并返回一个向量结果。 －3 M AT L A B利用辛普森规则执行递归的积分，计算误差为1 0 。 q u a d ( f c n , a , b , t o l ) 求g 的积分近似值，其相对误差由参数 t o l定义。否则，计算 过程同上。 quad(fcn,a b,tol,求g 的积分近似值，其相对误差由参数 t o l所定义。如果参数 p i c ) p i c 是非零值，则在图形中显示求值的点。 q u a d (. . ., t r a c e ) 如果t r a c e是非零值，则画出积分图形。 1 5 2 M A TLAB 5 手册 下载 q u a d 8 ( . . .) 可以与 q u a d一样用于相同的参数组合并返回相同的结果， 但使用更高精度的方法。因此，如果被积函数的导数在某一 1 q s i n x 区间内是不定的，例如： 0 dx ,使用此命令将会更好 一些。q u a d和quad8 都要求被积函数在整个区间里是有限的。 dblquad(f, min1,计算双变量函数f 的二重积分。函数中的第一个自变量用于