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

Спектр [pic], следовательно, должен быть суммой импульсов в этих точках.

Кроме того, поскольку возможно построить положительный полином, который

равен нулю в [pic] произвольно выбранных точках и нигде больше, то отсюда

следует, что [pic] имеет единственное спектральное представление в виде

суммы импульсов в общих нулях всех положительных полиномов [pic] так что

[pic].

В более широком смысле, теорема продолжимости совместно с теоремой

Каратеодори [16] показывает, что имеется по крайней мере одно спектральное

представление [pic] в виде суммы не более чем 2М импульсов.

Теорема представления: Если [pic], то существует [pic] и [pic], так

что

[pic] (4.7)

Доказательство теоремы представления можно найти в Приложении В. Это

представление и, таким образом, решение основной задачи оптимизации могут

быть не единственными. Дальнейшее обсуждение этой проблемы единственности

можно найти в Приложений С.

Если [pic] и местоположения импульсов в единственном решении [pic]

могут быть определены для данного [pic], то амплитуды импульсов могут быть

вычислены просто путем решения набора линейных уравнений. А сейчас мы

получим двойственную задачу оптимизации, которая дает [pic] и [pic], так

что [pic]. Тогда, если [pic] имеет единственное спектральное представление,

местоположения импульсов могут быть определены по нулям [pic]. Из теоремы

продолжимости следует

[pic] (4.8)

Так как [pic] и [pic], то отсюда следует, что [pic] и [pic] для всех [pic].

Кроме того, так как [pic] для некоторого [pic], то отсюда следует, что

[pic] (4.9а)

на множестве

[pic](4.9b)

и минимум достигается при [pic]. Решение этой двойственной задачи может не

быть единственным даже в случае временной последовательности, когда она

сводится к задаче собственного вектора, полученной Писаренко, и приводит к

интерпретации метода Писаренко в виде определения сглаживающего фильтра с

ограничениями по методу наименьших квадратов.

Пример 4.2 : Случай временной последовательности, [pic]. Как в примере

/3.1/

[pic].

Кроме того, если [pic] соответствует белому шуму единичной мощности,

[pic].

Таким образом, двойственная задача оптимизации сводится к нахождению

собственного вектора теплицевой матрицы, связанного с [pic],

соответствующего наименьшему собственному значению. Если имеется несколько

таких собственных векторов, импульсы располагаются в общих нулях

соответствующих полиномов. Любой нормированный собственный вектор,

соответствующий минимальному собственному значению, дает коэффициенты

сглаживающего фильтра, сумма квадратов величин которых ограничена единицей,

что дает наименьшую выходную мощность при наличии входного процесса,

корреляции которого описываются [pic][17].

1.4.2 Вычисление оценки Писаренко

При разработке алгоритмов вычисления оценки Писаренко можно

столкнуться с дискретной спектральной основой

[pic]

Для такой основы основная задача /4.4/ может быть переписана в виде

линейное программы стандартного вида

[pic] (4.11з)

так что для [pic]

[pic] (4.11b)

с N переменными и 2М ограничениями. Минимум равен [pic] и достигается для

[pic]. Основная теорема линейного программирования 18 эквивалентна теореме

представления в этом случае. При условии, что для этой линейной программы

существует решение, как показано в предыдущем разделе, основная теорема

гарантирует решение, в котором не более, чем 2М из [pic] не равны нулю, так

называемое, базовое решение.

Двойственная линейная программа [l5]

[pic] (4.12з)

так что для [pic]

[pic] (4.12b)

эквивалентная двойственной задаче /4.9/ для дискретной спектральной основы,

где ограничение

[pic] (4.13)

было использовано для исключения [pic] и где [pic]. Её минимум равен [pic]

и достигается при [pic].

Основная задача может быть решена при использовании симплекс-метода

[18]. Применение симплекс-метода к основной задаче приводит в результате к

существенно тому же результату /вычислительному алгоритму/, что и

применение, /одинарного/ метода замены к двойственной задаче [19]. Применив

соответствующий метод для избежания зацикливания [20], может быть получен

алгоритм, который гарантирует сходимость к оптимальному решению за конечное

