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

tuoni

我需要把三个mesh绘图俺的曲面放在一个坐标系内，并且将三个曲面分别设置为三个不同三个的纯色显示，然后在图中标记出曲面的交点，麻烦大神帮我看看怎么做？贴上源程序。非常感谢！
np=0:100;nd=0:100;ca=5;
[np,nd]=meshgrid(np,nd);
Pp2=1000+1.1.*nd+0.8.*(np-1);
! A6 i) D: c% W: u2 A# b, DPd2=1000+np+0.6.*(nd-1);
9 T) J. O2 e, E7 r; V! ]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;
6 N% ^( J$ j+ u
Ua1=Pp1.*np-ca.*(np+nd);
Ua2=Pd1.*nd-ca.*(np+nd);
% [( F+ O3 u* W8 d. N" UUa3=(Pp2-ca).*np+(Pd2-ca).*nd;
' ?0 o8 j& w. Z% xmesh(np,nd,Ua1,'FaceColor',[1,1,0]);
>> hold on;
>> mesh(np,nd,Ua2,'FaceColor',[1,0,0]);
>> hold on;
: [2 V& @5 W( ^  |4 B# x! O>>mesh(np,nd,Ua3,'FaceColor',[0,0,1]);
, c- o0 p# d+ s  p4 c
程序修改, mesh绘图, 图中标注交点
/ |% h/ j) l$ `2 c% k

我这样出来的图根本看不出哪个是哪个，

. z2 ~& i8 G5 a  ]0 Y. j7 A

Colbie

) h$ F2 s) B  A4 }
3 e1 \6 H8 f* D. n
, W( Z' D3 o! _& _/ p& F; ], J' Nh1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
9 Y2 O# l4 h! u6 f; h    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
9 D; M. e; N- m" j: _6 j% Nh2 = surf(np,nd,Ua2,'FaceColor',[1,0,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
h3 = surf(np,nd,Ua3,'FaceColor',[0,0,1],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
( I# t% W6 c7 t5 V) Y6 e    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8);
zdiff01 = Ua1 - Ua2;
/ T8 J' V4 }+ D* T- wC01 = contours(np,nd,zdiff01,[0 0]);
) E/ I+ k3 b1 fxline01 = C01(1,2:end);
yline01 = C01(2,2:end);
zline01 = interp2(np,nd,Ua1,xline01,yline01);
9 s2 _8 E7 O/ v! j2 x; Yzdiff02 = Ua1 - Ua3;
* c" n5 [% p0 l  S% qC02 = contours(np,nd,zdiff02,[0 0]);
' @* e* C  ^! N( C$ D/ e+ L  k* cxline02 = C02(1,2:end);
yline02 = C02(2,2:end);
zline02 = interp2(np,nd,Ua1,xline02,yline02);
zdiff03 = Ua2 - Ua3;
! V' |2 r- p' j4 AC03 = contours(np,nd,zdiff03,[0 0]);
6 v9 s' P  K( q+ q  E+ V* j3 w2 uxline03 = C03(1,2:end);
yline03 = C03(2,2:end);
0 i+ I+ y9 O3 D: }% Z) g% pzline03 = interp2(np,nd,Ua2,xline03,yline03);
0 D" p! z  Z) E" N: Jl1 = plot3(xline01,yline01,zline01,'g-','LineWidth',3);
& E: U3 ^- c* N. c5 ]l2 = plot3(xline02,yline02,zline02,'c-','LineWidth',3);
8 T: q, ?8 c, n0 Gl3 = plot3(xline03,yline03,zline03,'k-','LineWidth',3);
legend([h1 h2 h3 l1 l2 l3],{'Surface01','Surface02','Surface03','Intersection01','Intersection02','Intersection03'},'Location','SoutheastOutside');
( ]' `: q: G7 P% W9 k3 ?* ]: H
$ ]1 B& Y' l+ s) Y

tuoni

Colbie 发表于 2020-6-19 17:43
' i- X! D' H' _  }6 o: h* Dh1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','E ...

非常感谢大神。
* x2 }( s6 V5 G$ B4 Z5 J6 z