dspdata.psd

Спектральная плотность мощности

Синтаксис

Hpsd = dspdata.psd(Data) Hpsd = dspdata.psd(Data,Frequencies) Hpsd = dspdata.psd(...,'Fs',Fs) Hpsd = dspdata.psd(...,'SpectrumType',SpectrumType) Hpsd = dspdata.psd(...,'CenterDC',flag)

Описание

Примечание

Использование dspdata.psd не рекомендуется. Использовать pburg, pcov, periodogram, pmcov, pmtm, pwelch, или pyulear вместо этого.

Спектральная плотность мощности (PSD) предназначена для непрерывных спектров. Интеграл PSD в данной полосе частот вычисляет среднюю мощность в сигнале в этой полосе частот. В отличие от среднеквадратичного спектра, пики в этом спектре не отражают мощность на данной частоте. См. раздел avgpower способ dspdata для получения дополнительной информации.

Односторонняя PSD содержит суммарную мощность сигнала в частотном интервале от постоянного тока до половины скорости Найквиста. Двусторонняя PSD содержит общую мощность в интервале частот от постоянного тока до скорости Найквиста.

Hpsd = dspdata.psd(Data) использует данные спектральной плотности мощности, содержащиеся в Data, которая может быть в виде вектора или матрицы, где каждый столбец является отдельным набором данных. Значения по умолчанию для других свойств объекта показаны ниже:

Собственность	Значение по умолчанию	Описание
`Name`	`'Power Spectral Density'`	Вектор символов только для чтения
`Frequencies`	`[]` напечатать `double`	Вектор частот, на которых оценивается спектральная плотность мощности. Диапазон этого вектора зависит от `SpectrumType` значение. Для одностороннего диапазона по умолчанию - [0, δ) или [0, `Fs`/ 2) для нечетной длины и [0, δ] или [0 ,`Fs`/ 2] для четной длины, если`Fs` указывается. Для двухсторонних это значение равно [0, 2pi) или [0, `Fs`). Если не указать `Frequencies`создается вектор по умолчанию. Если выбран односторонний вариант, то целое число точек БПФ (nБПФ) для этого вектора предполагается четным. Если `onesided` выбрано и указано `Frequencies`, последняя частотная точка сравнивается со следующей-последней точкой и с δ (или `Fs`/ 2, если`Fs` указывается). Если последняя точка ближе к δ (или `Fs`/ 2), чем до предыдущей точки, nБПФ считается четным. Если он ближе к предыдущей точке, то nБПФ считается нечетным. Длина `Frequencies` вектор должен соответствовать длине столбцов `Data`.
`Fs`	`'Normalized'`	Частота выборки, которая `'Normalized'` если `NormalizedFrequency` является `true`. Если `NormalizedFrequency` является `false` `Fs` по умолчанию: `1`.
`SpectrumType`	`'Onesided'`	Интервал Найквиста, по которому вычисляется спектральная плотность мощности. Допустимые значения: `'Onesided'` и `'Twosided'`. Односторонняя PSD содержит суммарную мощность сигнала в половине интервала Найквиста. См. раздел `onesided` и `twosided` методы в `dspdata` для получения сведений об изменении этого свойства. Диапазон для половины интервала Найквиста составляет [0 pi) или [0 pi] в зависимости от количества точек БПФ. Для всего интервала Найквиста диапазон равен [0 2pi).
`NormalizedFrequency`	`true`	Нормализована ли частота (`true`) или нет (`false`). Это свойство устанавливается автоматически во время конструирования на основе `Fs`. Если `Fs` указывается, `NormalizedFrequency` имеет значение `false`. См. раздел `normalizefreq` метод в `dspdata` для получения сведений об изменении этого свойства.

Hpsd = dspdata.psd(Data,Frequencies) использует данные оценки спектральной плотности мощности, содержащиеся в Data и Frequencies векторы.

Hpsd = dspdata.psd(...,'Fs',Fs) использует частоту выборки Fs. Определение Fs использует набор линейных частот по умолчанию (в Hz) на основе Fs и наборы NormalizedFrequency кому false.

Hpsd = dspdata.psd(...,'SpectrumType',SpectrumType) определяет интервал, в течение которого вычисляется спектральная плотность мощности. Для данных, которые находятся в диапазоне от [0 δ) или [0 δ], установите значение SpectrumType кому onesided; для данных, которые находятся в диапазоне от [0 2δ), установите значение SpectrumType кому twosided.

Hpsd = dspdata.psd(...,'CenterDC',flag) использует значение flag для указания, является ли компонент с нулевой частотой (DC) центрированным. Если flag является true, это указывает на то, что компонент постоянного тока находится в центре двустороннего спектра. Установите flag кому false если компонент постоянного тока находится на левом краю спектра.

Методы

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

avgpower
centerdc
normalizefreq
onesided
plot
twosided

Например, для нормализации частоты и установки NormalizedFrequency значение true, используйте

Hpsd = normalizefreq(Hpsd)

Для получения подробной информации об использовании методов и печати спектра см. dspdata справочная страница.

Примеры

свернуть все

Разрешить компоненты сигнала

Открыть сценарий в реальном времени

Оценивают одностороннюю спектральную плотность мощности шумного синусоидального сигнала с двумя частотными составляющими.

Fs = 32e3;   
t = 0:1/Fs:2.96;
x = cos(2*pi*t*1.24e3)+ cos(2*pi*t*10e3)+ randn(size(t));
nfft = 2^nextpow2(length(x));
Pxx = abs(fft(x,nfft)).^2/length(x)/Fs;

Сохраните спектр в объекте данных PSD и постройте график результата.

Hpsd = dspdata.psd(Pxx(1:length(Pxx)/2),'Fs',Fs);  
plot(Hpsd)

Figure contains an axes. The axes with title Power Spectral Density contains an object of type line.

Создайте двусторонний спектр и постройте его график.

Hpsd = dspdata.psd(Pxx,'Fs',Fs,'SpectrumType','twosided');
plot(Hpsd)

Figure contains an axes. The axes with title Power Spectral Density contains an object of type line.

См. также

Представлен до R2006a

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