Оптимизации топологии датчика волнового фронта

Приемная антенна поста РТН представляла собой фазированную антенную решетку с параметрами: ширина диаграммы направленности в диапазоне частот, в горизонтальной плоскости ©0,5=600±100; ширина диаграммы направленности в диапазоне частот, в вертикальной плоскости: Ф0,5=200±100. На приемной антенне помещался малошумящий усилитель (Кш< 1,5 дБ, Ку> 25дБ). Длина кабеля снижения составляла 15м. Высота подъема приемных антенн и антенн СПИ была около 10 м над поверхностью земли.

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

Результаты эксперимента показывают то, что на основе АОИПС может быть построена РДС, принцип действия которой основан на измерении времени прихода по фронту радиоимпульса. Данная система может быть использована для решения задач радиоконтроля, охватывая зону контроля порядка 30 км на 30 км.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Логинов Н.А. Актуальные вопросы радиоконтроля в Российской Федерации. - М.: Радио и связь, 2000. - 240 .: ил.

2. Кондратьев В.С., Котов А.Ф., Марков Л.Н. Многопозиционные радиотехнические системы/; Под ред. Проф. В.В. Цветнова. - М.: Радио и связь, 1986. -264 с.

УДК 535.375:551.521

Д.А.Безуглов, И.А.Сахаров, И.В. Решетникова

ОПТИМИЗАЦИИ ТОПОЛОГИИ ДАТЧИКА ВОЛНОВОГО ФРОНТА

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

Одним из наиболее эффективных (иногда в сочетании с другими) способов ослабления возмущающего действия атмосферы на работу оптической системы является применение адаптивных методов и систем. Идеи, положенные в основу создания адаптивных систем, предложены сравнительно недавно [1].

При компенсации нестационарных волновых искажений, которые возникают при распространении излучения в оптически неоднородной среде, в качестве измерительного устройства в адаптивных оптических системах волнового сопряжения применяются датчики волнового фронта [3, 4, 5, 6, 7]. Они являются ключевыми элементами многих современных систем управления и коррекции излучения. При

создании адаптивных оптических систем волнового сопряжения проводят, как правило [5, 6, 7], косвенные измерения волнового распределения на апертуре адаптивной оптической системы. Затем эти измерения одним из численных методов пересчитывают в базис функций отклика гибкого зеркала. При этом, каждому из известных алгоритмов, как правило, присущи свои достоинства и недостатки, однако конечные характеристики разрабатываемых на их основе датчиков волнового фронта зависят от конкретной технической реализации последних. В связи с этим интерес к дальнейшему совершенствованию таких устройств не ослабевает. Хорошо известный гартмановский тест [1] предложенный первоначально для контроля телескопической оптики, был в дальнейшем использован для адаптивной оптики и является наиболее часто применяемым типом датчика волнового фронта. Изображение входного зрачка проецируется на матрицу линз. Все изображения формируются на одном фотоприемнике, обычно ПЗС-матрице. Когда приходящий волновой фронт плоский, все изображения расположены в правильной сетке, определенной геометрией матрицы линз. Как только волновой фронт искажается, изображения смещаются со своих номинальных положений.

Следует отметить, что универсальным разложением волнового фронта, отвечающим ряду условий оптимальности, является разложение Карунена—Лоэва [1]. Оно характеризуется следующими свойствами, обусловливающими его оптимальность: минимальной среднеквадратической ошибкой при удержании заданно -го числа членов в бесконечном ряде разложения, получением наибольшего по сравнению с любым другим разложением количества информации о представляемой усеченным рядом функции, какое бы число членов ряда ни удерживалось, а также некоррелированностью коэффициентов разложения, что упрощает дальнейшее использование результатов разложения и их анализ. Однако в силу того, что аналитическое такое разложение трудно представимо, для практических целей обычно [1] используют систему полиномов Цернике, достаточно близких к ним. Для аппроксимации функций отклика гибкого адаптивного зеркала используем систему полиномов Цернике, ортогональных (ортонормированных) внутри единичной окружности или окружности радиусом Я, представленных в полярных координатах г, 0 [1, 2, 8]:

Vп + 1Япт(г)л/2со8тв,т Ф 0

уіп + 1Япт

(г), т = о

где

2п 1

/ / (г, 0)2  (г=

0 0 7 7 ]

ит( ) (п-т)'2 (-1)5(п - з)!гп-2

