该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
close all;
reader_m;
topodir = 'Topo';
dapcname = 'Results_DAPC';
filterdir = 'Filter';
figdir = 'Figures';
for runs = 1:length(nn)
meanplot_dapc = [];
stdplot_dapc = [];
spplot_dapc = [];
mpplot_dapc = [];
for mindists = 1:length(mindist)
tempname = strcat(dapcname, '_MEAN_R', num2str(nn(runs)), '_D', num2str(mindist(mindists)));
tempname = strrep(tempname, '.', '_');
(['meanplot_dapc = [meanplot_dapc ' tempname ']']);
tempname = strcat(dapcname, '_STD_R', num2str(nn(runs)), '_D', num2str(mindist(mindists)));
tempname = strrep(tempname, '.', '_');
(['stdplot_dapc = [stdplot_dapc ' tempname ']']);
tempname = strcat(dapcname, '_STD_PERCENTAGE_R', num2str(nn(runs)), '_D', num2str(mindist(mindists)));
tempname = strrep(tempname, '.', '_');
(['spplot_dapc = [spplot_dapc ' tempname ']']);
tempname = strcat(dapcname, '_MEAN_PERCENTAGE_R', num2str(nn(runs)), '_D', num2str(mindist(mindists)));
tempname = strrep(tempname, '.', '_');
(['mpplot_dapc = [mpplot_dapc ' tempname ']']);
end
figure;
h = plot(mindist, stdplot_dapc ,'b-s', 'lineWidth', 1);
titlename = sprintf('Standard deviation of average detection range with %s nodes', num2str(nn(runs)));s
xlabel('minimum distance between readers (meters)');
ylabel('Standard deviation');
legend('DAPC');
set(gca,'XTick', mindist);
plotname = strcat(figdir, '/STD_PERC_FIX_SIZE_VAR_DENSITY_', num2str(nn(runs)), '.emf');
(['print -dmeta ' plotname]);
figure;
h = plot(mindist, mpplot_dapc ,'b-s', 'lineWidth', 2);
titlename = sprintf('Average percentage of time achieving desired range with %s nodes', num2str(nn(runs)));
title(titlename);
xlabel('minimum distance between readers (meters)');
ylabel('percentage of time achieving desired range');
legend('DAPC');
axis([3 14 min(min(mpplot_ppc), min(mpplot_dapc)) max(max(mpplot_ppc), max(mpplot_dapc))]);
set(gca,'XTick', mindist);
plotname = strcat(figdir, '/MEAN_PERC_FIX_SIZE_VAR_DENSITY_', num2str(nn(runs)), '.emf');
(['print -dmeta ' plotname]);
figure;
h = plot(mindist, spplot_dapc ,'b-s', 'lineWidth', 2);
titlename = sprintf('Standard deviation of percentage of time achieving desired range with %s nodes', num2str(nn(runs)));
title(titlename);
xlabel('minimum distance between readers (meters)');
ylabel('Standard deviation');
legend('DAPC');
set(gca,'XTick', nn);
plotname = strcat(figdir, '/STD_PERC_FIX_SIZE_VAR_DENSITY_', num2str(nn(runs)), '.emf');
(['print -dmeta ' plotname]);
%
close all;
end