filtr dolnoprzepsustowy(gotowy)

% filtr dolnoprzepsustowy częstotliwości do 317Hz fs=1798Hz
fs=1798;
f=500; %widmo większe od 317
N=1798;
%generwoanie sygnalu
t = 0:1/fs:(N-1)/fs;%wektor czasu
%t = (0:N-1)'/fs;
y=sin(2*pi*f*t);
subplot(311)
plot(t,y);
%widmo amplitudowe bez filtru
Y=fft(y);
x = (-N/2:(N-1)/2)*fs/N;
subplot(312)
%skala jest decybelowa
%oś x - czestotliwosc
% oś y - amplituda
%plot(x, 20*log10(abs(fftshift(Y))));
plot(x, 20*log10(abs((Y))));
%widmo z filtrem
%b - wspólczynniki
Nf=length(Y);
b=[ones(1,317) zeros(1,(Nf-317))];
Yf=filter(b,1,Y);
subplot(313)
N2=length(Yf);
x2=(-N2/2:(N2-1)/2)*fs/N2;
% oś x - częstotliwosc
% oś y - amplituda
plot(x2, 20*log10(abs(fftshift(Yf))));