Конспект лекций по курсу Лазерные измерительные системы (8 семестр)
Содержание
Характеристика современных ЛИС
Физические основы лазерных измерений
Простейшая математическая модель лазерного пучка (ЛП)
Измерения поляризации
Амплитудные измерения
Измерение поперечных размеров пучка
Измерения частоты
Рефрактометрия
ТЕНЕВОЙ МЕТОД
Измерения фазы
Локация и дальнометрия
Схема и принцип действия одного из типов оптического локатора для слежения за авиационными и космическими объектами
Эффект Доплера в оптике
ДОПЛЕР Христиан (30. XI 1803-17.III 1853) австрийский физик, математик и астроном, член Австрийской АН Родился в Зальцбурге.
Лазерная доплеровская анемометрия (ЛДА)
Эффект Доплера в лазерной анемометрии
Методы регистрации доплеровского сдвига частоты
Принцип фотосмешения
Основные типы оптических схем ЛДА
Геометрия волновых векторов в дифференциальной схеме ЛДА
Геометрия волновых векторов в инверсной дифференциальной схеме ЛДА
ФУНКЦИОНАЛЬНАЯ СХЕМА ЛДА
Требования к лазерам
Требования к СФЗИ
Требования к ПОС
Требования к рассеивающим частицам
Требования к ФП и процессору сигнала
Выполним сначала расчет сигнальной частоты по эффекту Доплера с учетом начального сдвига частот зондирующих пучков
Интерференционная модель работы ЛДА
Основные выводы:
Пространственные условия согласования оптических волн при фотосмешении
Эффективность фотосмешения
Поляризационное согласование
Амплитудное согласование
Амплитудное согласование двух цилиндрических пучков
Фазовое согласование
Фазовое согласование при фотосмешении двух плоских волн
Прямоугольная апертура
Формирование сигнала на разностной частоте при фотосмешении как результат нестационарной интерференции сигнальной и опорной
Устройства сдвига частоты лазерных пучков (Однополосные частотные модуляторы)
Устройство акустооптического частотного модулятора
Режимы работы АОМ
Ультразвуковой модулятор Рамана-Ната
К расчету частот дифрагированных пучков для ячейки Рамана-Ната на основе эффекта Доплера
Распределение мощности по дифрагированным пучкам разного порядка
Оценки эффективности использования излучения
Пример. Характеристики УЗМ Рамана-Ната, используемого в лабораторном практикуме лазерных измерительных систем
Ультразвуковой модулятор Брэгга
К расчету частоты дифрагированного пучка для ячейки Брэгга на основе эффекта Доплера
Двухчастотная схема с опорным пучком
Двухчастотная дифференциальная схема
Двухчастотная инверсно-дифференциальная схема
Вращающаяся фазовая или амплитудная решетка радиального типа
Вращающаяся фазовая или амплитудная решетка барабанного типа
Механический модулятор для сверхмалых сдвигов частоты
Частотный модулятор с вращающимся поляризатором
Принцип работы поляризационного модулятора
Рассеяние плоской монохроматической волны на сферической частице
Соотношение между комплексными амплитудами компонент падающего и рассеянного поля
Формулы Ми (G.Mie, 1908 г)
Основные положения теории Ми
Угловое распределение интенсивности рассеянного излучения
Дифференциальное сечение рассеяния
Индикатрисы рассеяния Ми
Таблица соответствия x и d=2r для разных l
Свойства индикатрис рассеяния Ми
Индикатрисы для больших частиц
Зависимости дифференциальных сечений от x для различных углов рассеяния
Рассеяние Рэлея (x<<1)
Рассеяние Рэлея
Рассеяние Рэлея-Ганса (|m-1|<<1; 2x|m-1|<<1)
Рассеяние Рэлея-Ганса
Приближение больших частиц (x>>1)
Рассеяние двух лазерных пучков микрочастицей
Расчет интенсивности рассеянного света
Фазовые характеристики рассеяния двух пучков
Взаимное дифференциальное сечение рассеяния
Двухпучковые индикатрисы рассеяния
Расчет мощности оптического сигнала ЛДА по дифференциальной схеме
Зависимость интегрального коэффициента согласования от угловой апертуры ПОС
Энергетический расчет ЛДА
Расчет интенсивности излучения на оси гауссова пучка
Расчет амплитудного значения мощности полезного оптического сигнала
Соотношения, использованные при выводе формулы для Psg0
Определение Y1,Y2
Учет потерь на оптических элементах
Потери при распространении излучения в рассеивающей среде
Обеспечение высокой эффективности фотосмешения рассеянных волн
Пример расчета сигнальной мощности
Расчет числа фотоотсчетов за пролет частицы
3090 отсчетов за пролет
82 отсчета за пролет
18 отсчетов за пролет
Применение ЛИС для измерения размеров микрочастиц
Метод спектральной прозрачности
Статистический метод флуктуаций
Броуновское уширение спектра рассеянного излучения
Измерение модуля дифференциального сечения рассеяния и глубины модуляции фототока в ЛДА
Фазовый метод в ЛДА
Пространственное разрешение ЛДА
Пробный объем ЛДА
К расчету весовой функции пробного объема
Весовая функция пробного объема
Влияние ПОС на пространственное разрешение ЛДА
Оптимизация диаметра диафрагмы поля зрения
Учет глубины резкости объектива ПОС
Работа ПОС при внеосевом приеме
Особенности ПОС в схемах с рассеянием назад
Ход лучей в ЛДА на обратном рассеянии со сдвигом частоты
Принципы построения многокомпонентных схем ЛДА
Трехпучковая двухкомпонентная схема для измерения двух ортогональных проекций скорости
Четырехпучковая двухкомпонентная схема для измерения двух ортогональных проекций скорости
Способы разделения полезных сигналов и подавления паразитных
Поляризационное разделение
Разделение по цвету (по спектральным линиям)
Частотное разделение
Частотное разделение
Примеры использования УЗМ
Временное (коммутационное) разделение
Разделение каналов по когерентности
Схемы для измерения осевой компоненты
б) схема с вычитанием сигналов:
Четырехпучковая схема для определения трех ортогональных проекций скорости
Временные и спектральные характеристики сигнала ЛДА
Частица пролетает вблизи центра измерительного объема
Частица пролетает на границе измерительного объема
Спектральный анализ сигнала ЛДА
Относительное уширение сигнальной составляющей
Для ЛДА с гауссовыми пучками и неограниченным полем зрения ПОС:
Специализированные процессоры сигнала ЛДА
Процессор сигнала ЛДА счетно-импульсного типа (каунтер)
Работа формирователя стробов
Принцип измерения длительности строба
Принцип проверки результатов
Факторы, определяющие погрешность оценки частоты.
Пример оценки погрешности, связанной с конечным временным разрешением таймера
Влияние шума на работу каунтера
Оптимальные условия работы каунтера
Шумы в ЛИС Отношение сигнал/шум
Расчет отношения сигнал/шум
Три основные составляющие шума на выходе фотоприемника
Приведение шума к эквивалентной оптической мощности на входе фотоприемника входу
Формула для расчета отношения сигнал/шум
При расчете (S/N)0:
Типичные характеристики используемых в ЛДА фотоприемников
Примечания
Основы аппаратурного спектрального анализа
Оценка энергетического спектра
Оценка энергетического спектра
Условия правильной работы последовательной схемы спектрального анализа
Основы аппаратурного спектрального анализа
Cредства аппаратурного спектрального анализа
Оценка энергетического спектра
Оценка энергетического спектра
Условия правильной работы последовательной схемы спектрального анализа
Влияние конечной ширины полосы УПФ
Возможные записи ДС при неправильном выборе параметров анализа
Дисперсионно-временные устройства спектрального анализа (одноканальные, параллельные, широкий диапазон 104-109 Гц)
Акустооптические анализаторы спектра
Панорамный анализатор спектра с компьютерной поддержкой
Работа ГУН в режиме ГКЧ
Три режима работы детектора и ОУ
Параллельный многоканальный спектроанализатор
Процессор сигнала ЛДА следящего типа
Принцип работы следящей системы
Графики сигналов
Графики сигналов в ЧД при f=f0
Особенности следящей системы
Тестирование следящей системы c помощью генератора ЧМ-сигналов
Особенности работы следящей системы
Парадоксы измерения спектра турбулентных пульсаций с помощью ЛДА со следящей системой
Оценка погрешности фазового шума
Результирующая погрешность
Результирующая погрешность
Погрешность калибровки
Электронные корреляторы
Линии задержки
Умножители
Осциллографический метод измерения двумерной ПВ
Цифровые методы
Оценка шума квантования АЦП
Фотонный коррелятор
Одночастичный сигнал ЛДА
Принцип действия простейшего фотонного коррелятора
Блок-схема однобитового фотонного коррелятора
Использование оптического спектрального анализа
Использование интерферометра Фабри-Перо, сканируемого давлением воздуха, в ЛДА с дифференциальной оптической схемой
Оптико-механическая сканирующая СФЗИ-ПОС (ЛДА для исследований сверхзвуковых потоков в аэродинамических трубах)
Система обработки сигнала
Некоторые технические характеристики ЛДА
Использование интерферометра Фабри-Перо для частотной демодуляции.
Использование интерферометра Фабри-Перо для визуализации поля скорости
Обобщенная структурная схема ЛИС
Классификация погрешностей
Классификация погрешностей ЛИС по характеру
Фундаментальная система уравнений измерений
Основы параметрической теории информации
Статистка погрешностей
Информационная энтропия
При стандартных условиях:
В случае, когда точность измерений меняется случайным образом в процессе эксперимента можно отбросить условие 4) и положить:
Свойства оценок
Несмещенность
Эффективность
Неравенство информации Рао-Крамера
Обобщение на случай измерения нескольких параметров
Пример анализа свойств оценок
Для определения эффективности используем неравенство информации
Оценка интенсивности турбулентных пульсаций
Эмпирическая дисперсия
Квантовые флуктуации лазерного излучения Распределение Пуассона
Задача Бернулли
Применение схемы Бернулли к фотоэффекту.
Интегральная интенсивность
Учет классических флуктуаций интенсивности
Формула Манделя
Избыточный шум
Параметр вырождения
Распределение числа фотоотсчетов для многомодового лазерного излучения
Сравнение распределений Пуассона (□) и Бозе-Эйнштейна (○) при <k> = 5
ПРИМЕРЫ ОЦЕНОК ПРЕДЕЛЬНОЙ ТОЧНОСТИ ЛАЗЕРНЫХ ИЗМЕРИТЕЛЬНЫХ СИСТЕМ
Квантовые пределы точности лазерной анемометрии
Квантовые пределы точности лазерной интерферометрии.
Методы получения оценок
Метод максимального правдоподобия (ММП)
Применение ММП при нормальном распределении погрешностей
Применение ММП при лапласовском распределении погрешностей
Байесовские оценки. Метод максимума апостериорной вероятности.
ДОВЕРИТЕЛЬНАЯ ОЦЕНКА РЕЗУЛЬТАТА СЕРИИ ОДНОРОДНЫХ ИЗМЕРЕНИИ
ПОГРЕШНОСТЬ КОСВЕННЫХ ИЗМЕРЕНИЙ
МНК ДЛЯ ЛИНЕЙНОЙ МОДЕЛИ
Проверка гипотез
Правила выбора решения
Исключение выбросов
Таблица критических значений отношения t для браковки “выскакивающих” значений (n – число приемлемых результатов, P –
Проверки нормальности
Анализ погрешностей лазерного интерферометра (Пример расчета точностных характеристик современной ЛИС)
Оптическая локация
Схема и принцип действия одного из типов оптического локатора для слежения за авиационными и космическими объектами
Области применения и характеристики
Отражательные характеристики целей
Эффективная поверхность рассеяния (ЭПР)
Дальность действия (L)
б) Статистические характеристики
Ослабление лазерного излучения в атмосфере
Рефракция и турбулентность атмосферы
Элементная база ИК-локаторов
Литература
АДАПТИВНАЯ ОПТИКА
Литература
7.69M
Категория: ФизикаФизика

Конспект лекций по курсу "Лазерные измерительные системы"

1. Конспект лекций по курсу Лазерные измерительные системы (8 семестр)

Лектор - профессор, доктор физико-математических наук
Смирнов Владимир Иванович
Учебный план
Лекции – 16 часа.
Практические занятия – 34 часов
Лабораторные работы – 16 часов.
Курсовой проект (56 часoв).
Защита курсового проекта.
Зачет по лабораторным работам.
• Экзамен.
1

2. Содержание


Лекция 1. Характеристика современных ЛИС. Физические основы
лазерных измерений. Эффект Доплера в лазерной анемометрии.
Основные типы оптических схем ЛДА. (5-40) Функциональная схема и
основные элементы ЛДА. Интерференционная модель ЛДА.
Пространственные условия согласования оптических волн при
фотосмешении. (41-64)
Лекция 2. Акустооптические частотные модуляторы Рамана-Ната.
Акустооптические частотные модуляторы Брегга. Модуляторы с
дифракционными и поляризационными устройствами. (65-85)
Рассеяние света. Рассеяние Рэлея. Рассеяние Рэлея-Ганса.
Рассеяние плоской световой волны сферической микрочастицей.
Теория Ми. (86-104)
Лекция 3. Рассеяние двух лазерных пучков микрочастицей.
Энергетический расчет ЛДА. Применения ЛИС для измерения
размеров микрочастиц. (105-131) Пространственное разрешение ЛДА.
Многокомпонентные ЛДА. Способы селекции сигналов. (132-159)
2

3.


Лекция 4. Временная и спектральная структура сигнала ЛДА.
Специализированный счетно-импульсный процессор ЛДА. Шумы в
ЛИС. Отношение сигнал/шум. (160-192) Основы аппаратурного
спектрального анализа. Дисперсионные и акустооптические
анализаторы спектра. Анализатор спектра последовательного типа с
компьютерной поддержкой. (193-214)
Лекция 5. Процессор ЛДА следящего типа. Корреляционные и
фотонкорреляционные методы обработки сигналов ЛИС. (215-239)
ЛДА со сканируемым интерферометром Фабри-Перо. Использование
интерферометра Фабри-Перо для частотной демодуляции.
Использование интерферометра Фабри-Перо для визуализации поля
скорости. (240-248)
Лекция 6. Анализ обобщенной структурной схемы ЛИС.
Классификация погрешностей. Фундаментальная система уравнений
измерений. Статистика погрешностей. (249-266) Свойства оценок.
Неравенство информации. Минимальная дисперсионная матричная
граница оценок параметров сигналов. Примеры анализа свойств
оценок на основе неравенства информации. (267-279)
3

4.


Лекция 7. Квантовые флуктуации лазерного излучения. Распределение
Пуассона. Статистика лазерного излучения. Формула Манделя.
Избыточный шум. Квантовые пределы точности лазерной анемометрии.
Квантовые пределы точности лазерной интерферометрии. (280-308)
Методы получения оценок. Методы МП, МАВ, НК, НК с взвешиванием,
НМ. Проверка гипотез. Матрица потерь. Исключение выбросов. Проверка
нормальности распределения. Доверительная оценка результата серии
однородных измерений. Погрешность косвенных измерений. Метод
наименьших квадратов для линейной модели. (309-330)
Лекция 8. Анализ погрешностей промышленных лазерных
интерферометров. (331-340) Лазерные локаторы и дальномеры.
Адаптивные лазерные измерительные системы. (341-263)
Литература. (364)
4

5. Характеристика современных ЛИС

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

6.

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

7. Физические основы лазерных измерений

• В основе работы ЛИС – решение обратной
задачи взаимодействия лазерного излучения
с изучаемой физической системой.
• По известным параметрам излучения на
входе системы и измеренным параметрам
излучения на выходе определяются
параметры системы, воздействующей на
излучение.
7

8. Простейшая математическая модель лазерного пучка (ЛП)

x2 y 2
E ( x, y, z; t ) e A( z ) exp 2
w ( z)
(t t0 ) 2
exp
exp i t kz ( x, y, z )
2
• При взаимодействии с физической системой могут меняться
следующие характеристики лазерного пучка:
• поляризация,
• амплитуда и связанные с ней интенсивность и мощность,
• поперечные размеры лазерного пучка,
• направление распространения,
• частота,
• фаза,
• время задержки.
8

9. Измерения поляризации

• Изменение состояния поляризации лазерного пучка может быть
обусловлено следующими физическими процессами:
• а) преломлением и отражением ЛП на границе раздела
двух сред;
• б) распространением пучка в анизотропных средах;
• в) рассеянием света в оптически неоднородной среде.
• Наиболее распространенными оптическими методами
исследования, в которых информативным параметром служит
поляризация, являются методы эллипсометрии и
поляриметрии, в частности поляризационно-оптический метод
исследования деформаций твердых тел.
• Методы эллипсометрии широко используются в производстве
полупроводниковых приборов для определения свойств
поверхности, показателя преломления и толщины тонких
пленок, нанесенных на образец.
9

10.

10

11. Амплитудные измерения

• Изменение амплитуды лазерного пучка и связанной с
ней мощности может быть обусловлено:
• а) преломлением и отражением на границе раздела
двух сред;
• б) поглощением световой энергии в среде;
• в) перераспределением плотности мощности
лазерного пучка за счет интерференции и
дифракции;
• г) рассеянием лазерного пучка в среде.
• На практике нашли распространение методы
фотометрии и радиометрии, методы скоростного
фотографирования, где информативными
параметрами являются плотность мощности и
мощность лазерного пучка.
11

12. Измерение поперечных размеров пучка

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

13. Измерения частоты

• Изменение частоты лазерного пучка может
быть обусловлено:
• а) распространением в нестационарных
средах;
• б) распространением в движущихся средах
(эффекты Доплера, Саньяка).
• На практике широкое распространение
получили методы лазерной доплеровской
анемометрии (ЛДА) и лазерной гироскопии,
которые применяются для измерения
линейных и угловых скоростей движения.
13

14.

14

15.

Внешний вид 3-компонентного ЛДА с аргоновым лазером
и с частотным сдвигом, работающего на обратном рассеянии
(проспект фирмы DISA)
15

16.

Схема 5-пучкового 3-компонентного ЛДА с аргоновым лазером
и с частотным сдвигом, работающего на обратном рассеянии
(проспект фирмы DISA)
1 – крышка и ослабитель; 2 – регулировка положения перетяжки; 3 – нейтральный делитель
пучков; 4 – блок модулятора Брега; 5 – модифицированный делитель пучков;
6 – устройство смещения; 7 – приемный блок обратного рассеяния; 8 – держатель;
9 – блок приемной оптики; 10 – модифицированный поляризационный делитель;
11 – интерференционный фильтр 514,5 нм; 12 - интерференционный светофильтр 488,0 нм;
13 – блок фотоумножителя; 14 – диафрагма поля зрения;
15 – устройство преобразования пучков; 16 – расширитель пучков.
16

17. Рефрактометрия

• Изменение направления распространения лазерного
пучка вызывается:
• а) отражением и преломлением на границе раздела
двух сред;
• б) распространением в оптически неоднородных
средах.
• На практике получили широкое распространение
теневые и рефрактометрические методы.
• Теневые методы широко используются в газовой
динамике при исследовании обтекания тел
различной формы дозвуковым и сверхзвуковым
газовым потоком.
17

18.

Рефрактометр ИРФ-470
Производитель: Казанский
оптико-мехалический завод
(КОМЗ)
Цена опт.: 4 900 р.
Цена: 5 200 р.
Портативный рефрактометр ИРФ-470 предназначен для экспрессного анализа состава и качества
пищевых, химических, нефтехимических продуктов, сырья, плодов, ягод, биологических проб, сточных
вод, сахаристости свеклы.
Показатель преломления 1,3 - 1,51 %
Диапазон измерений по сахарозе от 0 до 30 %
Специальное устройство с нанесенной шкалой и справочным устройством позволяет быстро определить:
процентное или массовое содержание сахара
содержание сухих веществ по сахарозе ("Brix")
концентрацию белка в сыворотке, крови и воде
содержание алкоголя, солей
долю сухих веществ в напитках, молоке, вине
качество меда, масел, лекарств
спелость плодов, ягод
сорт жидкого топлива
чистоту сточных вод
Рефрактометр используется на предприятиях пищевой, химической, фармацевтической,
нефтехимической промышленности, сельском хозяйстве. Прибор может использоваться в полевых
18
условиях

