%filtr
figure(2)
M=301;
h=fir1(M-1,0.6); % 0.06 wyszlo mi z proporcji
th= 0: 1/fpy:(M-1)/fpy;
subplot(2,2,1);
plot(th,h);
xlabel('czas[s]');
ylabel('h(t)');
title('odp imp filtru');

%modul filtru
subplot(2,2,2);
Nfh=2^9 % wieksze niz cos
Nf21h=Nfh/2+1;
vh=fft(h,Nfh);
wh=abs(vh);
fh= linspace(0,fpy/2,Nf21h);
plot(fh,wh(1:Nf21h));
xlabel('czest[Hz]');
ylabel('|H(f)|');
title('Modul widma odp imp');

% sygnal po filtracji

z=filter(h,1,y);
subplot(223);
plot(ty,z);
xlabel('czas[s]');
ylabel('z(t)');
title('sygnal po filtracji');

%modul widma po filtracji
vhy=fft(z,Nfy);
why=abs(vhy);
subplot(224);
plot(fy,why(1:Nf21y));
xlabel('czest [Hz]');
ylabel('|z(t)|');
title('modul widma po filtracji');