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

zophim

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

>> clear
& O7 d9 d3 a' q: Iap=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;
# X+ ?$ a( q) B" r6 j[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);
1 E! ^5 \) ^, E9 R! Qnd1=(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);
& P) R6 _5 ~$ f& U% uUap=Pp.*np1-ca.*(nd1+np1);
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);
" h# |, z0 k% f2 knp2=(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);
Uad=Pd.*nd2-ca.*(nd2+np2);
>> mesh(ap,ad,Uap)
9 g1 N+ R; v  [$ X4 [>> hold on
2 u# f$ m$ b/ z$ E# ]" i7 P6 t>> mesh(ap,ad,Uad)
; [" e0 w' W. u5 b" k7 m" M, ?>> hold off
>>
出来的图片如下：
) [7 U9 r& x2 H# q+ X: a
; V" T" o: k2 o1 E8 ]' {. c
+ j+ @! e, L8 H  ^0 S! U
: t2 U1 g8 |8 ^4 A( X" S2 _7 a

artic

你的Uap和Uad之间数量级差别太大，不变换z轴坐标的话，数量级较大的会“吃掉”较小的，做个log变换就好了。
( \0 O0 S, s& w  N  y& x# Jmesh(ap,ad,log(Uap + abs(min(min(Uap)))),'FaceColor','b','FaceAlpha',0.5,'MeshStyle','none')
hold on;
, k" K7 g, s& w& e+ S* Xmesh(ap,ad,log(Uad + abs(min(min(Uad)))),'FaceColor','r','FaceAlpha',0.5,'MeshStyle','none')