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

Убедитесь, что комплексная чирпа имеет действительную и мнимую части, которые равны, но имеют ${\mathrm{}}^{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)