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

Эта статья посвящена обработке сигналов и связи спектров с временными рядами. Для дальнейших применений в физических науках см. Spectrum (физические науки) .

«Спектральная плотность мощности» перенаправляется сюда. Не путать со Спектральной силой .

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

Форма голосового сигнала с течением времени (слева) имеет широкий спектр мощности звука (справа).

При обработке сигналов спектр мощности $S_{xx}(f)$ непрерывного времени сигнала $x(t)$ описывает распределение мощности на частотные составляющие $f$ составляя этот сигнал. Согласно анализу Фурье , любой физический сигнал можно разложить на ряд дискретных частот или спектр частот в непрерывном диапазоне. Среднее статистическое значение любого типа сигнала (включая шум ), анализируемое с точки зрения его частотного содержания, называется его спектром .

Когда энергия сигнала сосредоточена вокруг конечного интервала времени, особенно если его полная энергия конечна, можно вычислить спектральную плотность энергии . Чаще используется спектральная плотность мощности (PSD или просто спектр мощности ), которая применяется к сигналам, существующим в течение всего времени или в течение достаточно большого периода времени (особенно по отношению к продолжительности измерения), который также можно было бы иметь. прошло в течение бесконечного промежутка времени. Тогда PSD относится к спектральному распределению энергии, которое можно найти в единицу времени, поскольку полная энергия такого сигнала за все время обычно будет бесконечной. Суммирование или интегрирование спектральных компонентов дает полную мощность (для физического процесса) или дисперсию (в статистическом процессе), идентичную тому, что было бы получено путем интегрирования $x^{2}(t)$ во временной области, как это диктуется теоремой Парсеваля .

Спектр физического процесса $x(t)$ часто содержит важную информацию о природе $x$. Например, высота и тембр музыкального инструмента сразу определяются на основе спектрального анализа. Цвет . источника света определяется спектром электрического поля электромагнитной волны $E(t)$ поскольку он колеблется с чрезвычайно высокой частотой. Получение спектра из таких временных рядов включает преобразование Фурье и обобщения, основанные на анализе Фурье. Во многих случаях временная область на практике специально не используется, например, когда дисперсионная призма используется для получения спектра света в спектрографе или когда звук воспринимается посредством его воздействия на слуховые рецепторы внутреннего уха, каждый из которых из которых чувствителен к определенной частоте.

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

Единицы

См. Также: преобразование Фурье § Единицы измерения.

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

Например , когда сигнал определяется только напряжением , не существует уникальной мощности, связанной с указанной амплитудой. В этом случае «мощность» просто рассчитывается в квадрате сигнала, поскольку она всегда будет пропорциональна фактической мощности, передаваемой этим сигналом при заданном импедансе . Таким образом, можно использовать единицы V ² Гц ⁻¹ для PSD. Спектральная плотность энергии (ESD) будет иметь единицы В. ² с Гц ⁻¹ , поскольку энергия имеет единицы мощности, умноженные на время (например, ватт-час ).

В общем случае единицами PSD будут соотношение единиц дисперсии на единицу частоты; так, например, ряд значений смещения (в метрах) во времени (в секундах) будет иметь PSD в единицах квадратных метров на герц, м. ² /Гц. При анализе случайных колебаний единицы измерения г ² Гц ⁻¹ часто используются для PSD ускорения , где g обозначает перегрузку .

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

Односторонний против двустороннего

PSD может быть либо односторонней функцией только положительных частот, либо двусторонней функцией как положительных, так и отрицательных частот , но только с половиной амплитуды. PSD шума обычно односторонние в технике и двусторонние в физике.

Определение

Спектральная плотность энергии

«Спектральная плотность энергии» перенаправляется сюда. Не путать с энергетическим спектром .

Спектральная плотность энергии описывает, как энергия сигнала или временного ряда распределяется по частоте. Здесь термин энергия используется в обобщенном смысле обработки сигналов; то есть энергия $E$ сигнала $x(t)$ является: $$E\triangleq \int _{-\infty }^{\infty }\left|x(t)\right|^{2}\ dt.$$

