matlab trapz二重积分函数_matlab二重积分

一

使用两次一重积分

%%

二重积分

f= @(x,y)exp(sin(x))*ln(y)

，

y

从

5*x

积分到

x^2

，

x

从

10

积分到

20

1

后版本才有此函数

quad2d)

y1=quad2d(@(x,y) exp(sin(x)).*log(y),10,20,@(x)5*x,@(x)x.^2)

2

y2 =quadl(@(x) arrayfun(@(x) quadl(@(y)exp(sin(x)).*log(y),5*x,x.^2),x),10,20)

3

y3 = dblquad(@(x,y)exp(sin(x)).*log(y).*(y>=5*x & y<=x.^2),10,20,50,400)

详细请看吴鹏老师的文章二

使用

dblquad

函数

q = dblquad(fun,xmin,xmax,ymin,ymax,tol,method)

该函数求

f(x,y)

在

[a,b]×[c,d]

区域上的二重定积分。参数

tol

，

trace

的用法与函数

quad

完全相同。

例

8-5

计算二重定积分

(1)

建立一个函数文件：

function f=fxy(x,y)

global ki;

ki=ki+1; %ki

用于统计被积函数的调用次数

f=exp(-x.^2/2).*sin(x.^2+y);

(2)

调用

dblquad

函数求解。

global ki;ki=0;

I=dblquad('fxy',-2,2,-1,1)

ki

I =

ki = 1038