主要内容

整合领域内部的奇点

此示例显示如何拆分集成域以在边界上放置奇点。

使用匿名功能定义积分

积分的复合值积分

- 1 1 - 1 1 1 X + y D. X D. y

有一个奇点的时候x = y = 0并且一般来说,在线上奇异y = -x.

使用匿名功能定义此集成。

fun = @(x,y)((x + y)。^( -  1/2));

融合在广场上

整合乐趣在指定的方框上 - 1 ≤. X ≤. 1 - 1 ≤. y ≤. 1

格式q = Integral2(乐趣,-1,1,-1,1)
警告:非有限情况。整合不成功。奇点可能。
Q = Nan + Nani

如果集成区域内部存在奇异值,则集成无法收敛并返回警告。

将集成域拆分为两个三角形

您可以通过将集成域分成互补部分并将较小的集成拆分将集成域拆分为互补部分来重新定义。通过将奇点放置在域的边界上,避免集成错误和警告。在这种情况下,您可以将方形集成区域分成两个三角形沿着奇异线y = -x.并添加结果。

Q1 = Integral2(Fun,-1,1,-1,@(x)-x);Q2 = Integral2(Fun,-1,1,@(x)-x,1);Q = Q1 + Q2
q = 3.771236166328258  -  3.771236166328256i

当奇异值在边界上时,集成成功。

这种积分的确切值是

8. 2 3. 1 - 一世

8/3 * SQRT(2)*(1-I)
ans = 3.771236166328253  -  3.771236166328253i

也可以看看

||

相关话题