Генерация сигнала с частотой 600 Гц в течение 2 секунд. Сигнал состоит из чирпа с синусоидально изменяющимся частотным содержанием.

fs = 6e2;
t = 0:1/fs:2;
x = vco(sin(2*pi*t),[0.1 0.4]*fs,fs);

Сохраните сигнал в неформатированном массиве глубокого обучения. Вычислите кратковременное преобразование Фурье сигнала. Введите время выборки как duration скаляр. (В качестве альтернативы введите частоту дискретизации в виде числового скаляра.) Укажите, что входной массив находится в 'CTB' формат.

dlx = dlarray(x);
[yr,yi,f,t] = dlstft(dlx,seconds(1/fs),'DataFormat','CTB');

Преобразуйте выходные данные в числовые массивы. Вычислите величину кратковременного преобразования Фурье и отобразите ее как график водопада.

yr = extractdata(yr);
yi = extractdata(yi);
f = extractdata(f);
t = seconds(t);
ax = newplot;
waterfall(ax,f,t,squeeze(hypot(yr,yi))')
ax.XDir = 'reverse';
view(30,45)
ylabel('Time (s)')
xlabel('Frequency (Hz)')
zlabel('Magnitude')

Создайте $\mathrm{}\mathrm{}\mathrm{}$массив 3 × 160 (× 1), содержащий одну партию трехканального 160-образного синусоидального сигнала. Нормализованные синусоидальные частоты $\mathrm{составляют}\mathrm{}$δ/4 рад/образец$,\mathrm{}$λ/2 рад/образец $$\mathrm{и}\mathrm{}$$3λ/4 рад/образец. Сохранение сигнала в виде dlarray, указывая размеры по порядку. dlarray изменяет размеры массива на 'CBT' форма, ожидаемая сетью глубокого обучения. Отображение размеров массива.

x = dlarray(cos(pi.*(1:3)'/4*(0:159)),'CTB');
[nchan,nbtch,nsamp] = size(x)

nchan = 3

nbtch = 1

nsamp = 160

Вычислите короткое преобразование Фурье сигнала с глубоким обучением. Укажите прямоугольное окно из 64 образцов и длину БПФ 1024.

[re,im,f,t] = dlstft(x,'Window',rectwin(64),'FFTLength',1024);

dlstft вычисляет преобразование вдоль 'T' размерность. Выходные массивы находятся в 'SCBT' формат. 'S' размерность соответствует частоте в кратковременном преобразовании Фурье.

Извлеките данные из массивов глубокого обучения.

re = squeeze(extractdata(re));
im = squeeze(extractdata(im));
f = extractdata(f);
t = extractdata(t);

Вычислите величину кратковременного преобразования Фурье. Постройте график величины отдельно для каждого канала на участке водопада.

z = abs(re + 1j*im);

for kj = 1:nchan
    subplot(nchan,1,kj)
    waterfall(f/pi,t,squeeze(z(:,kj,:))')
    view(30,45)
end

xlabel('Frequency (\times\pi rad/sample)')
ylabel('Samples')