filip33Print

%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');