clear all;
N=12000;
fp=10000;

t=0:1/fp:(N-1)/fp;
x1= 4*chirp(t,100,t(end),450);
x2= 2*chirp(t,150,t(end),600);
x3= 4*sin(2*pi*t*350);
x4= 5*rand(1,N);
x=x1+x2+x3+x4;

figure(1);
%oczyt
subplot(221);
plot(t,x);
xlabel('czas[s]');
ylabel('x(t)');
title('sygnal');

%modul widma
Nf=2^14;
Nf21=Nf/2+1;
f=linspace(0,fp/2,Nf21);
v=fft(x,Nf);
w=abs(v);
subplot(222);
plot(f,w(1:Nf21));
xlabel('czest[Hz]');
ylabel('|x(f)|');
title('modul widma sygnalu');

%decymacja
dr=floor(fp/(2*600));
y=decimate(x,dr);
fpy=fp/dr;
Ny=length(y);
ty=0:1/fpy:(Ny-1)/fpy;
subplot(223);
plot(ty,y);
xlabel('czas[s]');
ylabel('y(t)');
title('sygnal po decymacji');

%modul widma decymacji
Nfy=2^11;
Nf21y=Nfy/2+1;
fy=linspace(0,fpy/2,Nf21y);
vy=fft(y,Nfy);
wy=abs(vy);
subplot(224);
plot(fy,wy(1:Nf21y));
xlabel('czest[Hz]');
ylabel('|y(f)|');
title('modul widma decymacji');

figure(2);
%odp imp filtru
M=301;
h=fir1(M-1,200/(fpy/2));
th=0:1/fpy:(M-1)/fpy;
subplot(221);
plot(th,h);
xlabel('czas[s]');
ylabel('h(t)');
title('odp imp filtru');
%modul widma odp im filtru
Nfh=2^nextpow2(M);
Nf21h=Nfh/2+1;
fh=linspace(0,fpy/2,Nf21h);
vh=fft(h,Nfh);
wh=abs(vh);
subplot(222);
plot(fh,wh(1:Nf21h));
xlabel('czest[Hz]');
ylabel('|h(f)|');
title('modul widma odp imp filtru');

%filtracja
z=filter(h,1,y);
subplot(223);
plot(ty,z);
xlabel('czas[s]');
ylabel('z(t)');

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

figure(3);

kmax=20;
rx=xcorr(z,z,kmax);
tr=-kmax/fpy:1/fpy:kmax/fpy;
subplot(121);
plot(tr,rx);


nbis=35;
subplot(122);
hist(z);



return