subplot中的图例

%% 根据云参数，生成云滴 ( 降水指标 )
clc;  clear;  close all;
DropN = 3000;   % 需要产生的云滴数目
Ex = [0.974738631575396,0.885884456986287,1.51764167435718,3.50383459588321,0.283601590079623;3.06679138513355,1.23197926358265,4.70566433440631,10.5534827517891,1.84722222222222;6.84029258893687,1.54887884409856,10.1471948371156,22.7961185132027,7.06362063214260];
En = [0.697350917852717,0.115364935532121,1.06267422001638,2.34988271863529,0.521206877380866;1.25783373460111,0.105633193505304,1.81384350090311,4.08087858713788,1.73879946997346;1.25783373460111,0.105633193505304,1.81384350090311,4.08087858713788,1.73879946997346];
He = [0.0697350917852717,0.0115364935532121,0.106267422001638,0.234988271863529,0.0521206877380866;0.125783373460111,0.0105633193505304,0.181384350090311,0.408087858713788,0.173879946997346;0.125783373460111,0.0105633193505304,0.181384350090311,0.408087858713788,0.173879946997346];
xrmin = [0, 0.5, 0, 0, 0];       
xrmax = [8, 2, 13, 27, 10];
txtx = [0, 0.5, 0, 0, 0];
txts = {('(a1)')  ('(a2)')  ('(a3)')  ('(a4)')  ('(a5)')};
pclr = hsv(3);  % 颜色条
for   i =1:3  % 画第1行
       subplot('Position',[0.1+0.3*(i-1) 0.7 0.25 0.25]) % 定位坐标位置及大小
       for  j =1:3  
             ex = Ex( j, i );
             en = En( j, i );
             he = He( j, i );
             Enp = normrnd(en, he,  DropN, 1);
             tmp = normrnd(ex, Enp, DropN, 1);
             miu = exp(-0.5*((tmp-ex)./Enp).^2);         
             if   j ==1
                  site = find(tmp<ex);
                  miu(site) = 1;
             end
             if  j == 3
                  site = find(tmp>ex);
                  miu(site) = 1;
             end
             set(gcf,'color',[1,1,1])
             scatter(tmp, miu, 3, pclr(j,:), 'filled','MarkerFaceAlpha',0.7);
             set(gca,'LineWidth',1.0,'Fontname','Times New Roman','FontWeight','Bold','FontSize',14)
             text(txtx(i), 1.08 ,txts(i) ,'color','k','FontName','Times New Roman','FontSize',12)
             xlim([xrmin(i), xrmax(i)]);        
             hold on; box on; 
       end   
end

for   i =4:5  % 画第2行
       subplot('Position',[0.1+0.3*(i-4) 0.37 0.25 0.25]) % 定位坐标位置及大小
       for  j =1:3  
             ex = Ex( j, i );
             en = En( j, i );
             he = He( j, i );
             Enp = normrnd(en, he,  DropN, 1);
             tmp = normrnd(ex, Enp, DropN, 1);
             miu = exp(-0.5*((tmp-ex)./Enp).^2);         
             if   j ==1
                  site = find(tmp<ex);
                  miu(site) = 1;
             end
             if  j == 3
                  site = find(tmp>ex);
                  miu(site) = 1;
             end
             set(gcf,'color',[1,1,1])
             h(j) =  scatter(tmp, miu, 3, pclr(j,:), 'filled','MarkerFaceAlpha',0.7);  % 赋值给句柄，方便后面加图例
             set(gca,'LineWidth',1.0,'Fontname','Times New Roman','FontWeight','Bold','FontSize',14)
             text(txtx(i), 1.08 ,txts(i) ,'color','k','FontName','Times New Roman','FontSize',12)
             xlim([xrmin(i), xrmax(i)]);        
             hold on; box on; 
       end   
end
legd = legend([h(1),h(2),h(3)],{'I','II','III'},'Box','off','Color','w','Position',[0.75 0.45 0.1 0.1],'FontSize',14); % 设置图例属性