акустооптический способ измерения смещений

Формула изобретения

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

Описание изобретения к патенту

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

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

Известен метод оптического зондирования поверхностных акустических волн (ПАВ) с опорной дифракционной решеткой, который может быть использован для контроля поперечных смещений объектов /1/. Метод, реализованный, например, на основе акустооптического (АО) измерителя координаты оптического пучка, заключается в использовании АО шкалы, выполненной на основе подложки с возбудителем ПАВ и опорной фазовой дифракционной решетки. Последняя установлена вблизи подложки на расстоянии , где - длина волны света, _пав - длина ПАВ. Для возбуждения ПАВ подводится сигнал с несущей частотой f_A, амплитудно-модулированный другим сигналом с частотой f_м.

Фаза сигнала изменяется при перемещении оптического пучка 1 относительно АО шкалы на величину -2 1/ _м, где _м - период огибающей модулированного сигнала на поверхности звукопровода, 1 - смещение пучка вдоль шкалы.

В устройстве с АО шкалой с периодом 100 мкм разрешающая способность достигала 2÷5 мкм.

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

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

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

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

Данный способ повышает разрешающую способность в направлениях, поперечных (ортогональных) лазерному лучу до значений _узв/1214 /4/, что при использовании в качестве светозвукопровода на основе воды (при f 8 МГц и _узв 200 мкм) составила значение 0,2 мкм. Как видно, этот способ не позволяет достичь разрешающей способности в направлениях, сонаправленных (коллинеарных) лазерному лучу до /1040 0,6 нм /5/ из-за того, что отношение _узв/ может достигать значений 10-350.

Также к недостаткам способа-прототипа следует отнести ограничения диапазона измерений размерами фотоприемника ( 1-2 мм) и функциональных возможностей из-за неоднозначности измерений при определении положения луча.

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

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

Сущность способа поясняется следующими чертежами.

Фиг.1 - блок-схема устройства измерения смещений, реализующего предлагаемый способ;

Фиг.2 - схема взаимодействия света и ультразвука в АО ячейке;

Фиг.3 - временные диаграммы работы устройства измерения смещений.

Устройство измерения смещений для реализации заявленного способа (Фиг.1) может включать: импульсный лазер 1, коллиматор 2, триппель-призма 3, закрепленная на каретке 4, интерферометр Маха-Цендера 6, состоящий из светоделителей 7 и 12, оптической схемы 8, отражателя 9, АО ячейки 10, диафрагмы 13, фотоприемника 14, импульсный интерполятор 15, схема синхронизации 16, состоящая из генератора 17, делителей частоты 18 и 21, формирователей импульсов 22 и 24, фаэовращателя 19, линии задержки 20, электронного ключа 23.

В момент светового импульса Т_св излучение импульсного лазера 1 преобразуется коллиматором 2 в коллимированный пучок и направляется на триппель-призму 3, которая закреплена на каретке 4, перемещающейся по контролируемой поверхности объекта 5. При перемещении каретки 4 вдоль оси Х триппель-призма 3 и отраженный от нее лазерный луч совершают вертикальные смещения вдоль оси Y 1_у, соответствующие отклонениям формы объекта 5. Отраженное от триппель-призмы 3 излучение освещает интерферометр Маха-Цендера 6, в котором светоделителем 7 оно делится на два пучка. Первый пучок (луч, поток) формирует опорное плечо 1 _оп интерферометра, следуя по пути: светоделитель 7 оптическая схема 8 отражатель 9 светоделитель 12 (пространственное совмещение с дифракционным порядком Е(+1)) диафрагма 13 фотоприемник 14.

