[x,y,z]=sphere(30);
X=a*x;   %%(a,b,c为椭球的半轴）
Y=b*y;
Z=c*z;
surf(X,Y,Z)

不好意思，看漏了极坐标。但好像surf函数只是支持笛卡尔坐标，那么可以先用pol2cart把极坐标转化成笛卡尔坐标就行了

谢谢阿，那怎么转化呢？

帮助里有：

    [X,Y,Z] = POL2CART(TH,R,Z) transforms corresponding elements of
    data stored in cylindrical coordinates (angle TH, radius R, height
    Z) to Cartesian coordinates X,Y,Z. The arrays TH, R, and Z must be
    the same size (or any of them can be scalar).  TH must be in radians.