麻烦大神帮我看看怎么做？

tuoni

6 y  \8 W$ O' }  k6 N我需要把三个mesh绘图俺的曲面放在一个坐标系内，并且将三个曲面分别设置为三个不同三个的纯色显示，然后在图中标记出曲面的交点，麻烦大神帮我看看怎么做？贴上源程序。非常感谢！
0 j$ b( F! c* |( h1 |" S6 Tnp=0:100;nd=0:100;ca=5;
[np,nd]=meshgrid(np,nd);
" B" m- m5 T) uPp2=1000+1.1.*nd+0.8.*(np-1);
$ e# O& l  T2 w$ X% mPd2=1000+np+0.6.*(nd-1);
Pp1=(2078.88+120+1106.6/6-(31/30).*np+0.62.*nd)./2;
Pd1=(2125.4+200+2991/55-0.775.*nd+(31/55).*np)./2;
1 Z" N8 d! [: C" |9 U# J7 l
8 g. @9 N/ Q4 W5 CUa1=Pp1.*np-ca.*(np+nd);
4 y0 C1 d3 d$ ?4 ]$ _3 e0 pUa2=Pd1.*nd-ca.*(np+nd);
0 [4 L9 @) f4 o, cUa3=(Pp2-ca).*np+(Pd2-ca).*nd;
6 [+ y5 a5 w( }  emesh(np,nd,Ua1,'FaceColor',[1,1,0]);
# h; K: n5 G) j* @0 h+ K>> hold on;
>> mesh(np,nd,Ua2,'FaceColor',[1,0,0]);
2 v2 l  H5 N. m; ?>> hold on;
7 ^5 L1 L8 ]* |# a9 N6 C>>mesh(np,nd,Ua3,'FaceColor',[0,0,1]);

8 p0 _, D* U4 B' g+ K( N程序修改, mesh绘图, 图中标注交点

8 z& b& v$ `: R7 ~9 Q+ H' Y$ ~5 l. o
我这样出来的图根本看不出哪个是哪个，

Colbie

$ Z' f& P* A; \  h3 C3 ?

h1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
h2 = surf(np,nd,Ua2,'FaceColor',[1,0,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
1 w& w8 z. Z! g" Z: Q6 o, j4 @    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
; o/ H& [' Y2 c& q, r& jh3 = surf(np,nd,Ua3,'FaceColor',[0,0,1],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8);
zdiff01 = Ua1 - Ua2;
C01 = contours(np,nd,zdiff01,[0 0]);
! V( W& @) E# y5 Txline01 = C01(1,2:end);
# X% |& n5 _% f0 R3 Wyline01 = C01(2,2:end);
5 I, C' y, }/ e: r; N, B2 \zline01 = interp2(np,nd,Ua1,xline01,yline01);
zdiff02 = Ua1 - Ua3;
" d# a" I2 o" ^0 {0 IC02 = contours(np,nd,zdiff02,[0 0]);
: g0 y; Q. P# [4 k! C4 Bxline02 = C02(1,2:end);
yline02 = C02(2,2:end);
/ f* Y9 K9 Z0 Ozline02 = interp2(np,nd,Ua1,xline02,yline02);
6 o% L5 D3 m+ Y0 |; Kzdiff03 = Ua2 - Ua3;
C03 = contours(np,nd,zdiff03,[0 0]);
xline03 = C03(1,2:end);
yline03 = C03(2,2:end);
: M2 }. E9 g$ d) U8 ozline03 = interp2(np,nd,Ua2,xline03,yline03);
+ l5 F* T8 [8 b* u$ Dl1 = plot3(xline01,yline01,zline01,'g-','LineWidth',3);
2 c0 R# M$ X+ xl2 = plot3(xline02,yline02,zline02,'c-','LineWidth',3);
l3 = plot3(xline03,yline03,zline03,'k-','LineWidth',3);
9 ^9 F- }* N5 k, a5 Y! Z2 f4 d. I& a8 Tlegend([h1 h2 h3 l1 l2 l3],{'Surface01','Surface02','Surface03','Intersection01','Intersection02','Intersection03'},'Location','SoutheastOutside');
) f; W- }. |& U/ c# r& N. O

tuoni

Colbie 发表于 2020-6-19 17:43
h1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','E ...

非常感谢大神。
" M1 T1 g8 F1 C9 _