19. ТЕНЕВОЙ МЕТОД

• -метод обнаружения оптич. неоднородностей в
прозрачных преломляющих средах и дефектов
отражающих поверхностей (напр., зеркал). Впервые
предложен в 1857 Л. Фуко (L. Foucault) для
отражающих поверхностей. В 1867 А. Теплером (A.
Toepier) этот метод был усовершенствован при
исследовании прозрачных преломляющих сред. Т. м.
наз. также шлирен-методом (от нем. Schliere - оптич.
неоднородность, свиль, шлир).
• Т. м. применяют при исследованиях распределения
плотности воздушных потоков, образующихся при
обтекании моделей в аэродинамических трубах,
используют для проекции на экран изображений
(получаемых в виде оптич. неоднородностей) в
пузырьковых камерах, в телевиз. системах проекции
на большой экран и др.,
19

20.

• В Т. м. пучок лучей от точечного или щелевого источника света
1 (рис.) линзой или системой линз и зеркал (2-2') направляется
через исследуемый объект (3 )и фокусируется на непрозрачной
преграде (5) с острой кромкой (на т. <н. н о ж е Ф у к о), так что
изображение источника проецируется на самом краю преграды.
• Если в исследуемом объекте нет оптич. неоднородностей, то
все идущие от него лучи задерживаются преградой.
• При наличии оптич. неоднородности (4 ) лучи будут
рассеиваться ею и часть их, отклонившись, пройдёт выше
преграды.
• Поставив за ней проекционный объектив (6) или окуляр, можно
на экране (7) получить изображение неоднородностей (8 )или
наблюдать их визуально.
20

21.


- фотография обтекания крыла самолета: а - теневой метод; б - метод
Теплера; в - интерференционный метод с использованием
интерферометра Маха - Цендера; 1 - источник евета, 2 - исследуемая
область течения; 3 - экран; 4 - линза, 5 - нож Фуко; 6 21
полупрозрачные зеркала, 7 - непрозрачные зеркала; 8 - компенсатор.

22.

22

23. Измерения фазы

• Изменение фазы лазерного пучка может быть
обусловлено следующими физическими процессами:
• а) изменением показателя преломления среды,
через которую распространяется пучок;
• б) изменением геометрического пути
распространения лазерного пучка.
• Изменение показателя преломления может быть
вызвано изменением температуры, давления,
плотности, химического состава, движением среды
(эффект Физо).
• На практике широкое распространение получили
различные интерференционные методы, в том числе
и голографические.
23

24.

Лазерный интерферометр ЛабоИнтерфер-1
Компактный лазерный интерферометр – это современный
высокоточный оптический метрологический инструмент. Он
предназначен для оценки и измерения параметров
отклонения формы плоских и сферических поверхностей
оптических деталей (небольшого диаметра).
Его преимущества: портативность и виброустойчивость.
Прибор может быть оснащён системами регистрации и
отображения информации различных конфигураций,
оригинальным программным обеспечением (по
согласованию с заказчиком).
Интерферометр может использоваться в метрологических
отделах, лабораториях и оптических производствах, а
также в учебных, научных и исследовательских целях.
·
·
·
·
·
Отличительные особенности:
высокая точность измерений
простота в эксплуатации
широкая гамма проекционных насадок (объективов)
наличие анти вибрационной платформы
небольшие линейные размеры
24

25. Локация и дальнометрия

• Изменение времени задержки лазерного
импульса может быть обусловлено
изменением оптического пути его
распространения.
• На практике нашли широкое
распространение методы лазерной локации и
дальнометрии.
25

26. Схема и принцип действия одного из типов оптического локатора для слежения за авиационными и космическими объектами

26

27. Эффект Доплера в оптике

• Эффект изменения частоты электромагнитной волны,
регистрируемой приемником, при движении источника и
приемника относительно друг друга.
• Особую известность эффект Доплера получил в связи с
открытием космологического красного смещения (эффект
Хаббла). Оно свидетельствует о расширении Вселенной.
ki
П
И
u
Используя преобразования Лоренца,
можно получить:
0 1 u / c
1 (u / c) cos
u c
0
0 u cos
c
0 kiu
2
2
27

28.

ki
а)
П
И
ki
б)
u
П
И
u
i kiu
i kiu
Источник и приемник
удаляются друг от друга
Источник и приемник
сближаются
28

29. ДОПЛЕР Христиан (30. XI 1803-17.III 1853) австрийский физик, математик и астроном, член Австрийской АН Родился в Зальцбурге.

• Профессор Венского университета и директор первого
в мире Физического института, организованного по его
инициативе.
• Физические работы в области оптики и акустики.
• В 1842 теоретически обосновал зависимость частоты
звуковых и световых колебаний, воспринимаемой
наблюдателем, от скорости движения наблюдателя и
источника колебаний (эффект Доплера).
• Исследования посвящены также аберрации света,
29
теории микроскопа, теории цветов.

30. Лазерная доплеровская анемометрия (ЛДА)

• - направление современной прикладной оптики,
связанное с разработкой когерентно-оптических
методов измерения скорости движения жидкостей,
газов и твердых тел.
• Эти методы основаны на регистрации доплеровского
сдвига частоты в свете, рассеянном исследуемой
средой.
• Наряду с термином ЛДА (LDA) в научной литературе
используется и другая терминология: ОДИС, ЛДИС
(LDV).
• Термин «анемометрия» образован от греческих слов:
anemos – ветер и metreo – измеряю.
30

31. Эффект Доплера в лазерной анемометрии


Особенность в том, что источник и приемник неподвижны, а движется
оптическая неоднородность среды, рассеивающая излучение по всем
направлениям.
u
И
u c
Ksi
2
Р
i
' i u cos 1 i kiu
c
ks
1
Z’ П
Ksi
ki
Z
(1)
'
s
s ' u cos 2 ' u cos 2 ' ksu
c
c
(2)
31

32.

D s i ksu kiu (ks ki )u;
(3)
u c
u
2 /
a/2
a/2
2 /
4 a
K si sin ; (7)
2
2
ks k ' ki , i ;
(ks ki ) K si ; (5)
Ksi
D K siu
(4)
(6)
4 a
D sin u cos
2
(8)
32

33.

D 2 a
D
sin u cos ;
2 2
uK
D
(9)
(10)
uK – проекция скорости на разностный волновой вектор Ksi;
– период виртуальной ИК при интерференции двух плоских волн с
волновыми векторами ks и ki.
2
;
2 sin( a / 2) K si
u K u cos ;
(11)
(12)
33

34.

D max
2uK
при a , ks ki ,
2
D min
0 при a 0, ks ki ,
u 1 м/с
(13)
(14)
34

35. Методы регистрации доплеровского сдвига частоты

• Прямой оптический спектральный анализ имеет
ограниченную область применения в связи с относительно
малыми сдвигами частоты. Разрешающая способность
оптических спектральных приборов не превышает 106-107.
• Основным, используемым на практике методом, является
метод фотосмешения.
• При фотосмешении на фотоприемник направляются две волны
и выделяется составляющая фототока на разностной частоте.
• В зависимости от того, какие волны направляются на
фотоприемник, различают три типа оптических схем ЛДА: схема
с опорным пучком, дифференциальная и обратная
дифференциальная схемы.
35

36. Принцип фотосмешения

E1 A1 sin( 1t )
E2 A2 sin( 2t )
2 1
2 1
0
2
E E1 E2 A1 sin( 1t ) A2 sin( 2t )
2
2
T
0
2 1
i (t ) B E 2
T
B E12
T
B E22
T
2 B E1E2
T
36

37.

1
sin( 1t ) sin( 2t ) [cos( ) cos( 2 0t )]
2
i (t ) BA12 sin 2 ( 1t )
BA1 A2 cos( t )
T
T
BA22 sin 2 ( 2t )
BA1 A2 cos( 2 0t )
T
T
B 2
i (t ) ( A1 A22 ) BA1 A2 cos[( 2 1 )t ]
2
37

38. Основные типы оптических схем ЛДА

• Простейшие реализации
основных оптических
схем ЛДА:
• а)-в) с опорным пучком,
• г) дифференциальная,
• д) инверснодифференциальная.
• 1,2 – объективы,
• 3 – фотоприемник,
• 4 – плоскопараллельная
пластина.
38

39. Геометрия волновых векторов в дифференциальной схеме ЛДА

D1 (ks k1)u
D 2 (ks k2 )u
sg D D 2 D1
sg Ku
K k1 k2
Частота сигнала не зависит от направления наблюдения. Можно
собирать рассеянное излучение в большом телесном угле, не
вызывая дополнительного уширения спектра сигнала.
39

40. Геометрия волновых векторов в инверсной дифференциальной схеме ЛДА

D1 (ks1 ki )u
D 2 (ks 2 ki )u
sg D D 2 D1
sg Ku
K ks1 ks 2
Частота сигнала не зависит от направления зондирующего пучка.
40

41. ФУНКЦИОНАЛЬНАЯ СХЕМА ЛДА

• Лазерный доплеровский анемометр (ЛДА) состоит из следующих
основных элементов:
• источника излучения - лазера (Л),
• системы формирования зондирующего излучения (СФЗИ),
• приемной оптической системы (ПОС), фотоприемника (ФП),
• электронной системы обработки сигнала - процессора ЛДА (П).
41

42. Требования к лазерам


В ЛДА используются в основном газовые лазеры непрерывного
действия, обладающие высокой пространственной и временной
когерентностью и имеющие наиболее стабильные характеристики.
При хорошем рассеянии можно использовать относительно дешевые
гелий-неоновые лазеры с излучением красного цвета мощностью 5 -50
мВт.
Если рассеяние слабое, то необходимы аргоновые лазеры с
излучением зеленого и синего цвета, мощностью в несколько ватт.
В последние годы широко используются полупроводниковые лазеры,
имеющие малые габариты, низковольтные блоки питания. Их
недостатками является малая длина когерентности, невысокая
стабильность и малый срок эксплуатации из-за деградации активной
среды. Как следствие, полупроводниковые лазеры не могут
использоваться для прецизионных измерений.
42

43. Требования к СФЗИ

• СФЗИ предназначена для формирования
зондирующих пучков, пересекающихся под
определенными углами в точке измерения и
имеющих заданные размеры, форму
пространственного распределения поля,
поляризацию и частотные сдвиги относительно
исходного излучения.
• Для измерения одной компоненты вектора скорости
СФЗИ формирует два лазерных пучка.
• Для измерения двух и трех компонент необходимы
многопучковые схемы ЛДА.
• СФЗИ может содержать однополосные модуляторы,
вращатели поляризации, интерференционные
фильтры и другие оптические элементы.
43

44. Требования к ПОС


Приемная оптическая система (ПОС), как правило, содержит
высококачественный объектив, собирающий рассеянное из области измерений
излучение и направляющий это излучение на фотоприемник.
В схеме прямого рассеяния ПОС располагается вдоль биссектрисы угла
пересечения пучков (оси СФЗИ) вперед по ходу лучей.
Если приемная система расположена вдоль биссектрисы, но в направлении
назад относительно падающих лучей, схема называется схемой обратного
рассеяния.
В ЛДА с обратным рассеянием СФЗИ и ПОС обычно удается совместить в
едином приемно-передающем блоке. При этом некоторые оптические элементы
выполняют двойные функции. Так, фокусирующая линза СФЗИ используется
для приема рассеянного излучения.
Существуют и другие варианты схем, в которых рассеянное излучение
собирается под некоторым углом к оси СФЗИ, - схемы внеосевого прямого и
обратного рассеяния.
В последние годы в СФЗИ и ПОС ЛДА все чаще используются световолоконные
элементы и световодные кабели для гибкого соединения
легких оптических головок с источником и приемником излучения, что
позволяет проще решать проблему сканирования точки измерения по потоку.
44

45. Требования к рассеивающим частицам

• Рассеяние лазерного излучения происходит на дискретных
неоднородностях потока - микрочастицах, естественным
образом присутствующих в потоке в виде инородных
включений, либо специально введенных в поток. Рассеяние на
этих частицах для большинства приложений ЛДА значительно
превышает рассеяние на флуктуациях показателя преломления
среды потока.
• Метрологические характеристики ЛДА существенно зависят от
оптических и гидродинамических параметров рассеивателей.
Поскольку параметры естественных неоднородностей трудно
поддаются контролю, целесообразно, когда это возможно,
очищать поток и вводить в него специально подобранные
частицы, например, с помощью генератора частиц.
• Однако в двухфазном потоке, содержащем, например, капли
воды в воздушной струе или пузырьки воздуха в воде,
рассеивающие частицы - неотъемлемая часть потока. Эти
частицы могут иметь разную форму и размеры.
45

46. Требования к ФП и процессору сигнала

• В качестве фотоприемника в ЛДА применяется ФЭУ,
редко -фотодиод, если мощность оптического
сигнала достаточно большая. Основное требование к
фотоприемнику - высокая чувствительность на
рабочей длине волны и быстродействие.
• Электронная система обработки сигнала (или
процессор) ЛДА предназначена для получения
оценок скорости и параметров турбулентности
потока в точке измерений на основании анализа
флуктуации фототока, возникающих в результате
детектирования рассеянного излучения.
46

47. Выполним сначала расчет сигнальной частоты по эффекту Доплера с учетом начального сдвига частот зондирующих пучков

Интерференционная модель работы ЛДА
Выполним сначала расчет сигнальной частоты по эффекту
Доплера с учетом начального сдвига частот зондирующих пучков
sg s1 s 2
D1 s1 1 (ks k1)u
D 2 s 2 2 (ks k2 )u
sg 12 Ku
12 1 2
K k1 k2
47

48.


Из полученных формул следует, что частота сигнала не зависит от
направления наблюдения и определяется составляющей скорости
вдоль разностного вектора K , поэтому K называют вектором
чувствительности ЛДА.
Для определенности дальнейшего частотного анализа сигналов
будем считать, что 1> 2. Это означает, что индексом 1 отмечается
тот из пучков, частота которого больше.
Начальный частотный сдвиг создается в ЛДА с помощью
однополосных модуляторов с частотами модуляции до нескольких
десятков МГц. Разность оптических частот по крайней мере на 5-6
порядков меньше 1 2.
48

49.

Направляем ось Х вдоль вектора K:
sg 12 Ku
f sg
ux u K
sg
2
ˆ
u cos(u , K )
ux
f sg f12
12
f12
2
2 sin( a / 2)
ux – проекция вектора скорости на вектор чувствительности ЛДА;
– период интерференционной картины в области пересечения
зондирующих пучков.
49

50. Интерференционная модель работы ЛДА


При интерференции волн с различными частотами образуется нестационарная
интерференционная картина с пространственным периодом .
Интерференционные полосы движутся вдоль X с фазовой скоростью u*= 12/K .
При движении рассеивающей частицы через такую интерференционную решетку в
результате пересечения темных и светлых полос происходит модуляция рассеянной
мощности, частота которой определяется скоростью движения частицы относительно
решетки.
Интерференционный подход дает объяснения всех особенностей графика зависимости
50
fsg(ux).

51. Основные выводы:

• ЛДА, использующий лазерные пучки с одинаковой частотой
излучения, не позволяет определять знак проекции вектора
скорости.
• Создание частотного сдвига зондирующих пучков во многих случаях
позволяет исключить такую неопределенность знака.
• Период ИК задает масштабный коэффициент перехода от частоты к
скорости в ЛДА.
• Период ИК , выраженный в мкм, численно дает величину изменения
скорости в м/с при изменении частоты сигнала на 1 МГц.
u x
f sg
51

52. Пространственные условия согласования оптических волн при фотосмешении

• Анализируем случай
фотосмешения линейно
поляризованных
когерентных волн
(сигнальной и опорной).
Es
(x,y)
S
Er
Es ( x, y; t ) es As ( x, y ) exp i st s ( x, y )
Er ( x, y; t ) er Ar ( x, y ) exp i r t r ( x, y )
r s
r s
2
52

53.

Суммарное поле на светочувствительной поверхности фотоприемника
E ( x, y; t ) Es ( x, y; t ) Er ( x, y; t )
*
* *
* * * *
EE ( Es Er )( Es Er ) Es Es Er Er Es Er Es Er
Фототок фотоприемника
*
e 1 0
i M
E ( x, y; t ) E ( x, y; t )dxdy
h 2 0 S
i
*
is ir irs exp( i rst ) irs exp( i rst )
rs r s
53

54.

Составляющие фототока при фотосмешении
e 1 0
B M
h 2 0
is B As2 ( x, y )dxdy
- коэффициент преобразования
- фототок от сигнальной волны
S
ir
2
B Ar ( x, y )dxdy
S
- фототок от опорной волны
Интерференционная компонента фототока:
irs B As ( x, y ) Ar ( x, y )es er exp i s ( x, y ) r ( x, y ) dxdy
S
54

55. Эффективность фотосмешения

irs
isir
es er As ( x, y ) Ar ( x, y ) exp i s ( x, y ) r ( x, y ) dxdy
S
1
2
2
2
A
(
x
,
y
)
dxdy
A
r ( x, y )dxdy
s
S
S
exp( i )
arg( ) arg(irs )
55

56.

i is ir irs exp[ i ( rst )]
*
irs
exp[ i( rst )]
*
irs irs
isir
exp[ i ( rs t )] exp[ i ( rs t )] 2 cos( rs t )
Основное уравнение фотосмешения когерентных волн
i is ir 2 isir cos( rst )
| | - эффективность фотосмешения
0 1
56

57. Поляризационное согласование

As ( x, y ) Ar ( x, y );
s ( x, y ) r ( x, y )
| es er | es er | cos |
es er
или
es er 1
es er 0
/ 4 2 / 2 0,71
57

58. Амплитудное согласование

es er 1;
s ( x, y ) r ( x, y )
As ( x, y ) Ar ( x, y )dxdy
S
1
2
2
2
A
(
x
,
y
)
dxdy
A
r ( x, y )dxdy
s
S
S
As ( x, y) const Ar ( x, y) 1
58

59. Амплитудное согласование двух цилиндрических пучков

As ( x, y ) As0circ (r / rs )
rs
rr
Ar ( x, y )
rs
rr
, rs rr и , rs rr
rr
rs
1
0
Ar circ (r / rr )
| |
rs/rr
1
59

60. Фазовое согласование

es er 1;
As ( x, y ) Ar ( x, y ) const
Эффективность фотосмешения:
1
S
exp{ i[ s ( x, y) r ( x, y)]}dxdy
S
60

61. Фазовое согласование при фотосмешении двух плоских волн

Z
kr
a/2
Es ( x, y; t ) As exp i st ks r
Er ( x, y; t ) Ar exp i r t kr r
K
a/2
ks
X
2
2
K 2 sin( a / 2)
2 sin( a / 2)
s r ks r kr r (ks kr )r
s r Kr
2
s ( x, y ) r ( x, y ) Kx
x
61

62. Прямоугольная апертура

lx l y
2 2
Y
1
exp( iKx )dxdy
l xl y l x l y
ly/2
lx/2
lx/2
ly/2
X
2
2
lx
2
1
exp( iKx )dx
lx l x
2
62

63.

1
iKl x
iKl x
iKl x
exp 2 exp 2
sin( Kl x / 2)
sinc (lx / )
Kl x / 2
sin( x)
sinc ( x)
x
sinc (lx / )
63

64. Формирование сигнала на разностной частоте при фотосмешении как результат нестационарной интерференции сигнальной и опорной

волн
I ( x; t ) I s I r 2 I s I r cos( rst Kx)
2 sin( a / 2)
u ph
rs
K
u ph f rs
64

