Calculating a surface integral over a regular shape

2019 11 22
1 个回答
回答已采纳
更新时间：2019 11 24
7 次查看（30 天）

0 个投票

let's suppose a function given in polar coordinate F(r,phi) and our purpose is to calculate the surface integral, say F(r,phi)dA over the region S defined by
S={ |Z|<b , |z-z0|>a } where Z=r*exp(1j*phi). it means the regions between the circles |Z|=b , |z-z0|=a . However, F(r,phi) has singularities inside the circle |z-z0|=a . therefore, we are not able to use integral(over circle |Z|<b )-integral(over circle |z-z0|<a) .
your help and consideration are much appreciated.

3 个评论
显示 1更早的评论 隐藏 1更早的评论

David Goodmanson
David Goodmanson 2019-11-24
Hi Oussama,
I am assuming that z and Z are basically the same thing, is that correct? Are z0 and 'a' such that the 'a' circle is totally contained in the b circle? Or the other way round? Is z0 real, or can it be complex?
Oussama GASSAB
Oussama GASSAB 2019-11-24
yes z is the same as Z (typing mistake). yes exactly the z0 and 'a' are chosen such that the circle |Z-z0|=a is totally inside the bigger circle |Z|=b.

请先登录，再进行评论。

请先登录，再回答此问题。

 采纳的回答

Oussama GASSAB
Oussama GASSAB 2019-11-24
编辑：Oussama GASSAB 2019-11-24

0 个投票

i have solved by using the following, where F(x,y)=F(r,phi). i have used x-y coordinate and quad2d.
the result was so accurate.
y1= @(x) -sqrt(b^2-x.^2) ; y2=@(x) sqrt(b^2-x.^2) ;
y3= @(x) -sqrt(a^2-(x-z0).^2) ; y4= @(x) sqrt(a^2-(x-z0).^2) ;
intg1 = quad2d(F,-b,z0-a,y1,y2) ;
intg2 = quad2d(F,z0-a,z0+a,y1,y3) ;
intg3 = quad2d(F,z0-a,z0+a,y4,y2) ;
intg4 = quad2d(F,z0+a,b,y1,y2) ;
integ_total=intg1+intg2+intg3+intg4 ;

更多回答（0 个）

类别

在 帮助中心 和 File Exchange 中查找有关 Numerical Integration and Differentiation 的更多信息

产品

版本

R2018b

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by