多重积分的MATLAB实现
Mg 分析实验二
实验二 多重积分
目的和要求
1. 学习MATLAB软件中与二重积分和三重积分相关的命令与编程
2. 利用MATLAB软件计算二重积分和三重积分
预备知识
1. 符号积分
利用int 命令进行二重积分、三重积分计算。
格式: int(int(f,y,y1(x),y2(x)),x,a,b)
int(int(int(f,z,z1(x,y),z2(x,y)),y,y1(x),y2(x)),x,a,b)
2. 数值积分
二重数值积分的指令:dblquad(fun,xmin,xmax,ymin,ymax)
三重(闭型)数值积分指令:triplequad(fun,xmin,xmax,ymin,
ymax,zmin,zmax)或用三个quadl函数
实验内容
1.
1.
11.. 二重积分的计算
2 2
例1 计算二重积分 (x + y − x)dxdy,其中 D是由直线y= 2,y= x及y= 2x
∫∫
D
所围成的闭区域.
2 y 2 2
解 该二重积分可以化为二次积分 dy (x + y − x)dx,输入命令
∫ ∫y
0
2
int(int(x^2+y^2-x,x,y/2,y),y,0,2)
得结果ans=13/6.
2π π sinx
例2 计算积分I= ∫π dy∫y−π x dx.
解 输入命令
int(int(sin(x)/x,x,y-pi,pi),y,pi,2*pi)
得结果ans =2.
Mg 分析实验二
若不借助于计算机计算,需要交换积分次序.
例3 求 2 2 .
∫∫sin(π(x + y ))dxdy
2 2
x +y ≤1
解 积分区域用不等式可以表示成 2 2 ,二重积分可
−1≤ x≤ 1,− 1− x ≤ y≤ 1− x
1 1−x2
2 2
化为二次积分 dx sin(π(x + y )dy,输入命令
∫ ∫ 2
−1 − 1−x
int(int(sin(pi*(x^2+y^2)),y,-sqrt(1-x^2),sqrt(1-x^2)),x,-1,1)
由输出结果可以看出,结果中仍带有int,表明matlab求不出这一积分的值.采用极坐标
2π 1
2
可化为二次积分 da rsin(πr )dr,输入命令∫ ∫
0 0
int(int(r*sin(pi*r^2),r,0,1),a,0,2*pi)
可得结果为ans=2.
注:二重积分也可以先进行符号运算,再得到数值解。
4 y
4
例44 ∫∫x dxdy,其中R={(x,y)| 0≤ x≤1,1≤ y≤ 2}。
R
syms xy; %定义两个符号变量
a=int(int(x^y,x,0,1),y,1,2) %积分
a =