65. Устройства сдвига частоты лазерных пучков (Однополосные частотные модуляторы)

• В лазерных измерительных системах часто
возникает необходимость сдвига частоты исходного
лазерного излучения:
0 1;
1 0 m 0 .
• Широкое распространение получили устройства
сдвига частоты на основе дифракции света на
ультразвуковых волнах – акустооптические
частотные модуляторы (АОМ)
65

66. Устройство акустооптического частотного модулятора

3
1- возбудитель УЗВ (пластинка из
монокристалла ниобата лития,
пьезокерамика – используется
пьезоэлектрический эффект).
a
2 – звукопровод (жидкость, стекло,
кристалл, пластмасса с хорошей
фотоупругостью n=n(p)).
k0, 0
3 – поглотитель УЗВ (выбирается в
соответствие с материалом
звукопровода). Если не ставить
поглотитель, возникнет стоячая волна
=> амплитудная модуляция.
k1, 1
2
1
66

67. Режимы работы АОМ

• Различают два предельных режима работы АОМ в
зависимости от коэффициента объемности Q:
Q 2 l sin / a
• l – длина взаимодействия луча с акустическим столбом;
a – длина ультразвуковой волны.
• Q<<1 (Q<0,3) => дифракция на плоской фазовой решетке
– режим Рамана-Ната (ячейка Рамана-Ната).
• Q>>1 (Q>4) => дифракция на пространственной фазовой
решетке – режим Брэгга (ячейка Брэгга).
67

68. Ультразвуковой модулятор Рамана-Ната

va
a
k0, 0
a sin p p
2
1
1
0
-1
X
-2
u=u0cos( mt+ 0)
p arcsin p p ,
a
a
( a )
p 0 p m ; p 0; 1; 2,
Для работы можно выбирать любые пары пучков.
Например, +1 и 1 => разность частот равна 2 m;
+2 и 3 => разность частот равна 5 m.
68

69. К расчету частот дифрагированных пучков для ячейки Рамана-Ната на основе эффекта Доплера

k0, 0
va
a
kp, p
p
0
u=u0cos( mt+ 0)
p 0 D ;
D (k p k0 )v a ;
4
2
D sin( p / 2)v a pv a ;
p p / a ;
va
a m
fa
;
a
2 2
p 0 p m
69

70. Распределение мощности по дифрагированным пучкам разного порядка

• Акустическая волна в звукопроводе может рассматриваться как
бегущая фазовая синусоидальная решетка. “Волна” показателя
преломления имеет вид:
2
n n0 n cos( mt K a x), K a
.
a
• Из теории дифракции на фазовой синусоидальной решетке известно,
что мощность дифрагированных пучков зависит от индекса фазовой
модуляции G0 и определяется выражением
Ip
2
I 0 J p (G0 ),
2
G0 l n.
• Jp – Функция Бесселя 1-го рода порядка р.
• Практически можно варьировать индекс модуляции, меняя
амплитуду модулирующего напряжения u0 , так как n~ u0.
70

71. Оценки эффективности использования излучения

• Максимально в первый
порядок уходит 34%
мощности при Г0=1,8 =>
Если использовать ±1
порядки эффективность
составит 68%.
• Если использовать 0 и 1
порядки, то максимальная
эффективность составляет
60% при Г0=1,4.
• Нулевой порядок исчезает
при Г0=2,4 и Г0=5,5.
• Первые порядки исчезают
при Г0=3,83.
71

72. Пример. Характеристики УЗМ Рамана-Ната, используемого в лабораторном практикуме лазерных измерительных систем


Возбудитель ультразвуковой волны – пьезокерамика ЦТС-19. (Ø=12
мм, f1р=1,9/d; частота выражается в МГц, толщина керамики - в мм).
Обычно работают на 3-ей резонансной частоте (fm~5МГц).
Звукопровод – дистиллированная вода, va=1500 м/с.
v a 1500
a
300 мкм;
fm
5
0,63 мкм 1
0,002.
a
=> При 2w=1мм необходимо расстояние около 1м, чтобы
дифрагированные пучки разделились.
72

73. Ультразвуковой модулятор Брэгга

• В режиме дифракции Брэгга коэффициент объемности
достаточно большой и необходимо рассматривать дифракцию на
объемной фазовой решетке.
• Дифрагированная волна возникает только при определенном
угле падения, который определяется условием Вульфа-Брэгга:
1
2
1’
2’
C
a
E
F
D
ED EF
ED EF a sin
2 a sin Б
73

74. К расчету частоты дифрагированного пучка для ячейки Брэгга на основе эффекта Доплера

1 0 D (k1 k0 )v a
va
3
k1, 1
2
2
a
k0, 0
1
m
4
D sin Бv a
4
va
D
v a 2
2 a
a
1 0 m
74

75. Двухчастотная схема с опорным пучком

• 1 – лазер; 2,4,5 – линзы; 3 – модулятор; 6 –
фотоприемник.
• Модулятор является не только частотно сдвигающим
устройством, но и выполняет функции лучевого
расщепителя.
75

76. Двухчастотная дифференциальная схема

• 1 – лазер; 2,4,5 – линзы; 3 – модулятор; 6 – фотоприемник.
• Модулятор является не только частотно сдвигающим
устройством, но и выполняет функции лучевого расщепителя.
• Особенность состоит в независимости ширины
интерференционной полосы в зондирующем поле от длины
волны лазерного излучения.
76

77. Двухчастотная инверсно-дифференциальная схема

Двухчастотная инверснодифференциальная схема
• 1 – лазер; 2,3,5 – линзы; 4 – диафрагма; 6 – модулятор; 7 –
фотоприемник.
• Модулятор является не только частотно сдвигающим
устройством, но и рекомбинационным элементом.
77

78. Вращающаяся фазовая или амплитудная решетка радиального типа

1
;
p
1 0 D (k1 k0 )v p ;
2
2
1 0 1v p
вр R;
p
1 0 вр N
Астигматизм дифрагированных пучков. Широкие пучки не годятся, так
как p~R.
Углы отклонения не зависят от частоты модуляции.
Максимальные сдвиги ограничены механической прочностью.
78

79. Вращающаяся фазовая или амплитудная решетка барабанного типа

1
;
p
1 0 D (k1 k0 )v p ;
2
2
1 0 1v p
вр R;
p
1 0 вр N
• Нет астигматизма пучков, так как R=const.
• Углы отклонения не зависят от частоты модуляции.
• Максимальные сдвиги ограничены механической прочностью.
79

80. Механический модулятор для сверхмалых сдвигов частоты

1 0 D (k1 k0 )v ;
v
k0, 0
k1, 1
4
1 0 v ;
v
1 0
.
/2
80

81. Частотный модулятор с вращающимся поляризатором

• 1,3 – четвертьволновые пластинки;
• 2 – вращающийся поляризатор;
• 4 – поляризационная призма;
81

82. Принцип работы поляризационного модулятора

cos 2
sin cos
P
;
2
sin
sin cos
(1)
Ex
1
E P j ;
y
врt; (2)
(3)
Ex cos 2 j sin cos ;
(4)
E y sin cos j sin 2 ;
(5)
82

83.

1
1
Ex (1 cos 2 ) j sin 2
2
2
1 1
(cos 2 j sin 2 )
2 2
1 1
exp( j 2 );
2 2
1
1
E y sin 2 j (1 cos 2 )
2
2
1
1
j j (cos 2 j sin 2 )
2
2
1
1
j j exp( j 2 );
2
2
(6)
(7)
83

84.

Ex 1 1 1
1
E 2 j 2 exp( j 2 врt ) j ;
y
(8)
• На выходе модулятора – две волны с круговой
поляризацией с правой и левой.
• Одна из волн сдвинута по частоте на 2 врt.
• Знак сдвига зависит от направления вращения.
84

85.

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

86. Рассеяние плоской монохроматической волны на сферической частице


О – центр рассеяния; M – точка наблюдения;
OM=R; – угол рассеяния; – азимутальный угол;
P – плоскость рассеяния;
ki – волновой вектор падающей волны;
ks – волновой вектор рассеянной волны.
86

87.

• Представление электрического вектора
падающей волны вектором Джонса:
E ( z , t ) A
E ( z , t ) A exp[ i ( t kz )];
II
II
A a exp(ia ), AII aII exp(iaII ).
• Представление электрического вектора
рассеянной волны вектором Джонса в точке
наблюдения:
Es (t ) As
E (t ) A exp( i st );
sII sII
As as exp(ia s ), AsII asII exp(ia sII ).
87

88. Соотношение между комплексными амплитудами компонент падающего и рассеянного поля

As exp(ikR) A
A ikR S A ,
II
sII
[S] - матрица рассеяния, которая в выбранной системе
координат имеет диагональный вид:
S
S
0
0
,
SII
S┴ и SII –безразмерные комплексные амплитудные функции
рассеяния для составляющих, перпендикулярной и
параллельной плоскости рассеяния.
88

89. Формулы Ми (G.Mie, 1908 г)

2n 1
S S (m, x, )
(an n bn n );
n 1n( n 1)
2n 1
SII SII (m, x, )
(bn n an n );
n 1n( n 1)
an,bn – коэффициенты, рассчитываемые с помощью
рекуррентных соотношений;
n, n – выражаются через полиномы Лежандра
и их производные.
89

90. Основные положения теории Ми

• Амплитудные функции рассеяния S┴ и SII зависят только от
трех параметров: m, x, .
• m – комплексный показатель преломления частицы
относительно окружающей среды:
m m1 / m2 ,
• m1 – показатель преломления частицы;
• m2 – показатель преломления среды.
• x – относительный размер частицы:
x kr 2 r / 2 rm2 / 0 ,
• r – радиус частицы.
• – угол рассеяния.
90

91. Угловое распределение интенсивности рассеянного излучения

as2 as2 as2II ;
2
2
2
as
a
2
k R
*
S S ;
2
2
asII
a II
k 2R
*
S S ;
2 II II
a a sin , aII a cos ;
I
*
2
*
2
I s 2 2 ( S S sin SIISII cos );
k R
91

92. Дифференциальное сечение рассеяния

1
( , ) 2 ( S S * sin 2 SII SII* cos2 );
k
I s ( R, , ) 2
( , )
R
I
• Дифференциальное сечение рассеяния не зависит от
интенсивности падающей волны и расстояния до точки
наблюдения, имеет размерность площади;
• численно дает:
• интенсивность рассеянной волны на расстоянии 1 метр при
единичной интенсивности падающей волны,
• мощность, рассеянную в единичном телесном угле в
направлении наблюдения.
92

93. Индикатрисы рассеяния Ми

Рис. 3.2.
с. 21 [2].
Y Y ( , ) lg[k 2 ( , )]
93

94. Таблица соответствия x и d=2r для разных l

Таблица соответствия x и d=2r для
разных
d
x kr ;
,мкм\ x 0,6
d x .
0,7
3
5
10
0,63
0,12
0,14
0,6
1,0
2,0
0,5
0,1
0,11
0,48
0,8
1,6
94

95. Свойства индикатрис рассеяния Ми

• С увеличением размера частиц растет число интерференционных
максимумов;
• индикатрисы вытягиваются в направлении рассеяния вперед;
• для больших частиц (диаметром несколько десятков мкм и более)
наблюдается интенсивный дифракционный максимум нулевого
порядка;
• Количество максимумов (или лепестков) индикатрисы в
диапазоне углов 0< < составляет примерно x;
• Угловая полуширина переднего лепестка приблизительно равна
/ .
• Другие лепестки примерно в два раза уже за исключением
заднего лепестка, который может быть частично или полностью
вырожденным.
• Для малых x индикатрисы при =0 отличаются от индикатрис при
= /2 отсутствием рассеяния в области углов 90°, что
соответствует релеевскому характеру рассеяния.
95

96. Индикатрисы для больших частиц

• Для больших частиц (диаметром несколько десятков мкм и
более) количество лепестков в индикатрисе велико,
• под некоторыми углами могут наблюдаться интенсивные
интерференционные пики, обусловливающие при освещении
белым светом появление радуг.
• В приближении больших частиц (x>>1) пользуются законами
геометрической оптики.
Рис. 3.3.
с. 22 [2].
96

97. Зависимости дифференциальных сечений от x для различных углов рассеяния

• Отмечается сложный, немонотонный, осциллирующий характер
зависимостей.
• Для всех кривых характерна тенденция увеличения рассеяния с
увеличением параметра размера частицы x, введенного Ми, и
относительного показателя преломления.
Рис. 3.4.
с. 23 [2].
97

98.

Полное сечение рассеяния:
T ( , )d .
4
Интегральное сечение рассеяния:
( , )d .
Коэффициент рассеяния:
Qs T /( r 2 ).
• При х Qs→2.=>
• Рассеяние на
больших
частицах имеет
дифракционную и
рефракционную
составляющие.
Рис. 3.5.
с. 24 [2].
98

99. Рассеяние Рэлея (x<<1)

Рассеяние Рэлея
(x<<1)
2
(
m
1)
3
S ix
;
2
(m 2)
2
(
m
1)
3
SII ix cos 2
;
(m 2)
x (m 1)
( , ) 2 2
(sin 2 cos 2 cos 2 ).
k (m 2)
6
2
8 x (m 1)
Qs
.
2
2
3(m 2)
6
2
2
• При φ=π/2 интенсивность рассеянного света одинакова по всем
направлениям.
• При φ=0 под углом θ=π/2 рассеянный свет отсутствует,
рассеяние вперед и назад одинаково.
99

100. Рассеяние Рэлея

2
2
2
2
, , x
sin cos cos
2 2
k m 2
6
x m 1
,
, 0.06
2
90
1201.5 .10 22 60
22
1 .10
150
23
5 .10
180
0
30
, 0 , 0.06
0
210
270
2
180
0
300
330
240
log k ( 0 , 0 , 0.7) 1.545
YMie
30
0
210
330
240
90
1201.5 .10 22 60
22
1 .10
150
23
5 .10
270
300
1
100

101. Рассеяние Рэлея-Ганса (|m-1|<<1; 2x|m-1|<<1)

Рассеяние Рэлея-Ганса
(|m 1|<<1; 2x|m 1|<<1)
S ix (m 1) 2 / U J 3 / 2 (U )
3
3
3
SII ix (m 1) 2 / U J 3 / 2 (U ) cos
3
3
3
U 2 x sin( / 2)
J 3 / 2 (U )
функция Бесселя степени 3/2.
• По мере возрастания размера частицы увеличивается
интенсивность света рассеянного вперед по сравнению с
рассеянием назад.
• При больших x индикатриса рассеяния имеет много лепестков.
101

102. Рассеяние Рэлея-Ганса

U x , 2 x sin
2
m 1.1
2
3
n , x x m 1
6
120
90
2
U x ,
2
3
p , x x m 1
120
0
270
300
60
p , 1
0
180
0.01 0.03
330
240
90
3
30
0.01 0.03
210
U x ,
JJ x , cos
2
150
30
180
2
6
60
150
n , 1
3
JJ x ,
2
3
JJ x , Jn , U x ,
2
210
330
240
270
300
102
2

103. Приближение больших частиц (x>>1)

Приближение больших частиц
(x>>1)
Рассеянный большой частицей свет можно представить
состоящим из двух частей.
Одна часть (в приближении геометрической оптики) учитывает
отражение и преломление на сферической поверхности.
Наиболее полно изучены характеристики рассеянного каплями
воды света.
Индикатриса имеет многочисленные максимумы и минимумы.
Усредненная индикатриса вытянута вперед.
В основном свет сосредоточен внутри конуса, ограниченного
углом рассеяния 700. В этой области находится 93,4%
мощности рассеянного света.
103

104.

Другая часть рассеянного света определяется дифракцией на
шаре.
Амплитудные функции дифрагированного света имеют вид:
xJ1 ( x sin )
S ( ) SII ( )
,
sin
J1 ( x sin ) - функция Бесселя 1-го рода 1-го порядка.
Дифракционный язык индикатрисы сосредоточен в конусе с
углом
3,83 / x.
Например, для частицы диаметром 50 мкм при =0,63 мкм он
составляет примерно 50’.
104

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

Рис. 3.6.
с. 25 [2].
• Направление наблюдения характеризуется углами ψ и χ.
105

106. Расчет интенсивности рассеянного света

I s I s1 I s 2 2 I s1I s 2 d cos( sg t d );
• Примеры расчетов дифференциального
коэффициента поляризационного согласования
Рис. 3.7.
с. 28 [2].
106

107. Фазовые характеристики рассеяния двух пучков

Рис. 3.8.
с. 29 [2].
107

108. Взаимное дифференциальное сечение рассеяния

I s R 2 I1 1 I 2 2 2 I1 1I 2 2 d cos( sg t d );
dP
IsR
d
2
dP
I1 1 I 2 2 2 I1I 2 12 cos( sg t d );
d
12 1 2 d ; d d exp(i d );
• σ12 – взаимное дифференциальное сечение рассеяния двух
лазерных пучков на одной частице.
• Модуль σ12 при I1=I2=0,5 Вт/м2 дает амплитудное значение
флуктуаций мощности рассеянного в единичном телесном угле
излучения на частоте сигнала Ωsg.
108

109. Двухпучковые индикатрисы рассеяния

(r )
2
Y12 ( ) lg[Re(k 12 )];
Рис. 3.9.
с. 30 [2].
109

110.

Рис. 3.10.
с. 31 [2].
110

111. Расчет мощности оптического сигнала ЛДА по дифференциальной схеме

P I1 1 I 2 2 2 I1I 2 12 cos( sg t );
1
1d ,
12
2
2 d ,
12 d 12 exp(i ),
1 2 d sin( d )d
arg( 12 ) arctg
.
cos(
)
d
1 2 d
d
111

112. Зависимость интегрального коэффициента согласования от угловой апертуры ПОС

0 0
Рис. 3.11.
с. 33 [2].
12
1 2
112

113. Энергетический расчет ЛДА

Psg (t ) 2 I1I 2 12 cos( sg t );
( Psg ) m 2 I1I 2 12 2 Ps1Ps 2
12
1 2
• На практике при проектировании ЛДА стремятся
обеспечить условие | |≈1. Тогда:
( Psg ) m 2 Ps1Ps 2 .
113

114. Расчет интенсивности излучения на оси гауссова пучка

2( x 2 y 2 )
P I 0 exp
dxdy;
2
w
2
exp(
x
)dx ;
w
P I0
;
2
2
2P
I0 2
w
114

115. Расчет амплитудного значения мощности полезного оптического сигнала

0
Psg
2
1
Y
3 PLTL 10
w
PL – мощность излучения лазера; – длина волны;
w – радиус перетяжек пучков в точке измерения;
– телесный угол приема рассеянного излучения.
TL T3 T1T2 ;
Т1, T2 – коэффициенты пропускания элементов СФЗИ и потока по
трассам 1-го и 2-го зондирующих пучков от выходного зеркала лазера до
их пересечения в точке измерения;
T3 – коэффициент пропускания потока и ПОС по трассе рассеянного
пучка от точки измерения до светочувствительной поверхности
фотоприемника.
115

116. Соотношения, использованные при выводе формулы для Psg0

0
Psg
2T3 I1I 2 12 ;
12 1 2 1 2
Y1
1 10 ;
2
Y2
2 10 ;
2
2T1PL
I1
;
w2
1
Y (Y1 Y2 );
2
2T2 PL
I2
;
2
w
116

117. Определение Y1,Y2

Y1 Y ( 1 ), Y2 Y ( 2 );
k1
2
a/2
1
a/2
0
k2
3
Схема осевого приема с рассеянием вперед.
1)
a
1 2 ;
2
Схема осевого приема с рассеянием назад.
2)
a
1 2 ;
2
Схема внеосевого приема.
3)
a
1 0 ,
2
a
2 0 .
2
117

118. Учет потерь на оптических элементах

