Как построить функцию сферы?

Question

Я пробовал рисовать сферу to look like this.

This earlier question содержит код, в котором я пробовал заговор.

Ответ на этот предыдущий вопрос, но теперь у меня возникли проблемы с построением функции сферы, как показано на рисунке, указанном ссылкой выше.

Редактировать Это функция сфера я использую:

function ph = sphereFN(x) 
    ph = sum(x.*x, 2); 

end 

Редактировать Результат должен выглядеть примерно так:

NB: Я изменил функцию из сферы на сферу FN, чтобы избежать конфликта с сферой матлаба.

Answer 1

Ваш sphereFN не может делать то, что вы ожидаете, что делать. Похоже, что вы пытаетесь осуществить тест sum of different powers на основе этой функции:

Однако, ваша реализация, как:

ph = sum(x.*x, 2); 

не похожа на эту функцию. Как об использовании этого в качестве отправной точки:

x = linspace(-1,1,25); 

for I=1:size(x,2) 
    for J=1:size(x,2) 
     s(I,J) = abs(x(I))^2 + abs(x(J))^3; 
    end 
end 

[xx,yy] = meshgrid(linspace(-1,1,25)); 

surfc(xx,yy,s) 

Как ваш, в сущности, описывающей 2D проблемы, эта грубая реализация функции будет достаточно. Вы можете поместить вложенный for -loops в функцию под названием sum_of_different_powers_2D (например) и вызвать ее, передав вектор x.

Редактировать

Вы можете получить форму, которая больше напоминает предполагаемому поверхность заменить команду внутри вложенной for -loops с:

s(I,J) = abs(x(I))^2 + abs(x(J))^2; 

Это будет напоминать эту функцию:

Редактировать 2

Важно понимать, что n в приведенном выше уравнении является измерением вашей проблемы. Что, как вы показали, это 2.

Редактировать 3

Я рекомендую использовать эту функцию:

function ph = sphereFN(x) 

for I=1:size(x,2) 
    for J=1:size(x,2) 
     ph(I,J) = abs(x(I))^2 + abs(x(J))^2; 
    end 
end 

end 

Answer 2

Matlab имеет встроенную функцию сферы, и в документации показано, как ее построить.

http://www.mathworks.com/help/matlab/ref/sphere.html