用matlab软件作锥面x^2+y^2=z^2和柱面(x-1)^2+y^2=1相交的图形，求大神帮忙写一下程序~~~

你看看程序，不懂得可以问我。

%程序

xgrid = -2:0.1:2

ygrid = -2:0.1:2

[x,y] = meshgrid(xgrid,ygrid)

z_p = sqrt(x.^2+y.^2)

z_n = -sqrt(x.^2+y.^2)

mesh(x,y,z_p)%锥面上半部分

hold on

mesh(x,y,z_n)%锥面下半部分

hold on

t=0:pi/20:2*pi

x1=sin(t)+1

y1=cos(t)

z1=linspace(-3,3,length(t))

X=meshgrid(x1)Y=meshgrid(y1)Z=[meshgrid(z1)]'

mesh(X,Y,Z)%圆柱面

xlabel('x')

ylabel('y')

zlabel('z')

axis equal

程序如下：

t=0:pi/20:2*pi

x=sin(t)*2

y=cos(t)*2

z=linspace(-5,5,length(t))

X=meshgrid(x)Y=meshgrid(y)Z=[meshgrid(z)]'

mesh(X,Y,Z)%第一个圆柱面

xlabel('x')

ylabel('y')

zlabel('z')

hold on

x1=sin(t)*2

z1=cos(t)*2

y1=linspace(-5,5,length(t))

X1=meshgrid(x1)Y1=[meshgrid(y1)]'Z1=meshgrid(z1)

mesh(X1,Y1,Z1)%第二个圆柱面

plot3(sin(t)*2,cos(t)*2,cos(t)*2,'k-*','LineWidth',2)%交线1

hold on

plot3(-sin(t)*2,-cos(t)*2,cos(t)*2,'k-*','LineWidth',2)%交线2

title('柱面x^2+y^2=4 和柱面x^2+z^2=4 相交的图形')

结果：

1、首先我们打开Matlab软件，默认d出窗口 Command Window，将已请确定的(x,y,z)值依次输入该窗口中。

2、清空工作空间与命令窗口后，输入x=[-pi:pi/100:2*pi]'y=sin(x)。

3、输入plot(x,y)绘制图形。

4、接着输入以下命令:x=0:0.1:5y=1+exp(x)。

5、输入plot(x,y,'-+b')绘制图形。

注意事项：

MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

---