&п (г) = 2 —т-----------------т-т----------------г.

5=0 з ![п + т)/2 - з](п - т)/2 - з]!

Величины п и т всегда целые и удовлетворяют условию п<ш, п - |т| = чётно. Индекс ] является порядковым номером моды и зависит от п и ш. Условие ортогональности в круге единичного радиуса имеет вид:

где Ш (г) =

1/п їбе г ^ 1

5]’ - символ Кронекера.

[0 1 бе Г > 1

Волновой фронт, измеренный датчиком, в базисе Цернике имеет следующий вид:

N

Ф

'ер (г,в) = Д(г,в)й/

где 1 - коэффициенты разложения волнового фронта по полиномам Цернике (сигналы с выхода датчика волнового фронта), число полиномов N в разложении определяют с помощью выражения [2]:

N =

-0,2944(Б г0)3 1п Бґ

2_ 5 \/з

где Б - диаметр апертуры, г0 - радиус когерентности, & - число Штреля.

При этом в известных устройствах подлежат измерению с помощью квадрантных фотоприемников [3, 4, 5, 6, 7] локальные наклоны волнового фронта в точках

апертуры пропорциональные величинам вида

дФ

где

дх ду

I = 1, Ь ; 1 = 1, К ; Ь х К -количество квадрантных фотоприемников датчика,

Ф(х:, у 1) -значения фазы на апертуре датчика. Затем эти измерения одним из 1 1

численных методов пересчитывают в значения фазы Ф(х:, у :) или коэффициен-

1 У

а.

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

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

вида:

: дФ(гі А )еді

дв

и радиальные наклоны пропорциональные соответственно

величинам вида,

дФ (гі А )ері . —

дг

і = 1, М . Для измерения этих величин предлага-

ется использовать двухплощадные фотоприемники, расположенные в точках апертуры на концентрических окружностях, при этом граница раздела двухплощадных фотоприемников совпадает с радиусом соответствующей окружности или перпендикулярна ему (рис.1 а,Ь).

Рис. 1. Расположение фотоприемников на апертуре датчика волнового фронта: а - тангенциальный датчик, Ь - радиальный датчик.

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

дф (г I $1 )ед1

нальные локальным наклонам волнового

дв

в точках с координатами

Г, $ , при этом координаты точек г, $ , могут быть выбраны, в принципе, произвольно. Для определенности рассмотрим, например, расположение точек, приведенное на рис.1а.

Для реконструкции фазы применим метод наименьших квадратов. Соответствующая квадратичная форма метода наименьших квадратов в этом случае будет иметь следующий вид:

J1 = 2

М ( дФ(г в )ер дФ(г в) ^

і=0

дв

дв

(

У

дФ

дв

дв

Т

где М - число точек измерений волнового фронта (фотоприемников).

Выражения для истинных значений градиента фазы можно представить в следующем виде:

дФ (г ,ві) д N

N Ж]

= — 2 а!^і (Гі ,ві) =2

дв дв ] =1 і і у =1 1 дв

Для вычисления в явном виде значений выражения (6) следует использовать выражение (1). Подставив (6) в (5) получим:

М - 2 і=0

Ґ

дФ(гі ,ві )ерІ дв

N д2і

2 Я: --------------

і=1 J дв

\ Ґ

У

дФ

----------------2 ак -------------

ч дв к =1 к дв у

Т

Значения коэффициентов а] найдем из М линейных уравнений, приравняв к нулю частные производные квадратичной формы Л по а]:

В • А = С

дJ1

да-

= 0

где В - матрица с коэффициентами

ь = М к,І і=0 дв дв

А - вектор-строка искомых коэффициентов полиномов Цернике а];

й = м »

п л - 7 1=0 дв дв

С - вектор-столбец правой части ' ^ ,

к, 7 = 1, N.

Решение системы (8) будет иметь следующий вид:

А = В-1 • С

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

Рассуждая аналогично, можно построить датчик волнового фронта на основе радиальных производных (рис. 1Ь).

Соответствующая квадратичная форма метода наименьших квадратов в этом случае будет иметь следующий вид:

м {дФ(гі,ві)ёф (г-,ві)^ (дФ(гі,ві)ёф

дг дг

дг дг

где М- число точек измерений волнового фронта.

Значения коэффициентов а] найдем из М линейных уравнений, приравняв к нулю частные производные квадратичной формы 12 по а] :

дЗ2

= 0

