想将两个mesh的三维图放在一个坐标系内，用了hold on函数但是还是只有第一个图出来

zophim

想将两个mesh的三维图放在一个坐标系内，用了hold on函数但是还是只有第一个图出来，麻烦大神帮忙指点一下，非常感谢！

; u) i1 @  i& ^1 N>> clear
ap=0:0.01:1.5;ad=0:0.01:1.5;ud=0;up=0;bp=0.60;bd=0.060;v=50;p=30;c=10;ca=5;m=50;n=50;
  G2 B' S0 O$ X  i! Q7 y+ Q[ap,ad]=meshgrid(ap,ad);
0 m3 a" }" ]+ o. oPp=-((ad.*ap - bd.*bp).*(ca.*(ad./(ad.*ap - bd.*bp) - bd./(ad.*ap - bd.*bp)) - (ap.*bd - bd.*bp + ap.*ud - bd.*up - bd.*c.*m + bd.*m.*p + ap.*n.*p - ap.*n.*v)./(ad.*ap - bd.*bp)))./(2.*bd);
& Z0 u$ B+ E6 u: w" \np1=-(Pp.*bd - ap.*bd + bd.*bp - ap.*ud + bd.*up + bd.*c.*m - bd.*m.*p - ap.*n.*p + ap.*n.*v)./(ad.*ap - bd.*bp);
3 W: n, d  e0 M" g: c+ B" d! Vnd1=(Pp.*ad + ad.*bp - bd.*bp + ad.*up - bp.*ud + ad.*c*m - ad.*m.*p - bp.*n.*p + bp.*n.*v)./(ad.*ap - bd.*bp);
Uap=Pp.*np1-ca.*(nd1+np1);
9 Z8 h* ?* V7 W. y" J) R4 d% @Pd =-((ad.*ap - bd.*bp).*(ca.*(ap./(ad.*ap - bd.*bp) - bp./(ad.*ap - bd.*bp)) - (ad.*bp - bd.*bp + ad.*up - bp.*ud + ad.*c.*m - ad.*m.*p - bp.*n.*p + bp.*n.*v)./(ad.*ap - bd.*bp)))./(2.*bp);
$ z# b+ B2 e& f( _np2=(Pd.*ap + ap.*bd - bd.*bp + ap.*ud - bd.*up - bd.*c.*m + bd.*m.*p + ap.*n.*p - ap.*n.*v)./(ad.*ap - bd.*bp);
nd2=-(Pd.*bp - ad.*bp + bd.*bp - ad.*up + bp.*ud - ad.*c.*m + ad.*m.*p + bp.*n.*p - bp.*n.*v)./(ad.*ap - bd.*bp);
; P+ J4 v: e) s) g0 Q5 `- wUad=Pd.*nd2-ca.*(nd2+np2);
) e5 l# A: O; V# x) p! Z( F>> mesh(ap,ad,Uap)
. l0 \) q: T- e! b) x" u0 J0 M>> hold on
/ o8 c' q  g7 z$ m+ B. N  H>> mesh(ap,ad,Uad)
>> hold off
6 A8 K7 Q  l( S5 E>>
) ~4 A" l3 a& H& i) x0 F出来的图片如下：

1 H( @: w; b" Z& @, P1 |4 P- d
6 @6 H% X3 L' ^; I3 k+ S
8 I. d% w; x+ F, M3 k

artic

你的Uap和Uad之间数量级差别太大，不变换z轴坐标的话，数量级较大的会“吃掉”较小的，做个log变换就好了。
0 R9 i7 w& l: X( mmesh(ap,ad,log(Uap + abs(min(min(Uap)))),'FaceColor','b','FaceAlpha',0.5,'MeshStyle','none')
hold on;
! v. a* E0 F3 n2 O& xmesh(ap,ad,log(Uad + abs(min(min(Uad)))),'FaceColor','r','FaceAlpha',0.5,'MeshStyle','none')