Matlab中任意区域的网格划分

Question

我需要在一个复杂的域中计算一个函数，并且仅在这个域中。它之外还有一些值，但我不想绘制它们。我不知道如何创建一个合适的meshgrid来做到这一点。

这里有一张我正在谈论的域名的图片：

​

​

Answer 1

这样如何：

[x, y]=meshgrid(linspace(0,10,1000),linspace(-5,5,1000));
r=sqrt(x.^2+y.^2);
theta=atan2(y,x);

a=4;
b=1;

I = r<2*a*cos(theta) & theta<pi/2 & theta>-pi/2 & r>b; % this is what you want

imagesc([0 10],[-5 5],I)

​

​

因此，根据问题中设置的逻辑条件，I是域中的1和外部的0的矩阵。我选择了a、b等的一些任意值来在图中演示它。