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

zophim

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

* v5 _4 T# P+ `7 s>> 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;
1 W' o/ R. a; H# n" u[ap,ad]=meshgrid(ap,ad);
Pp=-((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);
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);
. m3 g0 E8 o1 dnd1=(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);
, L" Y) ]* B( I$ p& TUap=Pp.*np1-ca.*(nd1+np1);
$ G% v) ?: R- Q: S6 {- D9 BPd =-((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);
5 R2 I- o1 U1 d5 l$ z; \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);
' \  v( v' S+ ?1 N( m. }/ W( UUad=Pd.*nd2-ca.*(nd2+np2);
) K& Y- T0 _$ U>> mesh(ap,ad,Uap)
>> hold on
4 m, M0 f* [6 J# o2 |' F>> mesh(ap,ad,Uad)
, S7 u* j8 R4 A0 q3 ?1 A: M$ f>> hold off
1 H+ I+ B8 b9 B- J0 [>>
出来的图片如下：
0 E' t6 u' A4 I; }
+ @; @. [+ G$ Q# A0 n

( U6 O+ g, E5 J, X

artic

你的Uap和Uad之间数量级差别太大，不变换z轴坐标的话，数量级较大的会“吃掉”较小的，做个log变换就好了。
" I+ S; [) S  l, W7 ]mesh(ap,ad,log(Uap + abs(min(min(Uap)))),'FaceColor','b','FaceAlpha',0.5,'MeshStyle','none')
5 ~7 l+ g: a, H# p6 w  M6 p6 [hold on;
$ M- g& F, o1 }mesh(ap,ad,log(Uad + abs(min(min(Uad)))),'FaceColor','r','FaceAlpha',0.5,'MeshStyle','none')