dsp.STFT

Кратковременный БПФ

Описание

dsp.STFT объект вычисляет кратковременное преобразование Фурье (STFT) входного сигнала временного интервала. Объект принимает системы координат данных временного интервала, буферизует их к желаемой длине окна и длине перекрытия, умножает выборки на окно, и затем выполняет БПФ на буферизированных окнах. Для получения дополнительной информации см. Алгоритмы.

Используйте STFT, чтобы анализировать содержимое частоты сигнала, который меняется в зависимости от времени.

Создание

Синтаксис

stf = dsp.STFT

stf = dsp.STFT(window)

stf = dsp.STFT(window,overlap)

stf = dsp.STFT(window,overlap,nfft)

stf = dsp.STFT(Name,Value)

Описание

stf = dsp.STFT возвращает объект, stf, это реализует кратковременный БПФ. Объектные процессы информационно-независимо через каждый вход образовывают канал в зависимости от времени.

stf = dsp.STFT(window) возвращает кратковременный объект FFT с набором свойств Окна к window.

stf = dsp.STFT(window,overlap) возвращает кратковременный объект FFT с Window набор свойств к window и набор свойств OverlapLength к overlap.

пример

stf = dsp.STFT(window,overlap,nfft) возвращает кратковременный объект FFT с Window набор свойств к window, OverlapLength набор свойств к overlap, и набор свойств FFTLength к nfft.

stf = dsp.STFT(Name,Value) возвращает кратковременный объект FFT с каждым заданным набором имени свойства к заданному значению. Можно задать дополнительные аргументы пары "имя-значение" в любом порядке.

Свойства

развернуть все

`Window` — Окно Analysis
`sqrt(hann(512,'periodic'))` (значение по умолчанию) | вектор

Окно Analysis в виде вектора из действительных элементов.

Объект буферизует вход в сегменты накладывающегося окна с помощью заданной длины окна и длины перекрытия, и затем умножает каждый перекрытый сегмент на окно.

Настраиваемый: да

Типы данных: single | double

`OverlapLength` — Перекройте длину
256 (значение по умолчанию) | положительное целое число

Количество выборок, которыми последовательные окна перекрываются в виде положительного целого числа. Перекрытие окна уменьшает артефакты на контурах данных.

`FFTLength` — Длина БПФ
512 (значение по умолчанию) | положительное целое число

Длина БПФ в виде положительного целого числа. Это свойство определяет длину STFT выход (количество строк). Длина БПФ должна быть больше или быть равна длине окна.

`FrequencyRange` — Частотный диапазон
`'twosided'` (значение по умолчанию) | `'onesided'`

Частотный диапазон, по которому кратковременный БПФ вычисляется в виде:

'twosided' – Кратковременный БПФ вычисляется для комплексных или действительных входных сигналов. Длина кратковременного БПФ равна значению, которое вы задаете в FFTLength свойство.
'onesided' – Односторонний кратковременный БПФ вычисляется для действительных входных сигналов только. Когда длина БПФ является четной, кратковременной длиной БПФ является FFTLength/2+1. Если длина БПФ является нечетной, длина кратковременного БПФ равна (FFTLength+1)/2.

Использование

Синтаксис

y = stf(x)

Описание

пример

y = stf(x) применяет кратковременный БПФ на вход x и возвращает частотный диапазон выход y.

Входные параметры

развернуть все

`x` — Входной сигнал
вектор | матрица

Входной сигнал временного интервала в виде вектора или матрицы. Если вход является матрицей, объект обрабатывает каждый столбец как независимый канал. Формат кадра (количество строк в x) должно быть равно или меньше, чем длина транзитного участка (длина окна − длина перекрытия).

Вход может быть сигналом переменного размера. Таким образом, формат кадра сигнала может изменить промежуточные вызовы объектного алгоритма, не вызывая release функция. Количество каналов должно остаться то же самое.

Если FrequencyRange свойство установлено в 'onesided', вход должен быть действительным. Если FrequencyRange свойство установлено в 'twosided', вход может быть действительным или комплексным.

Типы данных: single | double
Поддержка комплексного числа: Да

Выходные аргументы

развернуть все

`y` — STFT выход
вектор | матрица

Кратковременный FFT выходной параметр, возвращенный как вектор или матрица.

Если существует достаточно выборок (равно, чтобы скачкообразно переместить длину), чтобы сформировать STFT выход, y FFTLength- N матрица, где N является количеством входных каналов. Если существует недостаточно выборок, чтобы сформировать STFT выход, y isempty.

Тип данных выходных соответствий тот из входного сигнала.

Типы данных: single | double
Поддержка комплексного числа: Да

Функции объекта