Спектральная плотность энергии наиболее подходит для переходных процессов, то есть импульсных сигналов, имеющих конечную общую энергию. Конечная или нет, теорема Парсеваля (или теорема Планшереля) дает нам альтернативное выражение для энергии сигнала: $$\int _{-\infty }^{\infty }|x(t)|^{2}\,dt=\int _{-\infty }^{\infty }\left|{\hat {x}}(f)\right|^{2}\,df,$$ где: $${\hat {x}}(f)\triangleq \int _{-\infty }^{\infty }e^{-i2\pi ft}x(t)\ dt$$ значение Фурье преобразования $x(t)$ на частоте $f$ (в Гц ). Теорема справедлива и в случаях дискретного времени. Поскольку интеграл в левой части представляет собой энергию сигнала, значение $\left|{\hat {x}}(f)\right|^{2}df$ можно интерпретировать как функцию плотности, умноженную на бесконечно малый частотный интервал, описывающую энергию, содержащуюся в сигнале на частоте $f$ в частотном интервале $f+df$.

Следовательно, энергии спектральная плотность $x(t)$ определяется как:

${\bar {S}}_{xx}(f)\triangleq \left|{\hat {x}}(f)\right|^{2}$ ( Уравнение 1 )

Функция ${\bar {S}}_{xx}(f)$ и автокорреляция ​ $x(t)$ образуют пару преобразований Фурье, результат, также известный как теорема Винера – Хинчина (см. также Периодограмму ).

В качестве физического примера того, как можно измерить спектральную плотность энергии сигнала, предположим, $V(t)$ представляет собой потенциал (в вольтах ) электрического импульса, распространяющегося вдоль линии передачи с импедансом $Z$и предположим, что линия завершена согласованным резистором (так что вся энергия импульса доставляется на резистор и ни одна не отражается обратно). По закону Ома мощность, передаваемая на резистор за раз $t$ равно $V(t)^{2}/Z$, поэтому полная энергия находится путем интегрирования $V(t)^{2}/Z$ по времени по длительности импульса. Чтобы найти значение спектральной плотности энергии ${\bar {S}}_{xx}(f)$ на частоте $f$, можно было бы вставить между линией передачи и резистором полосовой фильтр , пропускающий только узкий диапазон частот ( $\Delta f$, скажем) вблизи интересующей частоты, а затем измерить полную энергию $E(f)$ рассеивается на резисторе. Значение спектральной плотности энергии при $f$ тогда оценивается как $E(f)/\Delta f$. В этом примере, поскольку мощность $V(t)^{2}/Z$ имеет единицы V ² Ой ⁻¹ , энергия $E(f)$ имеет единицы V ² ох ⁻¹ = J , а значит, и оценка $E(f)/\Delta f$ спектральной плотности энергии имеет единицы Дж Гц ⁻¹ , как требуется. Во многих ситуациях принято забывать этап деления на $Z$ так что спектральная плотность энергии вместо этого имеет единицы V ² Гц ⁻¹ .

Это определение напрямую обобщается на дискретный сигнал со счетным бесконечным числом значений. $x_{n}$ например, сигнал, дискретизированный в дискретные моменты времени $t_{n}=t_{0}+(n\,\Delta t)$: $${\bar {S}}_{xx}(f)=\lim _{N\to \infty }(\Delta t)^{2}\underbrace {\left|\sum _{n=-N}^{N}x_{n}e^{-i2\pi fn\,\Delta t}\right|^{2}} _{\left|{\hat {x}}_{d}(f)\right|^{2}},$$ где ${\hat {x}}_{d}(f)$ — Фурье с дискретным временем преобразование $x_{n}.$ Интервал выборки $\Delta t$ необходимо, чтобы сохранить правильные физические единицы и гарантировать, что мы восстановим непрерывный случай в пределе $\Delta t\to 0.$ Но в математических науках интервал часто устанавливается равным 1, что упрощает результаты в ущерб общности. (также см. нормализованную частоту )

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

Не путать со спектральным распределением мощности .

Спектр мощности измеренной анизотропии температуры космического микроволнового фонового излучения в угловом масштабе. Сплошная линия — теоретическая модель для сравнения.

