time-frequency analysis of multi-component chirp

1 view (last 30 days)
hi,
I am working on time-frequency analysis of multi-component chirp. I tried to use 'wsst' to separate the chirps but failed (see the figure below). Can anyone help me with this?
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;
wsst(y,Fs);

Accepted Answer

Sulaymon Eshkabilov
Sulaymon Eshkabilov on 19 Feb 2023
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])

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!