`step`	Запустите алгоритм Системного объекта
`release`	Высвободите средства и позвольте изменения в значениях свойств Системного объекта и введите характеристики
`reset`	Сбросьте внутренние состояния Системного объекта
`clone`	Создайте объект дублированной системы
`isLocked`	Определите, используется ли Системный объект
`getFrequencyVector`	Получите вектор из частот, на которых вычисляется кратковременный БПФ

Примеры

свернуть все

Кратковременное спектральное затухание

Скрипт Open Live Script

Кратковременное спектральное затухание достигается путем применения изменяющегося во времени затухания к кратковременному спектру сигнала с шумом. Усиление затухания определяется оценкой шумовой степени в каждом поддиапазоне спектра. Это усиление, когда применено шумный спектр, ослабляет поддиапазоны с более высокой шумовой степенью и подъемниками поддиапазоны с меньшей шумовой степенью.

Вот шаги, вовлеченные в выполнение кратковременного спектрального затухания:

Анализируйте шумный входной сигнал путем вычисления кратковременного преобразования Фурье (STFT).
Умножьте каждый поддиапазон преобразованного сигнала с действительным положительным усилением меньше чем 1.
Синтезируйте denoised поддиапазоны путем взятия обратного кратковременного преобразования Фурье (ISTFT). Сигнал resconstructed является denoised входным сигналом.

Используйте dsp.STFT и dsp.ISTFT объекты вычислить короткое время и обратные кратковременные преобразования Фурье, соответственно.

Шумный входной сигнал

Вход является звуковым сигналом, произведенным на уровне 22 050 Гц. dsp.AudioFileReader возразите читает этот сигнал в системах координат 512 выборок. Звуковой сигнал повреждается белым Гауссовым шумом, который имеет стандартное отклонение 0,05. Используйте audioDeviceWriter возразите, чтобы проигрывать шумный звуковой сигнал к аудио устройству вашего компьютера.

FrameLength = 512;
afr = dsp.AudioFileReader('speech_dft.wav',...
    'SamplesPerFrame',FrameLength);
adw = audioDeviceWriter('SampleRate',afr.SampleRate);

noiseStd = 0.05;
while ~isDone(afr)
    cleanAudio = afr();
    noisyAudio = cleanAudio + noiseStd * randn(FrameLength,1);
    adw(noisyAudio);
end
reset(afr)

Инициализируйте кратковременные и обратные кратковременные объекты преобразования Фурье

Инициализируйте dsp.STFT и dsp.ISTFT объекты. Установите длину окна, равную длине входного кадра и длине транзитного участка к 16. Длина перекрытия является различием между длиной окна и длиной транзитного участка, OL = WL – HL. Установите длину БПФ на 1 024.

WindowLength = FrameLength;
HopLength = 16;
numHopsPerFrame = FrameLength / 16;
FFTLength = 1024;

Окно использовалось для расчета STFT, и ISTFT является периодический hamming окно с длиной 512. ConjugateSymmetricInput флаг istf объект установлен в true, указание, что выход istf объект является сопряжено-симметричным сигналом.

win = hamming(WindowLength,'periodic');
stf = dsp.STFT(win,WindowLength-HopLength,FFTLength);
istf = dsp.ISTFT(win,WindowLength-HopLength,1,0);

Получите средство оценки

Следующий шаг должен задать параметры средства оценки усиления. Это усиление применяется к шумному спектру, чтобы ослабить поддиапазоны с более высокой шумовой степенью и подъемником поддиапазоны с меньшей шумовой степенью.

dec = 16;
alpha = 15;
stftNorm = (sum(win.*win) / dec).^2;

Спектральное затухание

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

while ~isDone(afr)
    cleanAudio =  afr();
    noisyAudio = cleanAudio + noiseStd * randn(FrameLength,1);
    y = zeros(FrameLength,1); % y holds the denoised audio frame
    
    % Feed audio to stft one hop-length at a time
    for index = 1:numHopsPerFrame        
        X = stf(noisyAudio((index-1)*HopLength+1:index*HopLength));        
        % Gain estimator
        Z = abs(X).^2 / (noiseStd^2 * alpha) / stftNorm;
        Z(Z<=1) = 1;
        Z = 1 - 1./Z;
        Z = sign(Z) .* sqrt(abs(Z));
        X = X .* Z;        
        % Convert back to time-domain
        y((index-1)*HopLength+1:index*HopLength) = istf(X);        
    end    
    % Listen to denoised audio:
    adw(y);
end

Совершенная реконструкция

Скрипт Open Live Script

Совершенная реконструкция состоит в том когда выход dsp.ISTFT совпадает с входом к dsp.STFT. Совершенная реконструкция получена если аналитическое окно, $g (n)$ , выполняется постоянное перекрытие - добавляют свойство (COLA) в размере транзитного участка R.