число шагов, хотя его воплощения обычно были медленными .

Задача чебышевской аппроксимации связана с вычислением оценки

Писаренко; она может быть сформулирована, как минимизация линейного

функционала на выпуклом пространстве, определенном ограничениями типа

линейных неравенств [l6]. Она также решалась с использованием симплекс-

метода /одинарная замена/. Однако для частной задачи чебышевской

аппроксимации непрерывных функций полиномами с одной переменной существует

вычислительный метод, который значительно быстрее симплекс-метода, это

метод многократной замены Ремеза. Хотя были сделаны попытки распространить

этот метод на более общие задачи [21], появившиеся в результате алгоритмы

не достаточно хорошо понятны; в частности, не доказана их сходимость.

И наконец, задачи недискретной оптимизации, включенные в вычисление

оценки Писаренко, /4.4/ к /4.9/, являются видом, известным, как

полубесконечные программы. Как теоретические, так и вычислительные аспекты

таких программ рассматриваются в сборнике статей, изданных Геттичем [22].

Резюме

Эта статья связана с тем, что вероятно является наиболее простой и

интересной задачей в обработке антенных решеток; оценкой спектра мощности с

известной основой при условии, что даны некоторые выборки его

корреляционной функции. Хотя и простая, эта задача сохраняет несколько

черт, которые являются общими для многих задач обработки решеток:

многомерные спектры, корреляционные выборки с неравномерными отчетами и

произвольные спектральные основы.

Исследование спектральных оценок, согласованных с корреляцией привели

к задаче продолжимости. Были даны две характеристики продолжаемости ста

задача, для случая временных последовательностей, известна как задача

тригонометрических моментов и ее решение включает рассмотрение

положительной определенности корреляционных выборок. Положительная

определенность может поэтому рассматриваться как специальный случай

продолжимости.

Базируясь на теоретической основе, разработанной при решении задачи

продолжаемости, метод Писаренко был распространен со случая временных

последовательностей на задачу обработки решетки. Было показано, что метод

Писаренко тесно .связан с задачек продолжимости. Было показано, что

вычисление оценки Писаренко включает решение линейной задачи оптимизации.

Было показало, что решение этой задачи не является единственным в общем

случае, хотя оно единственно для случая временной последовательности, где

задача линейном оптимизации сводится к задача собственных значений.

Хотя рассмотренная в этой статье задача спектральной оценки была

разработала для обработки решетки, теоретическая структура и результирующие

алгоритмы должна быть полезными в других многомерных задачах, например,

обработке изображений.

2.1 ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ ДЛЯ ОТКРЫТОГО РЕЗОНАТОРА С ОСЕСИММЕТРИЧНЫМ

ДИСКОМ

В § 9.3 было получено интегральное уравнение (9.39) для резонатора с

диэлектрическим телом в виде шара. Такая форма диэлектрика хороша для

анализа, но неудобна для практики.

Обычно приходится иметь дело с диэлектрическими образцами более

сложной формы, в частности с диэлектрическим диском. В такой ситуации

получить аналитическое выражение для ядра не удается, однако это не

является препятствием для нахождения решения задачи.

Действительно, ядро уравнения для резонатора с шаром (9.39) — это

сумма ядра для пустого резонатора и дополнительного члена, представляющего

собой поле, рассеянное шаром.

Запишем уравнение для резонатора с диском в аналогичном виде,

поскольку физическая картина явлении одна и та же:

[pic] (9.45)

Здесь [pic] - ядро пустого резонатора; Т — ядро, связанное с

рассеянием на диэлектрическом образце. Обсудим, что в сущности делается при

решении уравнения (9.39) методом Галеркина. Для определенности будем

считать, что в качестве базисных и весовых (см. приложение 2) взяты

собственные функции резонатора без шара, которые обозначим [pic] и будем

считать ортонормированными.

С первым слагаемым ядра все ясно, базисные функции являются его

собственными, и действие интегрального оператора с таким ядром эквивалентно

умножению на постоянную, являющуюся собственным значением пустого

резонатора:

[pic] (9.46)

Интегральный оператор со вторым слагаемым ядра представляет собой

