%Szûrõtervezés sorfejtéssel
n=100;
t=-n:n;
Wp=0.2;
b=Wp*sinc(Wp*t);a=1;		
%b=b.*hanning(2*n+1)';          %Karakterisztika simítása
%b=b.*sin(t/100*2*pi*15);
%mod=2*cos(2*pi*t/4);b=b.*mod;  %Karakterisztika eltolása
[H,w]=freqz(b,1);
figure(1);	hold on
plot(w/pi,20*log10(abs(H)+eps),'m'),grid on	%Absolutérték dB
%plot(w/pi,abs(H),'m'),grid on	%Absolutérték lin
figure(2)
plot(w/pi,unwrap(angle(H)))	%Fázis
t=1:1000;
x=sin(2*pi*t/4);
y=filter(b,a,x);
%figure(3);plot(y)