Numerical integration over a circular domain

Question

Kobye 2013-10-23

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/91217-numerical-integration-over-a-circular-domain

回答： Jan 2013-10-23

I have an analytic function defined over x=[0,100] and y=[0,100]. I do not know the precise workings of the function, just that it is smooth, takes in inputs (x,y) and returns a scalar.

It is quite straightforward to numerically integrate this function over a rectangular domain using the function integral2. However, I would like to integrate the function over a circle. For simplicity's sake let's have it centered at (x,y)=(50,50) with a radius of 10.

How would I be able to achieve this? Remember I cannot modify the function by substituting x = 50+10*cos(theta) and y=50+10*sin(theta).

An example function below.

func = @(x,y) x^3*y+x*y^2;
for i=1:100
    for j=1:100
        data(i,j)=func(i,j);
    end
end
surf(1:100,1:100,data)

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Jan 2013-10-23

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/91217-numerical-integration-over-a-circular-domain#answer_100608

I would try to transform the scalar field into polar coordinates and do numerical integration over a square window there. Of course you have to take care of the translation prior to coordinate transformation to make sure, that the circles center is located in the origin.