Приведенное выше определение спектральной плотности энергии подходит для переходных процессов (импульсных сигналов), энергия которых сосредоточена вокруг одного временного окна; тогда преобразования Фурье сигналов обычно существуют. Для непрерывных сигналов в течение всего времени необходимо скорее определить спектральную плотность мощности (PSD), которая существует для стационарных процессов ; это описывает, как мощность сигнала или временного ряда распределяется по частоте, как в простом примере, приведенном ранее. Здесь мощность может быть фактической физической мощностью или, что чаще, для удобства абстрактных сигналов, просто идентифицируется как квадрат значения сигнала. Например, статистики изучают изменение функции во времени. $x(t)$ (или по другой независимой переменной), и, используя аналогию с электрическими сигналами (среди других физических процессов), его принято называть спектром мощности, даже если здесь не задействована физическая мощность. Если бы нужно было создать физический источник напряжения , который последовал бы $x(t)$ и применил ее к выводам сопротивлением 1 Ом резистора , тогда действительно мгновенная мощность, рассеиваемая на этом резисторе, будет равна $x^{2}(t)$ ватты .

Средняя мощность $P$ сигнала $x(t)$ Таким образом, за все время определяется следующим средним временем, где период $T$ сосредоточено относительно некоторого произвольного времени $t=t_{0}$: $$P=\lim _{T\to \infty }{\frac {1}{T}}\int _{t_{0}-T/2}^{t_{0}+T/2}\left|x(t)\right|^{2}\,dt$$

Однако для понимания последующей математики удобнее иметь дело с ограничениями по времени в самом сигнале, а не с ограничениями по времени в границах интеграла. Таким образом, у нас есть альтернативное представление средней мощности, где $x_{T}(t)=x(t)w_{T}(t)$ и $w_{T}(t)$ есть единица в пределах произвольного периода и ноль в другом месте. $$P=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }\left|x_{T}(t)\right|^{2}\,dt.$$ Очевидно, что в случаях, когда приведенное выше выражение для P не равно нулю, интеграл должен неограниченно расти, поскольку T неограниченно растет. Именно по этой причине в таких случаях мы не можем использовать энергию сигнала, которая представляет собой расходящийся интеграл.

При анализе частотного содержания сигнала $x(t)$, можно было бы вычислить обычное преобразование Фурье ${\hat {x}}(f)$; однако для многих представляющих интерес сигналов преобразование Фурье формально не существует. Тем не менее, теорема Парсеваля говорит нам, что мы можем переписать среднюю степень следующим образом. $$P=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }|{\hat {x}}_{T}(f)|^{2}\,df$$

Тогда спектральная плотность мощности просто определяется как подынтегральная функция, указанная выше.

$S_{xx}(f)=\lim _{T\to \infty }{\frac {1}{T}}|{\hat {x}}_{T}(f)|^{2}\,$ ( Уравнение 2 )

Отсюда, благодаря теореме о свертке , мы также можем просмотреть $|{\hat {x}}_{T}(f)|^{2}$ как преобразование временной свертки Фурье $x_{T}^{*}(-t)$ и $x_{T}(t)$, где * представляет собой комплексно-сопряженное соединение. Принимая во внимание, что $${\begin{aligned}{\mathcal {F}}\left\{x_{T}^{*}(-t)\right\}&=\int _{-\infty }^{\infty }x_{T}^{*}(-t)e^{-i2\pi ft}dt\\&=\int _{-\infty }^{\infty }x_{T}^{*}(t)e^{i2\pi ft}dt\\&=\int _{-\infty }^{\infty }x_{T}^{*}(t)[e^{-i2\pi ft}]^{*}dt\\&=\left[\int _{-\infty }^{\infty }x_{T}(t)e^{-i2\pi ft}dt\right]^{*}\\&=\left[{\mathcal {F}}\left\{x_{T}(t)\right\}\right]^{*}\\&=\left[{\hat {x}}_{T}(f)\right]^{*}\end{aligned}}$$ и делая, $u(t)=x_{T}^{*}(-t)$, у нас есть: $${\begin{aligned}\left|{\hat {x}}_{T}(f)\right|^{2}&=[{\hat {x}}_{T}(f)]^{*}\cdot {\hat {x}}_{T}(f)\\&={\mathcal {F}}\left\{x_{T}^{*}(-t)\right\}\cdot {\mathcal {F}}\left\{x_{T}(t)\right\}\\&={\mathcal {F}}\left\{u(t)\right\}\cdot {\mathcal {F}}\left\{x_{T}(t)\right\}\\&={\mathcal {F}}\left\{u(t)\mathbin {\mathbf {*} } x_{T}(t)\right\}\\&=\int _{-\infty }^{\infty }\left[\int _{-\infty }^{\infty }u(\tau -t)x_{T}(t)dt\right]e^{-i2\pi f\tau }d\tau \\&=\int _{-\infty }^{\infty }\left[\int _{-\infty }^{\infty }x_{T}^{*}(t-\tau )x_{T}(t)dt\right]e^{-i2\pi f\tau }\ d\tau ,\end{aligned}}$$ где теорема свертки при переходе от 3-й к 4-й строке использовалась .

