instbw

Оцените мгновенную пропускную способность

свернуть все на странице

Синтаксис

ibw = instbw(x,fs)

ibw = instbw(x,t)

ibw = instbw(xt)

ibw = instbw(tfd,fd,td)

ibw = instbw(___,Name,Value)

[ibw,t] = instbw(___)

instbw(___)

Описание

пример

ibw = instbw(x,fs) оценивает мгновенную пропускную способность сигнала, x, произведенный на уровне fs. Если x матрица, затем функция оценивает мгновенную пропускную способность независимо для каждого столбца и возвращает результат в соответствующий столбец ibw.

пример

ibw = instbw(x,t) оценивает мгновенную пропускную способность x произведенный в то время, когда значения сохранены в t.

пример

ibw = instbw(xt) оценивает мгновенную пропускную способность сигнала, сохраненного в расписании MATLAB^® xt. Функция обрабатывает все переменные в расписании и все столбцы в каждой переменной независимо.

пример

ibw = instbw(tfd,fd,td) оценивает мгновенную пропускную способность сигнала чья плотность распределения времени, tfd, производится в значениях пропускной способности, сохраненных в fd и временные стоимости сохранены в td.

ibw = instbw(___,Name,Value) задает дополнительные опции для любого из предыдущих синтаксисов с помощью аргументов значения имени. Можно задать масштабный коэффициент или пределы частоты, используемые в расчете. Например, 'FrequencyLimits',[10 20] вычисляет мгновенную пропускную способность входа в диапазоне от 10 Гц до 20 Гц.

пример

[ibw,t] = instbw(___) также возвращает t, вектор из шагов расчета, соответствующих ibw.

пример

instbw(___) без выходных аргументов строит предполагаемую мгновенную пропускную способность.

Примеры

свернуть все

Мгновенная пропускная способность синусоидального щебета

Скрипт Open Live Script

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

fs = 6e2;

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

Вычислите спектрограмму сигнала и отобразите его как график водопада.

[p,f,t] = pspectrum(x,fs,'spectrogram');

