clear all;
close all;

% parametry osi czasu, generowanie osi czasu

N=1000;
fp=1000;
%t = 0:1/fp:(N-1)/fp;
t=linspace(0,(N-1)/fp,N);

% suma sinusow i szumu

x=6*sin(2*pi*250*t)+3*sin(2*pi*350*t)+2*randn(1,N);
subplot(211)
plot(t,x)
xlabel('czas[s]')
ylabel('sygnal')

% widmo

Nf=2^12;
v=fft(x,Nf);                % widmo zespolone
w=abs(v);                   % widmo amplitudowe

N21=Nf/2+1;
f=linspace(0,fp/2,N21);        % generowanie osi czestotliwosci
subplot(212)
plot(f,w(1:N21))
xlabel('czestotliwosc[Hz]')
ylabel('Modul widma')
%_________________________________________


fgr = 300;
fgrN = fgr / (fp/2);

M = 3001;

h = fir1(M-1, fgrN, 'low');

pof=filter(x,1,h);
pof_tra=fft(pof,Nf);
pof_tra_abs=abs(pof_tra);
figure(2);

grid on;
plot(f(1:N21),pof_tra_abs(1:N21))



