Оптическая схема 8 представляет собой последовательное оптическое соединение коллиматора и фокусирующей линзы. Коллиматор уменьшает весь диапазон поперечных смещений луча в m раз , а фокусирующая линза преобразует новый диапазон поперечных смещений луча в диапазон угловых отклонений пучка ± на входе фотоприемника 14. Так как угловые отклонения пучка ± при пространственном совмещении с Е(+1) пучком приводят к изменениям периода интерференционной картины _ик= /sin , то расстояние от оптической схемы 8 до фотоприемника 14 и, тем более, длина опорного плеча 1_оп должны быть такими, чтобы выполнялось условие _д< _ик/2,5 /6/, где _д - диаметр диафрагмы. С учетом оптической схемы (Фиг.1) и, принимая, что для малых углов tg sin , правомерно записать условие 1_оп>2,5 L_{у д}/2m . Тогда при =0,6 мкм, m=6, Lу=30 мм и _д=0,1 мм получается 1_оп>1 м. Проведенный расчет позволил получить условие, при котором на выходе фотоприемника 14 возможно получить сигнал с высоким соотношением сигнал/шум /6/.

Второй оптический луч после светоделителя 7 освещает АО ячейку 10 под углом падения _п. В каждом цикле через АО ячейку 10 пробегает ультразвуковой импульс 11 (движущаяся периодическая структура в прозрачной среде), являющийся как-бы бегущим «окном дифракции» для входного излучения. Если в момент светового импульса Т _св положения луча, освещающего АО ячейку 10, и «окна дифракции» пространственно пересекутся, то небольшой, локальный участок пучка продифрагирует и сформирует два порядка дифракции: первый Е(+1) и нулевой Е(0). При этом нулевой порядок дифракции Е(0) не изменяет первоначального направления движения, а первый порядок дифракции Е(+1) отклоняется на угол дифракции а _д и получает частотное приращение равное несущей частоте f₁. Значения _п и _д зависят от параметров режима дифракции, выбранного в АО ячейке 10. При этом началом сформированного первого порядка дифракции Е(+1) является центр АО взаимодействия, в качестве которого можно принять середину бегущего ультразвукового импульса 11.

Таким образом, при пространственном пересечении второго, оптического луча, образуемого после светоделителя 7, и «окна дифракции» излучение Е(+1) следует по следующему пути: «окно дифракции» АО ячейки светоделитель 12 диафрагма 13 фотоприемник 14. Весь путь излучения от светоделителя 7, через АО ячейку 10 и далее до фотоприемника 14 образует измерительное плечо 1_изм интерферометра Маха-Цендера 6 (Фиг.2).

