MATLAB绘制同时绘制两曲面，并显示交线，及交线在XOY平面投影

找个画图程序找半天，索性自己写了个 

%此处曲面以平面为例，和曲面相交，绘制交线，并显示XOY投影

% 定义x,y区间
x=-10:0.01:10;
y=x;
% 制作网格
[X,Y]=meshgrid(x,y);
% 平面Z1 = 4
Z1=4*ones(size(X));
% 平面Z2 = x^2/4 - x^2/9
Z2=X.^2/4-Y.^2/9;

%同时绘制Z1与Z2
mesh(X,Y,Z1);
hold on
mesh(X,Y,Z2);

% 寻找相交点，此处设为.65是为了效果更明显，严格为0
r0=(abs(Z1-Z2)<=.65);
% 记录交线坐标
zz=r0.*Z1;yy=r0.*Y;xx=r0.*X;

% 绘制两曲面交线
plot3(xx(r0~=0),yy(r0~=0),zz(r0~=0),'k*')

hold on 
% 绘制在XOY平面的投影
plot3(xx(r0~=0),yy(r0~=0),zeros(size(xx(r0~=0))),'k*');

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