магнитное поле тока на зеркалах, рассеянное шаром. Плотность тока задается

в виде [pic], а рассеянное поле рассчитывается на поверхности зеркала. При

решении (9.39) расчет рассеянного шаром поля проводится аналитически.

Однако ту же процедуру можно произвести численно, и тогда ограничения на

формулу диэлектрического образца в значительной степени снимаются.

Для расчета рассеянного поля будем применять интегральное уравнение

(3.85). Диэлектрический образец может быть произвольным телом вращения, в

частности диском.

После этих общих соображений рассмотрим процедуру решения (9.45)

последовательно. Функция U(x) ищется в виде

[pic] (9.47)

В соответствии с методом Галеркина (см. приложение 2) подставляем (9.47) в

(9.45), затем умножаем на [pic] и повторно интегрируем по образующей

зеркала. С учетом ортонормированности базисных функций имеет однородную

СЛАУ

[pic] (9.48)

где [pic] - собственные числа уравнения невозмущенного резонатора [см.

(9.46)].

Элементы матрицы СЛАУ выражаются интегралами

[pic] (9.49)

Последнюю формулу надо понимать как символическую. Она эквивалентна

процедуре расчета рассеянного поля, описанной выше. Остановимся на ней

подробнее.

Вначале необходимо найти поле на поверхности диэлектрического тела,

созданное током вида [pic] на зеркалах. Это можно было бы сделать с помощью

(3.8), (3.9), однако есть более простой путь, если ограничиться

рассмотрением тел небольших, на порядок меньших диаметра зеркал. Тогда

можно воспользоваться приближенным выражением для поля в резонаторе,

соответствующим приближенным функциям токов на зеркалах. На рис. 9.6

представлены графики распределения токов на зеркалах, соответствующие

низшему типу колебаний [pic] и колебанию, имеющему вариацию по радиусу

[pic]. Резонатор конфокальный с параметром [pic]. Вблизи оси плотность

тока, описываемая гиперсфероидальными функциями (кривые 1), практически не

отличаются от экспоненциальной функции, умноженной на полиномы Лагерра

(кривые 2), т. е. от гауссова пучка [68]. Радиальное распределение

отличается только масштабом по радиусу.

Таким образом, будем описывать поле в резонаторе вблизи его центра

приближенным .выражением в виде гауссова пучка

[pic] (9.50)

где

[pic];

R - радиус кривизны волнового фронта; W — радиус «освещенного пятна» в

пучке. Последняя величина определяется как радиус, на

[pic]

Рис. 9.6. Сравнение точных и приближенных кривых для

гиперсфероидальных функций:

1 - точные, 2 - приближенные кривые

котором интенсивность пучка спадает в е раз по отношению к центру пучка.

Характерной величиной для каждого пучка является наименьший радиус «пятна»

[pic]. Применительно к резонатору - это радиус «пятна» в центре, который

связан с длиной резонатора 1:

[pic] (9.51)

1 Как и ранее, все длины предполагаются умноженными на волновое число,

которое здесь соответствует действительной части собственной частоты

невозмущенного резонатора.

Величины W и R медленно меняются вдоль резонатора:

[pic] (9.52)

[pic] (9.53)

В центре резонатора [pic] Естественно в резонаторе существуют не

один, а два встречных гауссовых пучка, и вблизи центра поле основной моды в

приближении гауссова пучка имеет вид

[pic] (9.54)

На зеркале [pic] для конфокальной геометрии резонатора в соответствии

с (9.51)—(9.53) [pic], и распределение тока имеет вид1

[pic] (9.55);

Для следующего колебания «1, 0, q» поле в центре резонатора представляется

формулой

[pic] (9.56)

и на зеркалах

[pic](9.57)

Таким образом, поле в резонаторе без образца, соответствующее различным

модам, в приближении гауссова пучка нетрудно записать. Оно играет роль

первичного поля для задачи возбуждения диэлектрического образца.

Вычисляем эквивалентные токи на поверхности диэлектрика в предположении,

что основная поляризация поля [pic]. В обозначениях § 3.3 имеем:

1 Напомним, что в открытых резонаторах с круглыми зеркалами принята

