如何在Matlab中使用Surf绘制球面函数

Question

我正在尝试在下面绘制球面函数，但结果却是错误的

这是我正在使用的代码

x1 = [-10:1:10];
x2 = [-10:1:10];
y = zeros(1,21);

for i = 1:21
    y(i) = sphere([x1(i) x2(i)]);
end
Y = meshgrid(y);
surf(x1,x2,Y);
colormap hsv;

sphere.m

function [y] = sphere(x)
d = length(x);
sum = 0;
for i = 1:d
    sum = sum + x(i)^2;
end
y = sum;
end

Answer 1

sphere(10)

它是MatLab内置函数。

请负责任地享受。

如果需要查看源代码，请使用：当sphere函数不在路径help sphere时， edit sphere或help sphere 。

Answer 2

为了完整起见，您的代码无法正常工作，因为您只对[-10,10]中的一些\\ \\对（x，x）上的函数进行求值，因此不会涵盖整个域。 它可以这样工作：

x1 = [-10:1:10];
x2 = [-10:1:10];
y = zeros(1,21);

for i = 1:21
    for j=1:21
        Y(i,j) = sphere([x1(i) x2(j)]);
    end
end
surf(x1,x2,Y);
colormap hsv;

或更快（因为出于计算时间的原因，您应始终避免不必要的循环）：

x1 = meshgrid([-10:1:10]);
x2 = x1';

Y = x1.^2+x2.^2;

surf(x1,x2,Y)