如果求一般性的两个三维曲面的交线还是有一些难度的(尤其对于两个曲面都是隐函数的情况),但本题是两个单位柱面的交线,情况比较特殊,用参数方程比较容易。
不妨以x-y平面内的相角t为参变量,则满足第一个柱面方程的x和y分别为cos(t)和sin(t),而由第二个柱面方程不难解出z=sqrt(1-y^2)或z=-sqrt(1-y^2)。
参考代码:
[x,y,z]=cylinder(1,200);
z=2.2*z-1.1;
cla
a=0.3;
mesh(x,y,z,'edgecolor','none','facecolor','r','facealpha',a)
hold on
mesh(z,x,y,'edgecolor','none','facecolor','g','facealpha',a)
axis equal
t=linspace(0,2*pi,200);
x=cos(t);
y=sin(t);
z1=sqrt(1-y.^2);
z2=-sqrt(1-y.^2);
plot3(x,y,z1,x,y,z2,'linewidth',3)
view(60,10)
xlabel x; ylabel y; zlabel z