да  в • А = С

7 ; 1 ,

ь = М д21

к, 7 , о дг дг

где В - матрица с коэффициентами 1 _0 ;

А - вектор-строка искомых коэффициентов полиномов Цернике а];

п = м дЩ в )еу д17(г,в) С - вектор-столбец правой части ^ г=0 дг дг ,

к, 7 = 1, N.

Решение системы (11) будет иметь следующий вид:

А = в-1 • С

1

в—

Таким образом, для заданной конфигурации точек матрицы В-1 или 1 могут быть рассчитаны заранее, а алгоритм вычисления коэффициентов разложения в базисе Цернике вектора А сводится к вычислению вектора правой части С и мат-

В 1 в-

ричному умножению на матрицу В-1 или 1 .

Рис.2. Функциональная схема датчика тангенциального типа: 1-линза, 2— матрица двухплощадных фотоприемников, 3—дифференциальные усилители, 4— масштабные усилители блока суммирования, 5—сумматоры блока суммирования, 6—группа масштабных усилителей, 7—вторая группа сумматоров.

На базе рассмотренного метода восстановления волнового фронта может быть построена структура соответствующего датчика. Данная структура приведена на рис. 2 и реализует алгоритм следующим образом. Искаженный волновой фронт фокусируется матрицей линз 1 на матрицу двухплощадных фотоприемников 2. При этом в каждом локальном участке, ограниченном апертурой линзы, сфокусированное пятно смещается относительно оптической оси в зависимости от локального наклона волнового фронта, пропорционального величинам вида дФ(Гі ,ві)ё?1

дв

в точках с координатами

гг в

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

дф(г1 А )ё#

. Затем сигналы с выхода дифференциальных усилителей 3 поступают на блоки суммирования, где вычисляются сигналы, пропорциональные эле-

с /

ментам вектор-столбца правой части .

При этом коэффициенты усиления Ку масштабных усилителей 4 блока суммирования вычисляются в соответствии с выражением

Е6і =

дг,

дв

Выходные сигналы блоков суммирования 5, пропорциональные величине С) поступают на входы группы масштабных усилителей, коэффициенты усиления которых вычисляются следующим образом:

к, = ь;1-

KJ k,J .

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

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

При технической реализации датчика волнового фронта, при произвольном выборе мест расположения фотоприемников оказывается, что весовые коэффициенты Ку имеют небольшие значения по модулю. Это приводит к тому что в этих каналах отношение сигнал/шум значительно меньше чем в других, что обуславливает существенное ухудшения технических характеристик датчика волнового фронта. Исследования показали, что величина данных весовых коэффициентов зависит от расположения фотоприемников. Оптимальное расположение фотоприемников следует выбирать с учетом максимума следующего критерия. Для тангенциального датчика волнового фронта (рис.4):

N

J = max 2 —------------

3 7=1 дв

V 7 ,

для радиального датчика волнового фронта:

N (dZ7 (r,в) ^2

J4 = max 2 —-----------

4 7=1 dr

v 7 .

Рассмотрим пример оптимизации топологии тангенциального датчика волнового фронта.

По формуле (4) определяем число полиномов Цернике. Так, например Б=0,05м, L=5 км, максимальное значение M в соответствии с (4) не превышает 12^16. В нашем случае примем N=16. Для выбора координат оптимального расположения двухплощадных фотоприемников, расположенных в точках апертуры на концентрических окружностях радиуса r используем критерий (15). На рис.4. точки a, b, c, d являются оптимизированными полярными координатами расположения фотоприемников датчика на окружности r1=0,2 и точки e, f, g, h, i, j, k, l, m, n, которые являются координатами расположения фотоприемников датчика на окружности r4=0,8. Далее аналогично выбираются оптимизированные координаты для расположения квадрантных фотоприемников на окружностях r2=0,4, r3=0,6, r5=1. В итоге мы получаем оптимизированную топологию тангенциального датчика волнового фронта.

,е 1 9 • • % И Л • в 1 • ) \ к 1 Л т Л 4 п І

• •  • і ‘ • • • % • • * * • а 1 • < • ' • * . С

# %. % і 1 1 • / % • 1 • * і і • 1 *./

а ь с а

Рис. 3. Выбор оптимального расположения фотоприемников датчика.

Рассуждая аналогично, можно оптимизировать топологию радиального датчика волнового фронта используя критерий (16).

Линза

