Сгенерируйте сигнал, произведенный на уровне 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);

waterfall(f,t,squeeze(hypot(yr,yi))')

ax = gca;
ax.XDir = 'reverse';
view(30,45)

ylabel('Time (s)')
xlabel('Frequency (Hz)')
zlabel('Magnitude')

Сгенерируйте 3 160 (-by-1) массив, содержащий один пакет синусоидального сигнала с 160 выборками, с тремя каналами. Нормированные частоты синусоиды $\pi /4$ рад/отсчет, $\pi /2$ рад/отсчет, и $$3\pi /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 выборками и длину БПФ 1 024.

[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')