主要内容

Integral2.

数值计算二重积分

描述

例子

= integral2 (乐趣xminXmax.yminymax近似于功能的积分z =乐趣(x, y)在平面区域xminxXmax.ymin(x)yymax (x)

例子

= integral2 (乐趣xminXmax.yminymax名称,价值指定一个或多个附加选项名称,价值对论点。

例子

全部收缩

考虑该功能

f x y 1 x + y 1 + x + y 2

此函数未定义 x y 是零。Integral2.当奇点在积分边界上时性能最好。

创建匿名函数。

fun = @(x,y)1 ./(sqrt(x + y)。*(1 + x + y)。^ 2)
有趣的=function_handle与价值:@ (x, y) 1. /(√(x + y) * (1 + x + y) ^ 2)

集成在界限的三角形区域上 0 x 1 0 y 1 - x

Ymax = @(x) 1 - x;q = integral2(有趣,0 1 0,ymax)
q = 0.2854

定义功能

f θ. r r r cos θ. + r θ. 1 + r cos θ. + r θ. 2

fun = @(x,y)1 ./(sqrt(x + y)。*(1 + x + y)。^ 2); polarfun = @(theta,r) fun(r.*cos(theta),r.*sin(theta)).*r;

定义上限的函数 r

Rmax = @(theta) 1 /(sin(theta) + cos(theta));

对边界区域积分 0 θ. π. / 2 0 r r 一个 x

Q = Integral2(PolarFun,0,Pi / 2,0,Rmax)
q = 0.2854

创建匿名参数化函数 f x y 一个 x 2 + b y 2 有参数 一个 3. b 5

= 3;b = 5;(x,y) a*x。^ 2 + b * y ^ 2;

求这个区域的积分 0 x 5 - 5 y 0 .指定'迭代'方法和大约10位有效数字的准确性。

格式q = integral2(乐趣0 5 5 0,“方法”'迭代'...'ABSTOL',0,“RelTol”1、平台以及)
Q = 1.66666666666667e + 03

输入参数

全部收缩

被指定为函数句柄的Integrand,定义了在平面区域上集成的功能xminxXmax.yminx)≤yymaxx).这个函数乐趣必须接受两个相同大小的两个数组并返回一个相应值的数组。它必须执行元素明智的操作。

数据类型:function_handle

下限x,指定为有限或无限的实标量值。

数据类型:双倍的|

的上限x,指定为有限或无限的实标量值。

数据类型:双倍的|

下限y,指定为有限或无限的实标量值。您可以指定ymin作为函数的句柄(…的函数x)在整合非正天区域时。

数据类型:双倍的|function_handle|

的上限y,指定为有限或无限的实标量值。你也可以指定ymax作为函数的句柄(…的函数x)在整合非正天区域时。

数据类型:双倍的|function_handle|

名称-值参数

指定可选的逗号分隔的对名称,价值论点。的名字参数名和价值是相应的价值。的名字必须出现在引号内。您可以以任何顺序指定多个名称和值对参数name1,value1,...,namen,valuen

例子:'ABSTOL',1E-12将绝对误差容错设置为大约12位小数精度。

绝对误差容差,指定为逗号分隔对'ABSTOL'和一个非负实数量。Integral2.使用绝对误差公差来限制绝对误差|的估计- - - - - -|,是积分的计算值是(未知的)精确值。Integral2.如果降低绝对误差容忍度,可能会提供更多的小数点精度。默认值为1平台以及

请注意

abRELTOL.一起工作。Integral2.可能满足绝对误差容差或相对误差容错,但不一定是两者。有关使用这些公差的更多信息,请参阅提示部分。

例子:'ABSTOL',1E-12将绝对误差容错设置为大约12位小数精度。

数据类型:双倍的|

相对容错,指定为逗号分隔对组成“RelTol”和一个非负实数量。Integral2.使用相对误差公差来限制相对误差|的估计- - - - - -| / ||,是积分的计算值是(未知的)精确值。Integral2.如果减少相对误差容差,可能会提供更高的精度数字。默认值为1 e-6

请注意

RELTOL.ab一起工作。Integral2.可能满足相对误差容差或绝对误差容错,但不一定都是两者。有关使用这些公差的更多信息,请参阅提示部分。

例子:e-9 RelTol, 1将相对误差公差设置为大约9位有效数字。

数据类型:双倍的|

集成方法,指定为逗号分隔对组成“方法”和下面描述的方法之一。

集成方法 描述
“汽车” 对于大多数情况而言,Integral2.使用“瓦”方法。它使用了'迭代'当任何一个积分极限都是无限时。这是默认方法。
“瓦” Integral2.将积分区域转换为矩形形状,并根据需要将其细分为更小的矩形区域。积分极限必须是有限的。
'迭代' Integral2.调用积分执行迭代积分。外部积分进行评估xminxXmax..内部积分进行评估ymin(x)yymax (x).整合限制可能是无限的。

例子:“方法”,“瓦”指定瓷砖集成方法。

数据类型:char|字符串

提示

  • Integral2.功能尝试满足:

    ABS(Q-Q)<= MAX(ABSTOL,RELTOL * ABS(Q))
    在哪里是积分的计算值是(未知的)精确值。绝对公差和相对公差提供了一种权衡精度和计算时间的方法。通常,相对公差决定了积分的精度。然而,如果ABS(Q)足够小时,绝对公差决定了积分的精度。通常应该同时指定绝对公差和相对公差。

  • '迭代'当函数在积分区域内存在不连续时,该方法更有效。然而,最好的性能和精度出现在您在不连续点分割积分并对多个积分结果求和的时候。

  • 在整合非正天区域时,最好的性能和准确性发生yminymax,(或两者)都是函数句柄。避免将被积函数值设为零来在非矩形区域进行积分。如果您必须这样做,请指定'迭代'方法。

  • 使用'迭代'方法时yminymax(或两者)是无界功能。

  • 当副主义匿名函数时,请注意参数值持续函数句柄的寿命。例如,功能= @(x,y) x + y + a使用的价值一个当时乐趣是创造的。如果您稍后决定更改值一个,您必须将匿名功能重新定义为新值。

  • 如果你指定了积分的单精度极限,或者乐趣返回单精度结果,您可能需要指定更大的绝对和相对误差公差。

参考文献

[1] L.F. Shampine "MATLAB中的矢量自适应求积®,“计算和应用数学杂志, 211, 2008,页131 - 140。

[2] L.F. Shampine, "二维求积的MATLAB程序。应用数学与计算。第202卷,第1期,2008年,266-274页。

扩展功能

在R2012A介绍