function c=make3d(fun_str,is_cplxgrid,x_min,x_max,precision,in_a,in_b,in_c,diff_n,s_color,is_x_y_z,math_Mode,in_R,is_holdon)

%  函数的使用说明

% 直接显示出函数的图像，及区域极大值(导数为0)、区域极小值(导数为0)

% 输入的参数：

%                                fun_str:函数字符串

%                                                例如：'(2.*x-a)/power((x-1),2)' ，其中，指数都要写成‘power(a,b)’的形式， ‘*’应该写成‘.*’，‘/’应该写成‘./’

%                 is_cplxgrid: 是否构建一个极坐标的复数数据网格。

%                 一般在球状、椭球状模型时才需要构建复数数据网格(圆面)，否则一般构建矩形数据网格。

%                              0: 不构建；

%                              >0:产生(is_cplxgrid+1)*(2*is_cplxgrid+1)的极坐标下的复数数据网格，最大半径为1的圆面。

%                                x_max,x_min：需要在x轴的显示的区域极大值(导数为0)或区域极小值(导数为0)，也可以视为定义域

%                                precision ： 精度

%                                in_a,in_b,in_c：在函数需要用到的参数

%                                                例如：'(2.*x-a)/power((x-1),2)' ，其中，a=in_a

%                                diff_n：导数的阶，求n阶导数。0：不求导；>0: 求导；<0:求积分

%                                s_color ：线条颜色，例如‘r'、’b‘

% 线形：－实线 －. 点划线 －－长虚线 ：短虚线

% 符号 颜色   符号      线形

% b     蓝    .         点

% c     青    。        圈

% g     绿    ×      ×标记

% k     黑    －       实线

% m     紫红  *         星号

% r     红    ：       点线

% w     白    －.      点划线

% y     黄    －－     虚线

%                                is_x_y_z ： 判断以哪个轴进行拉伸，缺省是z轴：x:x轴；y：y轴；z：z轴

%                                is_holdon ：是否需要与以前的图像重叠，0:否；1：是

%                                c：返回值，暂时未定义

%               math_Mode:  数学模型，用字符串表示：'球面'、'椭圆锥面'、'椭球面'、'椭圆抛物面'、'双曲抛物面'、'单叶双曲面'、'双叶双曲面'、其他。

%               in_R： 半径

%   +eps 可以避免分母为0

% make3d('x1=sqrt((z+y.^2./b.^2).*a.^2);mesh(real(x1),y,z);hold on;mesh(-real(x1),y,z);title(''双曲抛物面'');',-10,10,0.06,3,1,1,-1,'r','x',-1,1,0);

% make3d('x1=sqrt((z.^2./a.^2+y.^2./b.^2+1).*c.^2);mesh(real(x1),y,z);hold on;mesh(-real(x1),y,z);title(''双叶双曲面'');',-10,10,0.06,3,1,1,-1,'r','x',-1,1,0);

% make3d('z1=sqrt((x.^2./a.^2+y.^2./b.^2-1).*c.^2);mesh(x,y,real(z1));hold on;mesh(x,y,-real(z1));title(''单叶双曲面'');',-10,10,0.06,1,1,1,-1,'r','z',-1,1,0);

if is_holdon

hold on;

x=0;y=0;

else

hold off;

clf;

end

format long;

a=in_a;

b=in_b;

c=in_c;

fn_min=0;fn_max=0;y=0;f2_min=0;tt='';

switch is_x_y_z

case 'x'

if is_cplxgrid==0

fn_min=x_min;fn_max=x_max;        % 设置函数的定义域[fn_min,fn_max]

y=(fn_min):precision:(fn_max*1);z=y;

[y,z]=meshgrid(y,z);        % 产生矩形数据网格

else

x=cplxgrid(is_cplxgrid);y=real(x);z=imag(x); % 产生(is_cplxgrid+1)*(2*is_cplxgrid+1)的极坐标下的复数数据网格，最大半径为1的圆面。

end

eval(fun_str);

case 'y'

if is_cplxgrid==0

fn_min=x_min;fn_max=x_max;        % 设置函数的定义域[fn_min,fn_max]

x=(fn_min):precision:(fn_max*1);z=x;

[x,z]=meshgrid(x,z);        % 产生矩形数据网格

else

y=cplxgrid(is_cplxgrid);x=real(y);z=imag(y); % 产生(is_cplxgrid+1)*(2*is_cplxgrid+1)的极坐标下的复数数据网格，最大半径为1的圆面。

end

eval(fun_str);

otherwise

if is_cplxgrid==0

fn_min=x_min;fn_max=x_max;        % 设置函数的定义域[fn_min,fn_max]

x=(fn_min):precision:(fn_max*1);y=x;

[x,y]=meshgrid(x,y);        % 产生矩形数据网格

else

z=cplxgrid(is_cplxgrid);x=real(z);y=imag(z); % 产生(is_cplxgrid+1)*(2*is_cplxgrid+1)的极坐标下的复数数据网格，最大半径为1的圆面。

x=x.*in_R;y=y.*in_R;

x=x.*a;y=y.*b;

end

switch math_Mode

case '球面'

%  a相当于x0, b相当于y0,c相当于z0,in_R相当于半径R

z1=sqrt(in_R.^2-(x-a).^2-(y-b).^2)+c;

surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('球面');

case '椭圆锥面'

z1=sqrt((x.^2./a.^2+y.^2./b.^2));

surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('椭圆锥面');

case '椭球面'

z1=sqrt((1-x.^2./a.^2-y.^2./b.^2).*c.^2);

surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('椭球面');

case '椭圆抛物面'

z1=x.^2./a.^2+y.^2./b.^2;

surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('椭圆抛物面');

case '双曲抛物面'

z1=x.^2./a.^2-y.^2./b.^2;

%                 surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('双曲抛物面');

surf(x,y,real(z1),'LineStyle','none');title('双曲抛物面');

case '单叶双曲面'

z1=sqrt((x.^2./a.^2+y.^2./b.^2-1).*c.^2);

surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('单叶双曲面');

%                 surf(x,y,real(z1),'LineStyle','none');

case '双叶双曲面'

z1=sqrt((x.^2./a.^2-y.^2./b.^2-1).*c.^2);

surf(x,y,real(z1),'LineStyle','none');hold on;surf(x,y,real(-z1),'LineStyle','none');title('双叶双曲面');

otherwise

title(math_Mode);

eval(fun_str);

end

end

grid on;

xlabel('x轴');ylabel('y轴');zlabel('z轴');

% axis([-10 10 -10 10 -10 10]);

% axis square;

axis equal

调用函数：

>> make3d('',100,-5,5,0.05,2,4,3,-1,'r','z','双叶双曲面',80,0);

>> make3d('',0,-5,5,0.05,2,4,3,-1,'r','z','双叶双曲面',3,0);

x^2/25+y^2/36  肯定是单叶双曲面啊！