for index=15:800
    temp = 0.0;
  for i = 1:15
    temp = temp + h1(i)*xr(index - (i-1));
  end
  yInterp(index) = temp;
end



 %a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb)
   %                       - a(2)*y(n-1) - ... - a(na+1)*y(n-na)