Сетка произвольной области в 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 и т. д., чтобы продемонстрировать это на графике ...