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

tuoni

我需要把三个mesh绘图俺的曲面放在一个坐标系内，并且将三个曲面分别设置为三个不同三个的纯色显示，然后在图中标记出曲面的交点，麻烦大神帮我看看怎么做？贴上源程序。非常感谢！
: D+ T, }/ O, v' W& snp=0:100;nd=0:100;ca=5;
3 W0 _& i5 }! Z0 A) t9 o- {( k$ E4 l! g[np,nd]=meshgrid(np,nd);
Pp2=1000+1.1.*nd+0.8.*(np-1);
Pd2=1000+np+0.6.*(nd-1);
3 e! H" L7 o  r) n  n9 H( CPp1=(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 ?' _, }  n) a; R% M4 x5 p( p5 ^
Ua1=Pp1.*np-ca.*(np+nd);
; k* ^3 y2 m. ~% V! bUa2=Pd1.*nd-ca.*(np+nd);
: O" N4 u6 |. R& Y( N# YUa3=(Pp2-ca).*np+(Pd2-ca).*nd;
* P2 Y, I. F; V- jmesh(np,nd,Ua1,'FaceColor',[1,1,0]);
>> hold on;
+ q' S3 D! p: J  |: }2 m% Y6 s>> mesh(np,nd,Ua2,'FaceColor',[1,0,0]);
>> hold on;
0 |3 I3 f" x2 l) k1 [( D>>mesh(np,nd,Ua3,'FaceColor',[0,0,1]);
$ J- P" _. Q) j2 ]' i
程序修改, mesh绘图, 图中标注交点

' D' w1 \8 e8 F) X5 h& H5 S
我这样出来的图根本看不出哪个是哪个，
8 X7 b6 [& ]) g/ r1 ]% ~9 O
2 z9 M( E+ z/ M$ d

Colbie

% b4 |3 t) s7 r4 @" g, F2 f
# A# A. u0 c5 E$ D. X/ H$ dh1 = 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;
; ]4 i0 j7 y5 k1 s( q7 Mh2 = surf(np,nd,Ua2,'FaceColor',[1,0,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
- }. ?# q! Z5 x- Y9 u$ B" E    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
  \! r" ~) L( ^! {: p# Qh3 = 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);
6 p; V/ s8 V" s2 B; f1 G; ?zdiff01 = Ua1 - Ua2;
& x6 Y$ d3 u/ d1 C2 T; }7 q. y# VC01 = contours(np,nd,zdiff01,[0 0]);
6 o' b$ B7 S# H. v5 Rxline01 = C01(1,2:end);
yline01 = C01(2,2:end);
* e" S- n; ~2 f- R9 i1 i4 [/ ~" Ozline01 = interp2(np,nd,Ua1,xline01,yline01);
zdiff02 = Ua1 - Ua3;
9 z% ]6 P+ m2 d: @8 z; NC02 = contours(np,nd,zdiff02,[0 0]);
xline02 = C02(1,2:end);
yline02 = C02(2,2:end);
zline02 = interp2(np,nd,Ua1,xline02,yline02);
2 |" F3 O5 J/ z! e8 y7 J& U' rzdiff03 = Ua2 - Ua3;
. n  o" }! J- F# nC03 = contours(np,nd,zdiff03,[0 0]);
0 C3 h+ J6 z+ b& Qxline03 = C03(1,2:end);
yline03 = C03(2,2:end);
. x6 t  N8 r9 qzline03 = interp2(np,nd,Ua2,xline03,yline03);
: j. W/ q8 @5 \/ Ql1 = plot3(xline01,yline01,zline01,'g-','LineWidth',3);
l2 = plot3(xline02,yline02,zline02,'c-','LineWidth',3);
l3 = plot3(xline03,yline03,zline03,'k-','LineWidth',3);
. p- M; p2 J7 y. B$ ]4 ilegend([h1 h2 h3 l1 l2 l3],{'Surface01','Surface02','Surface03','Intersection01','Intersection02','Intersection03'},'Location','SoutheastOutside');
; s- r. A% ~8 k7 r0 ^" n

tuoni

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

非常感谢大神。