0001 function stat_plot_buffer(model,nb_mob,export,format) 0002 0003 res.models = { 'street' 'city' 'hexagonal' }; 0004 res.model = char(res.models(model)); 0005 0006 res.mob = nb_mob; 0007 0008 res.tmax = 60; 0009 0010 for cov = 1:20 0011 res.coverages = 5:5:100; 0012 res.coverage = res.coverages(cov); 0013 0014 buffer_mcn_max_mean = load(sprintf('%s/%dusr/buffer-mcn-max-mean.dat',res.model,res.mob)); 0015 buffer_mcn_max_std = load(sprintf('%s/%dusr/buffer-mcn-max-std.dat',res.model,res.mob)); 0016 0017 buffer_mcn_mean_mean = load(sprintf('%s/%dusr/buffer-mcn-mean-mean.dat',res.model,res.mob)); 0018 buffer_mcn_mean_std = load(sprintf('%s/%dusr/buffer-mcn-mean-std.dat',res.model,res.mob)); 0019 0020 buffer_std_max_mean = load(sprintf('%s/%dusr/buffer-std-max-mean.dat',res.model,res.mob)); 0021 buffer_std_max_std = load(sprintf('%s/%dusr/buffer-std-max-std.dat',res.model,res.mob)); 0022 0023 buffer_std_mean_mean = load(sprintf('%s/%dusr/buffer-std-mean-mean.dat',res.model,res.mob)); 0024 buffer_std_mean_std = load(sprintf('%s/%dusr/buffer-std-mean-std.dat',res.model,res.mob)); 0025 0026 end 0027 0028 figure(7) 0029 plot([5:5:100],buffer_mcn_max_mean,'--b',[5:5:100],buffer_mcn_mean_mean,'-b',[5:5:100],buffer_std_max_mean,'--r',[5:5:100],buffer_std_mean_mean,'-r') 0030 axis([0 100 0 max(buffer_mcn_max_mean)*1.1]) 0031 grid ('on') 0032 xlabel('coverage in %') 0033 ylabel('number packets') 0034 legend('maximum for mcn','mean for mcn','maximum for std','mean for std') 0035 title(sprintf('buffer occupation of the terminals (%s model, %d users, %d min)',res.model,res.mob,res.tmax)) 0036 if export == 1 0037 exportfig(gcf,[res.model '-' num2str(res.mob) 'usr-buffer-max-mean' ],'color', 'cmyk','format',format); 0038 end 0039 0040 figure(8) 0041 hold on 0042 b1 = errorbar([5:5:100],buffer_mcn_max_mean,buffer_mcn_max_std,'b'); 0043 b2 = errorbar([5:5:100],buffer_std_max_mean,buffer_std_max_std,'r'); 0044 grid ('on') 0045 axis([0 100 0 max(buffer_mcn_max_mean)*1.5]) 0046 xlabel('coverage in %') 0047 ylabel('number packets') 0048 legend([b1,b2],'maximum for mcn','maximum for std') 0049 title(sprintf('maximum buffer occupation of the terminals (%s model, %d users, %d min)',res.model,res.mob,res.tmax)) 0050 hold off 0051 if export == 1 0052 exportfig(gcf,[res.model '-' num2str(res.mob) 'usr-buffer-max-error' ],'color', 'cmyk','format',format); 0053 end 0054 0055 figure(9) 0056 hold on 0057 b3 = errorbar([5:5:100],buffer_mcn_mean_mean,buffer_mcn_mean_std,'b'); 0058 b4 = errorbar([5:5:100],buffer_std_mean_mean,buffer_std_mean_std,'r'); 0059 grid ('on') 0060 axis([0 100 0 max(buffer_mcn_max_mean)]) 0061 xlabel('coverage in %') 0062 ylabel('number packets') 0063 legend([b3,b4],'mean for mcn','mean for std') 0064 title(sprintf('average buffer occupation of the terminals (%s model, %d users, %d min)',res.model,res.mob,res.tmax)) 0065 hold off 0066 if export == 1 0067 exportfig(gcf,[res.model '-' num2str(res.mob) 'usr-buffer-mean-error' ],'color', 'cmyk','format',format); 0068 end