Exponenta Banner

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 вектор, затем 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
Поддержка комплексного числа: Да

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

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

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

  • duration скаляр — временной интервал между последовательными выборками x.

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

свернуть все

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

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

Больше о

свернуть все

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

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

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

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

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

    σf2(t)=∫0∞(f−finst(t))2 P(t,f) df∫0∞P(t,f) df,

    где finst(t) мгновенная частота, возвращенная instfreq и оцененный

    finst(t)=∫0∞f P(t,f) df∫0∞P(t,f) df.

Ссылки

[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.

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

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

|

Введенный в R2021a

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

Поддержка