Френелевские потери:
2
n21 1 n 1,5 R 0.04 ( 4%);
Rn
; 21
n
n21 1
Например, при N=16
T16 (1 Rn )16 0,52 ( 52%).
Просветляющие интерференционные покрытия
снижают френелевские потери до <0,3%.
118

119. Потери при распространении излучения в рассеивающей среде

Закон Бугера:
I ( x) I ( x dx) осл ndxI ( x);
I(x)
I ' ( x) / I ( x) осл n;
I(x+dx)
I (l ) I 0 exp( осл nl )
осл рас погл T ;
T r Qs .
2
(см. рис. 3.5. на с. 24 [2])
Длина l пути луча в среде определяется геометрией
измерительной схемы;
Концентрация неоднородностей n определяется из расчета 1
частица в измерительном объеме;
119

120. Обеспечение высокой эффективности фотосмешения рассеянных волн

;
x
2
;
x
На практике диафрагменное число объектива:
D 1
при R 2 F
F 2
F
2;
D
1
0,063
16
120

121. Пример расчета сигнальной мощности

2
Дано : x 5; m 1,4; 0,01; 0,4;
x
0,63 10 6 ; w 50 10 6 ; Y 1 Y 2 Y 1;
h 6,6 10 34 ; c 3 108 ; TL 0,5; PL 5 10 3.
Psg0
2
1
Y
3 PLTL 10
w
0
Psg
1,28 10
9
Вт.
121

122. Расчет числа фотоотсчетов за пролет частицы

u 1; 0,1;
0
Psg
0
Psg ( w / u )
W
n
;
h
h(c / )
h(c / )
n 2 104.
1
0,7%.
n
122

123. 3090 отсчетов за пролет

123

124. 82 отсчета за пролет

124

125. 18 отсчетов за пролет

125

126. Применение ЛИС для измерения размеров микрочастиц

• Метод
асимметрии
индикатрисы
рассеяния
• Метод полной
индикатрисы –
многоканальная
регистрация под
разными углами
126

127. Метод спектральной прозрачности

• Используется закон Бугера и калибровочные
зависимости сечений рассеяния от длины волны
излучения.
I I 0 exp( nx ), ( ).
127

128. Статистический метод флуктуаций

• Измеряется число импульсов фототока в единицу
времени и их амплитуды при пролете частиц через
измерительный объем с калиброванным сечением.
• Возможно измерение концентрации и распределения
рассеивающих частиц по размерам
128

129. Броуновское уширение спектра рассеянного излучения

• В основе эффект Доплера.
• Диапазон размеров частиц от 1 нм (макромолекулы)
до 30 мкм.
• ДНК, бактерии, вирусы, клетки.
• Для измерений применяются методы лазерной
доплеровской спектроскопии, включая
фотонкорреляционные методы.
129

130. Измерение модуля дифференциального сечения рассеяния и глубины модуляции фототока в ЛДА

• Наглядное объяснение принципа измерений дает
интерференционная модель ЛДА
130

131. Фазовый метод в ЛДА

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

132. Пространственное разрешение ЛДА

• Пространственное разрешение R определяет число
независимых измерений скорости потока в единичном
объеме:
1
R
,
V pb
• Vpb - пробный (измерительный) объем.
• Пробный объем – объем той области пространства,
откуда поступает рассеянный свет, несущий
информацию о скорости.
132

133. Пробный объем ЛДА

• Информацию о скорости несет свет, рассеянный
частицей, пролетающей через интерференционную
область пересечения лазерных пучков.
133

134. К расчету весовой функции пробного объема

Z’
Z
a/2
a/2
X’’
X
k1
X’
Z’’
k2
134

135. Весовая функция пробного объема

Psg Psg ( x, y, z ) 2 I1 ( x, y, z ) I 2 ( x, y, z ) 12 ;
H 0 ( x, y , z )
Psg ( x, y, z )
Psg (0,0,0)
h1 ( x, y, z )h2 ( x, y, z );
( x' ) 2 ( y ' ) 2
h1 ( x' , y ' , z ' ) exp 2
;
2
w
( x' ' ) 2 ( y ' ' ) 2
h2 ( x' ' , y ' ' , z ' ' ) exp 2
;
2
w
135

136.

x' x cos( a / 2) z sin(a / 2);
x' ' x cos( a / 2) z sin(a / 2);
y ' y ' ' y;
[ x cos(a / 2) z sin(a / 2)]2 y 2
h1 ( x, y, z ) exp 2
;
2
w
[ x cos(a / 2) z sin(a / 2)]2 y 2
h2 ( x, y, z ) exp 2
;
2
w
x 2 cos 2 (a / 2) y 2 z 2 sin 2 (a / 2)
H 0 ( x, y, z ) exp 2
2
w
H 0 ( x, y, z ) H 0 (0,0,0) 1
136

137.

Условно определяют границы пробного объема
по уменьшению уровня сигнала в e2 раз. =>
x 2 cos 2 (a / 2) y 2 z 2 sin 2 (a / 2)
2
1.
2
2
w
w
w
Пробный объем представляет собой эллипсоид, оси которого
ориентированы вдоль осей координат и равны:
2w
x
;
cos(a / 2)
V pb
y 2w;
2w
z
;
sin(a / 2)
4 w3
3 sin a
137

138. Влияние ПОС на пространственное разрешение ЛДА

138

139.

1 – апертурная диафрагма;
2 – объектив ПОС;
3 – диафрагма поля зрения;
4 – фотоприемник.
• При правильной юстировке ПОС в плоскости диафрагмы поля
зрения P2 в рассеянном свете наблюдается сфокусированное
изображение рассеивателей, находящихся в плоскости P1,
проходящей через точку пересечения осей зондирующих
пучков.
• Плоскости P1 и P2 являются сопряженными и точка пересечения
осей пучков (центр измерительного объема) отображается в
центр отверстия диафрагмы поля зрения.
1 1 1
,
z1 z2 F
F – фокусное расстояние объектива.
139

140.

• В геометрическом приближении без учета аберраций,
дифракции и дефокусировок схема имеет четкие границы поля
зрения.
• Поле зрения определяется обратной проекцией отверстия
диафрагмы 3 из плоскости P2 в сопряженную плоскость P1 .
• При выходе частицы за границы поля зрения рассеянный свет
фокусируется вне отверстия диафрагмы 3 и не попадает на
фотоприемник.
• При достаточно большой глубине резкости объектива для
диафрагмы поля зрения с круглым отверстием диаметром d
область поля зрения представляет собой цилиндр с
диаметром основания
z1 d
d ' d ,
z2 G
• Г – коэффициент увеличения объектива в сторону
фотоприемника.
140

141. Оптимизация диаметра диафрагмы поля зрения

3w
d 'опт 1,5 x
,
cos(a / 2)
d опт
3w
.
G cos(a / 2)
• При d<<dопт эллипсоид пространственного
разрешения обрезается полем зрения ПОС, что
уменьшает эффективность использования лазерного
излучения.
• При d>>dопт ПОС собирает значительную часть
рассеянного излучения, не несущего информации о
скорости, что уменьшает отношение сигнал/шум.
141

142. Учет глубины резкости объектива ПОС

• В геометрическом приближении глубина резкости
объектива ПОС составляет
z резк
2d '
2d
.
G
• При z> zрезк диаметр изображения рассеивателя в
плоскости полевой диафрагмы превышает диаметр
последней, принимаемый полезный световой поток
уменьшается.
• С учетом глубины резкости объектива ПОС одним
из условий оптимальной работы ЛДА является
условие
z z резк .
142

143. Работа ПОС при внеосевом приеме


При малых углах пересечения пучков применение осевых схем
рассеяния становится неэффективным, так как линейный размер
измерительного объема вдоль оси СФЗИ z сильно возрастает.
Для ограничения размера измерительного объема вдоль оси СФЗИ
используют схемы внеосевого приема рассеянного излучения.
Если ось ПОС установлена под углом 0 к оси СФЗИ, линейный
размер пробного объема уменьшится до величины
d'
( z )'
.
sin 0
143

144. Особенности ПОС в схемах с рассеянием назад

d
d ' ;
G
F2
G ;
F1
Конструктивная реализация ПОС может быть очень разнообразной. В
ЛДА на обратном рассеянии объектив СФЗИ используется для сбора
рассеянного назад излучения. Приемная апертура ПОС ограничена
размерами призмы.
В общем случае ПОС может содержать набор линз и объективов,
поляризационные и согласующие элементы.
Диафрагма поля зрения может быть выполнена как в виде малого
отверстия в тонкой черненой металлической пластинке, так и в
световолоконном варианте, при этом ПОС и фотодетектор с блоком
144
питания могут быть разнесены на значительное расстояние.

145. Ход лучей в ЛДА на обратном рассеянии со сдвигом частоты

• Возможна установка дополнительных оптических элементов
между диафрагмой поля зрения и фотодетектором для
согласования рассеянного светового потока со
светочувствительной поверхностью и системой внутреннего
усиления фотоприемника.
145

146. Принципы построения многокомпонентных схем ЛДА

146

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

K1 k1 k3 ( рабочий);
K 2 k2 k3 ( рабочий);
K3 k1 k2 (паразитный);
a, , 1
K1 оси X ; K 2 оси Y ; K1 K 2
1 / a, 2 / ,
u x f D1 1 , u y f D 2 2 .
2
K3
2
K1
2
K2 ;
a ;
2
2
a 1 2 .
2
147

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

K1 k1 k2 ( рабочий);
K 2 k3 k4 ( рабочий);
K 3 k1 k 4 ;
K 5 k 2 k3 ;
K 4 k4 k2 ;
K 6 k3 k1;
(паразитныe);
K1 оси X ; K 2 оси Y ; K1 K 2
a, 1 1 / a, 2 / ;
u x f D1 1, u y f D 2 2 .
a 1 2 .
148

149. Способы разделения полезных сигналов и подавления паразитных

• Поляризационное разделение
• Разделение по цвету (по спектральным
линиям)
• Частотное разделение
• Временное (коммутационное) разделение
• Разделение по когерентности
149

150. Поляризационное разделение

• В ПОС используются
поляризационные
призмы.
150

151. Разделение по цвету (по спектральным линиям)

• Используются интерференционные зеркала и
призмы.
• Необходима компенсация хроматической аберрации
фронтальной оптики СФЗИ.
151

152. Частотное разделение

0+5 МГц
4
0 7 МГц
1
2
0+7 МГц
3
0 5 МГц
152

153. Частотное разделение

Пусть f D max 0,5 МГц;
S 2-4
1-3
2-3
3-4 1-4 1-2
f, МГц
2
10
12 14
153

154. Примеры использования УЗМ

154

155. Временное (коммутационное) разделение

• Синхронное переключение рабочих пар зондирующих пучков и
измерительных каналов.
• Нет потери информации, если t<< кор.
• Недостаток – коммутационные помехи.
155

156. Разделение каналов по когерентности

1
3
2
4
l12 l1 l2 0;
l34 l3 l4 0;
x
l13 l14 l23 l24 2 x
2 x >> lког
Паразитных сигналов нет.
156

157. Схемы для измерения осевой компоненты

K k1 k2 ;
• а) схема на встречных пучках:
k2
И
k1
Z
З
4
D uz ;
uz
fD
;
/2
П
/ 2;
Измеряется проекция вдоль оси Z.
157

158. б) схема с вычитанием сигналов:

D1 (k1 k2 )u ;
D 2 (k3 k1 )u ;
D (k3 k2 2k1 )u ;
4
2
2
D (cos 2a 1)u z u z sin a
158

159. Четырехпучковая схема для определения трех ортогональных проекций скорости

Условие ортогональности
и равенства модулей
векторов чувствительности:
a 2arctg (1/ 3) 60 .
0
Возможна реализация
в дифференциальном
варианте, с опорным
пучком и в инверсном
дифференциальном.
159

160. Временные и спектральные характеристики сигнала ЛДА

e
i P;
h
ia iM ;
i i1 i2 2 i1i2 cos( sg t );
160

161.

i i01H1 ( x, y, z ) i02 H 2 ( x, y, z )
2 i01i02 H ( x, y, z ) cos( sg t );
H i ( x, y, z ) hi ( x, y, z ) F ( x, y, z );
H ( x, y, z ) H 0 ( x, y, z ) F ( x, y, z );
H 0 ( x, y, z ) h1 ( x, y, z )h2 ( x, y, z );
161

162.

u c
x x(t ), y y (t ), z z (t ), i i (t );
x(t ) u (t t p ), y y0 , z z0 ;
i (t ) iНЧ (t ) iВЧ (t );
iНЧ (t ) i pd (t ), iВЧ (t ) isg (t );
162

163. Частица пролетает вблизи центра измерительного объема

163

164. Частица пролетает на границе измерительного объема

164

165. Спектральный анализ сигнала ЛДА

Комплексное преобразование Фурье
~
i ( f ) i (t ) exp( j 2 ft )df ;
Спектральная плотность сигнала
~
G ( f ) lim ( N / T ) | ip ( f ) |2 >;
T
165

166.

166

167. Относительное уширение сигнальной составляющей

f D
TD
k
k
;
fD
tT N ef
Спектральная разрешающая способность ЛДА
f D N ef
R
;
f D
k
• Nef – эффективное число интерференционных полос,
пересекаемых траекторией частицы.
• Коэффициент пропорциональности k зависит от
формы распределения интенсивности в пучках,
параметров ПОС и от выбора уровней для оценки fD
и Nef
167

168. Для ЛДА с гауссовыми пучками и неограниченным полем зрения ПОС:

f D 0,5
2 ln 2 u
cos(a / 2);
w
u u
f D 2 sin(a / 2);
f D 0,5
fD
ln 2
ctg (a / 2);
w
168

169.

N ef
x
2w 2 sin(a / 2) 4w
tg (a / 2);
cos(a / 2)
4w
ctg (a / 2)
;
N ef
f D 0,5 4 ln 2
fD
N ef
1
1,06
;
N ef
fD
R
0,94 N ef .
f D
169

170. Специализированные процессоры сигнала ЛДА

i(t)
Вход
П
u(t)
Выход
• П – процессор – электронная система обработки
сигнала.
• Назначение процессора – получение оценки
скорости.
170

171. Процессор сигнала ЛДА счетно-импульсного типа (каунтер)

Процессор сигнала ЛДА счетноимпульсного типа (каунтер)
• Принцип работы основан на подсчете числа n
периодов T0 эталонного высокочастотного
генератора импульсов за известное число N
периодов TD доплеровского сигнала ЛДА.
• Оценка частоты доплеровского сигнала
производится по формуле:
ˆf N .
D
nT0
• Соответственно оценки компоненты скорости:
uˆ x
fˆD
2 sin(a / 2)
171

172.

Блок схема каунтера
Входной
формирователь
4
8
Блок измерения
временных
интервалов
A
1
B
2
1
3
1
9
N/2
5
6
7
10
N В
11
N=2,4,8,16,32,64. Разрешение ±1 нс.
C
1 – усилители; 2 – фильтр верхних частот;
3 – фильтр нижних частот;
4 – регулировка усиления; 5 – буфер;
6 – триггер Шмита;
7 – формирователь стробов;
А – вход, В – фильтрованный выход.
12
13
D
8 – генератор тактовых импульсов;
9 – счетчик N/2; 10 – счетчик N;
11 – ЦАП; С – аналоговый выход;
12 – цифровой регистр;
D – цифровой выход.
172

173. Работа формирователя стробов

173

174. Принцип измерения длительности строба

T
T0
n 18
Tˆ nT0 ;
174

175. Принцип проверки результатов

2TˆN / 2
1

N
175

176. Факторы, определяющие погрешность оценки частоты.

• Быстродействие и надежность триггерной
системы.
• Временное разрешение при измерении
длительности стробов.
• Соответствие частоты сформированных
импульсов частоте сигнала (влияние шумов).
176

177. Пример оценки погрешности, связанной с конечным временным разрешением таймера

T0 1 нс;
f sg 200 МГц;
N 8;
N
TN
40 нс ;
f sg
n 40 1;
1 1
f
2,5%;
n 40
f sg 40 МГц f 0,5%;
177

178. Влияние шума на работу каунтера

i(t)
i
A
t
t
i(t ) Asin(2 f sg t );
2
S A
Df :
2;
N i
178

179.

i
di
i
t
dt
A
t
A
S
N
;
A2 f sg
t 0
2 A
;
Tsg
S
N
2 A
;
Tsg
t
Tsg
S
2
N
;
179

180.

Tn t2 t1;
Tn
2
t
1
2
t
2
Tn t2 t1;
2 t
1
Tsg
S
2
N
;
180

181.

Tn
Tsg ;
n
Tsg
f sg
Tsg
Tsg
Tsg
f sg
f sg
Tsg
S
2 n
N
1
S
2 n
N
Tsg
Tsg
;
;
1
S
2 n
N
;
181

182. Оптимальные условия работы каунтера

n N ef ;
f sg
S
S
2
;
N
N 0
1
S
N ef
N 0
;
182

183. Шумы в ЛИС Отношение сигнал/шум

183

184. Расчет отношения сигнал/шум

S
Df :
N
isa
e
MPsg ;
h
2
isa
2
ina
;
Psg 2 P1P2 ;
2
2
isa
e
2
4 M P1P2 ;
h
184

185.

Формула Шоттки: in2 > 2e i > f ;
e
i >
( P1 P2 Pb ) id ;
h
2
ina
> in2 > M 2 F iT2 >;
185

186. Три основные составляющие шума на выходе фотоприемника

• Дробовой шум детектируемого излучения
2 2
2
e
M F f
2
ins >
( P1 P2 Pb );
h
• Дробового шум темнового тока
2
ind
> 2eid M F f ;
• Тепловой шум
iT2
2
4kT f
>
;
RL
186

187. Приведение шума к эквивалентной оптической мощности на входе фотоприемника входу

e
ia
MPвх ;
h
Pвх.экв
h
eM
i 2na > ;
187

188.

• Приведенная ко входу фотоприемника эквивалентная
мощность дробового шума детектируемого излучения
2h F f
Pns
( P1 P2 Pb ) .
• Приведенная ко входу фотоприемника эквивалентная
мощность дробового шума темнового тока
Pnd
h
2id F f
.
e
• Приведенная ко входу фотоприемника эквивалентная
мощность теплового шума
2h
Pnt
eM
kT f
.
RL
188

189. Формула для расчета отношения сигнал/шум

S
N
2
Psg
2
2
Pns Pnd
2
Pnt
;
189

190. При расчете (S/N)0:

1;
Pb
0
Psg
P10 , P20 ;
2
Y
3 PLTL 10 ;
w
0
Pns
1
2h F f 0
( P1 P20 Pb ) ;
2
Y ( i )
0
Pi 3 PLTiT3 10
;
w
2
1
190

191. Типичные характеристики используемых в ЛДА фотоприемников

Тип
Квантовая
эффективность
Число темновых
фотоэлектронов в
секунду
Коэффициент
шума
ФЭУ
=20%
( =0,5 мкм)
=14%
( =0,6 мкм)
id/e ≈ 5·104 c-1 F ≈ 2,5
Фотодиоды
=72%
( =0,6 мкм)
id/e ≈ 1011 c-1
F≈1
191

192. Примечания

• В расчетах принимается стандартное значение
сопротивления нагрузки RL=50 Ом.
Предполагается, что постоянная времени =RCпар
не ограничивает временное разрешение.
• При использовании каунтера полоса f
принимается равной полосе пропускания
используемого ФЭУ f= fPM’.
• При использовании следящих систем обработки
полоса f принимается равной максимальной
ширине спектра доплеровского сигнала f≈ fPM/Nef.
• Коэффициент M внутреннего усиления ФЭУ
определяется из условия ограничения
максимального анодного тока imax<100 мкА.
192

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

• Оценка энергетического спектра сигнала
(распределения мощности сигнала по частотам)
1 ~
2
ˆ
G( f ) xT ( f , f ) .
T
• Алгоритм оценки: фильтрация, возведение в квадрат,
усреднение.
193

194. Оценка энергетического спектра

