Число действительных умножений при вычислении свертки двух N-точечных последовательностей
Недостатком этого метода является значительные ошибки округления, большой объем памяти, необходимый для хранения комплексных экспоненциальных коэффициентов, и все ещё значительный объем вычислений. 19. Вычисление линейной свертки с секционированием.
Существует два метода с секционированием свертки: метод перекрытия с суммированием и метод перекрытия с накоплением. Рассмотрим первый из них. Пусть более длинной, а в общем случае неограниченной является последовательность
где Линейная свертка последовательностей
Меняя порядок суммирования и учитывая то, что последовательности
Отсюда следует, что в самом начале вычисляется k -я частичная линейная свертка последовательностей Другой метод вычисления линейной свертки последовательностей, одна из которых значительно длиннее другой, также основан на секционировании более длинной последовательности. Однако в данном случае перекрываются смежные исходные секции. Ошибочные отсчеты круговых сверток отдельных секций отбрасываются. Остальные отсчеты накапливаются и из них формируется конечный результат. В этом методе неограниченная последовательность
![]()
Рис. 1.32. Вычисление свертки После этого вычисляются секционированные круговые свертки
Следует отметить, что в этом случае Представим последовательности
где Тогда частичные свертки
где Поэтому при формировании суммарной свертки последние ( Рис. 1.33. Вычисление свертки Таким образом, используя метод перекрытия с суммированием или метод перекрытия с накоплением, можно достаточно легко найти свертку короткой и очень длинной последовательностей. При этом требуемый результат получается в виде отдельных небольших секций, которые должны объединяться в одну последовательность.
20. Амплитудный спектр, спектр мощности. Определение и алгоритмы получения.
Основу спектрального анализа, как известно, составляет преобразование Фурье. При этом, для детерминированных периодических сигналов
где
– коэффициенты соответствующего ряда Фурье. Рассмотрим следующее определение спектральной плотности мощности, которое широко использовалось на практике до появления алгоритмов быстрого преобразования Фурье:
где x (t,w,Δw) – процесс (сигнал) на выходе полосового фильтра с центральной частотой w и полосой пропускания Δw. Запишем данное выражение в более привычном виде:
где В – ширина полосы пропускания узкополосного фильтра с центральной частотой f. Выясним физический смысл выражения для
– энергия, а величина
является средней мощностью сигнала. Тогда выражение
определяет среднюю мощность сигнала на выходе полосового фильтра с центральной частотой f и полосой пропускания В. Следовательно, спектральная плотность мощности
Для получения оценки спектра мощности в некотором частотном диапазоне достаточно иметь или параллельный набор полосовых фильтров или же один полосовой фильтр с изменяемой центральной частотой. В первом случае получают спектр с постоянным относительным разрешением, во втором – с постоянным абсолютным разрешением. Параллельный метод построения анализаторов спектра широко применяется и в настоящее время в спектральном анализе, звуковых измерениях и анализе шума. Однако, в последнее время наиболее часто используется определение спектральной плотности мощности, основанное на непосредственном преобразовании Фурье исследуемой реализации:
где М – оператор статистического усреднения. 21. Оценка спектра мощности на основе периодограммы. Свойства периодограммы. Методы получения состоятельных периодограммных оценок.
В последнее время наиболее часто используется определение спектральной плотности мощности, основанное на непосредственном преобразовании Фурье исследуемой реализации:
где М – оператор статистического усреднения. Из данного определения оценка спектральной плотности мощности может быть получена в следующем виде
где Здесь Основные свойства этой оценки (тут ошибка – см. в тетрадь(левого предела не должно быть)):
т. е. данная оценка является асимптотически несмещенной. Дисперсия данной величины
Это значит, что асимптотически несмещенная оценка Для получения эффективных оценок применяются методы сглаживания. В этом случае для получения правильных результатов при измерении
Это обстоятельство и лежит в основе получения состоятельных оценок спектральной плотности мощности, т. е. путем сглаживания (усреднения) оценки
0 M f i M f Рис. 1.36. Иллюстрация сглаживания по частоте Для этой же цели, кроме того, применяется и усреднение по коротким периодограммам. В этом случае исходная реализация исследуемого сигнала x(t) длительностью T делится на более короткие реализации xi(t) длительностью
Рис. 1.37. Иллюстрация сглаживания
Измерение (оценка) спектра мощности (часто называемого энергетическим спектром) дает возможность, например, получать информацию о динамических характеристиках линейных физических систем с постоянными параметрами, позволяет исследовать соотношения между процессами на входе и выходе таких систем, обнаруживать скрытые периодичности и т. д. В теоретических исследованиях принято чаще всего говорить об оценке спектра или спектральном оценивании. в настоящее время в спектральном анализе используются оценки спектральной плотности мощности, основанные на прямом преобразовании исходных данных и последующем их усреднении. Этот метод, как уже отмечалось, чаще называют методом периодограммной оценки спектра мощности. Для того, чтобы по отсчетам обрабатываемого сигнала можно было бы получить спектральные оценки в соответствующих единицах энергии или мощности, необходимо выражение для прямого ДПФ умножить, а для обратного ДПФ разделить на интервал дискретизации D t:
где В этом случае оценка спектральной плотности мощности будет определяться следующим образом:
где Эта оценка называется выборочным спектром, периодограммой Шустера или просто периодограммой. Данная оценка также не является состоятельной оценкой истинной спектральной плотности мощности (СПМ), так как дисперсия этой величины не стремится к нулю ни при каком сколь угодно большом значении N. Вследствие этого для получения состоятельных оценок требуется выполнение операции статистического усреднения. В этом случае будем иметь
Для расчетов используется выражение
которое называют исходной немодифицированной формой периодограммной оценки СПМ. Для сглаживания периодограммной оценки используются три основных метода: метод Даньелла (Даниелла), Бартлетта и Уэлча. В методе Даньелла осуществляется усреднение оценок, полученных по соседним частотам (усреднение по смежным частотам), Бартлетта – по ансамблю (по коротким временным последовательностям), а в методе Уэлча подход Бартлетта применяется к перекрывающимся реализациям для уменьшения смещения оценок из-за эффекта просачивания. Практическое использование этих трех процедур подтверждает их статистическую устойчивость для многих классов сигналов. Периодограмма Даньелла. Для сглаживания быстрых флуктуаций выборочного спектра в этом случае используется усреднение по соседним спектральным частотам. Если для вычисления выборочного спектра
Вычисление оценки по Даньеллу рекомендуется для случаев, когда анализируемое множество данных состоит из малого (100–500) или среднего (500–4000) числа выборок. Периодограмма Бартлетта. При этом подходе последовательность входных данных х (п)из N отсчетов делится на K неперекрывающихся сегментов по М отсчетов в каждом, так что Тогда i -ый сегмент будет определяться таким образом:
Затем на каждом из этих сегментов независимо вычисляется выборочный спектр:
Далее на каждой частоте, представляющей интерес, K отдельных немодифицированных периодограмм усредняются с тем, чтобы получить усредненную периодограмму Бартлетта.
Рис. 1.38. К иллюстрации периодограммы Бартлетта
Дисперсия рассмотренной оценки уменьшается с увеличением числа K, а величина смещения – увеличивается, так как при фиксированной выборке N с увеличением числа сегментов число выборок М в каждом из них уменьшается. Это приводит к ухудшению разрешающей способности спектрального анализа, так что приходиться находить компромиссное решение между значениями N и М. Данная оценка применяется при N > 2000. Периодограмма Уэлча. Уэлч модифицировал основную схему Бартлетта за счет использования перекрывающихся сегментов (рис. 1.39). Цель перекрытия– увеличить число усредняемых оценок спектральной плотности мощности при заданной длительности исходной реализации и тем самым уменьшить дисперсию результирующей оценки. На основе БПФ Уэлч разработал также эффективную вычислительную процедуру для реализации данного метода, что и сделало метод Уэлча самым популярным периодограммным методом спектрального оценивания. Если выборка из N отсчетов разбита на К сегментов по М отсчетов в каждом со сдвигом S отсчетов между соседними сегментами Например, 50 %-е перекрытие сегментов во многих случаях обеспечивает весьма эффективную реализацию данного метода на основе алгоритмов БПФ. Кроме того, в этом случае все данные используются дважды, за исключением М /2 отсчетов на каждом конце исходной N -точечной последовательности данных. Следует отметить, что на практике часто используется перекрытие до 70 %.
Рис. 1.39. Формирование периодограммы Уэлча
Также как и дисперсия периодограммы Бартлетта, дисперсия периодограммы Уэлча примерно обратно пропорциональна числу сегментов в предположении независимости сегментов (хотя перекрытие сегментов приводит к некоторой их взаимозависимости). Благодаря перекрытию по заданной выборке исходных данных можно сформировать большее число сегментов, чем в методе Бартлетта, что уменьшает величину дисперсии периодограммы Уэлча по сравнению с дисперсией периодограммы Бартлетта. 22. Основные проблемы цифрового спектрального анализа. Взвешивание. Свойства весовых функций. Модифицированные периодограммные оценки СПМ.
|