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

tuoni

我需要把三个mesh绘图俺的曲面放在一个坐标系内，并且将三个曲面分别设置为三个不同三个的纯色显示，然后在图中标记出曲面的交点，麻烦大神帮我看看怎么做？贴上源程序。非常感谢！
# {, o" I' h6 c; v2 \np=0:100;nd=0:100;ca=5;
[np,nd]=meshgrid(np,nd);
Pp2=1000+1.1.*nd+0.8.*(np-1);
5 }$ _7 B& y/ T7 F" l' APd2=1000+np+0.6.*(nd-1);
Pp1=(2078.88+120+1106.6/6-(31/30).*np+0.62.*nd)./2;
# ~. t) x* h, P. UPd1=(2125.4+200+2991/55-0.775.*nd+(31/55).*np)./2;

- \/ U+ c4 s# [Ua1=Pp1.*np-ca.*(np+nd);
$ m. w) [' q1 e. a0 }/ t7 R3 F% VUa2=Pd1.*nd-ca.*(np+nd);
7 U% e0 a/ y5 I* l$ \4 n5 EUa3=(Pp2-ca).*np+(Pd2-ca).*nd;
6 l/ O; }9 q0 |) A5 J6 Emesh(np,nd,Ua1,'FaceColor',[1,1,0]);
" H1 ?4 r  t' h) L8 Z>> hold on;
% J% c4 @+ p6 S0 f  m>> mesh(np,nd,Ua2,'FaceColor',[1,0,0]);
>> hold on;
, e3 B6 M8 ^9 ?; _9 C9 H>>mesh(np,nd,Ua3,'FaceColor',[0,0,1]);
6 j. C+ B0 z/ g9 c
+ f" H* `6 R1 x0 ]程序修改, mesh绘图, 图中标注交点
4 q$ W6 Z" X& }! M; v: L
" [; F5 H  W& C6 u
) b& n# R. ^' P) e, D, O, U我这样出来的图根本看不出哪个是哪个，
) R2 q3 k0 S4 R  f5 }( e- e

Colbie

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;
) o5 l  ^6 ~+ k) S/ |h2 = surf(np,nd,Ua2,'FaceColor',[1,0,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
; H+ B4 [6 h6 x2 X    '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,...
    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8);
+ Z" C, M% m6 c& Y% l, xzdiff01 = Ua1 - Ua2;
# ]. ^8 l: V  p- f4 d) |" jC01 = contours(np,nd,zdiff01,[0 0]);
5 W/ S: B) P$ Z) O: j9 }xline01 = C01(1,2:end);
. w; ?2 Q; y- {1 z: d1 Xyline01 = C01(2,2:end);
zline01 = interp2(np,nd,Ua1,xline01,yline01);
zdiff02 = Ua1 - Ua3;
4 F/ }6 O/ j& d3 v) W+ yC02 = contours(np,nd,zdiff02,[0 0]);
xline02 = C02(1,2:end);
yline02 = C02(2,2:end);
1 Y2 {" T* H0 ~1 Z0 P& b9 Qzline02 = interp2(np,nd,Ua1,xline02,yline02);
! o% u4 F9 W! D, k+ [zdiff03 = Ua2 - Ua3;
C03 = contours(np,nd,zdiff03,[0 0]);
xline03 = C03(1,2:end);
yline03 = C03(2,2:end);
zline03 = interp2(np,nd,Ua2,xline03,yline03);
+ ?1 I6 \" C* x5 E5 Ol1 = plot3(xline01,yline01,zline01,'g-','LineWidth',3);
l2 = plot3(xline02,yline02,zline02,'c-','LineWidth',3);
4 c% D2 E, s3 i" t" _8 `- d$ `l3 = plot3(xline03,yline03,zline03,'k-','LineWidth',3);
, M2 V# J0 R/ H, P5 a8 ^& o( U$ blegend([h1 h2 h3 l1 l2 l3],{'Surface01','Surface02','Surface03','Intersection01','Intersection02','Intersection03'},'Location','SoutheastOutside');
4 {# ^8 x  ^& a* f
- j2 H7 v/ m! E

tuoni

Colbie 发表于 2020-6-19 17:43
6 p7 i4 l+ H$ |3 [h1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','E ...

8 |- ?5 o6 D$ x+ U. F4 r非常感谢大神。