fs = 10000;
czas = 5;
t = 0: 1/fs:czas-1/fs;
y3 = chirp(t,0,5,1000,'quadratic',[ ],'concave');
hold on;
subplot(2,1,1);
plot(t, y3);
axis([0 1 -1 1]);
ylabel('Amplituda');
xlabel('Częstotliwość[Hz]');
title('Sygnał sinusoidalny o zmiennej częstotliwości');

subplot(2,1,2);
spectrogram(y3,window(@hamming, 256),0,256,fs,'yaxis');