matlab怎么给个图例标注,如何对同一坐标系下的三个曲面进行图例标注

说明：Question(1)我们知道在同一坐标系如果画两条不同的曲线，可以用语句：legende(........)注释图例，比如可以用普通线条注释第一条曲线，另外一条可用圈圈，或者“+”号来画第二条曲线(总之为了书面打印时能区分这两条不同曲线就好)，但是如果是三个曲面同时在一个坐标系出现，就不能像曲线那样一个用线标注，一个用其他的什么符号标注了。如果此时要对这三个曲面用legend(….)语句进行注释要怎样区别这三个曲面？(当然要在书面打印时能够区分哪个是第一个曲面，哪个是第二个曲面，哪个是第三个曲面)，还要求在图例中标注：a表达式表示：F(t,s,m)对m的偏导数  ，b表达式表示：F(t,s,m)对t的偏导数 c表达式表示：F(t,s,m)对s的偏导数 ，如何满足上述要求？请各位老师指教！该三个曲面的代码如下：  a=vectorize('(exp((3*m*(s - 1))/2) - cos((3^(1/2)*m*(s - 1))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - 1))/2))*((3*t*cos((3^(1/2)*m*t)/2))/2 - (3*t*exp(-(3*m*t)/2))/2 + (3^(1/2)*t*sin((3^(1/2)*m*t)/2))/2) - ((3*exp(-(3*m)/2))/2 + (3*cos((3^(1/2)*m)/2))/2 - (3^(1/2)*sin((3^(1/2)*m)/2))/2)*(cos((3^(1/2)*m*(s - t))/2) - exp((3*m*(s - t))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - t))/2)) + (s/2 - t/2)*(3*exp((3*m*(s - t))/2) - 3*cos((3^(1/2)*m*(s - t))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - t))/2))*(cos((3^(1/2)*m)/2) - exp(-(3*m)/2) + 3^(1/2)*sin((3^(1/2)*m)/2)) + (s/2 - 1/2)*(exp(-(3*m*t)/2) - cos((3^(1/2)*m*t)/2) + 3^(1/2)*sin((3^(1/2)*m*t)/2))*(3*exp((3*m*(s - 1))/2) + 3*cos((3^(1/2)*m*(s - 1))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - 1))/2))');a = eval(['@(t,s,m)',vectorize(a)]);  %a表示F(t,s,m)对m的偏导数[t,s,m] = meshgrid(0:.01:1,0:.01:1,0:.01:3.0167);  % 画图范围v = a(t,s,m);v(s>t) = NaN;  v(s==0)=NaN;  v(t==0)=NaN;h = patch(isosurface(s,t,m,v,0));isonormals(t,s,m,v,h)   set(h,'FaceColor','r','EdgeColor','none');xlabel('t');ylabel('s');zlabel('m');alpha(0.8)   grid on; view(3); axis equal; camlight; lighting gouraudhold onb=vectorize('(exp((3*m*(s - 1))/2) - cos((3^(1/2)*m*(s - 1))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - 1))/2))*((3*m*cos((3^(1/2)*m*t)/2))/2 - (3*m*exp(-(3*m*t)/2))/2 + (3^(1/2)*m*sin((3^(1/2)*m*t)/2))/2) - (m*(3*exp((3*m*(s - t))/2) - 3*cos((3^(1/2)*m*(s - t))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - t))/2))*(cos((3^(1/2)*m)/2) - exp(-(3*m)/2) + 3^(1/2)*sin((3^(1/2)*m)/2)))/2');b = eval(['@(t,s,m)',vectorize(b)]);  %b表示F(t,s,m)对t的偏导数[t,s,m] = meshgrid(0:.01:1,0:.01:1,0:.01:3.0167);  % 画图范围v = b(t,s,m);v(s>t) = NaN;  v(s==0)=NaN;  v(t==0)=NaN;h = patch(isosurface(t,s,m,v,0));isonormals(t,s,m,v,h)   set(h,'FaceColor','r','EdgeColor','none');xlabel('t');ylabel('s');zlabel('m');alpha(0.8)   grid on; view(3); axis equal; camlight; lighting gouraudhold onc=vectorize('(m*(exp(-(3*m*t)/2) - cos((3^(1/2)*m*t)/2) + 3^(1/2)*sin((3^(1/2)*m*t)/2))*(3*exp((3*m*(s - 1))/2) + 3*cos((3^(1/2)*m*(s - 1))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - 1))/2)))/2 + (m*(3*exp((3*m*(s - t))/2) - 3*cos((3^(1/2)*m*(s - t))/2) + 3^(1/2)*sin((3^(1/2)*m*(s - t))/2))*(cos((3^(1/2)*m)/2) - exp(-(3*m)/2) + 3^(1/2)*sin((3^(1/2)*m)/2)))/2');c= eval(['@(t,s,m)',vectorize(c)]);  %c表示F(t,s,m)对s的偏导数[t,s,m] = meshgrid(0:.01:1,0:.01:1,0:.01:3.0167);  % 画图范围v = c(t,s,m);v(s>t) = NaN;  v(s==0)=NaN;  v(t==0)=NaN;h = patch(isosurface(t,s,m,v,0));isonormals(t,s,m,v,h)   set(h,'FaceColor','r','EdgeColor','none');xlabel('t');ylabel('s');zlabel('m');alpha(0.8)   grid on; view(3); axis equal; camlight; lighting gouraudhold off