clear
%x = [3 4 6 8 11 10 9 7 5];
%load ecg;
x = sin([1:200]*0.01*pi);
del = 0.08;
x_hat = zeros(1,length(x));
x_hat_m1 = 0;
alpha = 0.95;

for n = 1:length(x),
   if n == 1
      x_tilde(n) = alpha*x_hat_m1
   else
      x_tilde(n) = alpha*x_hat(n-1);
   end
   d(n) = x(n) - x_tilde(n);
   d_hat(n) = quant(del, d(n));
   x_hat(n) = x_tilde(n) + d_hat(n);
end