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

tuoni

: |# [" E0 A! O9 h7 h3 K我需要把三个mesh绘图俺的曲面放在一个坐标系内，并且将三个曲面分别设置为三个不同三个的纯色显示，然后在图中标记出曲面的交点，麻烦大神帮我看看怎么做？贴上源程序。非常感谢！
np=0:100;nd=0:100;ca=5;
& v% {3 {- `; l0 `1 r' y$ [[np,nd]=meshgrid(np,nd);
8 e" S- f( H0 B" K, i0 MPp2=1000+1.1.*nd+0.8.*(np-1);
Pd2=1000+np+0.6.*(nd-1);
- x! N9 n9 X7 I8 _/ X7 t3 v2 EPp1=(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;

Ua1=Pp1.*np-ca.*(np+nd);
% g7 c$ Z) M6 e( M; [: ~! AUa2=Pd1.*nd-ca.*(np+nd);
. I6 K, u$ x) i. PUa3=(Pp2-ca).*np+(Pd2-ca).*nd;
( z7 y- c8 L, B( Z+ m$ Umesh(np,nd,Ua1,'FaceColor',[1,1,0]);
>> hold on;
7 X# U" z% G" H" A: U>> mesh(np,nd,Ua2,'FaceColor',[1,0,0]);
6 w5 z: l" e! d* K% {9 l) k1 y>> hold on;
3 d' J7 b# {  U4 I8 \>>mesh(np,nd,Ua3,'FaceColor',[0,0,1]);

, Y8 y) ^- c' {, E1 k( Q程序修改, mesh绘图, 图中标注交点
7 _2 d. C! l1 ~/ P

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

: X/ y1 R8 C& K, |4 }

Colbie

; O0 _+ S+ f7 x) f- d; Y; u
8 a. ]5 h) L/ _0 f
+ a: f5 o. h; _* d: ~' _0 _7 uh1 = surf(np,nd,Ua1,'FaceColor',[1,1,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
+ J! }& V0 g" n4 @( H    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
! A* n/ ~6 r- g5 d+ l! qh2 = surf(np,nd,Ua2,'FaceColor',[1,0,0],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
' }3 L) O4 R8 U7 R    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8); hold on;
) Z/ u9 o, ?0 D) L" o; lh3 = surf(np,nd,Ua3,'FaceColor',[0,0,1],'FaceAlpha',1,'FaceAlpha',0.7,'MeshStyle','default','EdgeColor','k','EdgeAlpha',0.6,...
6 z4 y: N; t3 L2 i; R- @; H    'AlignVertexCenters','on','LineStyle','-','LineWidth',0.8);
zdiff01 = Ua1 - Ua2;
C01 = contours(np,nd,zdiff01,[0 0]);
9 M5 I- F" [' }# e: M; q( y6 Gxline01 = C01(1,2:end);
yline01 = C01(2,2:end);
* C8 M+ x* U+ F6 D! r! \+ g# rzline01 = interp2(np,nd,Ua1,xline01,yline01);
1 [- v+ a5 ^2 L2 tzdiff02 = Ua1 - Ua3;
0 U3 L: l' p( d5 E- IC02 = contours(np,nd,zdiff02,[0 0]);
, M. I( C. \+ X- qxline02 = C02(1,2:end);
( y, d6 g, A: |, Z' C, G0 Jyline02 = C02(2,2:end);
zline02 = interp2(np,nd,Ua1,xline02,yline02);
zdiff03 = Ua2 - Ua3;
7 l4 q! i( |1 g9 fC03 = contours(np,nd,zdiff03,[0 0]);
xline03 = C03(1,2:end);
yline03 = C03(2,2:end);
- F! [8 V6 z$ B7 uzline03 = interp2(np,nd,Ua2,xline03,yline03);
, D5 ^4 D2 z6 Z$ r" |! rl1 = plot3(xline01,yline01,zline01,'g-','LineWidth',3);
l2 = plot3(xline02,yline02,zline02,'c-','LineWidth',3);
1 H% Z0 X# p" L; Ll3 = plot3(xline03,yline03,zline03,'k-','LineWidth',3);
' ]6 _: y+ U( o8 x, G' Vlegend([h1 h2 h3 l1 l2 l3],{'Surface01','Surface02','Surface03','Intersection01','Intersection02','Intersection03'},'Location','SoutheastOutside');
2 }" S: w, C8 b( v" V& y1 Q- i9 ]

tuoni

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

9 g' D) _+ F/ x非常感谢大神。
0 q2 d1 Y# t( j* ^% i