Here are the frequency and the time-frequency analyses of the Chirp signal:
Fs = 1000;
L = 2000;
T = 1/Fs;
t = (0:L-1)*T;
y1 = chirp(t,50,1,60);
y2 = chirp(t,50,1,70);
y3 = chirp(t,50,1,80);
y4 = chirp(t,50,1,90);
y5 = chirp(t,50,1,100);
y = y1+y2+y3+y4+y5;
figure('Name', 'Freq Analysis')
pspectrum(y, Fs) % Fre. Analysis
[p,f,t] = pspectrum(y,Fs,'spectrogram');
figure('Name', 'Time-Freq Analysis')
waterfall(f,t,p')
xlabel('Frequency (Hz)')
ylabel('Time (seconds)')
wtf = gca;
wtf.XDir = 'reverse';
view([30 45])
