Сохранять структуру meshgrid при подмножестве meshgrid

Question

Когда я подмножу массив, построенный из meshgrid, я не могу решить, как сохранить его структуру meshgrid. Таким образом, вы не можете использовать его при вызове mesh или surface. Я продемонстрирую это в моем примере построения единичной сферы.

Возможные альтернативные названия для этого вопроса:

1. Как сделать верхнюю половину meshgrid области с нуля в Matlab?
2. Как вы используете сетку для построения подмножества meshgrid?

Это мотивировано, используя следующую игрушку пример построения сферы единичного радиуса в Matlab с нуля, так что она подобна той, генерируемой:

[x, y, z] = sphere(100) 
mesh(x, y, z) 

Используя уравнение для сферы:

Определение meshgrid и z быть:

x = linspace(-1, 1, 201); 
y = linspace(-1, 1, 201); 
[x, y] = meshgrid(x, y); 
z = sqrt(1 - x.^2 - y.^2); 

Пока все хорошо, кроме z, принимает мнимые значения, где сфера не существует над плоскостью xy, то есть где-либо за пределами единичного круга.

Вызов mesh теперь возвращает ошибку:

>> mesh(x, y, z) 
Error using mesh (line 58) 
X, Y, Z, and C cannot be complex. 

Таким образом, логический шаг должен удалить все комплексные значения:

% get logical vector index where real z is 
LI = z == real(z) 

x = x(LI) 
y = y(LI) 
z = z(LI) 

Но теперь x, y и z больше не являются 3d-массивов, а также по вызову mesh дает еще одну ошибку:

>> mesh(x, y, z) 
Error using mesh (line 58) 
Z must be a matrix, not a scalar or vector. 

Итак, я вообще не знаю, как сохранить структуру meshgrid при подмножестве данных. Следовательно, я не могу сгенерировать верхнюю половину этой сферы с нуля.

Answer 1

В общем, вы можете «исключить» значения из графика при сохранении структуры матрицы, используя NaN значение. В вашем случае попробуйте следующее:

LI = z == real(z); 
z(~LI) = NaN; 
mesh(x,y,z);