Matlab中任意域的Meshgrid

Question

我需要在一个复杂的域中评估一个 function 并且只在这个域中。 它之外有一些值，但我不希望它们被绘制。 我不知道如何创建一个合适的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等的一些任意值来在 plot 中演示它......