x(t)
УПФ
y(t)
fф, fф
|y(t)|2
Квадратор
Усредняющее
устройство
z(t)
Т
1T
1
1 ~
2
2
2
z y (t ) dt
yT (t ) dt yT ( f ) df .
T0
T
T
~
~
~
yT ( f ) xT ( f )hф ( f );
2
1 ~
2~
z
xT ( f ) hф ( f ) df ;
T
194

195. Оценка энергетического спектра

G
f fф
~
;
hф ( f ) rect
f ф
fsg

f

f ф f sg
~
2
1 ~
2
z xT ( f ) hф ( f ) df Gˆ ( f ф ) f ф ;
T
z
ˆ
G( f f ф )
;
f ф
195

196.

• Чтобы измерить всю функцию G(f) нужно
просканировать fф по всему интересующему
диапазону частот.
• Практически удобнее сдвигать спектр сигнала
относительно фиксированной частоты настройки
УПФ.
fф, fф
x(t)
УПФ
M
квадратор
fгкч
ГКЧ
u(t)
ДС
uG (t)
1/ fф
усредняющее
устройство
196

197. Условия правильной работы последовательной схемы спектрального анализа

1) f ф f sg .
1
2) T >> кор
.
f ф
3) Rs
f ф
T
.
• Обеспечивается малая
величина ошибки смещения.
• Обеспечивается малая
величина случайной ошибки
(T/ кор n – число
независимых интервалов
оценивания).
• Ограничение на скорость
сканирования (время
смещения на fф не должно
быть меньше времени
усреднения).
197

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

• Спектральная плотность (энергетический
спектр)сигнала - распределения мощности
сигнала по частотам
1 ~
2
G( f ) lim
xT ( f ) .
T T
• Алгоритм оценки: фильтрация, возведение в
квадрат, усреднение.
198

199. Cредства аппаратурного спектрального анализа

• Панорамные анализаторы
спектра
• Акустооптические
• Дисперсионно-временные
• Электронные корреляторы
• Фотонные корреляторы
• Цифровые методы. БПФ
199

200. Оценка энергетического спектра

x(t)
УПФ
y(t)
fф, fф
|y(t)|2
Квадратор
Усредняющее
устройство
z(t)
Т
1T
1
1 ~
2
2
2
z y (t ) dt
yT (t ) dt yT ( f ) df .
T0
T
T
~
~
~
yT ( f ) xT ( f )hф ( f );
2
1 ~
2~
z
xT ( f ) hф ( f ) df ;
T
200

201. Оценка энергетического спектра

G
f fф
~
;
hф ( f ) rect
f ф
fsg

f

f ф f sg
~
2
1 ~
2
z xT ( f ) hф ( f ) df Gˆ ( f ф ) f ф ;
T
z
ˆ
G( f f ф )
;
f ф
201

202.

• Чтобы измерить всю функцию G(f) нужно
просканировать fф по всему интересующему
диапазону частот.
• Практически удобнее сдвигать спектр сигнала
относительно фиксированной частоты настройки
УПФ.
fф, fф
x(t)
УПФ
M
квадратор
fгкч
ГКЧ
u(t)
ДС
uG (t)
1/ fф
усредняющее
устройство
202

203. Условия правильной работы последовательной схемы спектрального анализа

1) f ф f sg .
1
2) T >> кор
.
f ф
3) Rs
f ф
T
.
• Обеспечивается малая
величина ошибки смещения.
• Обеспечивается малая
величина случайной ошибки
(T/ кор n – число
независимых интервалов
оценивания).
• Ограничение на скорость
сканирования (время
смещения на fф не должно
быть меньше времени
усреднения).
203

204. Влияние конечной ширины полосы УПФ

204

205. Возможные записи ДС при неправильном выборе параметров анализа

Время усреднения мало
Скорость сканирования велика
205

206. Дисперсионно-временные устройства спектрального анализа (одноканальные, параллельные, широкий диапазон 104-109 Гц)

x(t )
ДЛЗ
y (t )
y(t )
2
Детектор
Усредняющее
устройство
z (t ) y(t )
2
T
ДЛЗ – дисперсионная линия задержки
x(t )
t
При распространении в ДЛЗ происходит
размазывание относительно короткого
импульса за счет дисперсии скорости:
v v ( f ).
y (t )
t
Каждая спектральная компонента
имеет свой сдвиг:
( f ).
z (t )
=> Форма импульса z(t) подобна форме
t
энергетического спектра G(f) сигнала x(t).
Лит: Панорамные приемники и анализаторы спектра. Мартынов В.А.,Селихов Ю.И. 206
М: Сов. радио. 1980. Гл.7.

207. Акустооптические анализаторы спектра

va
x
F
a
0
2
2
1
1
1
0
x1
0
-1
X
0
x1 F
;
a
-1
-2
-2
I
va
a ;
f0
F 0
x1
f0 ;
va
u(t)=u0cos(2 f0 t+ 0)
x1 const f 0
207

208.

va
x
F
2
2
0
1
1
1
0
-1
X
x1
0
-1
-2
-2
I
u(t) – узкополосный случайный процесс
I1 ( x) G ( f )
xv a ;
f
F 0
208

209.

• Примечания:
• Параллельная фильтрация за счет угловой селекции,
возведение в квадрат и усреднение => форма распределения
интенсивности I1(x) подобна энергетическому спектру сигнала
u(t).
• Схема инверсная по отношению к спектральному оптическому
анализу с дифракционной решеткой.
• Если решетка регулярная, а излучение полихроматическое
дифракционная картина показывает спектр излучения.
• Если излучение регулярное (лазерное), а решетка случайнонеодородная, картина дифракции дает спектр нерегулярностей
решетки – энергетический спектр электрического сигнала
питающего УЗМ.
209

210. Панорамный анализатор спектра с компьютерной поддержкой

fsg
108,16 МГц
УПЧ-1
M1
8,16 МГц
M2
100 МГц
fгун
Г
ГУН
ПК

Д
T

uупр
УПЧ-2
ЭЛТ
ОУ
АЦС
210

211.

• Двойное преобразование частоты обеспечивает
широкий частотный диапазон и высокую
избирательность и чувствительность.
• Фактически это схема радиоприемника
супергетеродинного типа.
• Компьютерная поддержка существенно расширяет
возможности панорамного анализатора.
• Использование самых совершенных алгоритмов
обработки данных.
• Увеличение времени усреднения.
211

212. Работа ГУН в режиме ГКЧ

fгкч
Tразв=0,02-3 с
fдев
t
fдев =0,02-50 МГц
Tразв
212

213. Три режима работы детектора и ОУ

линейный
квадратичный
логарифмический
Полоса УПФ регулируется плавно в диапазоне от 3 до 70 кГц и
устанавливается фиксированной 300 кГц.
Время усреднения можно варьировать в диапазоне от 0,03 до 3 мс.
Динамический диапазон составляет 60 дб (1:1000).
Чувствительность 1 мкВ.
213

214. Параллельный многоканальный спектроанализатор

• У панорамного анализатора последовательного типа низкая
эффективность использования информации.
• При параллельном многоканальном анализе гребенкой
фильтров этот недостаток устраняется.
• Возникают трудности калибровки, учета или выравнивания
чувствительности и избирательности каналов. Дороговизна.
214

215. Процессор сигнала ЛДА следящего типа

Блок-схема
f1,Q1
ВФ
fsg
УПЧ-1
М
УО-1
f0

ГУН
f2,Q2

У
И
УПЧ-2
f0
УО-2
КФ
ЧД
215

216. Принцип работы следящей системы

• В результате действия обратной связи
управляющее напряжение uу меняется
таким образом, чтобы разность fг fsg
оставалась неизменной и равной f0.
• Таким образом, fг повторяет изменения
fsg и uу является аналогом пульсаций
скорости в точке измерения.
216

217. Графики сигналов

• На входе
t
• После ВФ
t
• После УПЧ-1
t
• После УО-1
t
217

218. Графики сигналов в ЧД при f=f0

Графики сигналов в ЧД при f f0
• После УО-1
• После УПЧ-2
• После УО-2
• На выходе КФ
+
-
+
-
+
-
+
-
uу=0
218

219.

Графики сигналов в ЧД при f f
/ 0
f
f
0
tg Q2
;
f0
2
f
а) f<f0
+
-
+
-
+
-
а) uу>0
+
-
b) f>f0
+
+
-
-
+
b) uу<0
+
219

220. Особенности следящей системы

• Два режима работы: поиск и слежение.
• В режиме поиска кольцо ОС разорвано и на ГУН подается
пилообразное напряжение. Когда уровень сигнала превысит
пороговое значение, происходит автоматический переход
системы в режим слежения.
• Предусмотрено запоминание на время пропадания сигнала.
• Предусмотрено согласование полосы пропускания УПЧ-2 с
полосой сигнала.
• Инерционность (постоянная времени, быстродействие) зависит
от коэффициента усиления и добротности УПЧ-1 и УПЧ-2.
• Чем больше усиление и меньше добротность, тем больше
быстродействие, но при этом возрастает шум и больше
вероятность срыва слежения. На практике приходится идти на
компромисс.
220

221. Тестирование следящей системы c помощью генератора ЧМ-сигналов

f дев
f sg f >
cos(2 f M 0 )
2
f
fдев
<f>
TM
1
TM
fM
t
221

222.

• Модуляционная
характеристика
U/U0
fдев=const
1
0,5
fM
fM*
• Девиационная
характеристика
U
fM=const
fдев
fдев*
222

223. Особенности работы следящей системы

• Работает, в основном, в многочастичном режиме.
• Следит за средней частотой сигнала, которая лишь
приближенно отражает скорость потока из-за
фазового шума.
• Фазовый шум связан со сменой ансамбля
рассеивателей за время пролета частицы через
измерительный объем.
• Чтобы уменьшить влияние фазового шума,
приходится увеличивать время усреднения, то есть
ограничивать быстродействие системы.
223

224. Парадоксы измерения спектра турбулентных пульсаций с помощью ЛДА со следящей системой

• 1 – результат
измерения;
Gu(f)
• 2 – истинный
энергетический
спектр
турбулентных
пульсаций скорости;
1
3
2
f
• 2 – энергетический
спектр фазового
шума.
224

225. Оценка погрешности фазового шума

f D 0,59
1
N эф T ( f D )0,5
;
• При увеличении времени усреднения Т погрешность
фазового шума становится сколь угодно малой, но
при этом уменьшается временное разрешение
системы, что при измерении скорости в
турбулентных потоках приводит к методическим
ошибкам.
225

226. Результирующая погрешность

f sg f12 f D ;
f sg
f sg
f sg
f sg
fD
1
n 2 S / N
;
fD
fD
fD
u
1
;
u
n 2 S / N
u
1
u
;
u
nD 2 S / N
n
;
nD
nD N эф
f12 0;
x
;
f f sg .
226

227. Результирующая погрешность

u
1
N эф
S / N 0
;
N эф >>1;
Total u 1% *u 0,3%;
S / N 0
1
2
u2 N эф
.
227

228. Погрешность калибровки

0
;
2m2 sin(a / 2)
2
2 2m 4 a
ctg 2 (a / 2) ;
0
2
228

229. Электронные корреляторы

i(t ) aI ( P, t );
P
ФП
i
S Sког ;
i1
i2
e M
a S
;
h
ФП ког ;
T
1
Bˆ i ( ) i1 (t )i2 (t )dt;
T0
ЛЗ
ДС
1/T
T
0

ког
.
Т

229

230. Линии задержки

Магнитофонная линия задержки (коррелограф)
x(t)
y(t)
y(t+ )
x(t)
z
z
u ленты
;
Частотный диапазон
20 20000 Гц
Акустическая линия задержки
x(t)
2
1
x(t+ )
l
l
;
u зв
uвозд 360 м/c; uводы 1500 м/c.
230

231. Умножители

Прямые методы
Использование ваттметра:
W (t ) u (t )i(t )
Использование эффекта Холла:
IB
const
b
1
I
b
B
2
Косвенные методы
Использование логарифмического усилителя:
log( xy) log x log y
Использование сумматоров и квадраторов:
1
1
2
2
xy [( x y) ( x y) ]; xy [( x y)2 x 2 y 2 ]
231
4
2

232. Осциллографический метод измерения двумерной ПВ


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

233. Цифровые методы

x(t)
АЦП
Дискретизация по времени
Квантование по уровню
xl l x; x
l(k)
1
tk k t; t кор
2
Теорема Котельникова о дискретном представлении процессов
T – длина реализации, B – полоса частот, N – число отсчетов.
Временное разрешение
1
t
2B
Спектральное разрешение
1
f
T
3
2B T
N
2BT
f t
4
5
233

234. Оценка шума квантования АЦП

2x
1
( x)2 x 0,29 x
12
При 8-битовом преобразовании (8-разрядный АЦП) =>
256 x
1000 60 дБ
0,29 x
x
0,1%
xmax
• После аналого-цифрового преобразования сигналов для оценки
статистических характеристик процессов, таких как ПВ, ХФ, КФ,
СП и др. можно использовать стандартные программы.
• В среде Mathcad удобно использовать программы быстрого
преобразования Фурье: CFFT, cfft, FFT, fft, ICFFT, icfft, IFFT, ifft.
234

235. Фотонный коррелятор

• Созданию фотонных корреляторов способствовали
два фактора:
• а) прогресс технологии в изготовлении ФЭУ с
высокой чувствительностью и быстродействием;
• б) достижения в технологии изготовления
интегральных схем (скорость обработки 109 бит/c).
Структура одноэлектронного импульса
на выходе ФЭУ
i
10-20 нс
t
100 нс
235

236. Одночастичный сигнал ЛДА

236

237. Принцип действия простейшего фотонного коррелятора

T
N
1
Bˆ ( ) i (t )i (t )dt Bˆ (l ) n (k )n (k l )
T0
k 1
(*)
T N t ; max L t ; l l t ; tk k t.
• Прореживание – счет двойками, тройками и т.д.
• Вычитание фона (n’→n).
• Более грубая однобитовая редукция используется
для сдвигового отсчета – клиппирование (привязка по
уровню):
N
1, n(k l ) > n0 ,
m( k l )
Bˆ (l ) n(k )m(k l )
0, n(k l ) n0 .
k 1
(**)
237

238. Блок-схема однобитового фотонного коррелятора

от генератора
тактовой частоты
фотоэлектронные
импульсы
К.Я.
1
0
0
1
1
0
логические ячейки
“И”
счетчики до 109
от 24 до 256 каналов
• Каждый импульс задающего генератора перебрасывает
содержимое (0 или 1) из клиппирующей ячейки (К.Я.) в 1-ый
регистр, содержимое 1-го регистра во 2-ой регистр и т.д.
• Импульсы проходят на счетчик, если в регистре 1.
238

239.

• Алгоритм (**) позволяет заменить операцию умножения на
логическую операцию “И”, что значительно увеличивает
быстродействие.
• Диапазон работоспособности по интенсивности 104-108 ф/c.
Снизу ограничение связано с шумовыми фотоотсчетами, сверху
– параметрами одноэлектронного импульса и конечным
временем его обработки (мертвое время).
• Области применения: звездная интерферометрия, лазерная
диагностика турбулентности, медико-биологические
исследования движения микроорганизмов и др.
Пример фотон-корреляционных измерений сигнала ЛДА
Поток с низкой степенью
турбулентности
Поток с высокой степенью
турбулентности
239

240. Использование оптического спектрального анализа

Гипотетическая схема
k0, 0
u
a
Л1
D
Э
ИФ-П
Л2
Л3
ks, s
F
s 0 (k s k 0 )u;
D (k s k 0 )u;
240

241.

k 2hn cos k
k 2hn sin k k
k
k
2hn cos k
1
;
2hn sin k 2hn sin k
tan k
k 1
k
1
;
k
rk
D 0
rk>0
rr rk
F2
rk
k
;
F
r
k k ;
F
;
rk
D 0
rk<0
rr rk
F2
;
241

242.

Интерферометр Ф-П как оптический фильтр
D1
Л1 ИФ-П
Л2
Используется центральное пятно с
максимальной концентрацией энергии.
D1
Pвых f ( ) Pвх ,
ФП
f( ) – частотная характеристика
интерферометра – коэффициент
пропускания оптического фильтра.
При высокой добротности:
f( )
f ( )
A2
2
2
B
(
)
m
m
m+1
m
m
m c
;
hn
m+2
c
;
hn
A
T 2c
2
2R h
;
B
(1 R 2 )c
2
;
2R h
T – коэффициент пропускания,
R – коэффициент отражения
242

243. Использование интерферометра Фабри-Перо, сканируемого давлением воздуха, в ЛДА с дифференциальной оптической схемой


ИФ-П помещается в барокамеру с повышенным давлением воздуха. Воздух
медленно стравливается. Давление плавно уменьшается по линейному
закону. Показатель преломления воздуха уменьшается по линейному закону.
Если в таком режиме на входе поступает световой поток от
дифференциального ЛДА с частотами s1 и s2 ( s1 - s2= D ), фототок на
выходе будет меняться в соответствии частотной характеристикой ИФ-П:
i(t) 1
2
1 2
1 2
Очевидно, что
t
tD
где
t
Поэтому, оценку доплеровского сдвига
частоты можно производить по простой формуле:
t D D
,
t
c
- область дисперсии.
2hn
c t D
D
243
2hn t

244. Оптико-механическая сканирующая СФЗИ-ПОС (ЛДА для исследований сверхзвуковых потоков в аэродинамических трубах)

З3
П1
П2
u
ПП
Л1
ЛГ-159
Л2
Р1
Р3
М
З1
Д1
поток
З2
Р2
244

245. Система обработки сигнала

Д1
Л1
БК
ИФ-П
Л2
Д2
ФП
СО
ДУ – двухканальный усилитель;
СО – светолучевой осциллограф
Д
Вид сигналов на экране СО
i(t)
ДУ
1 2
t
Особенности обработки:
синхронное детектирование,
раздельное наблюдение,
сверхразрешение.
245

246. Некоторые технические характеристики ЛДА

Диапазон по скорости: 5-5000 м/c
Диапазоны сканирования:
250 мм (вдоль оси ЛДА)
15 мм (перпендикулярно оси ЛДА)
Лазер ЛГ-159: гелий-неоновый, одночастотный,
мощность 5 мВт, длина волны 0,6328 мкм.
Частоты модуляции:
f1=459 Гц – 1ый луч,
f2=1065 Гц – 2ой луч.
246

247. Использование интерферометра Фабри-Перо для частотной демодуляции.

f( )
m
m+1
i( )
m+2
Выбираем рабочую точку
на склоне одного из пиков.
i(t)
t
s
i(t) воспроизводит s(t)
247
t

248. Использование интерферометра Фабри-Перо для визуализации поля скорости

Ki
Ki
u
u
K
K
Ks
Ks
Ki
Обычное фото вдоль Кs
Фото вдоль Кs через ИФ-П
(Выбираем рабочую точку на склоне одного
из пиков частотной характеристики ИФ-П) 248

249. Обобщенная структурная схема ЛИС

II
V
VI
7
III
I
1
2
3
4
5
6
10
8
9
11
IV
249

250.

• I
— измерительный генератор, формирующий
зондирующее излучение с заданными и
контролируемыми характеристиками;
• II — канал распространений и взаимодействия
излучения с объектом исследования;
• III —средство измерения параметров сигнального
излучения па выходе канала взаимодействия;
• IV — система обработки информации и управления
процессом измерений.
• К названным четырем структурным единицам при
анализе погрешностей следует добавить еще две,
непосредственно не входящие в систему, но
определяющие точность результата измерения: блок
V отображает исследователя, проводящего
измерения, а блок VI — внешние факторы, не
учтенные другими структурными единицами.
250

251.

