N= 20000; fp=10000; t=linspace(0,(N-1)/fp,N); %t=0:1/fp:(N-1)/fp %x=2*sin(2*pi*500*t)+5*sin(2*pi*2000*t); %x=randn(1,N); %1 os cz 2 fp 3 2 punkt kontrolny dl sygnalu f/fp x=chirp(t,1000,2,4000); subplot(2,1,1); plot(t,x); xlabel('czas[s]'); ylabel('syg'); Nf=2048; N21=Nf/2+1; v=fft(x,Nf); w=abs(v); f=linspace(0,fp/2,N21); %dzieli ekran w k odwolanie subplot(2,1,2); plot(f,w(1:N21)); xlabel('czestotliwosc[Hz]'); ylabel('Modul widma'); figure(2); spectrogram(x,256,128,256,fp);
niobe666