
- •Цифровой спектральный анализ
- •Оглавление
- •Введение
- •Благодарности
- •Обращение к читателю
- •1. Основные понятия, задачи и проблемы цифрового
- •1.1. Что такое цифровой спектральный анализ
- •1.2. Проблемы в области спектрального оценивания
- •1.2.1. Разрешение
- •1.2.2. Обнаружимость сигналов
- •1.2.4. Состоятельность оценки
- •1.4. Взаимная связь спектров континуальных и дискретных сигналов
- •2. Спектральный анализ и быстрое преобразование
- •2.2. Перестановка данных и двоичная инверсия
- •2.3. Алгоритм бпф с прореживанием по частоте
- •2.4. Вычисление обратного дпф с помощью алгоритма прямого дпф
- •2.5. Единый подход к алгоритмам бпф
- •2.6. Особенности цифровой реализации алгоритма бпф
- •3. Соотношения между бпф и цифровой фильтрацией
- •3.1. Некоторые характеристики спектрального анализа на основе бпф
- •3.2. Соотношение между "скачущим" бпф и гребенкой фильтров
- •4. Использование окон в спектральном анализе
- •4.1. Прямоугольное окно и бпф
- •4.2. Наиболее распространенные функции окна и их характеристики
- •5. Классические методы спектрального оценивания
- •5.1. Периодограммные методы спектрального анализа
- •5.1.1. Периодограмма Даньелла
- •5.1.2. Периодограмма Бартлетта
- •5.1.3. Периодограмма Уэлча
- •5.2. Корреляционный анализ и спектральная оценка
- •5.2.1. Метод получения оценки фвк
- •5.2.2. Коррелограммный метод оценки спм
- •5.4. Комбинированные периодограммно-коррелограммные оценки
- •6. Алгоритмы спектрального анализа в ограниченном
- •6.1. Алгоритм анализа с использованием бпф
- •6.2. Алгоритм Блюстейна
- •6.3. Алгоритм z-преобразования с использованием лчм-фильтрации
- •7. Спектральный анализ, основанный на линейном
- •7.1. Модели случайных процессов
- •7.2. Методы оценки параметров ар-моделей. Блочные алгоритмы
- •7.4.1. Геометрический алгоритм
- •7.4.2. Гармонический алгоритм Берга
- •7.8.1. Выбор порядка модели
- •7.8.2. Авторегрессионные процессы с шумом наблюдения
- •7.8.3. Оценивание мощности синусоидальных компонент
- •7.8.4. Частотное разрешение
- •7.9. Метод Прони
- •7.9.1 Исходный метод Прони
- •7.9.2. Метод наименьших квадратов Прони
- •Условные обозначения, принятые в книге
- •Литература
1.2. Проблемы в области спектрального оценивания
Интерес к цифровым методам спектрального анализа поддерживается тем улучшением характеристик, которое они обещают, а именно: высоким частотным разрешением, повышенной способностью к обнаружению слабых сигналов или же сохранением "достоверности" формы спектра при меньшем числе используемых параметров. Аналитически описать характеристики большинства методов в случае ограниченного временного анализа (т. е. в случае короткой реализации данных) весьма затруднительно, именно поэтому в литературе можно найти лишь очень малое количество эмпирических результатов. Это обусловило появление ряда проблем в области современного спектрального оценивания, некоторые из них кратко освещены ниже.
1.2.1. Разрешение
Спектральное разрешение относится к числу главных проблем современного спектрального оценивания, в особенности применительно к анализу коротких последовательностей данных. При этом то, что понимается под термином "разрешение", носит весьма субъективный характер. Одно из ранних определений принадлежит Релею, которое исходит из следующего рабочего определения для разрешения оптических телескопов с ограниченной пространственной апертурой: "Подобно тому, как оптическая сила телескопа измеряется близостью двойных звезд, которые он может разрешить, так и оптическую силу спектроскопа следует измерять близостью самых близких двойных линий в спектре, которые он может разрешить".
Рис. 1.2
В литературе принято характеризовать относительные величины разрешающей способности двух спектральных оценок на основе визуальных наблюдений. Рассмотрим две спектральные оценки, показанные на рис. 1.2. Спектральная оценка на рис. 1.2, а имеет единичный нормированный максимум с двумя слабо различимыми пиками, что должно указывать на присутствие двух спектральных составляющих. Если теперь нелинейным масштабированием сформировать новую спектральную функцию, которая показана на рис. 1.2, б, то она будет выглядеть, как спектр с "более высоким" разрешением. На самом же деле оба приведенных спектра несут в себе одну и ту же информацию, а это означает, что одного лишь визуального сравнения спектральных оценок, конечно же, недостаточно для суждения о характеристиках разрешения.
1.2.2. Обнаружимость сигналов
Использование
спектрального оценивания для обнаружения
присутствия сигналов приводит к
возникновению другой проблемы. Рассмотрим
периодограммы, показанные на рис. 1.3.
Обе они соответствуют одной и той же
1024-точечной реализации некоторого
процесса, состоящего из аддитивной
смеси одной синусоиды и белого шума при
отношении сигнал/шум равном -17 дБ.
Периодограмма на рис. 1.3, а
получена непосредственно по всем
1024 отсчетам, и в ней над СПМ шума отчетливо
видна спектральная линия синусоиды с
относительной частотой 0.18 . На рис. 1.3,
б
представлена сглаженная периодограмма,
полученная в результате разбиения
исходной последовательности на 32
сегмента по 32 отсчета в каждом и
последующего усреднения периодограмм
этих сегментов. На этой оценке более
очевидна равномерная спектральная
характеристика белого шума, но отклик,
обусловленный присутствием синусоиды,
здесь не обнаружим. Можно показать,
что обнаружимость сигнала по несглаженной
периодограмме превышает обнаружимость
этого сигнала по сглаженной периодограмме
в число раз, равное примерно корню
квадратному из числа сегментов, или
на 10lg(321/2)=7.5
дБ в рассматриваемом случае. С точки
зрения отображения полного спектра
процесса (и оценки его СПМ) лучше
сглаженная спектральная оценка,
показанная на рис. 1.3, б,
но с точки зрения обнаружимости сигнала
- лучше
Рис. 1.3
несглаженная оценка.
Таким образом, аргументы в пользу выбора той или иной спектральной оценки будут зависеть от того, интересует ли нас сглаженная оценка в пределах всего диапазона анализируемых частот или же нам важна более высокая степень обнаружимости сигнала на некоторых конкретных его участках. Поэтому алгоритмы спектрального оценивания вовсе не обязательно являются также и хорошими алгоритмами обнаружения.
1.2.3. О разграничении понятий оценивания параметров сигнала со
спектральным оцениванием
Задача спектрального оценивания подразумевает оценивание некоторой функции частоты. О характеристиках спектральной оценки судят по тому, насколько хорошо она согласуется с известным спектром тест-сигнала в некоторой непрерывной области частот. Имеются приложения, где основной интерес представляет лишь локальное поведение спектральной оценки в некоторых заданных диапазонах частот. Примером может служить оценивание частоты синусоиды в белом шуме, результаты которого представлены на рис. 1.3. Нетрудно видеть, что в качестве спектральной оценки лучше использовать оценку на рис. 1.3, б, так как она ближе к равномерному спектру белого шума. В качестве же оценки частоты синусоиды лучше использовать оценку на рис. 1.3, а, так как она позволяет более точно определить частоту этой синусоиды. К сожалению, эти оценки нередко путают в литературе, поскольку многие задачи оценивания параметров излагаются в контексте спектрального оценивания, даже если функция спектральной оценки не является искомой количественной величиной.
Задача оценивания параметров и задача спектрального оценивания требуют различного статистического подхода, хотя обе они могут перекрываться в некоторых частных случаях. Поэтому о характеристиках спектральной оценки нельзя полностью судить по ее характеристикам как оценки параметров. Так, например, сравнение оценок, представленных на рис. 1.3, показывает, что сглаженная периодограмма является хорошей спектральной оценкой, но плохой оценкой параметра, которым в данном случае является частота синусоиды.