clc; clear;
%%Ô²Öù¿ÇÌåģ̬ÕñÐÍ
m=3; n=2;                            %ģ̬²ÎÊý£ºÖáÏò3¸ö°ë²¨£¬ÖÜÏòÁ½¸öÕû²¨
L=4;                                 %¿ÇÌ峤¶ÈΪ4
r=2;                                 %¿ÇÌå°ë¾¶Îª2
A=1; B=1; C=1;                       %ÖáÏò£¬ÖÜÏò£¬¾¶Ïòģ̬Õñ·ù£¬¿É¾ùȡΪ1
z = (0:0.1:length(L))'; s = 0:pi/50:2*pi;
%Õñ¶¯³õʼµã
u=£¿;                               %ÖáÏò³õʼλÖà                
v=£¿;                               %ÖÜÏò³õʼλÖÃ
w=£¿;                               %¾¶Ïò³õʼλÖÃ
%Õñ¶¯Î»ÒÆ
X = A*cos(m*pi*z/L)*(cos(n*s));    %XÊÇÖáÏòλÒÆ
Y = B*sin(m*pi*z/L)*(sin(n*s));    %YÊÇÔ²ÖÜ·½ÏòµÄλÒÆ
R = C*sin(m*pi*z/L)*(cos(n*s));    %RÊǾ¶ÏòλÒÆ

[XX,YY,ZZ] = pol2cart(Y+v,R+w,X+u);
h=surf(XX,YY,ZZ);
set(h,'CData',YY);
axis equal ·µ»ØСľ³æ²é¿´¸ü¶à