$\sum_{m = - \infty}^{\infty} g (n - mR) = 1, \forall n \in Ζ (g \in COLA (R))$

Сигнал отлично восстановлен если выход dsp.ISTFT возразите совпадает с входом к dsp.STFT объект.

iscola Функция

iscola функционируйте проверяет, что заданное окно и перекрытие удовлетворяют ограничению COLA, чтобы гарантировать, что обратное кратковременное преобразование Фурье (ISTFT) приводит к совершенной реконструкции для немодифицированных спектров. Функция возвращает логический true если комбинация входных параметров совместима КОЛОЙ и логический false если нет. method аргумент функции установлен в 'ola' или 'wola' в зависимости от того, добавляют ли использование метода инверсии взвешенное перекрытие - (ВОЛЯ).

Проверяйте если hann() окном длины 120 выборок и продолжительность перекрытия 60 выборок является совместимый COLA.

winLen = 120;
overlapLen = 60;
win = hann(winLen,'periodic');
tf = iscola(win,overlapLen,'ola')

tf = logical
   1

Инициализация

Инициализируйте dsp.STFT и dsp.ISTFT Системные объекты с этим hann окно, которое является совместимым COLA. Установите длину БПФ равняться длине окна.

frameLen = winLen-overlapLen;
stf = dsp.STFT('Window',win,'OverlapLength',overlapLen,'FFTLength',winLen);
istf = dsp.ISTFT('Window',win,'OverlapLength',overlapLen,'WeightedOverlapAdd',0);

Восстановите данные

Вычислите STFT случайного сигнала. Установите длину входного сигнала равняться длине транзитного участка (длина окна – длина перекрытия). Поскольку окном является совместимый COLA, ISTFT, неизмененные спектры отлично восстанавливают исходный сигнал временной области.

Чтобы подтвердить, сравните вход, x к восстановленному выходу, y. Из-за задержки, введенной объектами, восстановленный выход смещен вовремя по сравнению с входом. Поэтому, чтобы выдержать сравнение, возьмите норму различия между восстановленным выходом, y и предыдущим входом, xprev. Норма очень мала, указывая, что выходной сигнал является отлично восстановленной версией входного сигнала.

n = zeros(1,100);
xprev = 0;
for i = 1:100
    x = randn(frameLen,1);
    X = stf(x);
    y = istf(X);
    n(1,i) = norm(y-xprev);
    xprev = x;
end       
max(abs(n))

ans = 1.6972e-13

ISTFT со взвешенным перекрытием - добавляют (ВОЛЯ)

В ВОЛИ второе окно вызвало окно синтеза, $f (n)$ , применяется после операции IFFT и перед перекрытием - добавляют. Окна синтеза и анализа обычно идентичны и обычно получаются путем пущения квадратного корня из окон, удовлетворяющих COLA (таким образом, гарантирующий совершенную реконструкцию).

iscola Функция

Проверяйте если sqrt(hann()) окном длины 120 выборок и продолжительность перекрытия 60 выборок является совместимая ВОЛЯ. Установите method аргумент iscola функционируйте к 'wola'. Выход iscola функция является 1 указанием, что этим окном является совместимая ВОЛЯ.

winWOLA = sqrt(hann(winLen,'periodic'));
tfWOLA = iscola(winWOLA,overlapLen,'wola')

tfWOLA = logical
   1

Восстановите данные с ВОЛЕЙ

Выпустите dsp.STFT и dsp.ISTFT Системные объекты и набор окно к sqrt(hann(winLen,'periodic')) окно. Чтобы использовать взвешенное перекрытие - прибавляют ISTFT сторона, набор 'WeightedOverlapAdd' к true.

release(stf);
release(istf);
stf.Window = winWOLA;
istf.Window = winWOLA;
istf.WeightedOverlapAdd = true;

n = zeros(1,100);
xprev = 0;
for i = 1:100
    x = randn(frameLen,1);
    X = stf(x);
    y = istf(X);
    n(1,i) = norm(y-xprev);
    xprev = x;
end       
max(abs(n))

ans = 4.6664e-15

norm из различия между входным сигналом и восстановленным сигналом очень мал, указывая, что сигнал был восстановлен отлично.

Больше о

развернуть все

Кратковременное преобразование Фурье

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

Сигнал области дискретного времени, который будет преобразован, разбит в короткие сегменты. Эти сегменты обычно перекрывают друг друга так, чтобы артефакты на контурах уменьшались. Каждым сегментом является Фурье, преобразованный, и комплексный результат добавляется к матрице, которая записывает величину и фазу для каждого момента времени и частоты.

STFT дают

$X_{m} (ω) = \sum_{n = - \infty}^{\infty} x (n) g (n - m R) e^{J ω n}$

