sm = [0];
mu = mean(ss.usdr.usdr,2);
a = 1/500;
for i = 2:length(mu)
    sm = [sm a*mu(i) + (1-a)*(sm(i-1))];
end
%%
std(sm(ss.usdr.bin > 1200))