waterfall(f,t,p')
ax = gca;
ax.XDir = 'reverse';
view(30,45)

Figure contains an axes. The axes contains an object of type patch.

Оцените и постройте мгновенную пропускную способность сигнала.

instbw(x,fs)
ylim([0 50])

Figure contains an axes. The axes with title Instantaneous Bandwidth Estimate contains an object of type line.

Мгновенная пропускная способность от спектрограммы

Скрипт Open Live Script

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

fs = 2000;
t = 0:1/fs:2-1/fs;

frq = [50:100:950]';

amp = (t > 4*(frq-frq(1))/fs);
x = sum(amp.*sin(2*pi*t.*frq).*exp(-3*t));

% To hear, type sound(x,fs)

plot(t,x)
xlabel('Time (s)')

Figure contains an axes. The axes contains an object of type line.

Вычислите и отобразите мгновенную пропускную способность сигнала.

[bw,bt] = instbw(x,t);

plot(bt,bw)
xlabel('Time (s)')
ylabel('Bandwidth (Hz)')

Figure contains an axes. The axes contains an object of type line.

Вычислите спектрограмму сигнала. Задайте разрешение времени 100 миллисекунд и 0 перекрытий между смежными сегментами. Используйте спектрограмму, чтобы оценить мгновенную частоту сигнала.

[p,ff,tt] = pspectrum(x,t,'spectrogram','TimeResolution',0.1,'OverlapPercent',0);

instfreq(p,ff,tt)

Figure contains an axes. The axes with title Instantaneous Frequency Estimate contains 2 objects of type image, line. This object represents Instantaneous Frequency.

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

instbw(p,ff,tt)

Figure contains an axes. The axes with title Instantaneous Bandwidth Estimate contains an object of type line.

Мгновенная пропускная способность расписания

Скрипт Open Live Script

Сгенерируйте сигнал, произведенный на уровне 14 кГц в течение 2 секунд. Частота сигнала варьируется как щебет, модулируемый Гауссовым. Сохраните сигнал как расписание MATLAB®.

fs = 14000;
t = (0:1/fs:2)';
s = vco(chirp(t+.1,0,t(end),3).*exp(-2*(t-1).^2),[0.1 0.4]*fs,fs);

sx = timetable(s,'SampleRate',fs);

Вычислите спектрограмму сигнала. Задайте утечку 0,2, разрешение времени 50 миллисекунд и 99% перекрытия между смежными сегментами. Отобразите спектрограмму.

opts = {'spectrogram','Leakage',0.2,'TimeResolution',0.05,'OverlapPercent',99};

[p,ff,tt] = pspectrum(sx,opts{:});

pspectrum(sx,opts{:})

Figure contains an axes. The axes with title Fres = 64.5333 Hz, Tres = 50 ms contains an object of type image.

Оцените и отобразите мгновенную пропускную способность сигнала.

instbw(p,ff,tt)

Figure contains an axes. The axes with title Instantaneous Bandwidth Estimate contains an object of type line.

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

свернуть все

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

Входной сигнал в виде вектора или матрицы. Если x вектор, затем instbw обработки это как один канал. Если x матрица, затем instbw вычисляет мгновенную пропускную способность независимо для каждого столбца и возвращает результат в соответствующий столбец ibw.

Пример: sin(2*pi*(0:127)/16)+randn(1,128)/100 задает шумную синусоиду.

Пример: [2 1].*sin(2*pi*(0:127)'./[16 64]) задает двухканальную синусоиду.

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

`fs` — Частота дискретизации
положительная скалярная величина

Частота дискретизации в виде положительной скалярной величины. Частота дискретизации является количеством выборок в единицу времени. Если модуль времени является секундами, то частота дискретизации находится в Гц.

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

`t` — Шаги расчета
вектор действительных чисел | `duration` скаляр | `duration` массив | `datetime` массив

Шаги расчета в виде вектора действительных чисел, a duration скаляр, a duration массив или a datetime массив.

duration скаляр — временной интервал между последовательными выборками x.
Вектор действительных чисел, duration массив, или datetime массив — момент времени, соответствующий каждому элементу x.

Пример: seconds(1) задает 1 вторую ошибку между последовательными измерениями сигнала.

Пример: seconds(0:8) указывает, что сигнал производится на уровне 1 Гц в течение 8 секунд.

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

`xt` — Введите расписание
расписание

Введите расписание. xt должен содержать увеличение, конечные времена строки.

Если расписание имеет пропавших без вести или дублирующиеся моменты времени, можно зафиксировать его с помощью советов в Чистом Расписании с Пропавшими без вести, Копией, или Неоднородные Времена.

Пример: timetable(seconds(0:4)',randn(5,1)) задает вероятностный процесс, произведенный на уровне 1 Гц в течение 4 секунд.

Пример: timetable(seconds(0:4)',randn(5,3),randn(5,4)) содержит вероятностный процесс с тремя каналами и вероятностный процесс с четырьмя каналами, оба произведенные на уровне 1 Гц в течение 4 секунд.

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

`tfd` — Плотность распределения времени
матрица

Плотность распределения времени в виде матрицы, произведенной на частотах, сохранена в fd и временные стоимости сохранены в td. Этот входной параметр поддерживается только когда 'Method' установлен в 'tfmoment'.

Пример: [p,f,t] = pspectrum(sin(2*pi*(0:511)/4),4,'spectrogram') задает плотность распределения времени синусоиды на 1 Гц, произведенной на уровне 4 Гц в течение 128 секунд, и также частот и времена, в которые она вычисляется.

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

`fd`, `td` — Частота и временные стоимости для плотности распределения времени
векторы

Частота и временные стоимости для плотности распределения времени в виде векторов. Эти входные параметры поддерживаются только когда 'Method' установлен в 'tfmoment'.

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

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'FrequencyLimits',[25 50] вычисляет мгновенную пропускную способность входа в диапазоне от 25 Гц до 50 Гц.

`'FrequencyLimits'` — Частотный диапазон
`[0 fs/2]` (значение по умолчанию для сигналов с действительным знаком) | `[-fs/2 fs/2]` (значение по умолчанию для сигналов с комплексным знаком) | двухэлементный вектор в Гц

Частотный диапазон в виде двухэлементного вектора в Гц. Если не заданный, 'FrequencyLimits' значения по умолчанию к [0 fs/2] для сигналов с действительным знаком и к [-fs/2 fs/2] для сигналов с комплексным знаком.

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

`'ScaleFactor'` — Масштабный коэффициент в течение спектрального момента
`sqrt(4*pi)` (значение по умолчанию) | действительный скаляр

Масштабный коэффициент в течение спектрального момента в виде действительного скаляра.

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

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

свернуть все

`ibw` — Мгновенная пропускная способность
вектор | матрица | расписание

Мгновенная пропускная способность, возвращенная как вектор, матрица или расписание с теми же размерностями как вход.

`t` — Времена оценок пропускной способности
вектор действительных чисел | `duration` массив | `datetime` массив

Времена оценок пропускной способности, возвращенных как вектор действительных чисел, a duration массив или a datetime массив.

Больше о

свернуть все

Мгновенная пропускная способность

instantaneous bandwidth неустановившегося сигнала является изменяющимся во времени параметром, который относится к распространению мгновенной частоты о ее среднем значении в установленный срок момент [1], [2].

instbw оценивает мгновенную пропускную способность как квадратный корень второго условного спектрального момента плотности распределения времени входного сигнала. Функция:

Вычисляет спектр мощности спектрограммы P (t, f) входа с помощью pspectrum функция и использование спектр как плотность распределения времени.
Оценивает мгновенное использование пропускной способности

$σ_{f}^{2} (t) = \frac{\int_{0}^{\infty} {(f - f_{inst} (t))}^{2} P (t, f) d f}{\int_{0}^{\infty} P (t, f) d f},$
где $f_{inst} (t)$ мгновенная частота, возвращенная instfreq и оцененный

$f_{inst} (t) = \frac{\int_{0}^{\infty} f P (t, f) d f}{\int_{0}^{\infty} P (t, f) d f} .$

Ссылки

[1] Boashash, B. “Оценивая и Интерпретируя Мгновенную Частоту Сигнала. I. Основные принципы”. Продолжения IEEE^® 80, № 4 (апрель 1992): 520–38. https://doi.org/10.1109/5.135376.

[2] Boashash, B. "Оценивая и Интерпретируя Мгновенную Частоту Сигнала. II. Алгоритмы и Приложения". Продолжения IEEE 80, № 4 (май 1992): 540–568. https://doi.org/10.1109/5.135378.

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

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

Указания и ограничения по применению:

Расписания не поддерживаются для генерации кода.

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

instfreq | pspectrum

Введенный в R2021a

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

instbw

Синтаксис

Описание

Примеры

Мгновенная пропускная способность синусоидального щебета

Мгновенная пропускная способность от спектрограммы

Мгновенная пропускная способность расписания

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

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

`fs` — Частота дискретизации
положительная скалярная величина

`t` — Шаги расчета
вектор действительных чисел | `duration` скаляр | `duration` массив | `datetime` массив

`xt` — Введите расписание
расписание

`tfd` — Плотность распределения времени
матрица

`fd`, `td` — Частота и временные стоимости для плотности распределения времени
векторы

Аргументы в виде пар имя-значение

`'ScaleFactor'` — Масштабный коэффициент в течение спектрального момента
`sqrt(4*pi)` (значение по умолчанию) | действительный скаляр

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

`ibw` — Мгновенная пропускная способность
вектор | матрица | расписание

`t` — Времена оценок пропускной способности
вектор действительных чисел | `duration` массив | `datetime` массив

Больше о

Мгновенная пропускная способность

Ссылки

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

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

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

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

Поддержка

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

instbw

Синтаксис

Описание

Примеры

Мгновенная пропускная способность синусоидального щебета

Мгновенная пропускная способность от спектрограммы

Мгновенная пропускная способность расписания

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

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

fs — Частота дискретизации положительная скалярная величина

t — Шаги расчета вектор действительных чисел | duration скаляр | duration массив | datetime массив

xt — Введите расписание расписание

tfd — Плотность распределения времени матрица

fd, td — Частота и временные стоимости для плотности распределения времени векторы

Аргументы в виде пар имя-значение

'ScaleFactor' — Масштабный коэффициент в течение спектрального момента sqrt(4*pi) (значение по умолчанию) | действительный скаляр

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

ibw — Мгновенная пропускная способность вектор | матрица | расписание

t — Времена оценок пропускной способности вектор действительных чисел | duration массив | datetime массив

Больше о

Мгновенная пропускная способность

Ссылки

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

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

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

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

Поддержка

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

`fs` — Частота дискретизации
положительная скалярная величина

`t` — Шаги расчета
вектор действительных чисел | `duration` скаляр | `duration` массив | `datetime` массив

`xt` — Введите расписание
расписание

`tfd` — Плотность распределения времени
матрица

`fd`, `td` — Частота и временные стоимости для плотности распределения времени
векторы

`'ScaleFactor'` — Масштабный коэффициент в течение спектрального момента
`sqrt(4*pi)` (значение по умолчанию) | действительный скаляр

`ibw` — Мгновенная пропускная способность
вектор | матрица | расписание

`t` — Времена оценок пропускной способности
вектор действительных чисел | `duration` массив | `datetime` массив

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