% filter of system output and input

function yf=filt1(X,fil,R,F)

R=inv(R');

for i=1:2*F
     R(i,:)=R(i,:)/(abs(R(i,:))*ones(2*F,1));
     end
Z=R*X';
yf=Z(fil);