Сгенерируйте щебет с линейным мгновенным отклонением частоты. Щебет отбирают с частотой дискретизации 1 кГц в течение 2 секунд. Мгновенная частота составляет 0 при t = 0 и пересекает 250 Гц при t = 1 секунде.

t = 0:1/1e3:2;
y = chirp(t,0,1,250);

Вычислите и постройте график спектрограммы щебета. Разделите сигнал на сегменты так, чтобы разрешение по времени составляло 0,1 секунды. Задайте 99% перекрытия между смежными сегментами и спектральной утечкой 0,85.

pspectrum(y,1e3,'spectrogram','TimeResolution',0.1, ...
    'OverlapPercent',99,'Leakage',0.85)

Сгенерируйте щебет с квадратичным мгновенным отклонением частоты. Щебет отбирают с частотой дискретизации 1 кГц в течение 2 секунд. Мгновенная частота составляет 100 Гц при t = 0 и пересекает 200 Гц при t = 1 секунде.

t = 0:1/1e3:2;
y = chirp(t,100,1,200,'quadratic');

Вычислите и постройте график спектрограммы щебета. Разделите сигнал на сегменты так, чтобы разрешение по времени составляло 0,1 секунды. Задайте 99% перекрытия между смежными сегментами и спектральной утечкой 0,85.

pspectrum(y,1e3,'spectrogram','TimeResolution',0.1, ...
    'OverlapPercent',99,'Leakage',0.85)

Сгенерируйте выпуклый квадратичный щебет, выбранный с частотой дискретизации 1 кГц в течение 2 секунд. Мгновенная частота составляет 400 Гц при t = 0 и пересекает 300 Гц при t = 1 секунде.

t = 0:1/1e3:2;
fo = 400;
f1 = 300;
y = chirp(t,fo,1,f1,'quadratic',[],'convex');

Вычислите и постройте график спектрограммы щебета. Разделите сигнал на сегменты так, чтобы разрешение по времени составляло 0,1 секунды. Задайте 99% перекрытия между смежными сегментами и спектральной утечкой 0,85.

pspectrum(y,1e3,'spectrogram','TimeResolution',0.1, ...
    'OverlapPercent',99,'Leakage',0.85)

Сгенерируйте вогнутый квадратичный щебет, выбранный с частотой дискретизации 1 кГц в течение 4 секунд. Задайте временной вектор так, чтобы мгновенная частота была симметричной относительно полуточечной точки интервала дискретизации с минимальной частотой 100 Гц и максимальной частотой 500 Гц.

t = -2:1/1e3:2;
fo = 100;
f1 = 200;
y = chirp(t,fo,1,f1,'quadratic',[],'concave');

Вычислите и постройте график спектрограммы щебета. Разделите сигнал на сегменты так, чтобы разрешение по времени составляло 0,1 секунды. Задайте 99% перекрытия между смежными сегментами и спектральной утечкой 0,85.

pspectrum(y,t,'spectrogram','TimeResolution',0.1, ...
    'OverlapPercent',99,'Leakage',0.85)

Сгенерируйте логарифмический щебет, выбранный с частотой дискретизации 1 кГц в течение 10 секунд. Мгновенная частота составляет 10 Гц первоначально и 400 Гц в конце.

t = 0:1/1e3:10;
fo = 10;
f1 = 400;
y = chirp(t,fo,10,f1,'logarithmic');

Вычислите и постройте график спектрограммы щебета. Разделите сигнал на сегменты так, чтобы разрешение по времени составляло 0,2 секунды. Задайте 99% перекрытия между смежными сегментами и спектральной утечкой 0,85.

pspectrum(y,t,'spectrogram','TimeResolution',0.2, ...
    'OverlapPercent',99,'Leakage',0.85)

Используйте логарифмическую шкалу для оси частоты. Спектрограмма становится линией с высокой неопределенностью на низких частотах.

ax = gca;
ax.YScale = 'log';

Сгенерируйте комплексный линейный щебет, выбранный с частотой дискретизации 1 кГц в течение 10 секунд. Мгновенная частота составляет -200 Гц первоначально и 300 Гц в конце. Начальная фаза равна нулю.

t = 0:1/1e3:10;
fo = -200;
f1 = 300;

y = chirp(t,fo,t(end),f1,'linear',0,'complex');

Вычислите и постройте график спектрограммы щебета. Разделите сигнал на сегменты так, чтобы разрешение по времени составляло 0,2 секунды. Задайте 99% перекрытия между смежными сегментами и спектральной утечкой 0,85.

pspectrum(y,t,'spectrogram','TimeResolution',0.2, ...
    'OverlapPercent',99,'Leakage',0.85)

Проверьте, что комплексный щебет имеет действительные и мнимые части, которые равны, но с ${90}^{\circ }$ различие фаз.

x = chirp(t,fo,t(end),f1,'linear',0) + 1j*chirp(t,fo,t(end),f1,'linear',-90);

pspectrum(x,t,'spectrogram','TimeResolution',0.2, ...
    'OverlapPercent',99,'Leakage',0.85)