Теперь, если мы разделим приведенную выше временную свертку на период $T$ и возьмем предел как $T\rightarrow \infty$, это становится автокорреляционной функцией неоконного сигнала $x(t)$, который обозначается как $R_{xx}(\tau )$, при условии, что $x(t)$ является эргодическим , что верно в большинстве, но не во всех практических случаях. $$\lim _{T\to \infty }{\frac {1}{T}}\left|{\hat {x}}_{T}(f)\right|^{2}=\int _{-\infty }^{\infty }\left[\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }x_{T}^{*}(t-\tau )x_{T}(t)dt\right]e^{-i2\pi f\tau }\ d\tau =\int _{-\infty }^{\infty }R_{xx}(\tau )e^{-i2\pi f\tau }d\tau$$

Отсюда мы видим, снова предполагая эргодичность $x(t)$, что спектральную плотность мощности можно найти как преобразование Фурье автокорреляционной функции ( теорема Винера–Хинчина ).

$S_{xx}(f)=\int _{-\infty }^{\infty }R_{xx}(\tau )e^{-i2\pi f\tau }\,d\tau ={\hat {R}}_{xx}(f)$ ( Уравнение 3 )

Многие авторы используют это равенство для фактического определения спектральной плотности мощности.

Мощность сигнала в заданном диапазоне частот $[f_{1},f_{2}]$, где $0<f_{1}<f_{2}$, можно рассчитать путем интегрирования по частоте. С $S_{xx}(-f)=S_{xx}(f)$, равное количество мощности можно отнести к положительным и отрицательным полосам частот, что составляет коэффициент 2 в следующей форме (такие тривиальные коэффициенты зависят от используемых соглашений): $$P_{\textsf {bandlimited}}=2\int _{f_{1}}^{f_{2}}S_{xx}(f)\,df$$ В более общем смысле, аналогичные методы могут использоваться для оценки изменяющейся во времени спектральной плотности. В этом случае интервал времени $T$ конечен, а не стремится к бесконечности. Это приводит к уменьшению спектрального охвата и разрешения, поскольку частоты менее $1/T$ не дискретизируются, и результаты получаются на частотах, которые не являются целым кратным $1/T$ не являются независимыми. При использовании одного такого временного ряда оцененный спектр мощности будет очень «зашумлен»; однако это можно облегчить, если можно оценить ожидаемое значение (в приведенном выше уравнении), используя большое (или бесконечное) количество краткосрочных спектров, соответствующих статистическим ансамблям реализаций $x(t)$ оценивается в течение указанного временного окна.

Как и в случае со спектральной плотностью энергии, определение спектральной плотности мощности можно обобщить на дискретные переменные времени. $x_{n}$. Как и раньше, мы можем рассмотреть окно $-N\leq n\leq N$ с сигналом, дискретизированным в дискретные моменты времени $t_{n}=t_{0}+(n\,\Delta t)$ за общий период измерения $T=(2N+1)\,\Delta t$. $$S_{xx}(f)=\lim _{N\to \infty }{\frac {(\Delta t)^{2}}{T}}\left|\sum _{n=-N}^{N}x_{n}e^{-i2\pi fn\,\Delta t}\right|^{2}$$ Обратите внимание, что единая оценка PSD может быть получена с помощью конечного числа выборок. Как и раньше, фактическая PSD достигается, когда $N$ (и таким образом $T$) приближается к бесконечности, и формально применяется ожидаемое значение. В реальных приложениях обычно усредняют PSD конечного измерения по многим испытаниям, чтобы получить более точную оценку теоретической PSD физического процесса, лежащего в основе отдельных измерений. Эту вычисленную PSD иногда называют периодограммой . Эта периодограмма сходится к истинной PSD как по количеству оценок, так и по временному интервалу усреднения. $T$ приближаться к бесконечности.

