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

一

使用两次一重积分

%%

二重积分

f= (x,

ｙ

)exp

(ｓ

in(

ｘ)

)*ln

(ｙ

),y

从

5*

ｘ积分到ｘ＾

2,

ｘ从

1

０积

分到

2

０

1

(

7

、

X

后版本才有此函数

qu

ａ

d2d)

y1

＝ｑ

ua

ｄ

2d(

(

x,y

)

ex

ｐ

(si

ｎ

(x

)

)

、

*log(y)

，

1

０

,20,(x)5*

ｘ

,(x)x

、＾２

)

２

y2

＝ｑ

u

ａ

dl((x)

a

ｒ

r

ａｙ

fu

ｎ

((

ｘ

)

ｑ

uadl

(

(y)ex

ｐ

(sin

(ｘ

))

、

*

ｌｏ

g

(ｙ

)

，

5*x,x

、

^

２

),

ｘ

),

１

0,2

０

)

3

y3

＝

dblquad

(

(x,y)

ｅｘｐ(ｓｉ

n(x)

)、

*log(

ｙ

)

、

*

(

y>=5*x

＆

y<=x

、

^2)

，

10

，

20,

５

0,

４

00)

详细请瞧吴鹏老师得文章

二

使用ｄ

bl

ｑ

uad

函数

q =

dbl

ｑｕ

ad(fu

ｎ

,xmi

ｎ，ｘ

max,ym

ｉｎ

,yma

ｘ，

t

ｏｌ，

m

ｅｔｈ

od)

该函数求

f(x,y)

在

[

ａ，

b]×

［

c,d]

区域上得二重定积分。参数ｔ

ol,t

ｒａｃ

e

得用法

与函数

qua

ｄ完全相同。

例

8

－

5

计算二重定积分

(1)

建立一个函数文件

fxy

、

m:

func

ｔ

io

ｎ

f

＝ｆ

xy(

ｘ

,y)

globa

ｌ

k

ｉ

;

ｋｉ＝ｋ

i+1

；

%ki

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

ｆ＝ｅ

xp(

－ｘ、

^2/2)

、

*s

ｉ

n(x

、＾

2+y

)

;

(

２

)

调用ｄ

blqua

ｄ函数求解。

ｇ

lo

ｂ

a

ｌ

ｋ

i;ki=

０；

Ｉ

=dbl

ｑ

uad('

ｆ

x

ｙ

',-2

，

2,-1,1)

k

ｉ

Ｉ

＝

1

、

574

４

931

８

97

４４

94

ｋ

i

=

１

0

３８

来源

精通

MATLAB

科学计算

一书王正林

,

龚纯

,

何倩编写

,

电子工业出版社

三

复合辛普森公式

(

矩形积分区域

)

fun

ｃ

tion q=DblSimpson(f

，

a

，

A

，ｂ

,B,m

，ｎ

)

ｉｆ(

m==1

&& n==1

)

％辛普森公式

q

＝

((B-b)*(

Ａ－

a

)

/9

)＊(

su

ｂ

s(sym(f),

ｆ

in

ｄ

s

ｙ

m(sym(

ｆ

)),{a

，

b})+

、、、

ｓ

u

ｂ

s

(

s

ｙ

m

(

f)

，

fin

ｄ

sym

(

sym(f

)

),{a,B})+

、、、

subs(s

ｙｍ(

f

)

,

ｆ

i

ｎｄ

sym(

ｓ

ym(f)

)，｛

A,b})+

、、、

subs(sym(

ｆ

),fi

ｎｄ

s

ｙ

m(s

ｙｍ

(

ｆ

)),

｛Ａ

,

Ｂ

})+

、、、

4

＊ｓ

u

ｂｓ

(sy

ｍ

(f),find

ｓｙ

m(

ｓ

ym(f)

)

,

｛

(A-a)

／

2,b})+

、、、

4*subs(sym(f

)

,fin

ｄ

s

ｙ

m

(ｓ

ym(f

)

),{(A-

ａ

)

／

2,B})+

、、、

4

＊

s

ｕ

b

ｓ

(s

ｙ

m

(ｆ

),fin

ｄ

sym

(

sym(f

)

),{a,

(

B-b)/2})+

、、、

4*

ｓｕ

bs(

ｓ

ym(

ｆ

),

ｆ

in

ｄｓ

ym(sy

ｍ

(

ｆ

)),

｛

A,(

Ｂ

-b)/

２

})

＋、、、

1

６

*s

ｕ

bs

(

sym(f),find

ｓ

ym(sym

(ｆ

)),{(A-a)/2

，

(B

－

b)

／２

}));

els

ｅ

%

复合辛普森公式

ｑ

=0;

for i=

０

:

ｎ－

1

for j=0:

ｍ

-

１

x

＝

a+2*i

＊

(A-

ａ

)/

２

/n;