Gdu

N=4000;
fp=1000;
t=0:1/fp:(N-1)/fp;

%sygnał - sinusoida o zmiennej częstotliwości
x=chirp(t, 100, 4.0, 400);
subplot(3,2,1);
plot(t,x);
xlabel('czas[s]');
ylabel('sygnal');

%moduł widma
Nf=2^12;
N21=Nf/2+1;
v=fft(x,Nf);
wx=abs(v);
f=linspace(0,fp/2,N21);
subplot(3,2,2);
plot(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

%odpowiedź impulsowa
M=201;
h=fir1(M-1,0.5);
th=0:1/fp:(M-1)/fp;
subplot(3,2,3);
plot(th,h);
xlabel('czas[s]');
ylabel('odp. impulsowa');

%moduł transmitancji
v=fft(h,Nf);
wh=abs(v);
subplot(3,2,4);
plot(f,wh(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul transmitancji');

%sygnał po filtracji
y=filter(h,1,x);
subplot(3,2,5);
plot(t,y);
xlabel('czas[s]');
ylabel('sygnal po filtracji');

%widmo amplitudy po filtracji
v=fft(y,Nf);
wy=abs(v);
subplot(3,2,6);
plot(f,wy(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('widmo aplitudy po filtracji');



%część 2:

%[x,fpx]=wavread('/home/c/mbi04czap.wav');
[x,fpx]=audioread('/mbi04czap.wav');
Nx=length(x);
tx=0:1/fpx:(Nx-1)/fpx;
subplot(2,2,1);
plot(tx,x);
xlabel('czas[s]');
ylabel('sygnal');

Nfx=2^15;
N21x=Nfx/2+1;
v=fft(x,Nfx);
wx=abs(v);
fx=linspace(0,fpx/2,N21x);
subplot(2,2,2);
plot(fx,wx(1:N21x));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

dr=8;
y=decimate(x,dr);
Ny=length(y);
fpy=fpx/dr;
ty=0:1/fpy:(Ny-1)/fpy;
subplot(2,2,3);
plot(ty,y);
xlabel('czas[s]')
ylabel('sygnal po decymacji');

Nfy=2^13;
N21y=Nfy/2+1;
v=fft(y,Nfy);
wy=abs(v);
fy=linspace(0,fpy/2,N21y);
subplot(2,2,4);
plot(fy,wy(1:N21y));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma sygnalu po decymacji');



%część 3:

%sinusoida
N=4000;
fp=10000;
t=0:1/fp:(N-1)/fp;
x=sin(2*pi*740*t);
subplot(331);
plot(t,x);
xlabel('czas[s]');
ylabel('sinus');

%histogram x
subplot(332);
nbins=51;
hist(x,nbins);
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');

% autokorelacja x
subplot(333);
kmax=1000;
rx=xcorr(x,x,kmax); %(sygnał, sygnał, "jak daleko w przeszłość i przyszłość mamy patrzeć żeby szukać korelacji")
tr=-kmax/fp:1/fp:kmax/fp;
plot(tr,rx);
xlabel('przesuniecie [s]');
ylabel('autokorelacja');

%szum o rozkładzie równomiernym
subplot(334);
y=rand(1,N);
plot(t,y);
xlabel('czas[s]');
ylabel('szum rownomierny');

%histogram y
subplot(335);
nbins=51;
hist(y,nbins);
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');

%autokorelacja y
y=y-mean(y); % e
subplot(336);
ry=xcorr(y,y,kmax);
plot(tr,ry);
xlabel('przesuniecie [s]');
ylabel('autokorelacja');

%szum o rozkładzie normalnym, szum Gaussowski
subplot(337);
z=randn(1,N);
plot(t,z);
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');

%histogram z
subplot(338);
nbins=51;
hist(z,nbins);
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');

%autokorelacja z
subplot(339);
rz=xcorr(z,z,kmax);
plot(tr,rz);
xlabel('wartosc probki');
ylabel('liczba probek w przedziale');