• Измерительный генератор I состоит из
метрологического лазера 1, системы формирования
зондирующего излучения 2 и аппаратуры контроля и
управления 3 параметрами метрологического
лазера и оптической системы 2.
• Канал взаимодействия II включает в себя
исследуемый облучаемый объект или объект,
оптические характеристики которого функционально
связаны с измеряемым параметром исследуемого
объекта 5, среды распространения зондирующего 4 и
сигнального излучений 6, посторонние облучаемые
объекты 7, источники фонового излучения 8,
аппаратуру контроля и управления 9 параметрами
канала взаимодействия.
251

252.

• Средство измерения параметров сигнального
излучения III состоит из входной согласующей
оптики 10 и блока первичных и вторичных
измерительных преобразователей 11.
• Первичный измерительный преобразователь
трансформирует оптический сигнал а электрический
и может функционировать в двух основных режимах:
в режиме прямого фотодетектирования, когда на
вход воздействует только сигнальное излучение, и в
режиме фотогетеродинирования, когда
осуществляется фотосмешение сигнального
излучения с опорным излучением измерительного
генератора.
• В режиме фотогетеродинирования блок 10 служит
для согласования сигнальной и опорной волн.
• Вторичные измерительные преобразователи служат
для согласования первичных сигналов с системой
обработки измерений.
252

253. Классификация погрешностей

• Классификация погрешностей ЛИС по
происхождению
Личные
погрешности
Внешние
погрешности
Инструментальные
погрешности
Погрешности
модели
(методические)
Погрешности
классификации
253

254.


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

255.

• К методическим относятся погрешности, обусловленные
приближенностью принятого метода нахождения конечного
результата, используемого в алгоритме, и имеющие место при
отсутствии помех; погрешности измерения первичной
информации из-за нестабильности характеристик первичных
преобразователей, шумов усилителей, наводок в каналах связи;
погрешности, обусловленные переходом к цифровой форме
обработки (динамические погрешности АЦП, ошибки
квантования сигнала по времени и по уровню, нестабильность
шага квантования и т. д.); ошибки округления при вычислениях,
аппроксимации, используемые в алгоритмах вычислений,
погрешности, связанные с конечностью разрядной сетки ЭВМ,
методом вычисления стандартных функции.
• Грубая погрешность модели, связанная с отнесением
измерении параметров постороннего объекта к исследуемому
объекту, называется ошибкой классификации. Например, в
лазерной доплеровской анемометрии процессоры иногда
срабатывают не по доплеровскому сигналу, а по различным
узкополосным наводкам, возникающим, если измерительные
сети плохо экранированы. При этом возможна ошибка
классификации, когда за доплеровский сдвиг может быть
принята несущая частота одной из радиостанций.
255

256. Классификация погрешностей ЛИС по характеру

Систематические
погрешности
Случайные
погрешности
Грубые
погрешности
256

257.

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

258.

• Случайной погрешностью называют погрешность,
которая изменяется случайным образом при повторных
измерениях одной и той же величины. Она отражает
менее существенные связи и ее невозможно в точности
воспроизвести, создавая те или иные условия
наблюдений.
• Частный вид случайной ошибки, когда эта ошибка
намного превосходит номинальные характеристики
погрешности измерительной системы, называют грубой
ошибкой. Обычно грубая ошибка связана с поломкой,
неисправностью, резким нарушением условий
эксперимента, грубым просчетом обслуживающего
персонала. Грубые погрешности следует исключать из
результатов измерений. Для этого существуют
специальные статистические критерии объективной
браковки данных.
258

259.

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

260. Фундаментальная система уравнений измерений

1( y1, y2 ,..., ym ) 1 x1
2 ( y1, y2 ,..., ym ) 2 x2
.......... .......... .......... .......... .
n ( y1, y2 ,..., ym ) n xn
• Представленная система уравнений называется фундаментальной
системой уравнений измерений. Она является основным объектом
исследования в параметрической теории информации и может быть
компактно записана в символической векторной форме:
( y) x
260

261. Основы параметрической теории информации

x ( x1, x2 ,..., xn )
y ( y1, y2 ,..., ym )
( 1, 2 ,..., n )
( 1, 2 ,..., n )
— n-мерный вектор результатов измерений;
— m-мерный вектор определяемых параметров;
— n-мерный вектор моделей;
— n-мерный вектор погрешностей.
• Типичная задача лазерной диагностики состоит в том, чтобы
оценить параметры электрических сигналов на выходе
фотоэлектрического преобразователя лазерной измерительной
системы, которые в свою очередь непосредственно или
косвенным образом связаны с параметрами исследуемого
физического процесса.
261

262.

• Фундаментальная система уравнений измерений всегда
является недоопределенной. Она содержит n уравнений, в
которых имеется n+m неизвестных.
• Не имеет решения в классическом математическом понимании.
• Решается только в статистическом смысле. Подбираются
наилучшие по определенным критериям оценки неизвестных
параметров.
• Статистическое решение определяется видом моделей и
статистикой погрешностей.
• Примеры:
• измерение параметров лазерного пучка;
• измерение параметров турбулентности;
• измерение параметров интерференционной картины.
262

263. Статистка погрешностей

• Центральная предельная теорема теории
вероятностей указывает на то, что во многих
практических ситуациях распределение вероятности
погрешностей является нормальным и описывается
гауссовой функцией:
( x x0 ) 2
.
w( x)
exp
2
2
2
1
263

264. Информационная энтропия

• Законы распределения погрешностей можно
обосновывать и с помощью принципа максимума
информационной энтропии (неопределенности)
результатов в условиях конкретного эксперимента.
Информационная энтропия
S w( x) ln[ w( x)]dx.
• Решается вариационная задача нахождения
экстремума (максимума) функционала S при заданных
условиях, накладываемых на w(x).
264

265. При стандартных условиях:

1) w( x) 0
2)
w( x)dx 1
3)
xw( x)dx x0
4)
2
2
(
x
x
)
w
(
x
)
dx
0
экстремум информационной энтропии достигается в
случае нормального распределения
( x x0 ) 2
wн ( x)
exp
2
2
2
1
265

266. В случае, когда точность измерений меняется случайным образом в процессе эксперимента можно отбросить условие 4) и положить:

> 0 .
В этом случае принцип максимума информационной
энтропии дает лапласовское распределение
погрешности результата измерения х
x x0
1
wл ( x)
exp
2
2
,
, 0
и релеевское распределение для
2
2
w( ) 2 exp
266

267. Свойства оценок

• Оценкой может быть названа любая функция
измерений, принимаемая в качестве подходящего
значения определяемого параметра.
• Как функция измерений, оценка является случайной
величиной.
• Качество оценки характеризуется следующими
основными свойствами:
• 1) несмещенность;
• 2) эффективность;
• 3) состоятельность.
267

268. Несмещенность


Смещением оценки называется разность
b yˆ > y.
Если среднее значение оценки совпадает с истинным значением y, b=0
и оценка называется несмещенной.
Источниками смещения оценок являются систематические
погрешности.
b – смещение;
р2 – рассеяние;
2 – дисперсия.
2р b 2 .
268

269. Эффективность


Для класса несмещенных оценок можно указать нижнюю границу
дисперсии оценки (минимальную дисперсионную границу) , которая
соответствует полному использованию информации при обработке
данных.
Оценка, обладающая минимальной из возможных дисперсий,
называется эффективной.
Эффективность любой произвольной оценки, обладающей некоторой
дисперсией можно характеризовать относительной эффективностью
E
2
y min
.
2
y
Эффективность может меняться от 0 до 1. При обработке данных
следует стремиться к тому, чтобы E было близким к единице.
269

270. Неравенство информации Рао-Крамера

b 0; yˆ > y 0; ( yˆ > y)w( x / y )dx 0.
X
( yˆ > y )w( x / y )dx 0.
y X
( yˆ > y) y w( x / y)dx w( x / y)dx 0.
X
X
=1
270

271.

( yˆ > y) y w( x / y)dx 1.
X
w( x / y ) ln w( x / y ) w( x / y ).
y
y
( yˆ > y) y ln w( x / y) w( x / y)dx 1.
X
A
AB > 1;
B
AB > 1.
2
271

272.

Воспользуемся неравенством Коши-Буняковского
2
A
( yˆ y)
2y
2
B
2
2
AB .
ln w( x / y)
y
1
ln w( x / y )
y
2
1.
2y
.
min
2
272

273.

ln w( x / y)
y
2
I информация по Фишеру.
2y
min
2
y
1
.
I
1
.
I
273

274. Обобщение на случай измерения нескольких параметров

I I информационная матрица.
ln w( x / y)
I i, j
yi
T
ln w( x / y)
y
j
;
1
D I
минимальна я дисперсионная
матричная граница ( МДМГ )
Диагональные элементы МДМГ дают минимальные дисперсии оценок
параметров, недиагональные элементы – корреляции этих оценок.
274

275. Пример анализа свойств оценок

• Измерение средней скорости турбулентного потока каунтером
ЛДА.
• Измерения независимые, равноточные с 0 несмещенные
<xi>=y с нормальной статистикой.
• Принимаем в качестве оценки среднее арифметическое
единичных измерений скорости.
• Исследуем свойства такой оценки.
1 n
yˆ xi ;
n i 1
n
w( x / y )
i 1
x y 2
exp i 2 ;
2 0
2 0
1
1 n
1
yˆ xi ny y b 0;
n i 1
n
2y
1 n
1 n
02
1
2
D xi
D xi 2 n 0 ;
2
n
n
n i 1
n i 1
275

276. Для определения эффективности используем неравенство информации

I ln w( x / y)
y
xi y
I
y i 1 2 02
2
n
I
2
y min
1
04
2
n
1
xi y
04 i 1
xi y
2
I 1 0 2y
n
min
;
n
2
i 1
2
1
n 02
04
2y E 1;
n
02
2
;
;
lim 2y 0.
n
Таким образом, выбранная оценка является несмещенной,
эффективной и состоятельной.
276

277. Оценка интенсивности турбулентных пульсаций

z ( xi y ) ;
2
(1)
1 n
ˆz ( xi yˆ )2 ;
n i 1
(2)
1 n
yˆ xi ;
n i 1
n
n n
n n
1 n
1 n
1
2
1
zˆ ( xi x j )2 xi2 2 xi x j 2 xi x j ;
n i 1
n j 1
n i 1
n i 1 j 1
n i 1 j 1
n n
1
zˆ x 2 2 xi x j ;
n i 1 j 1
277

278.

n
n
n
n
n
2
2
2
2
x
x
x
x
x
n
x
(
n
1
)
x
;
i j i i j
i 1 j 1
i 1
i 1 j 1(i j )
2
n
n
n 2
1 2
2
2
2
zˆ x 2 x
x
(
1
)(
x
x
);
2
n
n
n
1
zˆ z z
n
(3)
• Из (3) следует, что оценка (2) является смещенной при конечном
n. Смещение равно
1
bzˆ z;
n
lim bzˆ 0.
n
• Оценка (2) является асимптотически несмещенной.
278

279. Эмпирическая дисперсия

• Из проведенного анализа следует, что несмещенная
оценка z определяется по формуле:
1 n
2

( xi yˆ ) ;
n 1 i 1
(4)
• Оценка дисперсии (4) при любых n несмещенная и
асимптотически эффективная при сформулированных
условиях. Оценка (4) утверждена ГОСТом и
называется эмпирической дисперсией.
279

280. Квантовые флуктуации лазерного излучения Распределение Пуассона

• Анализ статистики фотоотсчетов можно выполнить,
основываясь на следующих двух постулатах.
• Среднее число фотособытий на площади фоточувствительной
поверхности, малой по сравнению с площадью когерентности
падающего света, за время, много меньшее времени
когерентности света, пропорционально интенсивности
падающей волны I(x,y,t), длине интервала времени t и
рассматриваемой площади фоточувствительной поверхности S.
• Вероятность каждого фотособытия не зависит от других
фотособытий и, следовательно, числа фотособытий,
происходящих в любых двух неперекрывающихся
пространственно-временных интервалах, статистически
независимы.
• При такой постановке задача о статистике фотоотсчетов
сводится к классической задаче математической теории
вероятности — задаче Бернулли.
280

281. Задача Бернулли

• В схеме независимых испытаний Бернулли некоторое
событие при каждом испытании может произойти с
вероятностью р.
• Вероятность того, что это событие не произойдет,
равна q = 1-р. Тогда, вероятность того, что в серии из
n независимых испытаний событие произойдет ровно
k раз, равна
k k
Pn (k ) Cn p
где
Cnk
n!
(n k )! k!
1 p
n k
,
— число сочетаний из n по k.
281

282. Применение схемы Бернулли к фотоэффекту.

• Пусть — среднее число эмитированных фотоэлектронов в
единицу времени, причем m постоянно на некотором интервале t.
• Разобьем интервал t на n одинаковых интервалов t = t/n.
Величина р = t имеет смысл среднего числа фотоэлектронов,
эмитируемых за t.
• При n→∞, для любого конечного t, p→0 так, что при некотором
достаточно большом n будет выполнено неравенство р <<1.
• В этом случае р можно рассматривать как вероятность появления
одного электрона на интервале t, так как вероятностями
появления двух и более электронов на t можно пренебречь, как
малыми более высокого порядка по р.
• При р <<1 на каждом малом интервале t эмитируется либо один
электрон с вероятностью р, либо эмиссии вообще нет с
вероятностью q =1-р.
• Тем самым задача свелась к схеме независимых испытаний
Бернулли.
282

283.

• Вероятность эмиссии k электронов за время t при этих
условиях можно определить, подставив в (1) р = a/n. Получим
a
Pn (k ) Cnk
n
k
a
1
n
n k
,
где a = t — среднее число эмитированных электронов за t.
Переходя в (1.2) к пределу при n→∞, получаем:
k
n!
a a
lim Pn (k ) lim
1
n
n
n (n k )! k! n
n k
n
a
k 1
n
n(n 1)...( n k 1) a n
ak
a
lim
lim 1 .
k
k!
k ! n n
n
n a
1
n
283

284.

• Отсюда, учитывая, что lim n ln 1
n
a
a , получаем
n
ak
lim Pn (k ) P(k )
exp( a).
k!
n
• Распределение P(k) называется распределением Пуассона.
Используя угловые скобки для обозначения усреднения, для
первых двух моментов распределения Пуассона имеем
k a;
2k k k
2
a;
• Таким образом, для распределения Пуассона абсолютная
среднеквадратичная флуктуация возрастает пропорционально
корню квадратному из среднего значения числа частиц, а
относительная среднеквадратичная флуктуация уменьшается
обратно пропорционально корню квадратному из среднего числа
частиц:
k
k a;
k / k 1/
k 1/ a.
284

285.

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

286.

• Примеры выборочных функций (реализаций) пуассоновского
процесса и соответствующих распределений вероятности
показаны ниже на рисунках. Рассмотрены три случая при <k> =
10, 100 и 1000.
• Относительный уровень шума для пуассоновского процесса
составляет 1/ и при увеличении среднего числа фотособытий
асимптотически стремится к нулю.
286

287. Интегральная интенсивность

• При фотоэлектрической регистрации излучения средняя
скорость счета фотособытий определяется энергией света W,
достигающей фоточувствительной поверхности приемника за
время регистрации, средней энергией кванта h < > (h —
постоянная Планка, < > — средняя оптическая частота
излучения) и квантовым выходом (средним числом
фотособытий, вызываемых одним упавшим фотоном, ≤1):
a = ( W)/(h<n>).
• Энергию W называют интегральной интенсивностью. Если
интенсивность (плотность мощности) падающего излучения I не
меняется на интервале регистрации t и постоянна по
светочувствительной поверхности фотоприемника S, то
W = ISt.
287

288.

• Используя представление об интегральной
интенсивности, распределение Пуассона можно
записать в виде
W
h
P(k )
k!
k
exp W
h
.
• Для типичных лазерных измерительных систем,
например из области лазерной анемометрии, при
временном разрешении t ≈ 10-8 с и квантовой
эффективности ≈ 0.1 среднее число фотоотсчетов
для излучения лазера с длиной волны =0,63 мкм и
мощности оптического сигнала 10-9 Вт составит 3,2.
• Такая количественная оценка показывает значимость
дробового квантового шума для лазерных
измерений.
288

289. Учет классических флуктуаций интенсивности

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

290. Формула Манделя

P(k ) P(k / W ) w(W )dW ;
0
W
h
P(k )
k!
0
k
exp W
h
w(W )dW ;
290

291. Избыточный шум

• Для среднего значения и дисперсии случайной величины k из
формулы Манделя можно найти:
k
W
h
;
2k
W
h
2
2 W
h2
2
.
• Существенное значение имеет тот факт, что дисперсия числа
фотоотсчетов состоит из двух слагаемых.
• Первое слагаемое, пропорциональное полной энергии,
определяется квантовыми шумами.
• Второе слагаемое, пропорциональное дисперсии флуктуаций
классической интенсивности, часто рассматривают как
избыточный шум, добавляемый к чисто пуассоновским
флуктуациям из-за случайной модуляции скорости счета
фотонов вследствие флуктуаций классической интенсивности.
291

292. Параметр вырождения


В общем случае отношение дисперсий классических флуктуаций и флуктуаций
квантового (дробового) шума равно так называемому параметру вырождения
фотоотсчетов dс = dw/ , где dw — волновой параметр вырождения.
Волновой параметр вырождения dw определяют как среднее число фотонов,
приходящихся на один интервал когерентности падающего излучения.
Используя параметр вырождения, дисперсию числа фотоотсчетов можно
представить в универсальной для рассмотренной полуклассической теории
форме
k2= <k>(1+ dс).
Видно, что при dс <<1 дробовой шум преобладает над классическим шумом.
При dс >>1 интенсивность классических флуктуаций значительно превосходит
интенсивность квантовых шумов.
Исключение составляет излучение хорошо стабилизированного одномодового
лазера, для которого интенсивность I = const , классические флуктуации просто
отсутствуют и распределение фотоотсчетов является пуассоновским, хотя
параметр вырождения может быть очень большим.
292

293. Распределение числа фотоотсчетов для многомодового лазерного излучения

• В этом случае компоненты напряженности
электрического поля лазерного излучения имеют
нормальное статистическое распределение и
соответственно для интегральной интенсивности
можно написать
W
1
;
w(W )
exp
W
W
• Формула Манделя в этом случае дает распределение
Бозе-Эйнштейна:
k
1
k>
P(k )
.
1 k > 1 k >
293

294. Сравнение распределений Пуассона (□) и Бозе-Эйнштейна (○) при <k> = 5

Сравнение распределений Пуассона (□)
и Бозе-Эйнштейна (○) при <k> = 5
Дисперсия распределения Бозе-Эйнштейна для любого
<k> может быть выражена через среднее значение числа
частиц в виде
k2= <k> + <k>2.
• Первое слагаемое определяется квантовыми флуктуациями,
второе слагаемое связано с флуктуациями классической
интенсивности.
• Относительный уровень шума в этом случае составляет более
100%. Как следствие, многомодовое излучение не обеспечивает
294
предельной точности измерений.

295. ПРИМЕРЫ ОЦЕНОК ПРЕДЕЛЬНОЙ ТОЧНОСТИ ЛАЗЕРНЫХ ИЗМЕРИТЕЛЬНЫХ СИСТЕМ


Прикладные задачи лазерной диагностики связаны с измерением
параметров лазерных пучков, параметров пространственной и
временной модуляции лазерного излучения при его рефракции,
дифракции и рассеянии на исследуемых объектах.
Особое место занимает лазерная интерферометрия, в которой
измеряются параметры интерференционных полей. К
интерференционным задачам можно отнести не только обработку
сигналов амплитудных интерферометров и интерферометров
интенсивностей, но и сигналов Фурье-спектрометров, фотонных
корреляторов, лазерных анемометров.
Интерференционный эксперимент напрямую связан с когерентностью
света. Основные статистические характеристики излучения, такие как
степень когерентности, взаимная когерентность, взаимная
спектральная плотность, непосредственно выражаются через
характеристики интерференционного поля.
295

296.

