имеется ли единственная согласованная с. корреляцией спектральная оценка и, если нет, как выбрать нужную ? Действительно, единственная оценка не существует, за исключением весьма специальных случаев; задача метода спектральной оценки состоит в выборе одного из ансамбля спектров, удовлетворяющего согласованию корреляции, положительности и ограничениям спектральной основы. Раздел IУ касается метода Писаренко [5] , который включает моделирование корреляционных измерений в виде суммы двух компонентов. Один, шумовой компонент известной спектральной формы, но неизвестной амилитуды, делается настолько большим, насколько это возможно без превращения второго компонента в непродолжаемый. Показано, что спектральная оценка по методу Писаренко решает линейную задачу оптимизации. Решение этой задачи оптимизации будет всегда существовать, если корреляционные измерения являются продолжаемыми. Показало, что тактически метод Писаренко тесно связан с вопросом продолжаемости и алгоритм вычисления оценки Писаренко будет также служить в качестве теста продолжаемости. Показано, что оценка Писаренко не является всегда единственной в общем случае, хотя она единственна для случая временной последовательности, где задача линейной оптимизации сводится к задаче на собственные значения.

1.1. Задача обработки решетки

Вообразим многомерную однородную среду, поддерживающую волновое поле с комплексными значениями u(x, t) и содержащую решетку датчиков. Волновое паче будет предполагаться однородным и стационарным, так что его статистики второго порядка описываются корреляционной санкцией r , или эквивалентно, спектром мощности [6].

(2.1)

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

В инженерной литературе более принято представлять спектр мощности посредством положительной функции спектральной плотности . В этом представлении

(2.2)

где - некоторая фиксированная мера, которая позволяет интерпретировать выражение /2.2/ в виде многомерной поверхности или объемного интеграла, возможно взвешенного, над частотно-волновым векторным пространством.

Если дана 'функция спектральной плотности мощности , то возможно определить соответствующую положительную меру путем требования, чтобы мера подмножества В частотно-волнового векторного пространства равнялась интегралу функции спектральной плотности по В:

(2.3)

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

Каждый из ПИП производит временную функцию, которая является волновым полем U, подвергнутым выборке в точке пространства. Совокупность временных функций, образуемых всеми ПИП, выход или отклик решетки, должна быть обработана с тем, чтобы обеспечить оценку спектра мощности частоты-волнового вектора. Стохастический характер волнового поля неизменно приводит к случайным 'изменениям любой спектральной оценки, основанной на выходе решетки. Чтобы противодействовать этому эффекту, спектральные оценки часто базируются на устойчивых статистиках, получаемых с выхода решетки. Обычным примером такой статистики является корреляционная оценка, вычисляемая посредством умножения выхода одного ПИП на задержанный во времени выход второго ПИП с усреднением по времени. Эта обработка дает в результате оценку корреляционной функции с временной задержкой, соответствующей запаздыванию во времени и пространственным разделением, которое является вектором расстояния между ПИП. Процесс усреднения обеспечивает статистически стабильные оценки корреляции, что дает в результате статистическую стабильность спектральной оценки, основанной на этих корреляционных оценках. Важно отметить, что оценки корреляций доступны только для конечного множества междатчиковых расстояний и временных задержек [8]. Тема ошибок корреляционных оценок не будет затрагиваться. Эта. статья касается скорее свойств множеств истинных корреляционных выборок и спектральных оценок, основанных на корреляционных выборка.

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

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

вещественные значения. К счастью, обычный анализ Фурье является часто удовлетворительным, когда данные избыточны, а также неявным при узкополосном характере многих датчиков. Там, где ограниченные данные во временное измерении делают упомянутый выше подход не практичным, а доступными являются широкополосные решетки датчиков, полная задача может трактоваться посредством включения временных переменных и в векторы и k. Тогда будет описывать разделение как в пространстве, так и во времени, a k волновой вектор пространства-времени. Будем полагать, что принят один из этих двух подходов; следовательно временные переменные и будут опущены.