Рис. 4 . Оптимизированное расположение тангенциальных датчиков на апертуре.

При проведении оценки эффективности адаптивных оптических систем на фоне гауссовских шумов следует учитывать расположение датчиков в точках апертуры на концентрических окружностях датчиков (рис. 1 а, Ь) и при оптимизированном расположении датчиков (рис. 4).

Также была исследована зависимость дисперсии ошибки восстановления волнового фронта Dв от дисперсии шумов Dш в каналах датчика для топологии рис.1 a, Ь и оптимизированная топология рис.4 при N=16. Анализ результатов показал, что предложенный вариант оптимизации топологии датчика волнового фронта позволяет существенно повысить точность восстановления волнового фронта.

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

дФ(х,у) дФ(rsini,rcosd) дФ . „ дФ ---^----------------’-----L =---sini +--cos i

дг дг дх ду

дФ( х, у) дФ( r sini, r cosi) дФ дФ

------------ =-------------------------- =------r cos i----------r cos i

дв дr дх ду

В выражениях (17), (18) учтена связь декартовой и полярной системы координат.

Выводы

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

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Воронцов М.А., Шмальгаузен В.И. Принципы адаптивной оптики. Москва: Наука, 1985. 336 с.

2. Лукьянов Д.П., Корниенко А.А., Рудницкий Б.Е. Оптические адаптивные системы. / Под ред. Д.П. Лукьянова. М.: "Радио и связь", 1989

3. Безуглов Д.А., Мищенко Е.Н. Метод сплайн-аппроксимации в задаче восстановления волнового фронта.// Известия РАН. Серия физическая, 12, 1992, с.156-160.

4. Безуглов Д.А., Мищенко Е.Н., Мищенко С.Е. Адаптивные оптические системы. Методы восстановления волнового фронта, разработка структур систем и новой элементарной базы. //Оптика атмосферы и океана, 3, 1995 (обзор), с.364-380.

5. А.с. 1647496 СССР, МКИ5 G 02 B 27/00. Датчик волнового фронта / Безуглов Д.А., Мищенко Е.Н, Тюриков В.Л. Опубл. в БИ. 1991. №17

6. Дегтярев Г.Л., Маханько А.В., Чернявский С.М., Чернявский А.С. Модальный датчик волнового фронта.// «Оптика атмосферы и океана», 15, № 12 2002г. с.1078-1083

7. А.с. 1720051 СССР, МКИ5 G 02 B 26/06. Датчик волнового фронта / Безуглов Д.А., Мищенко Е.Н, Крымский М.И., Серпенинов О.В. Опубл. в БИ. 1992. №.10.

8. Noll R.J. Zemike polynomials and atmospheric turbulence. J. Opt. Soc. Amer. 1976. V66. P. 207-211.

УДК 681.3.067:621.396.2

Д.М. Голубчиков

СТРУКТУРА И ПРИНЦИПЫ ФУНКЦИОНИРОВАНИЯ СИСТЕМЫ КВАНТОВОГО РАСПРЕДЕЛЕНИЯ КЛЮЧЕЙ ID 3000 CLAVIS

В мире существует всего три компании, предлагающие системы квантовой криптографии, предназначенные для применения в коммерческих приложениях. Две находятся в Европе - Id Quantique и Smart Quantum, одна в США - MagiQ Technology. Производимое ими оборудование строится по схожим схемам и незначительно отличается по своим характеристикам. В данном обзоре будет детально рассмотрена система квантового распределения ключей Id 3000 Clavis, производства компании Id Quantique. Она предназначена для проведения исследований в области квантовой криптографии и предоставляет пользователю широкие возможности по настройке и оценке параметров и характеристик квантового канала и оборудования предназначенного для его формирования[1].

Вначале будут кратко рассмотрены принципы функционирования системы. Далее проведено детальное описание блоков и модулей, входящих в ее состав.

Описание системы квантового распределения ключей

Квантовая система распределения ключей Id 3000 Clavis состоит из двух устройств размещенных в 19-ти дюймовых корпусах и пакета программного обеспечения для управления устройствами. Первое устройство является приемопередающим модулем и носит кодовое название QKDS-B или Bob, второе устройство является кодирующим модулем, не содержит приемо-передающей аппаратуры, работающей с квантовыми состояниями, и носит кодовое название QKDS-A или Alice[2].

Рис. 1. Внешний вид системы квантового распределения ключей Id 3000 Clavis.