Параллельная схема регистрации с линейкой приемных элементов:
1 — исследуемое оптическое поле; 2 — линейка или матрица фотоприемников
Каждый элемент накапливает фотоотсчеты и в конце общего периода
счета регистрируется n-компонентный вектор k.
Каждая составляющая вектора k есть число фотособытий ki,
зарегистрированных отдельным i-м элементом многоэлементного
фотоприемника за время регистрации.
296

297.

Последовательная схема регистрации с рассеивающей частицей:
1 — исследуемое оптическое поле, 2 — рассеивающая частица,
движущаяся в исследуемом поле, 3 — приемник рассеянного излучения.
Для стационарных полей при перемещении одного фотоприемника в
сечении исследуемого поля излучения можно получить эквивалентный
n-компонентный вектор фотоотсчетов.В модификации такой схемы
используется сканирование оптического поля относительно
неподвижного фотоприемника, например с помощью вращающегося
зеркала.
Временную развертку пространственного распределения
интенсивности и эквивалентный вектор фотоотсчетов k можно также
получить, регистрируя неподвижной приемной оптической системой
излучение, рассеянное движущейся в этом поле микроскопической 297
частицей.

298. Квантовые пределы точности лазерной анемометрии

• Рассмотрим три наиболее характерные схемы лазерных
измерений скорости.
• В простейшей времяпролетной схеме ЛВА-1 для облучения
может использоваться один лазерный пучок и измеряться
длительность Т1 сигнала, возникающего при пересечении пучка
частицей в области перетяжки радиусом r = l1.
• В схеме ЛВА-2 для облучения используются два пучка и
измеряется время задержки Т2 между двумя сигнальными
импульсами, возникающими при последовательном пролете
частицей двух пучков, разнесенных на расстояние l =l2.
• В доплеровской схеме ЛДА измерение скорости производится
по измерению периода Т3 сигнала, возникающего при
пересечении частицей интерференционного поля с
пространственным периодом L = l3 , образованного в результате
пересечения двух когерентных лазерных пучков.
298

299.

• Оценка скорости во всех трех описанных схемах производится по
формуле
ui li / Ti ,
• где индекс i принимает значение 1, 2, 3 в зависимости от типа
схемы (ЛВА-1, ЛВА-2 или ЛДА соответственно).
• Математические модели сигналов имеют следующий вид:
1 (t ) 01 exp[ (t t0 ) 2 / 2T12 ];
2 (t ) 02 {exp[ (t t0 T2 / 2) 2 / 2 22 ] exp[ (t t0 T2 / 2) 2 / 2 22 ]};
3 (t ) 03 exp[ (t t0 ) 2 / 2 32 ]{1 m cos[ 2 (t t0 ) / T3 o ]}.
Используются следующие обозначения: r — радиус лазерного пучка по
уровню интенсивности e-1/2; L = /[2sin(a/2)], — длина волны
лазерного излучения, a — угол пересечения пучков в схеме ЛДА; u —
вектор скорости частицы; m — глубина модуляции, 0 — начальная
299
фаза сигнала.

300.

Картины облучения частицы в лазерных анемометрах ЛВА-2 и ЛДА
Примеры регистрации последовательности отсчетов в ЛВА-2 и ЛДА
300

301.

Таким образом, сигнал ЛВА1 описывается
трехпараметрической моделью; сигнал ЛВА2 —
четырехпараметрической моделью; сигнал ЛДА —
шестипараметрической моделью.
При измерении скорости информативными являются
параметры T1, T2, T3, но поскольку другие параметры в условиях
реальных измерений заранее неизвестны, речь должна идти о
совместном оценивании, что требует расчетов минимальных
дисперсионных матричных границ оценок параметров сигналов.
301

302.

• Излучение лазера будем предполагать идеально когерентным
с постоянной интенсивностью и пуассоновской статистикой
фотонов. В этих условиях найдем информационные для
сигналов ЛВА1, ЛВА2 и ЛДА. Обращением этих матриц
получим минимальные дисперсионные матричные границы.
• Рассмотрим предельный случай ограничения точности
фотонным шумом сигнала, считая фотоприемник идеальным с
квантовой эффективностью =1 линейным преобразователем
числа квантов в единицу времени n' на входе в фототок i = en'
на выходе.
• Анализ результатов показывает, что минимальные
относительные среднеквадратичные погрешности измерения
скорости при пуассоновской статистике фотоотсчетов во всех
случаях обратно пропорциональны корню квадратному из n,
где n — общее среднее число квантов излучения,
регистрируемых фотоприемником от всего сигнального
импульса.
1
n
302

303.

– коэффициент пропорциональности, зависящий от формы сигнала,
или форм-фактор.
Для сигнала ЛВА-1
1
.
2
Для сигнала ЛВА-2
2r
.
l
Для сигнала ЛДА
2
cot a
2 r (m)
1
.
N (m)
303

304.

• Предположим, что энергетика оптических схем ЛВА1, ЛВА2 и
ЛДА одинакова, и при сравнении потенциальных погрешностей
n можно считать одинаковым.
• Тогда из сопоставления форм-факторов следует, что
времяпролетный метод с двумя пучками имеет преимущество
по точности перед времяпролетным методом с одним пучком.
Вместе с тем видно, что выигрыш по точности достигается за
счет увеличения размеров измерительного объема, поскольку
возможности уменьшения размеров пучков во всех схемах
одинаково ограничены аберрациями оптических систем.
• Анализ также показывает, что при одинаковой локальности
измерений и большой глубине модуляции потенциальная
точность доплеровского метода примерно на два порядка выше,
чем у времяпролетного.
304

305.

• Полученные соотношения для потенциальных погрешностей
нетрудно обобщить на случай фотодетектора с < 1 и
внутренним усилением, характеризуемым фактором шума F, не
нарушающим пуассоновской статистики фотоотсчетов. Для
этого достаточно вместо числа фотонов n подставить в
полученные соотношения величину h /F.
• Если взять типичные значения F = 2,5, h = 0,1, Nэф = 10, то для
обеспечения составляющей погрешности, вызванной лишь
квантовым шумом излучения лазера, не более 0,1%
потребуется, чтобы на светочувствительную поверхность
фотоприемника упало в среднем не менее nmin = 2,5∙104
фотонов за полное время пролета рассеивающей частицы
через облучаемую область.
• Знание nmin позволяет рассчитывать в зависимости от
конкретной постановки задачи предельные дальность, скорость,
параметры рассеивающей частицы или необходимую мощность
лазера.
305

306. Квантовые пределы точности лазерной интерферометрии.

Квантовые пределы точности лазерной
интерферометрии
.
• Многие задачи лазерной интерферометрии связаны с
измерением параметров интерференционного поля,
имеющего конечные размеры.
• Пусть центр распределения расположен в начале координат и
нормаль к полосам задает направление оси X.
• Огибающая распределения интенсивности по x имеет
гауссову форму с полушириной r, определенной по уровню 1/e
от максимального значения, ширина полос и видность V
постоянны вдоль х.
• Параметрическая модель такого поля имеет вид
x 2
2 x
I ( x) I 0 exp 1 V cos
.
r 2
306

307.

• Если уравнение движения рассеивающей частицы задать в
виде x = ut (линейное сканирование с постоянной скоростью
вдоль оси x), модель скоростной функции можно представить в
следующей обобщенной форме:
(t ) 0 exp( b 2t 2 ) 1 V cos( 2 Nbt ) ,
где b = u/r; N = r/ - эффективное число интерференционных
полос.
• Модель содержит пять параметров: y = ( 0,b,V,N, ). Часть из
них в зависимости от конкретного содержания эксперимента
можно считать известными константами.
307

308.


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

309. Методы получения оценок

ˆ ; ˆ - вектор невязок;
( y) x; x ( yˆ )
( ˆ ) - целевая функция;
Формальная процедура:
ищется экстремум целевой функции по
ŷ.
Примеры целевых функций:
n
2
ˆ
ˆ
( ) i
- метод наименьших квадратов (МНК);
i 1
n
( ˆ ) ˆ i
i 1
- метод наименьших модулей (МНМ);
309

310. Метод максимального правдоподобия (ММП)

Функция правдоподобия – условная плотность
вероятности как функция вектора параметров:
w( x / y )
или
L( x / y) ln w( x / y).
Система уравнений правдоподобия:
L( x / y) 0.
y
310

311. Применение ММП при нормальном распределении погрешностей

Равноточные измерения
1 n
2
L( x / y )
(
x
y
)
n ln( 2 0 ).
i
2 02 i 1
ММП при нормальном распределении дает МНК
L( x / y )
y
n
( xi y) 0;
i 1
1 n
( xi y).
02 i 1
1 n
yˆ xi ;
n i 1
02

.
n
311

312.

Применение ММП при нормальном
распределении погрешностей
Неравноточные измерения
n
( xi y )
i 1
i2
wi
0;
1
i2
;
n
ˆ ( xi y ) 2 wi ;.
i 1
ММП дает МНК со взвешиванием.
n
xi wi
yˆ i 1
n
wi
i 1
;

1
n
;
wi
i 1
312

313. Применение ММП при лапласовском распределении погрешностей

1 n
L( x / y ) xi y n ln( 2 );
i 1
• В этом случае целевая функция
соответствует методу наименьших модулей.
• Уравнение правдоподобия решить нельзя,
так как нет непрерывности производной
функции правдоподобия по у.
313

314. Байесовские оценки. Метод максимума апостериорной вероятности.

Формула Байеса
P( A B) P( A) P( B / A);
P( A B) P( B) P( A / B);
P( B) P( A / B)
P( B / A)
;
P( A)
w( y) w( x / y)
w( y / x )
;
w( x )
ММАВ в отличие от ММП позволяет учесть априорную вероятность
оценки. При увеличении числа измерений ММАВ оценки
совпадают с ММП оценками, так как информация, получаемая из
эксперимента становится больше любой априорной информации.
314

315. ДОВЕРИТЕЛЬНАЯ ОЦЕНКА РЕЗУЛЬТАТА СЕРИИ ОДНОРОДНЫХ ИЗМЕРЕНИИ


Пусть x1, x2, ..., xn — единичные результаты серии из n независимых
измерений величины y. Точечная оценка дисперсии этих результатов
может быть представлена в виде эмпирической дисперсии:
2эмп
1 n
2
ˆ
(
x
y
)
,
i
n 1 i 1
1 n
yˆ xi .
n i 1
Предположим, что xi имеет нормальное распределение со
средним y и дисперсией 2, которые априори неизвестны.
Тогда случайная величина
t
yˆ y
2эмп / n
.
описывается распределением Стьюдента Sk(t), где k — параметр
распределения, называемый числом степеней свободы. Для
рассматриваемого случая k = n—1.
315

316.

• Поскольку
lim 2эмп 2 ,
n
распределение Sk(t) при k сходится к стандартному
нормальному. При конечных k оно отличается от нормального.
Примерный вид Sk(t) при различных k показан на рисунке
ниже.
0.5
0.4
dt ( x , 2)
0.3
dt ( x , 4)
dnorm ( x , 0 , 1) 0.2
0.1
0
0
4
4
3
2
1
0
x
1
2
3
4
4
316

317.

• Вероятность того, что
равна
P
t t
t
Sk (t )dt.
t
• Можно решить и обратную задачу нахождения доверительного
интервала [–t*,t*] при заданной надежности Р. Доверительная
оценка может быть выражена через граничное значение t*,
являющееся функцией числа степеней свободы k и надежности
Р следующим образом:
2
y yˆ t (k , P ) эмп .
n
• Коэффициенты t*(k, P) находятся из таблиц процентных точек
распределения Стьюдента. В таблице на следующем слайде
приведены коэффициенты t*(k, P) для наиболее типичных
ситуаций в практике учебных лабораторных измерений.
317

318.

Распределение Стьюдента
Значение t*=t*(k, Р)
P
0,90
0,95
0,99
4
2,13
2,78
4,60
5
2,02
2,57
4,03
6
1,94
2,45
3,71
7
1,90
2,37
3,50
8
1,86
2,31
3,36
9
1,83
2,26
3,25
10
1,81
2,23
3,17
20
1,73
2,09
2,85
1,65
1,96
2,58
k
318

319. ПОГРЕШНОСТЬ КОСВЕННЫХ ИЗМЕРЕНИЙ

На практике часто возникает потребность дать оценку величины y, которая связана с
некоторыми другими величинами x1, x2, ..., xn известной функциональной
зависимостью
y f ( x1, x2 ,..., xn )
При заданных доверительных оценках величин
xi xˆi xi , i 1,2,...n
доверительную оценку y можно представить в виде
у уˆ у,
yˆ f ( xˆ1, xˆ2 ,..., xˆn ).
В случае некоррелированных погрешностей x1, x2, ..., xn
2
2
2
f
f
2 f
2
xn 2 ,
x1
x2 ...
у
x1
x2
xn
В любом случае верхняя граница погрешности y определяется неравенством
у
f
f
f
x1
x2 ...
xn .
x1
x2
xn
319

320.

• Пример. Для расчета скорости в ЛДА используется
формула
u
f
,
2 sin( a / 2)
где — длина волны лазерного излучения; f — частота
доплеровского сигнала; a — угол пересечения пучков.
В результате измерений получено f =1,25 0,01 МГц,
a = (209 1) 10-3 рад, для красной линии гелий-неонового лазера
взята по справочным данным с точностью до четырех знаков
= (6328 1) 10-10м.
Рассчитывая погрешность u по формуле, представленной на
предыдущем слайде, после деления полученного результата на u
для относительной погрешности измерения скорости получим
2
2
u
f 1 2
2
cot
(
a
/
2
)
a
u
f 4
2,5 10 8 0,64 10 4 0,23 10 4 0,01.
320

321. МНК ДЛЯ ЛИНЕЙНОЙ МОДЕЛИ

Пусть эксперимент описывается линейной по некоторому
фактору t моделью
(t ) A Bt
и получена серия результатов
xi x(ti ) (i 1, 2 , ..., n).
При нормальной статистике погрешностей и независимости
результатов измерений использование метода максимального
правдоподобия для нахождения параметров А и В приводит к
МНК и следующим оценкам:
2
t
i xi ti ti xi
ˆ
A
;
2
2
n ti ti
n ti xi ti xi
ˆ
B
.
2
n ti2 ti
321

322.

Погрешности оценок можно рассчитать по формуле
погрешности косвенных измерений. После несложных расчетов
получим
2A
2B
2 ti2
n ti2 ti
2
n 2
n ti2 ti
2
;
,
где 2 — дисперсия единичного результата измерения.
322

323.

Иногда для приближенного оценивания пользуются
полученными формулами линейного МНК без анализа
статистики погрешностей (формальный линейный МНК), сводя
нелинейные двухпараметрические модели к линейному виду
заменой переменной и функциональным преобразованием.
Например, экспоненциальная модель
(t ) A exp Bt
после логарифмирования приводится к линейной
1(t ) A1 B1t,
где
1(t ) ln (t ), A1 ln A, B1 B.
Сначала определяют Â1 и B̂1 по формулам линейного МНК,
где в качестве результатов берут
ln x1, ln x2 ,..., ln xn .
Затем рассчитывают
Aˆ exp Aˆ1 ;
Bˆ Bˆ .
323

324. Проверка гипотез

• Проверяется простая гипотеза H0 против простой
альтернативы H1.
• Два рода ошибок
• a - вероятность отвергнуть правильную гипотезу H0
(ошибка 1-го рода) – уровень значимости.
• - вероятность принять ложную гипотезу H1 (ошибка 2го рода – ложная тревога).
• Два рода правильных решений
• (1- a) - вероятность принять правильную гипотезу H0.
• (1- ) - вероятность отвергнуть ложную гипотезу H1 –
мощность правила выбора решения.
324

325. Правила выбора решения

• В простейшем случае задаются уровнем значимости (обычно
выбирают a=5%). Следует иметь в виду, что с уменьшением a
возрастает .
• В более сложном варианте задаются потери (риски) при
принятии того или иного решения.
Матрица потерь: по главной диагонали
расположены расходы на правильные
решения, а по побочной – платы (потери)
за ошибочные решения.
П00
П
П10
П01
П11
Определяется среднее значение потерь (или среднего риска)
R qП00 pП11 q( П01 П00 )a p( П10 П11 )(1 ),
где p и q априорные вероятности истинности гипотез H0 и H1. Правило выбора
состоит в том, что принимается та гипотеза, средний риск для которой меньше.
325

326.

• Как задавать риски?
• На этот вопрос данная теория не отвечает.
• Бывают случаи, когда риски можно задать в
денежном выражении, исходя из
экономических рыночных ситуаций, но это
крайне редко.
• Как правило, риски задаются интуитивно на
основании предшествующего опыта.
326

327. Исключение выбросов

Если величина среднеквадратической ошибки заранее
неизвестна, то она оценивается приближенно по результатам
измерений, т. е. вместо нее применяют эмпирический стандарт
1
2
ˆ
э
xi x ,
n 1
ˆx 1 xi .
n
При этом абсолютную величину разности между
«выскакивающим» и средним значением остальных n
(приемлемых) результатов делят на эмпирический стандарт и
полученное отношение
t xв xˆ э ,
сравнивают с критическими значениями из таблицы (см.
следующий слайд). Если при данном числе n приемлемых
результатов отношение t оказывается между двумя критическими
значениями при надежностях P1 и P2, (P2 > P1), то с надежностью
вывода, большей P1, можно считать, что «выскакивающее»
значение содержит грубую ошибку, и исключить его из
дальнейшей обработки результатов.
327

328. Таблица критических значений отношения t для браковки “выскакивающих” значений (n – число приемлемых результатов, P –

надежность вывода)
n\P
0,95
0,98
0,99
0,999
n\P
0,95
0,98
0,99
0,999
5
3,04
4,11
5,04
9,43
20
2,145
2,602
2,932
3,979
6
2,78
3,64
4,36
7,41
25
2,105
2,541
2,852
3,819
7
2,62
3,36
3,96
6,37
30
2,079
2,503
2,802
3,719
8
2,51
3,18
3,71
5,73
35
2,061
2,476
2,768
3,652
9
2,43
3,05
3,54
5,31
40
2,048
2,456
2,742
3,602
10
2,37
2,96
3,41
5,01
45
2,038
2,441
2,722
3,565
11
2,33
2,89
3,31
4,79
50
2,030
2,429
2,707
3,532
12
2,29
2,83
3,23
4,62
60
2,018
2,411
2,683
3,492
13
2,26
2,78
3,17
4,48
70
2,009
2,399
2,667
3,462
14
2,24
2,74
3,12
4,37
80
2,003
2,389
2,655
3,439
15
2,22
2,71
3,08
4,28
90
1,998
2,382
2,646
3,423
16
2,20
2;68
3,04
4,20
100
1,994
2,377
2,639
3,409
17
2,18
2,66
3,01
4,13
1,960
2,326
2,576
3,291
18
2,17
2,64
2,98
4,07
328

329. Проверки нормальности

В качестве приближенного метода проверки нормальности
распределения применяют метод, связанный с оценками центральных
моментов третьего и четвертого порядков. В случае нормального
распределения случайных ошибок эти моменты равны соответственно
3 0
и
4 3 4 .
Оценками этих моментов по результатам эксперимента служат
эмпирические центральные моменты. Для удобства сравнения
подсчитывают безразмерные характеристики, связанные с этими
моментами: эмпирические коэффициенты асимметрии и эксцесса
ˆ1
3
ˆ
x
x
i
(n 2) 3э
и
где x̂ — оценка среднего значения и
ошибки рассчитываются по формулам:
1
xˆ xi ;
n
ˆ 2
э
4
ˆ
x
x
i
3,
(n 4) 4э
—эмпирический стандарт
1
2
ˆ
э
xi x .
n 1
329

330.

