matlab多重积分编程,多重积分的MATLAB实现.pdf

多重积分的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 =２．

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 =