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

Question

Кто-нибудь знает, возможно ли найти спектральную плотность мощности сигнала с пробелами в нем? Например (в синтаксисе matlab это то, с чем я знаком)

ta=1:1000; 
tb=1200:3000;

t=[ta tb];    % this is the timebase
signal=randn(size(t));  this is a signal

figure(101)

plot(t,signal,'.')

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

Любые мысли приветствуются.

Answer 1

Периодограмма Lomb-Scargle обычно используется для выполнения анализа неравномерно распределенных данных (выбранных в произвольные моменты времени) или когда часть данных отсутствует.

Вотпара реализаций MATLAB:

• lombscargle.m (FEX)
• Lomb (Lomb-Scargle) периодограмма (FEX)
• lomb.m - инструменты ЭКГ от Гари Клиффорда

Answer 2

Отсутствие сегментов базисных векторов Фурье приводит к точно такой же FT, то есть PSD, что и использование полного базиса, но умножение на нули в пределах заполнения нулями в любых сигнальных «пробелах».

Answer 3

Я нашел это Неравномерное БПФ , но я не уверен, что это именно то, что мне нужно, поскольку это может быть действительно для данных, которые в основном отбираются на неравномерной временной основе, а не для равномерно распределенных данных с значительные пробелы. Я попробую!