
- •Цифровой спектральный анализ
- •Оглавление
- •Введение
- •Благодарности
- •Обращение к читателю
- •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. Метод наименьших квадратов Прони
- •Условные обозначения, принятые в книге
- •Литература
5. Классические методы спектрального оценивания
Существуют два формальных, но эквивалентных типа методов определения спектральной плотности мощности (СПМ). Прямые методы основаны на вычислении квадрата модуля преобразования Фурье для бесконечной последовательности данных с использованием соответствующего статистического усреднения и исторически получили название периодограммных методов. Косвенные методы основаны на использовании бесконечной последовательности значений данных для расчета автокорреляционной последовательности (АКП), преобразование Фурье которой дает искомую оценку СПМ. Они называются коррелограммными методами спектральной оценки. Эти методы основаны на связи АКП и СПМ дискретным преобразованием Винера-Хинчина.
Оба типа методов применяются сравнительно давно и поэтому получили название классических методов спектрального оценивания. Все эти методы обычно используются применительно к дискретно-временным последовательностям отсчетов данных.
Мы говорим об оценке СПМ, так как на практике в спектральном анализе используют конечные последовательности данных вместо бесконечных, что обусловливает невозможность получения точного спектра. Оценка характеризуется следующими параметрами.
Смещением
оценки
параметра
называется разность между истинным
значением этого параметра и математическим
ожиданием этой оценки:
Несмещенной
оценкой
называется такая, для которой
.
Дисперсией оценки называется средний квадрат отклонения оценки от ее среднего значения
Говорят, что оценка состоятельна, если с увеличением числа наблюдений смещение и дисперсия стремятся к нулю.
Классические методы спектрального оценивания относятся к числу наиболее устойчивых (робастных) методов спектрального оценивания. Они применимы почти ко всем классам стационарных сигналов и шумов, тогда как методы высокого разрешения (основанные на линейном моделировании) оказываются рабастными только в случае ограниченного класса стационарных сигналов.
В классических методах, как правило, используются алгоритмы БПФ, поэтому они оказываются наиболее вычислительно эффективными. В получаемых с их помощью оценках СПМ (в отличие от методов высокого разрешения) высота спектральных пиков линейно зависит от мощности синусоид, присутствующих в данных.
Основной недостаток классических методов спектрального оценивания обусловлен искажающим воздействием просачивания по боковым лепесткам из-за использования конечных последовательностей данных. Взвешивание позволяет ослабить влияние боковых лепестков, но лишь за счет ухудшения спектрального разрешения.
Разрешение, обеспечиваемое классическими методами, не может превышать величины, обратной длине записи данных, и не зависит от характеристик анализируемых данных. Для улучшения статистической устойчивости спектральных оценок может быть использовано псевдоусреднение по ансамблю за счет сегментации данных, но это также ухудшает спектральное разрешение. Обработка конечных записей данных требует принятия компромиссов относительно разрешения, устойчивости (минимизации дисперсии оценки) и подавления просачивания.
И периодограммные, и коррелограммные методы дают оценки, которые обладают аналогичными статистическими характеристиками и в целом выглядят примерно одинаково, если не считать небольших визуальных различий в тонких деталях формы спектров. На практике предпочтение отдается наиболее эффективным в вычислительном отношении методам. До широкого распространеня ЭВМ в основном применялся корреляционный метод. С появлением алгоритмов БПФ и спецпроцессоров ЦОС предпочтение было отдано периодограммным методам.