В результате пространственного совмещения с малым углом и интерференции импульсных разночастотных оптических потоков опорного 1_оп и измерительного 1_изм плеч интерферометра Маха-Цендера 6 на выходе фотоприемника 14 формируется импульсный частотный измерительный сигнал U_фп(t) rect(t/Т_св)·sin(2 f₁t+( о), где _о=2 (1_оп-1_из)/ , - начальная разность фаз в сбалансированном интерферометре Маха-Цендера при 1_у=0. Далее сформированный сигнал U_фп (t) следует на первый вход импульсного интерполятора 15.

Формирование электрических сигналов в устройстве происходит следующим образом (Фиг.3). Генератор 17 формирует электрический сигнал U₁(t)=U_1msin2 f₁t, поступающий на первый вход фазовращателя 19 и вход делителя частоты 18. Делители частоты 18 и 21 используются в устройстве для деления частоты входных сигналов на коэффициент деления n, который одинаков для обоих делителей, зависит от параметров устройства и для малых углов падения _п и дифракции _д может находиться в диапазоне [20; 40]. Делитель частоты 18 создает на своем выходе сигнал U₂(t)=U _2msin2 f₂t, где f₂=f_i/n, поступающий на второй вход импульсного интерполятора 15 и вход формирователя импульсов 24. Последний создает на выходе сигнал в виде коротких электрических импульсов U₈(t)=rect(t/T_CB ), приводящих к генерации импульсным лазером 1 импульсов света Т_св.

На второй вход фазовращателя 19 от импульсного интерполятора 15 подается код N_др, управляющий фазовым сдвигом (временной задержкой) сигнала U ₁(t). Этот блок позволяет осуществлять плавное дробное изменение фазового сдвига ( _др)в пределах от 0 до 2 л - ( _др=N_др· ), где - минимальная дискрета фазового сдвига. Выходной сигнал фазовращателя 19 U₃(t)=U_3msin(2 f₁t+( _др) поступает на первый вход электронного коммутатора 23 и на вход линии задержки 20.

Линия задержки 20 представляет собой цепочку последовательно включенных элементов задержки, коммутируемых цифровым кодом N_ц, приходящим от импульсного интерполятора 15. При этом минимальная дискрета задержки равна длительности периода частоты f1-T_min =1/f₁, а полное время задержки - T_max=N _ц·T_min=N_ц/f1. Соответственно, вводимый задержкой фазовый сдвиг пропорционален целым периодам 2 - _ц=2 N_u. Последовательное включение фазовращателя 19 и линии задержки 20 позволяет изменять в широких пределах общий суммарный фазовый сдвиг (временную задержку) сигнала U ₄(t)=U_4msin(2 f₁t+ ), где = _ц+ _др=2 N_u+N_др· =N_вых , где N_u и N_др - цифровые коды составляющие код N_вых.

Сигнал U ₄(t) приводит к созданию на выходе делителя частоты 21 импульсного сигнала U₅(t), с частотой f₂ =f₁/n, где n - коэффициент деления. Затем импульсный сигнал U₅ (t) проходит через формирователь импульсов 21 и преобразуется в сигнал U₆(t), состоящий из электрических импульсов. Этот сигнал следует на второй вход электронного коммутатора 23, в результате чего на его выходе формируется сигнал U ₇(t) rect(t/Т_узв)(sin(2 f₁t+ ) в виде электрических импульсов Т_узв, заполненных короткими импульсами несущей частоты f₁. Импульсы Т_узв, поступая на АО ячейку 10, возбуждают в ней соответствующие бегущие ультразвуковые импульсы (движущиеся периодические структуры) 11, синхронизированные со световыми импульсами Т_св . В описании принимается, что в процессе измерений несущая частота неизменна f_t=const, а длительности световых и ультразвуковых импульсов находятся в соотношении: Т_узв 4Т_св.

В результате пространственного совмещения и интерференции импульсных разночастотных оптических потоков, как описано ранее, на выходе фотоприемника 14 формируется импульсный сигнал U_фп(t), который является измерительным сигналом и, как сигнал обратной связи для всей схемы, следует на первый вход импульсного интерполятора 15. На второй вход последнего приходит импульсный сигнал U₁(t) от генератора 17, который является опорным, а сигнал U₂(t), поступающий на его третий вход, является стробирующим сигналом «запуска измерений».

По приходу стробирующего импульса U₂(t) импульсный интерполятор 15 начинает преобразование разности фаз между сигналами U_фп(t) и U₁ (t) в цифровой код N_вых, старшие разряды которого являются кодом N_ц, а младшие - N_др. Цифровой код N_вых является выходным и одновременно используется в качестве сигналов обратной связи N_u и N_др . В зависимости от значений двух кодовых компонент N_ц и N_др в фазовращателе 19 и линии задержки 20 вводятся пропорциональные временные задержки (фазовые сдвиги), которые, изменяя в совокупности суммарный фазовый сдвиг , смещают «окно дифракции» и компенсируют фазовый сдвиг разбаланса интерферометра, возникающий из-за поперечных смещений луча 1у.

Описание работы устройства составлено с учетом использования режима дифракции Брэгга. Однако допустимо использовать дифкрацию света и в режиме Рамана-Ната. Так же принимается, что при малой длительности световых импульсов Т_св смещения бегущего ультразвукового импульса 11 минимальны и не приводят к значительному «размазыванию», расширению первого порядка дифракции Е(+1).

Значение коэффициентов деления n делителей частоты 18 и 21, находящееся в диапазоне от 20 до 40, выбрано из расчета малых углов падения _п и дифракции _д. Однако, в зависимости от различных параметров устройства, в первую очередь от используемой АО ячейки 10, значение n может отличаться от выбранных ранее.

На Фиг.1 и Фиг.2 изображены схемы интерферометра Маха-Цендера 6, в которых поток опорного плеча 1_оп, формируемый светоделителем 7, пространственно совмещается и интерферирует с первым порядком дифракции Е(+1). Несмотря на это, в предлагаемом способе возможно использование в качестве опорного потока нулевого пучка дифракции Е(0).

В заявленном способе используется режим синхронизации освещения (импульсами света Т_св) и возбуждения (импульсами ультразвука Т_узв) АО ячейки. За счет стробоскопического эффекта для каждого импульса света смещения ультразвукового импульса минимальны, поэтому принимается, что в каждом цикле освещения формируется неподвижное «окно дифракции» с собственной небольшой локальной областью и центром АО взаимодействия (фиг.2). Поперечные смещения луча 1_у относительно «окна дифракции» приводят к смещению центра АО взаимодействия и, в итоге, к сонаправленным (коллинеарным) смещениям оптического потока 1" и 1' в измерительном плече интерферометра. Появившийся сигнал разбаланса в виде дополнительного фазового сдвига световых волн _св преобразуется в компенсирующее фазовое рассогласование режима синхронизации, смещающее «окно дифракции», чтобы уравновесить изменения _св (т.е. реализуется измерительный метод уравновешивания):

_св= .

Выражение для _св ((св записывается в виде:

,

где и - фазовые сдвиги, возникающие от изменения оптических путей 1" и 1', соответственно, 1' - изменения оптического пути до центра АО взаимодействия (от светоделителя 7 до бегущего ультразвукового импульса 11), 1" - изменения оптического пути после центра АО взаимодействия (от бегущего ультразвукового импульса 11 до светоделителя 12).

Исходя из приведенных формул и с учетом оптической схемы (Фиг.2), получаем:

N_вых= / =2 ( 1" и 1')/ · =2 ( 1_у[sec( _п+ _д)-ctg( _п+ _д]/ · ,

где _п+ _д - углы падения и дифракции в АО ячейке, соответственно, - минимальная дискрета компенсирующего фазового рассогласования .

Как видно из данной формулы, полученная зависимость реализует «абсолютный отсчет» для поперечных смещений луча _у, а полный диапазон фазового рассогласования и кода N_вых соответствуют смещениям луча в пределах диапазона L_у. Этим расширяются функциональные возможности.

При использовании АО ячейки с _п+ _д=45° формула преобразуется к виду: N _вых=2 1_у/ · , из которой для одного кванта (N_вых=1), равного =1°, разрешающая способность составляет ~1,7 нм. Как видно, по разрешающей способности предлагаемый способ более чем 100 раз превышает прототип. Дальнейшее уменьшение суммы _п+ _д содержит резерв повышения разрешающей способности.

Как следует из алгоритма предлагаемого способа первый порядок дифракции Е(+1), формируемый в АО ячейке, пространственно стабилизирован при любом положении освещающего пучка. Следовательно, весь диапазон поперечных смещений луча 1_у возрастает от размеров фотоприемника 1-2 мм (как в прототипе) до размеров, определяемых выражением (Фиг.2):

L_у=h_аом-d_п-l_аом tg _п,

где h_аом и l _аом - высота и ширина АО ячейки, соответственно, d _n - диаметр оптического пучка.

Исходя из приведенных уравнений для увеличения диапазона измерений L_у и повышения разрешающей способности необходимо подбирать такой режим дифракции в АО ячейке, чтобы уменьшать _п+ _д. Как следует из /7/, наилучшими для этих требований подходят АО ячейки на основе одно- и двуосных кристаллов, работающих в режиме анизотропной дифракции. Для для таких АО ячеек /8/ при малых углах _п и _д диапазон измерений L_у может достичь 10÷30 мм, что в 5-15 раз больше, чем у прототипа.

Согласно алгоритму работы устройства значение Т_св не должно быть больше 1-2 тактов несущей частоты f₁-Т_св<(1-2)/f ₁. При использовании АО ячейки на основе кристалла парателлурита (Те0₂), как в /8/, с f₁=25 МГц, v_узв 600 м/с и l_аом (1 см) временные параметры измерений будут следующими:

Т_св<40÷80 нс, цикл измерения (для малых углов _п, _д (п, (д) Т_ц>l_аомtg· _п/v_узв 1,5·10^-6 с. Полученные значения приемлемы при использовании наиболее распространенных импульсных лазеров, например, как в /9/.

Максимальная «поперечная» скорость v_{y max} при условии двух измерений за период (длину волны К) по критерию Найквиста, с учетом ранее определенных параметров, равна v_{y max}= /2Т_ц 0,2 м/с. Как видно, полученное значение соразмерно максимальным «продольным» скоростям v_{x max}=0,2-0,3 м/с, например, как у интерферометров НР-5525А (5526А) или ИПЛ-30 (30К1, 30К2) /6/.

Как следует из алгоритма работы устройства время измерений импульсного интерполятора не должно превышать длительность светового импульса - Т_св<40÷80 нс. Практическая реализация такого импульсного интерполятора возможна при введении режима синхронизации, стробирования в конструкциях интерполяторов, широко используемых в интерферометрах с фазовым выходом, например, как в /10/.

Введение фазового сдвига может быть осуществлено, например, за счет использования цифрового формирователя фазового сдвига с положительным или с отрицательным знаком, в ручном или автоматическом режимах. Для текущего уровня развития схемотехники фазометрии возможно для 8-ми управляющих двоичных разрядов достичь минимального дискретного значения 1,41°. Для сдвига фазы сигнала в предлагаемом устройстве также допустимо совместное использование цифроаналогового преобразователя с управляемым напряжением фазовращателями, например, Merrimac модели PSFM-4-30 с точностью 1°÷5° /11/ или применение систем фазовой автоподстройки частоты.

Предлагаемый акустооптический способ измерения смещений позволяет повысить точность измерений перемещений объектов в направлении, ортогональном когерентному световому потоку (лазерному лучу), и довести разрешающую способность до величины ~1-5 нм, а также расширить диапазон измеряемых поперечных смещений до ~10-30 мм, что позволяет обойтись без применения альтернативных методов измерений для повышения достоверности измерений. При этом реализацию предлагаемого способа можно осуществить на базе устройств, известных из уровня техники.

Источники информации

1. Комоцкий В.А., Котюков М.В. Акустооптический измеритель координаты оптического пучка. // Автометрия. - 1991, - № 5. - с.110-113 (аналог).

2. А. с. № 1610252 МПК G01B 11/00. Способ измерения пространственных перемещений объекта // Телешевский В.И., Яковлев Н А., Игнатов С.А. Опубл. в БИ № 44, 1990 (аналог).

3. А.с. № 1765691 МПК G01B 21/00. Способ измерения смещений объекта. Телешевский В.И., Яковлев Н.А. // Опубл. в БИ 1992, № 36 (прототип).

4. Яковлев Н.А. Построение лазерных систем для измерения перемещения по трем координатам на основе акустооптического преобразования измерительной информации. Автореф. дисс. на соиск. уч.степ, к.т.н. 05.11.16. М., 1991.

5. Игнатов С.А. Повышение разрешающей способности лазерных измерительных систем для контроля оборудования ГПС методом акустооптоэлектронной обработки информации. Автореф. дисс. на соиск. уч.степ, к.т.н. 05.11.16. М., 1987.

6. Иванов В.А. Привалов В.Е. Использование лазеров в приборах точной механики. СПб, Политехника, 1993.

7. Балакший В.И., Мартынова М.В., Румянцев А.А. Дифракция света на акустическом импульсе. / Юптика и спектроскопия, 1998, т.84, № 5, с.860-866.

8. Балакший В.И., Казарьян А.В. Акустооптическая стабилизация направления лазерного пучка. // Квантовая электроника, 1998, т.25, № 11, с.988-992.

9. Акустооптический сканер на основе бегущей акустической ЛЧМ-линзы // Вовк Ю.М., Затолокин В Н., Рудаков И.Б. и др. Автометрия, № 1, с.54-62.

10. Патент РФ № 2016381 МПК G01B 21/00. Способ автоматической интерполяции порядка интерференции и устройство для его осуществления // Михальченко Е.П., Рюмин А.В., Яковлев Н.А. Опубл. в БИ № 13, 1994.

11. Эрлих М., Филипс Л., Вагнер Дж. Управляемый напряжением акустооптический фазовращатель / Приборы для научных исследований. Русский перевод "Review Scientific Instrument" 1988, № 11, с.57-59.

Ehrlich M.J., Philips L.C., Wagner J.W. Voltage-controlled acousto-optic shifter / 7 Review Scientific Instrument 59, 1988, No.11, pp.2390-2392.