0001
0002
0003 function [PltPolForTimeInter, U, L] = getRobustPltPolForTimeInter(ALL_Family_plRec, t_start, t_end, own_force_id, enemy_cmp_id, pltNo, alternativeNumber)
0004
0005
0006
0007
0008 [noOfParticles, noOfstructs] = size(ALL_Family_plRec{1}{1,1});
0009 U = []; L = []; PltPolForTimeInter = [];
0010 for t = t_start:t_end
0011 curr_alt = [];
0012 for p = 1:noOfParticles
0013 curr_alt = [curr_alt; ALL_Family_plRec{t}{own_force_id, enemy_cmp_id}{p, pltNo+1}.PltPol(alternativeNumber)];
0014 end;
0015 PltPolForTimeInter = [PltPolForTimeInter; mean(curr_alt)];
0016 U = [U; max(curr_alt) - mean(curr_alt)];
0017 L = [L; mean(curr_alt) - min(curr_alt)];
0018 end;
0019 T = [t_start:t_end]';
0020 errorbar(T, PltPolForTimeInter, L, U)