clc;
clear;

fs=48000;                       %vzorkovaci frekvence
f=500;                          %pozadovana frekvence
Q=2*pi*f/fs;
B=2*cos(Q);
A=sin(Q);

A=floor(A*32000)/32000;         %ukazka vlivu kvantizace koeficientu
B=floor(B*32000)/32000;

N=1000;

y4=zeros(1,N);                  %vystupni pole 1x1000 
y1=A;
y2=0;
y3=0;

for n=1:N
  y2=y1;
  y1=y3;
  y3=y1*B-y2;
  y4(n)=-y3;
end
    
axis([0 1000 -1.2 1.2]);
hold;
plot(y4);