где,

x (n) – Входной сигнал во время n.
Функция окна g (n) - Length M.
X _m (ω) – DTFT оконных данных, сосредоточенных во время mR.
R Скачкообразно переместите размер, в выборках, между последовательным DTFTs.

Если окно g (n) имеет постоянное перекрытие - добавляет свойство (COLA) в размере транзитного участка R, то есть, если

$\sum_{m = - \infty}^{\infty} g (n - m R) = 1, \forall n \in Z (g \in C O L A (R)),$

затем сумма последовательного DTFTs в зависимости от времени равняется DTFT целого X сигнала (ω):

$\begin{array}{l} \sum_{m = - \infty}^{\infty} X_{m} (ω) ≜ \sum_{m = - \infty}^{\infty} \sum_{n = - \infty}^{\infty} x (n) g (n - m R) e^{J ω n}, \\ = \sum_{n = - \infty}^{\infty} x (n) e^{J ω n} \cdot \sum_{m = - \infty}^{\infty} g (n - m R), \\ = \sum_{n = - \infty}^{\infty} x (n) e^{J ω n} \cdot 1, \\ ≜ {DTFT}_{ω} (x) = X (ω) . \end{array}$

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

Величина, в квадрате из STFT, уступает spectrogram представление Степени Спектральная Плотность функции.

Алгоритмы

Вот эскиз того, как алгоритм реализован:

Входной сигнал временного интервала буферизуется на основе заданной пользователями длины окна (WL) и длина перекрытия (OL). Размер транзитного участка, R, задан как R = WL – OL. Буферизированные окна умножаются на заданное пользователями окно длины WL. STFT выход является БПФ этого продукта. Количеством выборок временного интервала, требуемых сформировать новый FFT выход, является R.

Вот рисунок того, как случайный сигнал похож в исходном временном интервале после умножения с накладывающимися окнами, и после применения БПФ на умноженных окнах:

Ссылки

[1] Аллен, J.B., и Л. Р. Рэбинер. "Объединенный Подход к Кратковременному анализу Фурье и Синтезу'', Продолжения IEEE, Издания 65, стр 1558–1564, ноябрь 1977.

Документация

dsp.STFT

Описание

Создание

Синтаксис

Описание

Свойства

`Window` — Окно Analysis
`sqrt(hann(512,'periodic'))` (значение по умолчанию) | вектор

`OverlapLength` — Перекройте длину
256 (значение по умолчанию) | положительное целое число

`FFTLength` — Длина БПФ
512 (значение по умолчанию) | положительное целое число

`FrequencyRange` — Частотный диапазон
`'twosided'` (значение по умолчанию) | `'onesided'`

Использование

Синтаксис

Описание

Входные параметры

`x` — Входной сигнал
вектор | матрица

Выходные аргументы

`y` — STFT выход
вектор | матрица

Функции объекта

Примеры

Кратковременное спектральное затухание

Совершенная реконструкция

Больше о

Кратковременное преобразование Фурье

Алгоритмы

Ссылки

Расширенные возможности

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

Объекты

Блоки

Документация DSP System Toolbox

Поддержка

Документация

dsp.STFT

Описание

Создание

Синтаксис

Описание

Свойства

Window — Окно Analysis sqrt(hann(512,'periodic')) (значение по умолчанию) | вектор

OverlapLength — Перекройте длину256 (значение по умолчанию) | положительное целое число

FFTLength — Длина БПФ512 (значение по умолчанию) | положительное целое число

FrequencyRange — Частотный диапазон 'twosided' (значение по умолчанию) | 'onesided'

Использование

Синтаксис

Описание

Входные параметры

x — Входной сигнал вектор | матрица

Выходные аргументы

y — STFT выход вектор | матрица

Функции объекта

Примеры

Кратковременное спектральное затухание

Совершенная реконструкция

Больше о

Кратковременное преобразование Фурье

Алгоритмы

Ссылки

Расширенные возможности

Генерация кода C/C++ Генерация кода C и C++ с помощью MATLAB® Coder™.

Смотрите также

Объекты

Блоки

Документация DSP System Toolbox

Поддержка

`Window` — Окно Analysis
`sqrt(hann(512,'periodic'))` (значение по умолчанию) | вектор

`OverlapLength` — Перекройте длину
256 (значение по умолчанию) | положительное целое число

`FFTLength` — Длина БПФ
512 (значение по умолчанию) | положительное целое число

`FrequencyRange` — Частотный диапазон
`'twosided'` (значение по умолчанию) | `'onesided'`

`x` — Входной сигнал
вектор | матрица

`y` — STFT выход
вектор | матрица

Генерация кода C/C++
Генерация кода C и C++ с помощью MATLAB® Coder™.