Сгенерируйте щебет с линейным мгновенным отклонением частоты. Щебет производится на уровне 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)