积分和微分方程组.docVIP

下载 第11章 积分和微分方程组 在有效的 MATLAB命令帮助下，可以求解出定积分和普通微分方程的数字解并绘制出其 图形。 11.1 积分 在MATLAB中能求解如下形式的定积分并给出数字解： 有许多方法都可以能够解决积分问题 (又叫做求面积 )。如果要用 MATLAB监控整个计算 过程，可以使用quad命令。同样能计算出被积函数 g的值，并且让 MATLAB使用梯形规则和 trapz命令计算出积分。当只有离散的数据点和被积函数的数学表达式为未知时，这种方法 是非常有效的。 命令集107 定积分计算  trapz(x,y) trapz(y) trapz(x,A) Z=trapz(x,A,dim) cumtrapz(A,dim) quad(fcn,a,b)  计算出函数x的积分并将结果返回到y。向量x和y有相同的长度， (xi, yi)代表曲线上的一点。曲线上点的距离不一定相等，x值也 不一定有序。然而，负值间距和子区间被认为是负值积分。 计算方法同上，但x值间隔为1。 将A中每列的值带入 x的函数算出其积分，并返回一组包含 积分结果的向量。 A的列向量必须和向量 x的长度相同。 在矩阵im指定的维内进行数据积分。如果给定向量 x， 则x的长度必须与size(A, dim)相同。 返回大小和 A相同的数组，包含的是将矩阵 A进行梯形积分 的累积值。如果dim已给定，则在dim维内进行计算。 返回在区间[a, b]上g的积分近似值。字符串fcn包含一个与g相 对应的MATLAB函数名，也就是预定义函数或者是 M文件。 这个函数接收一个向量参数，并返回一个向量结果。 MATLAB利用辛普森规则执行递归的积分，计算误差为10－3。 quad(fcn,a,b,tol) 求g的积分近似值，其相对误差由参数 tol定义。否则，计算 过程同上。 quad(fcn,a b,tol 求g的积分近似值，其相对误差由参数 tol所定义。如果参数 pic) quad(...,trace) pic是非零值，则在图形中显示求值的点。 如果trace是非零值，则画出积分图形。 152  MATLAB 5 手册  下载  quad8(...)  可以与quad一样用于相同的参数组合并返回相同的结果， 但使用更高精度的方法。因此，如果被积函数的导数在某一 1 区间内是不定的，例如： q???∫0 s i n x ,使用此命令将会更好 一些。quad和quad8 都要求被积函数在整个区间里是有限的。 dblquad(f, min1,计算双变量函数f的二重积分。函数中的第一个自变量用于 max1,min2,max2, tol,trace,order) 内层积分。内层积分在 min1 和max1之间进行，外层积分在 min2和max2之间进行。变量 tol指定相对误差。 trace的使用 方法与quad相同。根据字符串 order，对于相同的访问， dblquad能选择使用quad、quad8和许多用户定义的积分 方法，并返回与quad相同的变量。 输入quaddemo可以看到一个演示实例。 ■ 例11.1 下面用不同的方法来计算下列积分： (a) 使用trapz命令。首先创建一个有 x值的向量。用5和10两个值进行计算： 然后创建x的函数y: 现在计算出积分值： 返回 (b) 使用quad命令。首先在M文件中创建函数。此文件 integrand.m包含函数，如下： 首先以标准误差计算积分，然后再以指定误差计算积分。 下载  第11章 积分和微分方程组  153  给出 (c) 使用quad8命令：使用在(b)中创建的M文件,然后输入： 这是MATLAB所能给出的最精确的结果。 (d) 使用cumtrapz命令能很容易地计算出不同区间的积分。 (e) 计算二重积分： 如图11-1。首先创建一个包含函数 M文件:integrand2.m: 然后用quad命令计算对于固定的 x值在y方向的一些积分值： 现在已计算出在y方向的15个积分值。trapz命令能使用这些值来计算二重积分： 输入下列语句可以得到一个积分区域的图形： 154  MATLAB 5 手册  下载  结果如图11-1所示。命令mesh和view定义在13.5节中。 图11-1 函数e-x2-y2在区间[0，1]×[0，1]的上的图形 ■  不定积分  f(t)dt 不能使用上面的命令来计算。  M AT L A B 中的数学符号工具箱和 MATLAB的编辑器能提供处理这些积分的命令。 11.2 常微分方程组 下面来研究常微分方程系统 ODE, 该系统处理的是初始值已知的一阶微分方程。在本节中 主要讨论这种类型的微分方程，同时也会举出两个有关边界值问题的