function [ eps fs] = epsilon( x )
xf=reshape(x,19200,[]);
[m n]=size(xf);
eps=zeros(n,1);
fs=zeros(n,1);
for i=1:n
[fs(i) eps(i)]=find_pss(xf(:,i),0,false);
end