View difference between Paste ID: uAN2VcYb and gzzxHiqJ

-%część 2:
+%część 3:
 
-%[x,fpx]=wavread('/home/c/mbi04czap.wav');
+%sinusoida
-[x,fpx]=audioread('/mbi04czap.wav');
+N=4000;
-Nx=length(x);
+fp=10000;
-tx=0:1/fpx:(Nx-1)/fpx;
+t=0:1/fp:(N-1)/fp;
-subplot(2,2,1);
+x=sin(2*pi*740*t);
-plot(tx,x);
+subplot(331);
 plot(t,x);
-ylabel('sygnal');
+
 ylabel('sinus');
 
 
 
 
 %histogram x
-Nfx=2^15;
+subplot(332);
-N21x=Nfx/2+1;
+nbins=51;
-v=fft(x,Nfx);
+hist(x,nbins);
-wx=abs(v);
+xlabel('wartosc probki');
-fx=linspace(0,fpx/2,N21x);
+ylabel('liczba probek w przedziale');
-subplot(2,2,2);
+
-plot(fx,wx(1:N21x));
+
-xlabel('czestotliwosc[Hz]');
+
-ylabel('modul widma');
+
 
 % 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")
-dr=8;
+tr=-kmax/fp:1/fp:kmax/fp;
-y=decimate(x,dr);
+plot(tr,rx);
-Ny=length(y);
+xlabel('przesuniecie [s]');
-fpy=fpx/dr;
+ylabel('autokorelacja');
-ty=0:1/fpy:(Ny-1)/fpy;
+
-subplot(2,2,3);
+
-plot(ty,y);
+
-xlabel('czas[s]')
+
-ylabel('sygnal po decymacji');
+
 
 
 
 
 
 %szum o rozkładzie równomiernym
 subplot(334);
-Nfy=2^13;
+y=rand(1,N);
-N21y=Nfy/2+1;
+plot(t,y);
-v=fft(y,Nfy);
+
-wy=abs(v);
+ylabel('szum rownomierny');
-fy=linspace(0,fpy/2,N21y);
+
-subplot(2,2,4);
+
-plot(fy,wy(1:N21y));
+
-xlabel('czestotliwosc[Hz]');
+
-ylabel('modul widma sygnalu po decymacji');
+
 
 %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');