следующая индексация мод : первый индекс - число вариаций по R, второй -

число вариаций по [pic], а третий - число вариаций по [pic]

[pic] (9.58)

Теперь необходимо возвратиться к азимутальным гармоникам вида [pic],

поскольку ЭВМ — программы для диэлектрических тел вращения сделаны

применительно к ним. Первичные токи представляют собой сумму первой и минус

первой гармоник. Каждую из них можно выделить, используя формулу Эйлера. В

результате решения задачи возбуждения диэлектрического тела, а конкретно

диска, получаем значения эквивалентных токов в дискретных и достаточно

часто расположенных точках образующей. Зависимость от [pic] этих токов

известная. Если объединить токи первой и минус первой гармоник, она будет

такой же, как и у первичных токов (9.58).

Следующий этап — вычисление рассеянных диском полей на зеркалах. Для

этого используются формулы (3.8), (3.9). Выра жения для элементов тензорной

функции Грина следует упрос тить, как и при выводе уравнений (9.5)—(9.8),

т. е. положить [pic], а для функции [pic] использовать асимптотическую

формулу (9.22). Последняя содержит множитель, учитывающий набег фазы на

половине размера резонатора (расстояние от образца до одного из зеркал).

Такой же набег фаз имеется в первичном для диэлектрического образца поле.

Этот сдвиг присутствует также в (9.56) и (9.57). Все это позволяет вынести

за знак интеграла множитель [pic], такой же, как и из основного ядра. Этот

множитель, как и ранее, дает основную частотную зависимость. Ядра без него

от частоты зависят слабо, и в них частота полагается равной действительной

части собственной частоты пустого генератора.

Теперь уже можно вычислить элементы матрицы (9.48). Для определения

элемента [pic] берется рассеянное поле, возбужденное нулевой модой пустого

резонатора, т. е. [pic], затем оно в соответствии с (9.49) домножается на

(9.55) и интегрируется. При этом необходимо помнить, что базисные функции

предполагались нормированными. Поэтому функцию (9.55) необходимо

предварительно пронормировать. В силу осевой симметрии системы

поверхностный интеграл (9.49) можно представить в координатах вращения.

Интеграл по [pic] берется аналитическим, а по радиальной координате [pic] -

численно. Остальные элементы [pic] отыскиваются точно так же.

Далее решается задача на собственные значения, а затем с помощью

формул (9.40) и (9.41) находятся изменения добротности и сдвиг частоты.

2.2 ИНТЕГРАЛЬНОЕ УРАВНЕНИЕ ОТКРЫТОГО РЕЗОНАТОРА С ДИЭЛЕКТРИЧЕСКИМ

ДИСКОМ, НЕСООСНЫМ С ЗЕРКАЛАМИ [72]

При проведении измерений параметров диэлектрика образец в виде диска

часто удобнее расположить несоосно с зеркалами и, в частности, так, чтобы

оси резонатора и диска были перпендикулярны (рис. 9.7). Такое расположение

диска нарушает осевую симметрию задачи. В общем случае отход от осевой

симметрии очень -сильно усложняет решение, поскольку теряется основное

преимущество систем вращения — независимость отдельных азимутальных

гармоник полей.

[pic]

Рис. 9.7. Геометрия открытого резонатора с несоосными зеркалом и диском

Однако в рассматриваемой задаче анализа полей в высокодобротном

открытом резонаторе несоосность вносит технические, но не принципиальные

затруднения. Действительно, для измерений параметров диэлектрический

образец берется небольшим по сравнению с размерами резонатора. Поэтому его

внесение в резонатор не приводит к переходу к другой моде, а лишь несколько

меняет добротность и резонансную частоту той моды, которая существовала без

диэлектрика. Таким образом, за счет фильтрующих свойств резонатора новых

азимутальных гармоник не появляется и основная трудность в несоосных

системах вращения снимается. Надо лишь следить за тем, чтобы на других

азимутальных гармониках у пустого резонатора не было поблизости от частоты

рабочей моды других высокодобротных мод.

Метод решения задачи остается в общих чертах тем же, что и в

Страницы: 1, 2, 3, 4