Matlab，錯誤使用衝浪，X，Y，Z和C不能複雜

Question

我運行下面的代碼，但得到X，Y，Z和C不能是複雜的錯誤，任何想法是什麼錯誤？

k=1; 
u = linspace(0,2*pi,72); 
v = [-3:.2:-1,1:.2:3]; 
[U,V] = meshgrid(u,v); 
r=sqrt((4*V.^-k)./(cos(U).^2+k*sin(U).^2)); 
X = r.*cos(U); 
Y = r.*sin(U); 
Z = V; 

這是我想要得到的圖像：你得到的是複雜的錯誤，因爲r是複數

function simple_math_functions_animation1 
clc, close all, clear all 

hf1=figure(1);hold on,grid on,axis equal, view([1 -1 1]) 
set(hf1,'Color','w');set(hf1,'Position',[300, 600, 500, 400]); 
xlabel('x');ylabel('y'),zlabel('z'); 

k=1; 
u = linspace(0,2*pi,72); 
v = [-3:.2:-1,1:.2:3]; 
[U,V] = meshgrid(u,v); 
r=sqrt((4*V.^-k)./(cos(U).^2+k*sin(U).^2)); 
X = r.*cos(U); 
Y = r.*sin(U); 
Z = V; 

surf(X,Y,Z,'EdgeColor',[0.5 1. 0.2],'FaceColor',[1 0.2 0.8],'FaceAlpha',0.6); 

XYZ=[reshape(X,1,prod(size(X))); 
    reshape(Y,1,prod(size(Y))); 
    reshape(Z,1,prod(size(Z))); 
    ones(1,prod(size(Z)))]; 
phi=[0 : pi/20 : 50*pi]; 
h=[]; axis([-20 20 -20 20 -20 20]); 

for beta=phi % animation loop ***************** 

    T=[cos(beta) -sin(beta) 0 0;  % rotation matrix 
     sin(beta) cos(beta) 0 0; 
     0   0  1 0; 
     0   0  0 1]; 

    XYZ1=T*XYZ; % coordinates changing 
    X1=reshape(XYZ1(1,:),size(X));Y1=reshape(XYZ1(2,:),size(Y));Z1=reshape(XYZ1(3,:),size(Z)); 
    pause(0.1);if ~isempty(h),delete(h);end 
    h=surf(X1,Y1,Z1,'EdgeColor',[0.5 1. 0.2],'FaceColor',[0.2 0.2 0.8],'FaceAlpha',0.6); 

end  % ****************************************** 

end 

Answer 1

：

http://adasu.info/plates.png

完整的代碼。 r用於X和Y，所以當對這些輸入使用surf時，最終會出現該錯誤。這是有道理的，因爲你的V範圍有負值，當你設置k=1這個表達式：

r=sqrt((4*V.^-k)./(cos(U).^2+k*sin(U).^2)); 

您有效地試圖利用價值的平方根V，其中一些是否定的，因此r是複雜的價值。如果你看看你上傳的實際圖片，你會錯過的V的力量。因此：

r=sqrt((4*V.^2-k)./(cos(U).^2+k*sin(U).^2)); 

當我這樣做，那麼嘗試運行你的代碼，我得到這個：