Если оба сигнала обладают спектральной плотностью мощности, то кросс-спектральную плотность аналогичным образом можно рассчитать ; Поскольку PSD связана с автокорреляцией, то же самое относится и кросс-спектральная плотность к взаимной корреляции .

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

Некоторые свойства PSD включают в себя:

• Спектр мощности всегда действителен и неотрицательен, а спектр действительнозначного процесса также является четной функцией частоты: $S_{xx}(-f)=S_{xx}(f)$.
• Для непрерывного случайного процесса x(t) автокорреляционная функция R _xx ( t ) может быть восстановлена ​​по ее спектру мощности S _xx (f) с помощью обратного преобразования Фурье
• Используя теорему Парсеваля , можно вычислить дисперсию (среднюю мощность) процесса путем интегрирования спектра мощности по всей частоте: $$P=\operatorname {Var} (x)=\int _{-\infty }^{\infty }\!S_{xx}(f)\,df$$
• Для реального процесса x ( t ) со спектральной плотностью мощности $S_{xx}(f)$, можно вычислить интегральный спектр или спектральное распределение мощности $F(f)$, который определяет среднюю мощность в ограниченной полосе частот , содержащуюся в частотах от постоянного тока до f, используя: $$F(f)=2\int _{0}^{f}S_{xx}(f')\,df'.$$ Обратите внимание, что предыдущее выражение для полной мощности (дисперсии сигнала) представляет собой особый случай, когда f → ∞ .

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

См. Также: Когерентность (обработка сигналов)

Учитывая два сигнала $x(t)$ и $y(t)$, каждый из которых обладает спектральной плотностью мощности $S_{xx}(f)$ и $S_{yy}(f)$, можно определить перекрестную спектральную плотность мощности ( CPSD ) или перекрестную спектральную плотность ( CSD ). Для начала рассмотрим среднюю мощность такого комбинированного сигнала. $${\begin{aligned}P&=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }\left[x_{T}(t)+y_{T}(t)\right]^{*}\left[x_{T}(t)+y_{T}(t)\right]dt\\&=\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }|x_{T}(t)|^{2}+x_{T}^{*}(t)y_{T}(t)+y_{T}^{*}(t)x_{T}(t)+|y_{T}(t)|^{2}dt\\\end{aligned}}$$

