CC BY
16
Геннадий Васильевич Подгузов родился в 1941 г., МВТУ им. Н.Э. Баумана окончил в 1970 г. Канд. техн. наук, доцент кафедры "Физика" МГТУ им. Н.Э. Баумана. Автор 20 научных работ в области исследований свойств оптических материалов.
G.V. Podguzov (b. 1941) graduated from the Bauman Moscow Higher Technical School in 1970. Ph. D. (Eng.), assoc. professor of "Physics" department of the Bauman Moscow State Technical University. Author of 20 publications in the field of research of properties of optical materials.
УДК 681.785.574
А. Н. Морозов, С. И. Светличный, С. Е. Табалин
ВЛИЯНИЕ НЕРАВНОМЕРНОСТИ СКОРОСТИ ПЕРЕМЕЩЕНИЯ ЗЕРКАЛА ФУРЬЕ-СПЕКТРОРАДИОМЕТРА НА РЕГИСТРИРУЕМЫЕ СПЕКТРЫ
Исследованы фазовые ошибки, возникающие вследствие неравномерности скорости движения зеркала фурье-спектрорадиометра. Показано, что гармонические воздействия на подвижное зеркало фурье-спектрорадиометра приводят к существенным искажениям регистрируемого спектра.
Фурье-спектрорадиометры нашли широкое применение для решения задач экологического мониторинга природных и техногенных процессов [1, 2]. Для идентификации и анализа различного рода загрязнений в открытой атмосфере необходимо верное воспроизведение спектрального состава регистрируемого излучения при требуемом спектральном разрешении не более 4 см-1 [3]. Методика получения спектров с помощью фурье-спектрорадиометров основывается на преобразовании Фурье регистрируемых при движении зеркала интерфе-рограмм [4]. Неравномерность указанного движения зеркала приводит к фазовым искажениям, которые, в свою очередь, вызывают ошибки при расчете регистрируемых спектров.
По этой причине одной из основных проблем, возникающих при разработке фурье-спектрорадиометров, является создание системы управления подвижным зеркалом, обеспечивающей сохранение постоянной скорости движения V подвижного зеркала интерферометра. При этом требования, налагаемые на допустимую нестабильность скорости движения зеркала, могут быть достаточно жесткими.
Кроме поддержания равномерности скорости движения, система управления подвижным зеркалом фурье-спектрорадиометра должна обеспечивать проведение одного измерения за время £ < 0,3 с. Спектральное разрешение не более 4 см-1 может быть реализовано только в случае, если величина рабочего перемещения зеркала интерферометра составляет не менее ±1, 25 мм от точки нулевой разности оптического хода.
Для решения задачи поддержания постоянной скорости движения зеркала в оптико-механическом блоке фурье-спектрорадиометра обычно применяется лазерный референтный канал динамического контроля положения зеркала. Он предназначен для выработки сигналов системы управления и коммутации перемещения в пространстве зеркала интерферометра, а также сигнала запуска аналого-цифрового преобразования интерферограммы, снимаемой с ИК-канала [5, 6].
Вследствие искажения сигнала ИК-канала и его запаздывания на величину т0 при прохождении через усилительный тракт возникают амплитудные и фазовые искажения регистрируемой интерферограммы, а следовательно, и измеряемых спектров [7]. Аналогичный результат получается таюке вследствие аппаратной задержки между моментом выработки сигнала считывания значений с ИК-приемника и моментом их оцифровки аналого-цифровым преобразователем (АЦП).
Если считать усилительный тракт фурье-спектрорадиометра идеальным (т0 = 0), то неравномерность хода зеркала не будет влиять на точность воспроизведения регистрируемых спектров. Более того, если возникшая фазовая ошибка в интерферограмме постоянна или монотонно зависит от волнового числа, то это не вызывает фатальных ошибок при восстановлении спектра. Решаемая ниже задача связана, прежде всего, с амплитудными и фазовыми ошибками интерферограм-мы, которые, в конечном счете, уже нелинейно зависят от волнового числа. Подобные задачи возникают при создании мобильных приборов с применением микрокриогенных систем (МКС) охлаждения ИК-фотоприемника, которые влияют на систему управления линейного двигателя.
На рис. 1 приведен пример двух спектров одного и того же абсолютно черного тела, полученных с помощью фурье-спектрорадиометра, нестабильность скорости подвижного зеркала которого достигала 30.. .50%. Указанная нестабильность была вызвана вибрациями, связанными с работой компрессора МКС. Хорошо видно, что периодические (гармонические) изменения скорости подвижного зеркала во время скана вызывают искажения, достигающие на краях спектра двукратной величины .При этом не происходит изменения величин интегралов под спектральными кривыми, поскольку значения максимумов
В(!?),отн. ed. 120-
О
700 900 1100 1300 ¡),ctf}
Рис. 1. Вид спектров излучения Б(р) абсолютно черного тела при нестабильности скорости подвижного зеркала, равной 30... 50 %
интерферограмм в точке нулевой разности оптического хода (17803 и 17422 квантов АЦП соответственно) практически совпадают.
Наличие указанной ранее задержки то приводит как к амплитудной фильтрации регистрируемой интерферограммы, так и к ошибкам в определении разности хода лучей в интерферометре:
что вызывает фазовую ошибку, зависящую от скорости движения зеркала и от знака скорости (при движении зеркала в разных направлениях эта величина имеет противоположные знаки). Коэффициент 2 в этой формуле указывает на удвоение разности хода лучей в интерферометре Майкельсона по сравнению с физическим перемещением зеркала.
ри анализе влияния нестабильности скорости зеркала на вид получаемого спектра будем считать, что иска ения и запаздывания сигналов можно описать в виде двух функций: функции амплитудных искажений а(У) и функции фазового запаздывания <^(У). Конкретный вид этих функций зависит от схем усилительного тракта и системы выработки сигнала считывания. В простейшем случае, если регистрируемый фотоприемником сигнал проходит через инерционное звено (ЯС-фильтр) с постоянной времени то, функции а(У) и <^(У) можно представить в виде
Дж = 2toV,
a(V)
1
(1)
' 1 + (4nvT0V )2' <^(V) = 4nvT0V,
(2)
(р,рад
OA
О -0,0
_g fl |____._
' 700 900 ГШ то ¡>,см~1
Рис. 2. График зависимости фазы от волнового числа при движении зеркала фурье-спектрорадиометра в противоположных направлениях: 1 — в прямом направлении; 2 — в обратном
где v — волновое число. В формулах (1) и (2) учтено то, что изменение разности хода лучей в два раза превышает физическое смещение подвижного зеркала интерферометра.
Заметим, что фаза р линейно зависит от волнового числа v и меняет свой знак в зависимости от направления движения зеркала. На рис. 2 приведены зависимости фазы р от волнового числа v при движении зеркала фурье-спектрорадиометра в противоположны« направлениях. На этом рисунке отчетливо просматривается линейный характер поведения функции р = <^(v), что соответствует выражению (2).
Далее будем полагать, что скорость движения зеркала V мало изменяется при его перемещении x:
V (x) = Vo + SV (x), (3)
где V0 = const, а изменения скорости SV(x) малы: |SV(x)| ^ V0. Тогда функции амплитудных и фазовых искажений можно представить в виде функций от перемещения:
a(V (x)) = a(x), V(V (x)) = <P(x).
Анализ условий, которым должна удовлетворять нестабильность скорости дви ения зеркала интерферометра, для достаточно иро-кого класса частотных характеристик приемного тракта выполнен в работе [7]. Здесь рассмотрим только частный случай этого анализа,
позволяющего сделать оценки искажений, реально наблюдающихся в фурье-спектрорадиометре среднего класса спектрального разрешения.
Полученную при измерениях спектральную плотность В^) с учетом амплитудных и фазовых искажений можно представить в виде
сю
В^) = J а(х)1(х) соэ(2пга + (4)
—с
где интерферограмма I(х) связана с реальным спектром В^') регистрируемого излучения следующим соотношением:
сю
I(х) = J В(^) сов(2п^ж)^'. (5)
—с
Отметим, что функция В(V') является четной: В^') = В(—V').
Считая значения функций амплитудных искажений и фазового запаздывания малыми: |а(х) — 1| ^ 1 и |<^(ж)| ^ 1, преобразуем в первом приближении выражение (4) к виду
сю сю
¿(V а(х)1 (х) — J <^(ж)1 (х) (6)
—с —с
Подставляя в (6) выражение (5) и изменяя порядок интегрирования, получим
сю сю
¿(V )=У ¿(V') ^У а(х) соэ(2^'ж)
—с —с
сю сю
— J В(V'^ ^ У соэ(2^'ж) ¿V',
—сю —сю
или
BB(v ) =
сю сю
= У B(v')Q J a(x)(cos(2n(v' + v)x) + cos(2n(v' - v)x))dx^dv'-
— ю —ю
сю сю
- J B(v') Q J ^(x)(sin(2n(v' + v)ж) - sin(2n(v' - v)x))dx^ dv'.
—с —с
(7)
С учетом малости изменения скорости (см. формулу (3)) представим функции а(ж) и <^(ж) в первом приближении в виде
a(x) = a(V)
^(x) = ^(V)
da
v=Vo dV
+
v=Vo dV
¿V (x) = a0 + av ¿V (x), (8)
v=vo
¿V (x) = ^0 + ¿V (x). (9)
V=V0
Если функции а(ж) и <^(ж) имеют соответственно вид (1) и (2), то коэффициента, введенные в выражениях (8) и (9), принимают форму
1
a0 =
av =
1 + (4nvroV))2' 32n2v 2T02Vo
(10) (11)
(1 + (4nvToVo)2)2' ^0 = 4nvToVo, ^v = 4nvTo .
Подставляя соотношения (8) и (9) в формулу (7), получаем следую -щее выражение:
ß(v ) =
сю сю
= ao jB(v') Q j(cos(2n(v' + V)x) + cos(2n(v' - v)x))dx^ dv'+
— с —с
сю сю
+ av ß(v') ( 1 ¿V(x)(cos(2n(v'+v)x)+ cos(2n(v' - v)x))dx j dv'-
—с —с
сю сю
1
- Jß(v') Q J(sin(2n(v' + V)x) - sin(2n(v' - v)x))dx^dv'-
—с —с
сю сю
- ^v Jß(v')Q/¿V(x)(sin(2n(v'+v)x) — sin(2n(v'-v)x))dx^dv'.
—с —с
(12)
Далее воспользуемся следующими формулами [8]:
сю
/ cos(2nz>x)dx = ¿(¿>), (13)
J зт(2пг/ж)^г = 0, (14)
—с
где ¿(¿>) — дельта-функция, и введем обозначения
сю
ССов(г>) = У ¿V(х) соз(2пг>ж)^ж, (15)
—с сю
Свт(^) = У ¿V(ж) зт(2пг>ж)^г. (16)
—с
Тогда выражение (12) примет более простой вид:
B(v) = ao J B(v')(i(v' + v) + ¿(v' - v))dv'+
—с
сю
+ av [ B(v')(GCos(v' + v) + Gcos(v' - v))dv'-
— ^ у в+ V) — — V
—с
Выполняя интегрирование в первом слагаемом этого выражения и учитывая четность функции В (V'), получим
¿(V) = аоВ(V) + у ВИ(Ссов(^ + V) + Ссов^' — V—
—с
сю
— ^У В(V')(Сзт^' + V) — Сзт^' — V(17)
—с
Проведем анализ полученного выражения. Если считать, что процесс ¿V(х) в интервале волновых чисел vmax] описывается модель белого ума, для которого
Ссов(^) = Свт(£) = ау
vmax vmin
где ау — величина средних квадратичных флуктуаций скорости, то функцию ¿(V) можно оценить следующим образом:
¿(V) = (ао + ау ау )B(v).
Считая, что ауоУ ^ а0, уровень относительного изменения спектра, связанного с флуктуациями скорости движения зеркала, можно оценить по формуле
¿В аУ пол
= — оУ' (18)
В ао
или с учетом (10) и (11) окончательно имеем
¿ß 32n2v2To2Vo
ß 1 + (4nvToVo)
^v, (19)
где в выражении для ау не учитывается знак минус .
Из формулы (19) следует, что ограничения на относительные флуктуации скорости для рассматриваемого случая определяются выражением
оу = 1 + (4путоуО)2 ¿В (20)
V) 2(4п^тоК))2 В ' ()
Если считать параметры в формуле (20) следующими: V = 1000 см-1, т0 = 2 • 10-5 с, Р)=1 см/с, а предельные относительные изменения спектральной плотности — ¿В/В = 10-2, то относительные флуктуации скорости в этом случае не должны превышать 8,5 %. Аналогичный результат получен в работе [7]. При повышении допустимого уровня ошибки определения спектра допустимые относительные флуктуации скорости таюке повышаются. Напомним, что сделанные оценки справедливы только в случае, если флуктуации скорости малы.
Таким образом, если изменения скорости движения зеркала описываются широкополосным случайным процессом, то требования по равномерности хода зеркала существенно зависят как от допустимой ошибки восстановления спектра, так и от постоянной времени отклика электронного тракта т0. Из выражения (19) такке следует, что в случае идеального электронного тракта (т0 = 0) ошибки восстановления спектра не будет.
Рассмотрим другой пример изменения скорости движения зеркала. Будем считать, что вариации скорости носят гармонический характер:
¿V(ж) = ¿V) вт^п^ж + а0), (21)
где ^ = //(2Р0), / — частота гармонических изменений скорости движения зеркала, а0 — фаза колебаний. Для случая, когда / = 100 Гц, а V) = 1 см/с, имеем ^ = 50 см-1.
Гармонические изменения скорости, описываемые формулой (21), могут возникнуть вследствие вибрационных воздействий на подвижное зеркало фурье-спектрорадиометра. Причиной таких воздействий,
2
кроме вибрации компрессора МКС, могут быть и внешние периодические воздействия на фурье-спектрорадиометр.
Подставляя соотношение (21) в выражение (15), получим
Gcosi^) = ÍV J cos(2nz/x) sin(2nv0x + a0)dx =
—ю
сю
= ÍV0 cos a0 J cos(2ni>x) sin(2nv0x)dx+
-с
сю
+ ÍV0 sin a0 / cos(2nz>x) cos(2nv0x)dx. (22)
После вычисления интегралов в формуле (22) с учетом соотношений (13) и (14) имеем
„ /_ч ÍV0 cos a0 ~\\ ■ fn f ~\W7 Gcos(i/) =-2- (sin(2n(v0 + г/)ж) + sin(2n(^0 - г/)ж))аж+
-с
сю
ÍV0 sin a0 /" _ _
+--2- (cos(2n(v0 + + cos(2n(v0 — =
— сю
= (*(*, + P) + í(vb - P)). (23)
Аналогично, подставляя выражение (21) в формулу (16), получим
сю
Gsin(P) = ÍVo J sin(2nPx) sin(2nv0x + a0)dx =
—ю
сю
= ÍV0 cos a0 J sin(2nPx) sin(2nv0x)dx+
—с сю
+ ÍV0 sin a0 J sin(2nPx) cos(2nv0x)dx =
—с
сю
ÍV0 cos a0
2
— СЮ
(cos(2n(v0 — — cos(2n(v0 +
8Vo sin ao _ч ч . /г, / w j
+---- (sin(2n(vo + v)x) — sin(2n(vo — v)x))dx =
—ю
8Vo cos ao/t-z t-/ ^„ч
^^-0 (8(vo — v) — 8(vo + is)). (24)
Подставив (23) и (24) в выражение (17), приведем его к окончательному виду:
B(v) = ao B(v) +
сю
av SVo sin ao f /wcv , ч ,, , ч + -4- B(v')(8(vo + v' + v) + 8(vo — v' — v) +
-с
+ 8(vo + v' — v) + 8(vo — v' + v ))dv'—
сю
yv8Vo cos ao
4
B(v')(8(v0 - v' - v) - 8(v0 + v' + v)- 8(v0 - v' + v) + 8(v0 + v' - v))dv' =
av 8Vo sin ao,~, , = aoB(v) +-2-(B(v + vo) + B (v — vo)) +
, Уv8Vocos ao (тз( ч , Л
+--2-(B (v + vo) — B (v — vo)). (25)
Анализ полученного выражения (25) показывает, что в зависимости от значения фазы колебаний ao второе слагаемое приводит к увеличению или уменьшению значений спектральной плотности во всем регистрируемом спектре, а третье — к его перекосу. Для несколько завышенной оценки искажений спектра можно воспользоваться следующей формулой:
8B av sin a^ yv cos ao
=-8V +--8V. (26)
B ao 2ao
При получении этого выражения предполагалось, что значение суммы спектральных плотностей B(v + vo) + B(v — vo) близко к 2B(v), а их разность B(v + vo) — B(v — vo) не превышает B(v).
Если в формуле (26) положить ao = п/2, то она становится аналогичной выражению (18). Поэтому сделанные выше оценки для случая ирокополосного ума справедливы и для определения амплитудных искажений, вызванных неравномерностью скорости движения зеркала при гармонических воздействиях.
Рассмотрим влияние второго слагаемого, связанного с запаздыванием фазы сигнала, на ошибки определения спектра. Используя оценки,
B(i)),omn.si. 120
SO
40
0
1
У 2 \
/
'700 900 1100 1300 и, см~
Рис. 3. Расчетные спектры: 1 — для а0 = 0; 2 — для а0 = п
сделанные по второму слагаемому формулы (26) (при а0 = 0), для принятых выше значений параметров фурье-спектрорадиометра получим
8Б 6У
= 0,13 • —.
Б Уо
Как следует из этого выражения, если гармонические изменения скорости достаточно большие, достигающие десятков и более процентов от средней скорости, то и перекосы спектра оказываются существенными. На рис. 3 приведены спектры, полученные путем расчета по формуле (25) для следующих параметров фурье-спектрорадиометра: У0 = 1 см/с, 8У = 0,3 см/с, и0 = 50 см-1, т0 = 2 • 10-5 с. Расчет был проведен для двух случаев: при а0 = 0 (кривая 1) и а0 = п (кривая 2).
Сравнивая графики, приведенные на рис. 3 и рис. 1, видно хорошее качественное совпадение построенной теории возникновения перекосов спектров и реальных экспериментальных данных.
Из выражения (25) следует, что величина перекоса оказывается пропорциональной величине нестабильности скорости движения подвижного зеркала. При уменьшении нестабильности скорости движения отличие спектров уменьшается и их можно достаточно хорошо скорректировать путем умно ения на линейну функци вида
/ (V) = Ь0 + М,
где коэффициенты Ь0 и Ь1 могут подбираться по критерию наименьших квадратов.
Таким образом, проведенный анализ влияния неравномерности скорости движения зеркала на ошибки определения регистрируемого спектра показывает, что существенные искажения спектра могут
возникнуть при неравномерностях скорости, достигающих десяти или более процентов.
СПИСОК ЛИТЕРАТУРЫ
1. Дворук С. К., Корниенко В. Н., Кочиков И. В., Морозов А. Н. и др. Применение фурье-спектрорадиометрии для анализа необратимых химических и физических процессов в природе // Вестник МГТУ им. Н.Э. Баумана. Сер. "Естественные науки". - 2001. - № 1. - С. 87-101.
2. B e i l A., D a u m R., M a t z G. et al. Remote sensing of atmospheric pollutants by passive FTIR spectrometry // Proc. of SPIE. - 1998. - V. 3493. - P. 32-43.
3. Flanigan F. Prediction of the limit of detection of hazardous vapors by passive infrared with the use of MODTRAN // Appl. Optics. - 1996. - V. 35. - P. 6090-6098.
4. Белл Р. Д ж. Введение в фурье-спектроскопию: Пер. с англ. - М.: Мир, 1975. - 382 с.
5. Кочиков И. В., Мазничко А. А., Морозов А. Н. и др. Фурье-спектрометр-радиометр для контроля выбросов промышленных предприятий // Наука - производству. - 1998. - № 8. - С. 46-51.
6. Горчаковский С. Н., Дроздов М. С., Ивлев О. А. и др. Дистанционный газовый анализ атмосферы при помощи пассивного фурье-спектрометра. Конструкция и тестовые испытания // Изв. РАН. Энергетика. - 1999. - № 2. -С. 111-119.
7. Вагин В. А., Гершун М. А., Ж и ж и н Г. Н. и др. Светосильные спектральные приборы / Под ред. К.И. Тарасова. - М.: Наука, 1988. - 264 с.
8. Бендат Дж., Пирсол А. Применение корреляционного и спектрального анализа: Пер. с англ. - М.: Мир, 1983. - 312 с.
Статья поступила в редакцию 21.01.2005 Андрей Николаевич Морозов родился в 1959 г., окончил в 1981 г. МВТУ им. Н.Э. Баумана Д-р физ.-мат. наук, профессор, зав. кафедрой "Физика" МГТУ им. Н. Э. Баумана. Автор более 100 научных работ в области прецизионных измерений и физической кинетики.
A.N. Morozov (b. 1959) graduated from the Bauman Moscow Higher Technical School in 1981. D. Sc. (Phys.-Math.), professor, head of "Physics" department of the Bauman Moscow State Technical University. Author of more than 100 publications in the field of high precision measuring and physical kinetics.
Сергей Иванович Светличный родился в 1950 г., окончил в 1973 г. Московский физико-технический институт. Канд. физ.-мат. наук, директор филиала Института энергетических проблем химической физики РАН. Автор более 60 научных работ в области химической физики, спектроскопии.
S.I. Svetlichny (b. 1950) graduated from Moscow Physical and Technical Institute in 1973. Ph. D. (Phys.-Math.), director of Affiliation of the Institute for Energy Problems of Chemical Physics of the Russian Academy of Sciences. Author of over 60 publications in the field of chemical physics, spectroscopy.
Сергей Егорович Табалин родился в 1961 г., окончил в 1984 г. МВТУ им. Н.Э. Баумана. Научный сотрудник НИИ ПММ МГТУ им. Н.Э. Баумана. Автор более 30 научных работ в области физической оптики и прецизионных измерений.
S.Ye. Tabalin (b. 1961) graduated from the Bauman Moscow Higher Technical School in 1984. Researcher of "PMM" research institute of the Bauman Moscow State Technical University. Author of over 30 publications in the field of physical optics and high precision measuring.
