matlab画曲面

xx=linspace(-1,1,50);

yy=linspace(-2,2,100);

[x,y]=meshgrid(xx,yy);

z=x.^2+2*y.2;

surf(x,y,z)

结果

x=-20:0.5:20;
y=-20:0.5:20;
[X,Y]=meshgrid(x,y); %x，y平面bai网格化
Z=-X.^2+X-X.*Y;
mesh(X,Y,Z) %绘制三维网格曲面（du图1）
figure(2)
surf(X,Y,Z) % 绘制三维彩色曲面（图2）

x=-20:0.5:20;
这是你要求的区间，从-20到20，0.5是步长，即-20，-19.5，-19，-18.5，-18，…，19.5，20

1、先定义bai参数t,再给出方程,最后画图即可,比如：du
t=0:1:100;
x=f(t);
y=g(t);
plot(x,y);
2、原理与1中相同,不同的是先把zhi极坐标化为直角dao坐标系中的坐标即可,比如：
a=0:1:100;
x=cos(a);
y=sin(a);
p=x+y;
plot(a,p);
此时p=cos(θ)+sin(θ),其他情况下依次类推即可.

参数方程画曲面：

h=k/40;
v = (ucos(0.1847)40+sin(0.1847)(ucos(k)+65sin(k)-200))/sin(k);
x0=u;
y0=vsin(0.1847)-65;
z0=v*cos(0.1847);
x1=-cos(k)*cos(h)*x0+sin(k)*cos(h)*y0-sin(h)z0+200cos(h);
y1=cos(k)*sin(h)*x0-sin(k)*sin(h)*y0-cos(h)z0-200sin(h);
z1=-sin(k)*x0-cos(k)*y0;

其中143.544<=u<=157.921
0.0402pi<=k<=0.2176pi

代码如下
[u,k]=meshgrid(143.544:0.01:157.921, (0.0402:0.01:0.2176)pi);
h=k/40;
v=(ucos(0.1847)40+sin(0.1847)(u.cos(k)+65sin(k)-200))./sin(k);
x0=u;
y0=vsin(0.1847)-65;
z0=vcos(0.1847);
x1=-cos(k).*cos(h).*x0+sin(k).*cos(h).*y0-sin(h).z0+200cos(h);
y1= cos(k).*sin(h).*x0-sin(k).*sin(h).*y0-cos(h).z0-200sin(h);
z1=-sin(k).*x0-cos(k).*y0;
surf(x1,y1,z1)
shading interp

%[u,k]=meshgrid(143.544:0.01:157.921, (0.0402:0.01:0.2176)*pi);
%t1 = -pi : 0.05 : pi
%t2 = -pi : 0.05 : pi

e = 0
r1 = 3
r = 10
q = pi/6

t1 =meshgrid( 0: 0.01 : pi/2)
t2 =meshgrid( 0 : 0.01 /2: pi/4)

xm1 = e + r * cos(t1) - r1 * cos(t1) .* (cos(q) - cos(t2))
ym1 = r1 * (sin(q) + sin(t2))
zm1 = r * sin(t1) - r1 * sin(t1) .* (cos(q) - cos(t2))
surf(xm1, ym1, zm1)

e = 0
r1 = 3
r = 10
q = meshgrid( 0: 0.04 : 2 * pi)

t1 =meshgrid( 0: 0.01 : pi/2)
t2 =meshgrid( 0 : 0.01 /2: pi/4)

xm1 = e + r * cos(t1) - r1 * cos(t1) .* (cos(q) - cos(t2))
ym1 = r1 * (sin(q) + sin(t2))
zm1 = r * sin(t1) - r1 * sin(t1) .* (cos(q) - cos(t2))
surf(xm1, ym1, zm1)

e = 0
r1 = 3
r = 10
q = meshgrid( 0: 2pi/1000 : 2pi)

t1 =meshgrid( 0: pi/2/1000 : pi/2)
t2 =meshgrid( 0 : pi/4/1000: pi/4)

xm1 = e + r .* cos(t1) - r1 * cos(t1) .* (cos(q) - cos(t2))
ym1 = r1 * (sin(q) + sin(t2))
zm1 = r .* sin(t1) - r1 * sin(t1) .* (cos(q) - cos(t2))
surf(xm1, ym1, zm1)

e = 0
r1 = 3
r2 = 3
r = 10
q = pi/4
s = 1
t1 =meshgrid( 0: 0.01 : pi/2)
t2 =meshgrid( 0 : 0.01 /2: pi/4)
t3 = meshgrid( 0 : 0.01 /2: pi/4)

x = e +r * cos(t1) - r1 * cos(t1) * (cos(q) - cos(s)) + r2 * cos(t1) * (cos(s) + cos(t3))
y = r1 * (sin(q) + sin(s)) + r2 * (sin(s) + sin(t3))
z = r * sin(t1) - r1 * sin(t1) * (cos(q) - cos(s)) + r2 * sin(t1) * (cos(s) + cos(t3))

surf(x,y,z)

,