Используя те же обозначения и методы, которые использовались для вывода спектральной плотности мощности, мы используем теорему Парсеваля и получаем $${\begin{aligned}S_{xy}(f)&=\lim _{T\to \infty }{\frac {1}{T}}\left[{\hat {x}}_{T}^{*}(f){\hat {y}}_{T}(f)\right]&S_{yx}(f)&=\lim _{T\to \infty }{\frac {1}{T}}\left[{\hat {y}}_{T}^{*}(f){\hat {x}}_{T}(f)\right]\end{aligned}}$$ где опять же вклады $S_{xx}(f)$ и $S_{yy}(f)$ уже понятны. Обратите внимание, что $S_{xy}^{*}(f)=S_{yx}(f)$, поэтому полный вклад в перекрестную мощность, как правило, в два раза превышает реальную часть любого отдельного CPSD . Как и прежде, отсюда мы преобразуем эти произведения в преобразование Фурье временной свертки, которое при делении на период и доведении до предела $T\to \infty$ становится преобразованием Фурье функции взаимной корреляции . $${\begin{aligned}S_{xy}(f)&=\int _{-\infty }^{\infty }\left[\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }x_{T}^{*}(t-\tau )y_{T}(t)dt\right]e^{-i2\pi f\tau }d\tau =\int _{-\infty }^{\infty }R_{xy}(\tau )e^{-i2\pi f\tau }d\tau \\S_{yx}(f)&=\int _{-\infty }^{\infty }\left[\lim _{T\to \infty }{\frac {1}{T}}\int _{-\infty }^{\infty }y_{T}^{*}(t-\tau )x_{T}(t)dt\right]e^{-i2\pi f\tau }d\tau =\int _{-\infty }^{\infty }R_{yx}(\tau )e^{-i2\pi f\tau }d\tau ,\end{aligned}}$$ где $R_{xy}(\tau )$ это взаимная корреляция $x(t)$ с $y(t)$ и $R_{yx}(\tau )$ это взаимная корреляция $y(t)$ с $x(t)$. В свете этого PSD рассматривается как частный случай CSD для $x(t)=y(t)$. Если $x(t)$ и $y(t)$ являются действительными сигналами (например, напряжением или током), их преобразованиями Фурье ${\hat {x}}(f)$ и ${\hat {y}}(f)$ по соглашению обычно ограничиваются положительными частотами. Следовательно, при типичной обработке сигналов полный CPSD — это всего лишь один из CPSD , масштабированный в два раза. $$\operatorname {CPSD} _{\text{Full}}=2S_{xy}(f)=2S_{yx}(f)$$

Для дискретных сигналов x _n и y _n взаимосвязь между кросс-спектральной плотностью и кросс-ковариацией равна $$S_{xy}(f)=\sum _{n=-\infty }^{\infty }R_{xy}(\tau _{n})e^{-i2\pi f\tau _{n}}\,\Delta \tau$$

Оценка

Основная статья: Оценка спектральной плотности

Целью оценки спектральной плотности является оценка спектральной плотности случайного сигнала из последовательности временных выборок. В зависимости от того, что известно о сигнале, методы оценки могут включать параметрические или непараметрические подходы и могут быть основаны на анализе во временной или частотной области. Например, общий параметрический метод включает в себя подгонку наблюдений к авторегрессионной модели . Распространенным непараметрическим методом является периодограмма .

Спектральная плотность обычно оценивается с использованием методов преобразования Фурье (таких как метод Уэлча другие методы, такие как метод максимальной энтропии ), но также можно использовать и .

Связанные понятия

Не путать со спектральной плотностью (физика) .

• Спектральный центр сигнала — это середина его функции спектральной плотности, т. е. частота, которая делит распределение на две равные части.
• Граничная частота спектра ( SEF ), обычно выражаемая как «SEF x », представляет собой частоту, ниже которой находится x процентов от общей мощности данного сигнала; обычно x находится в диапазоне от 75 до 95. В частности, это популярный показатель, используемый при мониторинге ЭЭГ , и в этом случае SEF по-разному используется для оценки глубины анестезии и стадий сна .
• Спектральная огибающая — это огибающая спектральной плотности. Он описывает один момент времени (точнее, одно окно). Например, при дистанционном зондировании с использованием спектрометра спектральная огибающая объекта является границей его спектральных свойств, определяемой диапазоном уровней яркости в каждом из спектральных диапазонов . интересующих
• Спектральная плотность является функцией частоты, а не временем. Однако можно вычислить спектральную плотность небольшого окна более длинного сигнала и построить график зависимости от времени, связанного с окном. Такой график называется спектрограммой . Это основа ряда методов спектрального анализа, таких как кратковременное преобразование Фурье и вейвлеты .
• «Спектр» обычно означает спектральную плотность мощности, как обсуждалось выше, которая отражает распределение содержания сигнала по частоте. Для передаточных функций (например, график Боде , чирп ) полная частотная характеристика может быть представлена ​​в виде двух частей: зависимость мощности от частоты и фаза от частоты — фазовая спектральная плотность , фазовый спектр или спектральная фаза . Реже этими двумя частями могут быть действительная и мнимая части передаточной функции. Это не следует путать с частотной характеристикой передаточной функции, которая также включает фазу (или, что то же самое, действительную и мнимую часть) как функцию частоты. во временной области Импульсная характеристика $h(t)$ как правило, не может быть однозначно восстановлено только по спектральной плотности мощности без фазовой части. Хотя это также пары преобразований Фурье, здесь нет симметрии (как в случае автокорреляции ), заставляющей преобразование Фурье быть действительным. См. Ультракороткий импульс#Спектральная фаза , фазовый шум , групповая задержка .
• Иногда можно встретить амплитудную спектральную плотность ( ASD ), которая является квадратным корнем PSD; ASD сигнала напряжения имеет единицы измерения В Гц. ^−1/2 . Это полезно, когда форма спектра довольно постоянна, поскольку изменения ASD будут тогда пропорциональны изменениям самого уровня напряжения сигнала. Но с математической точки зрения предпочтительнее использовать PSD, поскольку только в этом случае площадь под кривой имеет смысл с точки зрения фактической мощности на всей частоте или в заданной полосе пропускания.