О малости этих характеристик обычно судят по сравнению их с
их средними квадратическими ошибками, соответственно равными
6 n 1
n 1 n 3
для ˆ1,
24n n 2 n 3
n 1 2 n 3 n 5
для ˆ 2 ,
где n — количество измерений.
Если хотя бы одна из указанных характеристик по абсолютной
величине значительно (в 2—3 раза) превосходит свою среднюю
квадратическую ошибку, то нормальность закона распределения
следует подвергнуть сомнению и провести более тщательный анализ
результатов эксперимента. В противном случае для такого сомнения
нет оснований.
330

331. Анализ погрешностей лазерного интерферометра (Пример расчета точностных характеристик современной ЛИС)

x0
3
4
1
2
x1
L
x2
Упрощенная схема лазерного интерферометра для измерения длины [4]
Лазерный интерферометр состоит из лазера 1, светоделительного кубика 2,
двух отражателей 3, 4. Отражатель 4 связан с объектом и перемещается
на искомый отрезок L.
331

332.

При движении отражателя со скоростью v(t) выходной сигнал
фотоприемника имеет следующий вид:
t
u (t ) A0 A cos[4 ( x0 x1 ) / (4 / ) v (t )dt].
(1)
0
где A0 — интерференционный фон; A — амплитуда переменной
составляющей; (x0 – x1) — некомпенсированная разность хода;
— длина волны лазерного излучения в воздухе.
Фаза переменной части сигнала связана с величиной перемещения
следующим образом:
t
(4 / ) v (t )dt] (4 / ) L.
0
Пусть /2 N1 - число полос ИК, отсчитываемое устройством.
С другой стороны, N1=2L/ и, значит, L=N1 /2 .
332

333.

Формула для расчета длины L, приведенной к нормальным
условиям, может быть представлена в виде
L N1k 0 ,
(2)
где N1 — число длин волн излучения лазера, укладывающихся в
2L; 0 — длина волны в вакууме; k — коэффициент пересчета,
равный
k [1 a1 (T1 20)] / 2nв ,
где ai — коэффициент линейного расширения измеряемого
объекта; T1 — его температура во время измерения; nв —
показатель преломления воздуха.
333

334.

• Анализ формулы (1) показывает, что точность измерения длины
принципиально ограничена неопределенностью значения
длины волны лазера в вакууме, дискретностью отсчета
результата измерения, погрешностью определения показателя
преломления воздуха.
• Точность определения приведённой длины ограничивается
также погрешностями измерения температуры образца и
неопределенностью температурного коэффициента a1.
334

335.

• Неопределенность значения 0 приводит к погрешности
измерения длины
1 L R2 d2э ,
где R — воспроизводимость длины волны лазера; dэ —
погрешность аттестации 0 по первичному эталону. В
соответствии с принятой поверочной схемой для аттестации ho
в интерферометрах допускается использование рабочих
эталонов, аттестованных по первичному эталону. Таким
рабочим эталоном служит обычно гелий-неоновый лазер,
стабилизированный по насыщенному поглощению в иоде—127.
• Методика аттестации основана на фотосмешении излучений
поверяемого и образцового лазеров и выделении сигнала
биений при сканировании длины резонатора образцового
лазера. Погрешность определения длины волны поверяемого
лазера при этом обычно составляет не более 6 10-9.
335

336.

• Воспроизводимость аттестованного значения длины волны
лазера R зависит от выбранного метода стабилизации частоты.
Чаще используется метод стабилизации по провалу Лэмба. При
этом опытным путем установлено, что воспроизводимость
составляет 10-8—10-7 и определяет в основном погрешность 1.
Фирма «Хьюлетт—Паккард»—наиболее известный
изготовитель лазерных интерферометров — проводит
аттестацию длины волны каждого изготовляемого лазера и
указывает значение длины волны с погрешностью ±10-7, откуда
1 = 10-7L.
336

337.

• Дискретность отсчета N1 принципиально ограничена
дробовыми флуктуациями электрического сигнала. Используя
модель (1) и методику оценки потенциальной точности на
основе неравенства информации, можно получить, что при
помощи лазера мощностью 0,1 мВт, потерях в интерферометре
10% и скорости движения отражателя 20 м/мин погрешность
определения дробной доли порядка интерференции составит
0,25 10-8. Последняя цифра отражает высокие возможности
лазерной интерферометрии. Для их реализации необходимо
поддерживать стабильность параметров экспериментальной
установки на таком же уровне.
• Практически вибрации, тряски, температурные колебания,
турбулентные шумы и другие факторы приводят к «дрожанию»
интерференционной картины с амплитудой (2—3) /8. Поэтому
для большинства применений дискретность отсчета N1
выбирают равной 1/8. Оптимальность такого выбора
подтверждается практикой. Таким образом, составляющая
погрешности измерения длины, связанная с N1 может быть
оценена как 2 = /8.
337

338.

• Показатель преломления воздуха определяется по известным
значениям давления p, температуры t и влажности e. Для
0=0,63299140 зависимость nв nв(p,t,e) с учетом дисперсионной
формулы Эдлена имеет следующий вид:
(nв 1) 108 38,39 p /(1 0,003671 t ) 5,61e.
• При небольших отклонениях от нормальных условий (p = 760 мм
рт. ст., t = 20°С, e = 10 мм рт. ст.) используется линейное
приближение
n (nв nв0 )[a(t 20) b( p 760) c(e 10)] 10 8.
где a = —93 1/°C; b = 36 1/мм рт. ст.; c = —6 1/мм рт. ст.; nв0 —
показатель преломления воздуха при нормальных условиях.
Таким образом, погрешность измерения длины, связанная с
неопределенностью значения nв, ограничивается точностью
измерений p,t,e. При p = e =1 мм рт. ст. и t = 0,1°C она
составляет 3 = 1,4 10-7L.
338

339.

• Температурные деформации детали приводят к погрешностям
измерения длины, которые можно оценить по формуле
4 L [da(t 20)]2 (adt )2 .
При da=10-7, a=11,5 10-6 1/°C, t = 0,01oC, t=(20±1)°C получим
4=1,4 10-7L.
• Существенными для лазерных интерферометров являются
также погрешности, вызванные дрейфом нуля, неточностями
юстировки и установки образца. Дрейф нуля обусловлен
тепловыми деформациями опорного плеча и нескомпенсированного расстояния (x0 – x1). Наиболее эффективным способом
уменьшения дрейфа нуля является конструктивное отделение
интерферометра от лазера и приближение его к объекту
(выносная конструкция). Суммарную погрешность дрейфа
нуля, юстировки и выставления для типичных серийных
образцов можно оценить как 5= 0,3 мкм.
339

340.

• При измерениях на больших трассах следует принимать во
внимание влияние турбулентности воздуха. В [4] приведены
данные одного из экспериментов, в котором на трассе длиной
80 м были зарегистрированы флуктуации оптического пути 3 10-7
м в течение нескольких минут. Движение воздуха является
одной из основных причин, ограничивающих надежную работу
лазерных интерферометров на расстояниях 60—100 м.
Наиболее простой способ уменьшения влияния этого фактора
— ограждение луча в измерительном плече интерферометра.
• Подводя итог анализу, можно сделать вывод, что суммарная
погрешность лазерного интерферометра при указанных
условиях составляет примерно 0,5 мкм/м. Именно такая цифра
приводится обычно в проспектах зарубежных фирм,
изготавливающих лазерные интерферометры, в качестве
гарантированной величины погрешности.
340

341. Оптическая локация

• - совокупность методов обнаружения, измерения
координат, а также распознавания формы удалённых
объектов с помощью электромагнитных волн
оптического диапазона — от ультрафиолетовых до
дальних инфракрасных.
• О. л. позволяет с высокой точностью (до нескольких
десятков см) производить картографирование
земной поверхности, поверхности Луны, определять
расстояние до облаков, самолётов, космических,
надводных и подводных (используя зелёный участок
спектра) объектов, исследовать распределение
инверсионных и аэрозольных слоев в атмосфере.
341

342.

• Практически создание оптических локаторов с
большой дальностью действия, высокими точностью
и разрешающей способностью стало возможным
только с появлением таких мощных источников
когерентного излучения, как оптические квантовые
генераторы — лазеры.
• В О. л. используются те же принципы определения
координат, что и в радиолокации: оптический локатор
облучает объект с помощью передатчика и
принимает отражённое от него излучение при
помощи приёмника.
• Электрический сигнал на выходе приёмника
содержит информацию о параметрах лоцируемого
объекта; характеристики этого сигнала в среднем
пропорциональны координатам объекта.
342

343.

• Методы обнаружения объектов оптическим
локатором и определения их угловых координат в
основном такие же, как в теплопеленгации, а методы
определения дальности такие же, как в
радиолокации.
• Вследствие квантового характера взаимодействия
лазерного излучения с детектором приёмника и
когерентности лазерного излучения методы
обработки сигнала в О. л. являются статистическими.
• Если оптический локатор определяет только
расстояние до объектов, он называется
электрооптическим дальномером.
343

344. Схема и принцип действия одного из типов оптического локатора для слежения за авиационными и космическими объектами

344

345. Области применения и характеристики

• Обнаружение целей и другие задачи военного
назначения.
• Получение изображения местности.
• Определение параметров спутников.
• Стыковки космических объектов.
• Высокое разрешение по дальности, углу, скорости.
• ИК-диапазон:
• CO2-10,6 мкм;
• ИАГ-Nd3+-1,06 мкм;
• п/п GaAs-0,9 мкм.
345

346. Отражательные характеристики целей

• а) Энергетические характеристики
• Модель диффузного рассеивателя по закону Ламберта.
P0
P P0 cos
P
P0
- мощность излучения,
отраженного перпендикулярно
поверхности объекта.
- угол наблюдения.
346

347. Эффективная поверхность рассеяния (ЭПР)


ЭПР служит для сравнения отражательных свойств любых целей
независимо от их формы, материала, геометрически размеров. ЭПР
численно равна площади такой плоской поверхности, которая при
изотропном расеянии создает на приемной апертуре лазерной
локационной системы интенсивность излучения, равную фактической.
P I .
• I - интенсивность лазерного излучения вблизи цели;
• - ЭПР;
• P - эквивалентная мощность, рассеиваемая по всем направлениям.
Интенсивность на приемной апертуре:
• L - дальность до цели.
P
I
I
.
2
2
4 L 4 L
Мощность на приемной апертуре:
IS
Pпр I S
.
2
4 L
• S - площадь приемной апертуры.
347

348. Дальность действия (L)

Pп
I
;
2
( L )
Pп- мощность передатчика;
Pпр N
Pпр
Pп S
;
2 4 2
(2 ) L
2 - угол расходимости передатчика;
Pп S
L 4
;
2
2
(2 ) N
Pп S
;
2
2
(2 ) N
- потенциал локатора (размерность – м2 ).
L 4 ;
С учетом потерь в оптико-механическом тракте, затухания в атмосфере
и деполяризации при отражении от цели:
L 4 ;
348

349. б) Статистические характеристики

• Если рассеиватель диффузный, то в силу ЦПТ имеем гауссово
нормальное распределение плотности вероятности для
компоненты напряженности, релеевское распределение для
амплитуды напряженности и экспоненциальное распределение
для интенсивности.
• Поскольку для гетеродинного приема is~A , а для прямого
детектирования is~I , статистические характеристики сигнала в
этих двух случаях оказываются разными.
A2
A
w( A) 2 exp 2 ;
a
2a
1
I
w( I ) exp ;
I
I
349

350.

• При гетеродинном приеме нужно учитывать частичную
пространственную когерентность. Согласно теореме ВанЦиттерта-Цернике
L
,
r
• где - радиус когерентности, r - размер цели. Если <<d (d -
диаметр приемной поверхности), то режим одномодовый и
статистика шума релеевская.
• При наличии зеркальной компоненты, статистика усложняется.
Если имеется лишь один зеркальный блик, получаем
распределение Райса
2
2
A
A
A
з AAз
w( A) 2 exp
I
,
2 0
2
a
2a a
• Аз – амплитуда зеркальной компоненты; I0 – модифицированная
функция Бесселя нулевого порядка.
• Зная статистику шума, можно рассчитать вероятность
правильного обнаружения цели и вероятность ложной тревоги.
350

351. Ослабление лазерного излучения в атмосфере

• Ослабление связано с рассеянием и поглощением.
• Резонансное поглощение CO2 , H2O, озон (O3).
• ИАГ-Nd3+-1,06 мкм – исключительно слабое
поглощение в чистой атмосфере;
• CO2-10,6 мкм – заметное поглощение, совпадение с
линиями резонансного поглощения CO2 и H2O.
• Вместе с тем, с учетом выпадения осадков и
появлением примесей в целом длина волны 10,6 мкм
оказывается предпочтительней.
351

352. Рефракция и турбулентность атмосферы

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

353. Элементная база ИК-локаторов


Передающие устройства: CO2;ИАГ-Nd3+; п/п GaAs.
Сканирующие устройства.
Механические (преломляющие, отражающие, вращающиеся
зеркальные барабаны, пирамиды).Надежность и большие углы
отклонения. Инерционность, трудно обеспечить режим быстрого
сканирования, сопровождения цели.
Пьезоэлектрические дефлекторы. Повышенное быстродействие,
малые углы сканирования.
Приемные устройства.
Лавинные фотодиоды. Малые габариты. Питание менее 100 В.
Чувствительность, сравнимая с ФЭУ; полоса до 100 ГГц. Для
автоматического сопровождения цели используются квадрантные
фотодетекторы с коммутацией сигналов с четырех площадок.
Модулирующие устройства.
Акустооптические и электрооптические (ячейки Керра, Фарадея).
Устройства оптической развязки.
Применяются для разделения отраженного и зондирующего сигналов в
локаторах с единым приемо-передающим трактом. (Часто
используется поляризационная развязка.)
353

354. Литература

• Криксунов Л. 3., Усольцев И. Ф., Инфракрасные
системы обнаружения, пеленгации и
автоматического сопровождения движущихся
объектов, М., 1968;
• Волохатюк В. А., Кочетков В. М., Красовский P. P.,
Вопросы оптической локации, М., 1971;
• Курикша А. А., Квантовая оптика и оптическая
локация, М., 1973.
• Протопопов В.В., Устинов Н.Д. Инфракрасные
лазерные локационные системы. М.:Воениздат, 1987.
354

355. АДАПТИВНАЯ ОПТИКА

• АДАПТИВНАЯ ОПТИКА, в астрономии —
автоматическая оптико-механическая система,
предназначенная для исправления в реальном
времени атмосферных искажений изображения,
которое дает телескоп.
• Системы адаптивной оптики применяются
в оптических и инфракрасных телескопах наземного
базирования для повышения четкости изображения.
Они необходимы также для работы астрономических
интерферометров, используемых для измерения
размеров звезд и поиска их близких спутников,
особенно планет.
355

356.

• Системы адаптивной оптики имеют
и неастрономические приложения: например, когда
требуется наблюдать форму искусственных
спутников Земли с целью их опознания.
• Разработка систем адаптивной оптики началась
в 1970-е годы и приобрела особый размах в 1980-е
в связи с программой "звездных войн", включавшей
разработку лазерного противоспутникового оружия
наземного базирования.
• Первые штатные системы активной оптики начали
работать на крупных астрономических телескопах
около 2000 года.
356

357.

• Развитие фото- и видеотехники позволило быстро
фиксировать изображение объекта в режиме
киносъемки с последующим отбором наиболее
удачных изображений.
• Были разработаны и более тонкие методы
апостериорного анализа изображений, например,
методы спекл-интерферометрии, позволяющие
выявлять в размытом атмосферой пятне
расположение и яркость объектов с заранее
известными свойствами, таких как "точечные" звезды.
• Математические методы восстановления изображений
также позволяют повышать контраст и выявлять
мелкие детали. Но указанные методы неприменимы
в процессе наблюдения.
357

358.

• Задача адаптивной оптики состоит в нейтрализации
в реальном времени искажений, вносимых
атмосферой в изображение космического объекта.
• Обычно адаптивная система работает совместно
с системой активной оптики, поддерживающей
конструкцию и оптические элементы телескопа
в "идеальном" состоянии.
• Действуя совместно, системы активной и адаптивной
оптики приближают качество изображения
к предельно высокому, определяемому
принципиальными физическими эффектами
(в основном — аберрацией света на объективе
телескопа).
358

359.

• Для выработки сигналов, управляющих формой
корректирующего зеркала, обычно анализируется
мгновенное изображение яркой одиночной звезды.
• В качестве приемника используется анализатор
волнового фронта, размещенный у выходного зрачка
телескопа.
• Через матрицу из множества небольших линз свет
звезды попадает на ПЗС-матрицу, сигналы которой
оцифровываются и анализируются компьютером.
• Управляющая программа, изменяя форму
корректирующего зеркала, добивается того, чтобы
изображение звезды имело идеально "точечный"
вид.
359

360.

• Искусственная звезда. Для быстрого анализа
изображения в системе адаптивной оптики
используется опорная звезда, которая должна быть
весьма яркой, поскольку ее свет делится
анализатором волнового фронта на сотни каналов
и в каждом из них регистрируется с частотой около
1 кГц. К тому же яркая опорная звезда должна
располагаться на небе вблизи изучаемого объекта.
• Однако в поле зрения телескопа далеко не всегда
встречаются подходящие звезды: ярких звезд
на небе не так много, поэтому до недавних пор
системам адаптивной оптики были доступны
наблюдения лишь 1% небосвода. Чтобы снять это
ограничение, было предложено использовать
"искусственный маячок", который располагался
бы вблизи изучаемого объекта и помогал
зондировать атмосферу.
360

361.

• Эксперименты показали, что для работы активной
оптики очень удобно при помощи специального
лазера создавать в верхних слоях атмосферы
"искусственную звезду" (LGS = Laser Guide Star) —
маленькое яркое пятно, постоянно присутствующее
в поле зрения телескопа. Как правило, для этого
используется лазер непрерывного действия
с выходной мощностью в несколько ватт,
настроенный на частоту резонансной линии натрия
(например, на линию D2 Na).
• Его луч фокусируется в атмосфере на высоте около
90 км, там, где присутствует естественный слой
воздуха, обогащенный натрием, свечение которого
как раз и возбуждается лазерным лучом. Физический
размер светящейся области составляет около
1 м, что с расстояния в 100 км воспринимается как
объект с угловым диаметром около 1’’.
361

362.

362

363.


Адаптивная оптика. М., 1980
Токовинин А.А. Звездные интерферометры. М., 1988
Уокер Г. Астрономические наблюдения. М., 1990
363

364. Литература

Основная
1. Информационная оптика:. // Н.Н. Евтихиев, О.А. Евтихиева, И.Н. Компанец,
А.Е.Краснов, Ю.Н. Кульчин, С.Б. Одиноков, Б.С. Ринкевичюс // Под. Ред. Н.Н.
Евтихиева. – М.: Изд. МЭИ, 2000. – 612 с.
2. Б.С. Ринкевичюс, В.И.Смирнов. Вопросы метрологии лазерных
измерительных систем. М.: МЭИ, 1989.
3. В.И.Смирнов, Г.М.Янина. Расчет и проектирование лазерных анемометров.
М.: МЭИ, 1996.
4. В.И.Смирнов. Фундаментальные пределы точности лазерных измерений.
М.: МЭИ, 2006.
Дополнительная
1. Б.С. Ринкевичюс. Лазерная диагностика потоков. М.: МЭИ, 1990.
2. Ю.Н. Дубнищев, Б.С. Ринкевичюс. Методы лазерной доплеровской
анемометрии. М.: Наука, 1982.
3. Дж. Бендат, А. Пирсол. Прикладной анализ случайных данных. М.: Мир,
1989.
4. Дж. Гудмен Статистическая оптика. М.: Мир. 1988. 528 с.
5. В. И. Мудров, В. Л. Кушко. Методы обработки измерений. М.: Сов. радио.
1983. 192 с.
6. Й. Бард. Нелинейное оценивание параметров. М.: Статистика. 1979. 349 с.
7. В.И. Смирнов Когерентность излучения. Задачи и расчетное задание. М.:
Издательство МЭИ, 2000. 28 с.
364
English     Русский Правила