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 по данному диапазону частот вычисляет среднюю степень в сигнале по тому диапазону частот. В отличие от среднеквадратического спектра, peaks в этом спектры не отражают степень на данной частоте. Смотрите avgpower метод dspdata для получения дополнительной информации.

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

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, вектор по умолчанию создается. Если односторонний выбран, то целое число точек БПФ (nFFT) для этого вектора принято, чтобы быть четным.

Если onesided выбран и вы задаете Frequencies, последняя точка частоты сравнивается с предпоследней точкой и к π (или Fs/2, если Fs задан). Если последняя точка ближе к π (или Fs/2) чем он к предыдущей точке, nFFT принят, чтобы быть четным. Если это ближе к предыдущей точке, nFFT принят, чтобы быть нечетным.

Длина Frequencies вектор должен совпадать с длиной столбцов Data.

Fs

'Normalized'

Частота дискретизации, которая является 'Normalized' если NormalizedFrequency true. Если NormalizedFrequency false Fs значения по умолчанию к 1.

SpectrumType

'Onesided'

Интервал Найквиста, на котором вычисляется спектральная плотность мощности. Допустимыми значениями является 'Onesided' и 'Twosided'. Односторонний PSD содержит общую степень сигнала в половине интервала Найквиста. Смотрите onesided и twosided методы в dspdata для получения информации об изменении этого свойства.

Область значений для половины интервала Найквиста является [0 пи), или [0 пи] в зависимости от количества точек БПФ. Для целого интервала Найквиста область значений [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, это указывает, что компонент DC находится в центре двухстороннего спектра. Установите flag к false если компонент DC находится на левом крае спектра.

Методы

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

  • avgpower

  • centerdc

  • normalizefreq

  • onesided

  • plot

  • twosided

Например, чтобы нормировать частоту и установить NormalizedFrequency параметр к истинному, использовать

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 object. The axes object 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 object. The axes object with title Power Spectral Density contains an object of type line.

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

| | | | | |

Представлено до R2006a
Для просмотра документации необходимо авторизоваться на сайте