Приложения

Дополнительная информация: Спектр.

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

Электротехника

Спектрограмма FM-радиосигнала с частотой по горизонтальной оси и временем, увеличивающимся вверх по вертикальной оси.

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

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

Космология

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

См. также

• Биспектр
• Яркостная температура
• Цвета шума
• Спектральный анализ методом наименьших квадратов
• Спектральная плотность шума
• Оценка спектральной плотности
• Спектральная эффективность
• Спектральная утечка
• Спектральное распределение мощности
• Вероятность Уиттла
• Функция окна

Ссылки

• Биролини, Алессандро (2007). Инженерия надежности . Берлин ; Нью-Йорк: Springer Science & Business Media. ISBN  978-3-540-49388-4 .
• Браун, Роберт Гровер; Хван, Патрик Ю.К. (1997). Введение в случайные сигналы и прикладную фильтрацию Калмана с помощью упражнений и решений Matlab . Нью-Йорк: Вили-Лисс. ISBN  978-0-471-12839-7 .
• Давенпорт, Уилбур Б. (младший); Рут, Уильям Л. (1987). Введение в теорию случайных сигналов и шума . Нью-Йорк: Wiley-IEEE Press. ISBN  978-0-87942-235-6 .
• Имтиаз, Сайед Анас; Родригес-Вильегас, Эстер (2014). «Алгоритм с низкими вычислительными затратами для обнаружения быстрого сна с использованием одноканальной ЭЭГ» . Анналы биомедицинской инженерии . 42 (11): 2344–59. дои : 10.1007/s10439-014-1085-6 . ПМК   4204008 . ПМИД   25113231 .
• Иранманеш, Саам; Родригес-Вильегас, Эстер (2017). «Система обнаружения спящего шпинделя со сверхнизким энергопотреблением на кристалле». Транзакции IEEE в биомедицинских схемах и системах . 11 (4): 858–866. дои : 10.1109/TBCAS.2017.2690908 . hdl : 10044/1/46059 . ПМИД   28541914 . S2CID   206608057 .
• Марал, Жерар (2004). VSAT-сети . Западный Суссекс, Англия ; Хобокен, Нью-Джерси: Уайли. ISBN  978-0-470-86684-9 .
• Миллер, Скотт; Чайлдерс, Дональд (2012). Вероятность и случайные процессы . Бостон, Массачусетс: Академическая пресса. ISBN  978-0-12-386981-4 . ОСЛК   696092052 .
• Нортон, член парламента; Карчуб, Д.Г. (2003). Основы анализа шума и вибрации для инженеров . Кембридж: Издательство Кембриджского университета. ISBN  978-0-521-49913-2 .
• Оппенгейм, Алан В.; Вергезе, Джордж К. (2016). Сигналы, системы и вывод . Бостон: Пирсон. ISBN  978-0-13-394328-3 .
• Рискен, Ханнес; Франк, Тилль (1996). Уравнение Фоккера-Планка . Нью-Йорк: Springer Science & Business Media. ISBN  978-3-540-61530-9 .
• Стоун, Джонатан Ю. (2000). Цифровая обработка сигналов . Нью-Йорк Вайнхайм: Wiley-Interscience. ISBN  978-0-471-29546-4 .

Внешние ссылки

• Скрипты Matlab для определения спектральной плотности мощности