般区域重重积分MATLAB计算方法.docVIP

[【原创】] 一般区域二重、三重积分MATLAB计算方法 这里讨论的计算方法指的是利用现有的MATLAB函数来求解，而不是根据具体的数值计算方法来编写相应程序。目前最新版的2009a有关于一般区域二重积分的计算函数quad2d（详细介绍见/viewthread.php?tid=873479），但没有一般区域三重积分的计算函数，而NIT工具箱似乎也没有一般区域三重积分的计算函数。本贴的目的是介绍一种在7.X版本MATLAB（不一定是2009a）里求解一般区域二重三重积分的思路方法。需要说明的是，上述链接里已经讨论了一种求解一般区域二重三重积分的思路方法，就是将被积函数“延拓”到矩形或者长方体区域，但是这种方法不可避免引入很多乘0运算浪费时间。因此，新的思路将避免这些。由于是调用已有的MATLAB函数求解，在求一般区域二重积分时，效率和2009a的quad2d相比有一些差距，但是相对于延拓函数的做法，效率大大提高了。下面结合一些简单例子说明下计算方法。譬如二元函数f(x,y) = x*y，y从sin(x)积分到cos(x)，x从1积分到2，这个积分可以很容易用符号积分算出结果 1. syms x y 2. int(int(x*y,y,sin(x),cos(x)),1,2) ]结果是 -1/2*cos(1)*sin(1)-1/4*cos(1)^2+cos(2)*si