CC BY
124. Alphonse G. A. Broad-band acousto-optic deflectors using sonic gratings for first-order beam steering // RCA Rev. 1972. Vol. 33, № 3. P. 543-594.
5. Магдич Л. Н., Молчанов В. Я. Акустооптические устройства и их применение. М.: Радио и связь, 1978. 112 с.
S. S. Shibaev, D. P. Volik, V. V. Rozdobudko
Taganrog institute of technology of Southern federal university
Acousto-optic receiver-cymometer based on deflector with anthiphased ultrasound excitation
The acousto-optic SHF cymometer structure, containing deflector with antiphased transducers grating and two laser sources with different light waves, is analyzed. Unlike known acousto-optic receiver configurations the offered one is characterized with twice as much instantaneous analysis band.
Acousto-optic SHF deflector, antiphased transducer, laser source, analysis band
Статья поступила в редакцию 29 февраля 2008 г.
УДК 621.37/39
Б. Н. Вольфовский, С. С. Шибаев, В. В. Роздобудько
Технологический институт Южного федерального университета
I Алгоритмы вычисления частоты в акустооптических измерителях параметров радиосигналов
Предложены и проанализированы два эвристических алгоритма вычисления частоты как абсциссы оси симметрии видеоимпульса, считываемого с фотоприемника в акустооптических измерителях параметров радиосигналов. Приведены теоретические и экспериментальные оценки погрешностей. Экспериментально оцененная погрешность измерения частоты в расширенном динамическом диапазоне уровней входных сигналов не превышает 10 % от шага сетки частот измерителя.
Алгоритм, погрешность, акустооптический измеритель, распределение интенсивности
У различных по назначению измерителей параметров радиосигналов информация об измеряемом параметре часто характеризуется положением на оси абсцисс видеоимпульса, формируемого на выходе измерителя. При этом алгоритм работы устройства сводится к вычислению положения оси симметрии видеоимпульса на фотоприемнике и постановке ему в соответствие измеряемого параметра. Такие схемные и алгоритмические решения используются, например в измерителях, предназначенных для вычисления азимута или угла места цели, положения импульса на оси времени, положения спектральной функции на оси частот или положения светового сигнала на фотоприемнике в акустооптических измерителях частоты [1].
В настоящей статье рассмотрены два алгоритма решения подобных задач. Первый из них, в котором искомая абсцисса оси симметрии видеоимпульса ищется как абсцисса точки пересечения двух прямых, пригоден для линейного режима, а второй, в котором используется свойство симметрии видеоимпульса, - для нелинейного режима (режима ограничения) работы измерителя.
38
© Вольфовский Б. Н., Шибаев С. С., Роздобудько В. В., 2008
У УС
fB fA Fc fC fD f
fB fA Fc fC fD
б
f
a
Рис. 1
Алгоритмы апробированы в эксперименте при обработке выходной информации в макете акустооптического измерителя частоты (АОИЧ) [2]. В связи с тем, что дальнейшее изложение ведется применительно к АОИЧ, сделаем необходимые пояснения.
В АОИЧ выполняется прямое преобразование Фурье акустического аналога анализируемого радиосигнала и его результат формируется на входе дискретного фотоприемника (линейке ПЗС) в виде распределения интенсивности светового сигнала (РИСС). Продетек-тированные значения РИСС, которые можно рассматривать как ординаты видеоимпульса, считываются с фотоприемника и оцифровываются. Указанное распределение интенсивности соответствует спектру мощности анализируемого радиосигнала. Оно изменяет свое положение (местонахождение) на фотоприемнике в зависимости от частоты, поэтому определение положения РИСС на фотоприемнике, в частности нахождение абсциссы его оси симметрии (в случае симметричного РИСС) равносильно измерению частоты. При действии на входе АОИЧ гармонического сигнала оцифрованное РИСС является аппаратной функцией АОИЧ, которой и определяются его метрологические параметры [1], [2].
Определение частоты как абсциссы точки пересечения двух прямых. Пусть РИСС, формируемые на фотоприемнике, симметричны относительно частоты /с, принимаемой за истинную частоту сигнала (рис. 1). Для реализации алгоритма должны быть известны частоты /а , /в, /с, /б , соответствующие четырем фотодиодам (ФД) А, В, С, В, расположенным
вблизи оси симметрии РИСС, и уровни сигналов, снимаемых с этих ФД, уа , Ув, Ус, Ув .
Пусть РИСС представлена гауссоидой
где а - коэффициент формы. Шаг расположения ФД, пересчитанный в частотную область, постоянен и составляет А/ .
На рис. 1, а представлен случай, когда истинный максимум РИСС располагается вблизи ФД А с некоторым смещением А/ в сторону увеличения частоты. Ломаная ВАСБ представляет собой кусочно-линейную аппроксимацию РИСС, построенную по снятым с ФД сигналам. Уравнения составляющих ломаной у = gl/ + у = g2/ + С2, где gl и g2 - угловые коэффициенты, а с и С2 - свободные члены, вычисляемые, соответственно, по формулам
(1)
& = (УА - Ув )/(fA - fB ) = (УА - Ув )/AF;
(2)
g2 = (Уб - Ус )/(fD - fe) = (Уб - Ус )/AF; (3)
c1 = (yBfA - yAfB )/( fA - fB ) = (yBfA - yAfB )/AF i (4)
c2 = (yefü - yüfe )/(fD - fe) = (yefü - yüfe )/AF • (5)
Абсцисса f 0 точки пересечения прямых BA и üC вычисляется по формуле [3]:
f0 = (С1- c2 )/(g2 - gl) • (6) С учетом формы РИСС (1) сигналы, снимаемые с ФД A-D, могут быть записаны как
yA = exP(-afA) = exP(-aAf 2); yB = exP L-a (Af+AF )2 J; ye = exP L-a (AF -Af )2 J;
2
yü - exp -a (2AF - Af) . После подстановки этих зависимостей в (2)-(5) абсцисса точки пересечения прямых BA и DC определится как
f0 = -Af -AF [2exp ( 2aAfAF ) exp (-aAF 2 ) - exp ( 4aAf AF ) exp ( -4aAF 2 ) -1]/Z , (7)
где
Z = exp (4aAf AF) exp (-4aAF2 ) - exp (-aAF2 ) [exp (2aAf AF) - exp (-2aAf AF)] -1. (8) Абсцисса fo принята за оценку частоты сигнала Fc. Методическая погрешность вычисления частоты 5f составляет 5f = Fc - fo. Из (7) и (8) следует, что 8f является функцией трех величин: Af, a и AF . Две последние - коэффициент формы a и частотный интервал между диодами AF - постоянны и являются параметрами измерителя частоты. В отличие от них частотный интервал Af меняется с изменением частоты Fc входного радиосигнала; он определяет ориентацию РИСС относительно ФД, и от него зависят уровни сигналов на ФД.
Математическое моделирование показывает, что параметры a и AF влияют на параметры измерителя частоты и погрешность 5f не порознь, а совместно. Этот результат согласуется с теоретическими ожиданиями, поскольку и от a , и от AF зависит распределение уровней сигналов на ФД.
На рис. 2 показана расчетная зависимость максимальной относительной погрешности измерения частоты Sfmax/AF от произведения a(AF)2. Анализ этой зависимости показывает, что наименьшая максимальная погрешность вычисления частоты Sfmaxo описанным методом достигается при выборе a и AF из условия
a (AF )2 = 0.624. (9)
При выполнении (9) максимальная погрешность вычисления частоты составляет 1.18 % от AF .
В качестве иллюстрации на рис. 3 представлена методическая погрешность алгоритма f (Af) для коэффициента формы a = 2.32 МГц-2 и частотного интервала между ФД AF = 518.2 кГц. На этом рисунке частота Af измеряемого сигнала варьируется в пределах AF . Как следует из рис. 3, при a(AF) = 0.624 погрешность 5f внутри частотной дискреты AF изменяется в пределах ±6.1 кГц. 40
Рис. 2 Рис. 3
Алгоритм определения значения частоты состоит в следующем:
• определяются номер к ФД, сигнал с которого имеет максимальную величину, и соответствующая ему частота;
• измеряются сигналы Ук-1, Ук, Ук+1 на трех ФД, средний из которых имеет максимальный сигнал;
• при выполнении неравенства Ук+1 > Ук-1 (рис. 1, а) вводятся обозначения: уа = Ук;
Ув = Ук-1; ус = Ук+1; измеряется Ук+2, к°т°рый обозначается yD = Ук+2;
• при выполнении неравенства Ук+1 < Ук-1 (рис. 1, б) вводятся обозначения: ус = Ук;
Уа = Ук-1; yD = Ук+1; измеряется Ук - 2, который обозначается ув = Ук -2;
• по формулам (2)-(5) вычисляются угловые коэффициенты и свободные члены, а по формуле (6) оценивается частота fo .
Вычисления по формулам (4)-(6) можно упростить, перенеся начало отсчета частот на ФД, дающий максимальный сигнал. Для этого необходимо выполнить следующее:
• частоту, соответствующую ФД, дающему максимальный сигнал, положить равной нулю ( fA = 0 для рис. 1, а или fC = 0 для рис. 1, б);
• частоты соответствующие соседним ФД, обозначить: для рис. 1, а - fB = -AF; fe =AF ; fD = 2AF ; для рис. 1, б - fB =-2AF; fA =-AF ; fD = AF ;
• с учетом этих обозначений по формулам (2)-(5) вычислить угловые коэффициенты &1 , &2 и свободные члены q, C2;
• по формуле (6) вычислить значение частоты f0Q ;
• истинное значение частоты сигнала определить с учетом абсолютных значений fA или
fe как fo = foo + fA (для рис. 1, а) или fo = foo + fe (для Рис. 1, б).
Экспериментальная оценка точности определения частоты выполнялась в макете АОИЧ с гауссовским РИСС, проецируемым на фотоприемник. АОИЧ функционировал в диапазоне частот 1.5...2.o ГГц с частотным интервалом между ФД AF = 518.2 кГц . Измерения выполнялись на 5o частотах в полосе AFa « 6AF. Частота сигнала Fc в указанной полосе изменялась с дискретностью ÁF/8 « 65 кГц . Зависимости оценки точности от частоты сигнала для средней и для верхней частей динамического диапазона АОИЧ приве-
А/, кГц
А/, кГц
- 125 -
- 160
- 55 -
- 90 -
- 20
- 40
0
1749 1750 1751 1752 Fc, МГц Рис. 4
1749 1750 1751 1752 Fc, МГц
Рис. 5
f, кГц
/, кГц
а (AF )2 = 0.7
9
0
- 9
- 18
0 130 260 390 А/, кГц
0 130 260 390 А/, кГц
Рис. 7
Рис. 6
дены на рис. 4 и рис. 5 соответственно. Из данных зависимостей следуют ожидаемая периодичность изменения погрешности, а также ее смещение преимущественно в область отрицательных значений. Это смещение (систематическая составляющая) связано с дрейфом частотной шкалы АОИЧ во время измерений. Скорость дрейфа составляла 2...3 кГц/мин, а время измерений для кривых на рис. 4 - около 5 мин, а для кривых на рис. 5 - около 3 мин.
После центрирования графиков приблизительное значение погрешности составило ±60 кГц (рис. 4) и ±28 кГц (рис. 5). Погрешность включала представленную на рис. 3 погрешность метода и погрешность, связанную с имевшей место в экспериментах асимметрией РИСС.
Расчеты показали, что в эксперименте значение а(АР)2 изменялось в пределах 0.58...0.7, причем сигналу в середине динамического диапазона (рис. 4) соответствовали
а = 2.3...2.607 МГц-2 и а(ЛР)2 = 0.618...0.7, а более сильному сигналу (рис. 5) -
а = 2.16...2.403 МГц-2 и а(ЛР)2 = 0.58...0.645. Графики методических погрешностей
вычисления частоты для крайних значений параметра (0.58 и 0.7) а(АР) представлены на рис. 6 и 7 соответственно. Анализ этих графиков и сопоставление их с приведенными ранее экспериментальными оценками погрешностей позволили утверждать, что при работе АОИЧ в середине динамического диапазона методическая погрешность составляет около ±37 кГц (из общей погрешности ±60 кГц), а при работе в верхней динамического диапазона части - около ±17 кГц (из общей погрешности ±28 кГц ).
Определение частоты в режиме сильного сигнала. С увеличением мощности исследуемого сигнала и, соответственно, освещенности линейки ФД АОИЧ, ФД, на которые
а б
Рис. 8
приходится центральная часть РИСС, могут перейти в режим насыщения. В этом случае выходной сигнал АОИЧ оказывается ограниченным (рис. 8, интервалы Ж).
Рассмотрим алгоритм определения частоты сигнала Рс, считая, как и прежде, расстояния между частотами, приходящимися на любые соседние ФД, одинаковыми и равными ЛР , а РИСС симметричным относительно частоты Рс.
При смещении оси симметрии РИСС Рс вправо относительно ближайшего ФД (рис. 8, а) необходимо определить частоты, соответствующие точкам С и О (/с и /о). При этом оценка /о частоты Рс определится выражением
/о = (/с + /о )/2. (10)
Частота /с известна как соответствующая одному из ФД АОИЧ. Точка О лежит на прямой АВ на уровне сигнала ус. Ее абсцисса /о вычисляется как
/о = /А-А/, (11)
где А/ - частотный интервал между точками О и Е. Определим его на основании пропорции АН/ВН = АЕ/ОЕ, вытекающей из подобия треугольников ВАН и ОАЕ. Тогда ОЕ = ВН (АЕ/АН). С учетом, что ВН = АР , АЕ = уа - Ус , АН = уа - Ув, получим
А/ = АР (уа - Ус )/(УА - УВ ). (12)
Для определения положения оси симметрии РИСС следует подставить (11) в (10):
/о = (/с + /А -А/)/2, (13)
а в последнее выражение - величину А/, полученную из (12).
Формула (13) справедлива, если уа > Ус > Ув (рис. 8, а). Если же центр симметрии РИСС смещен относительно ближайшего ФД влево (рис. 8, б), (т. е. ус > Уа > Уб ), формулы для вычисления частоты будут иметь вид:
А/ = АР (ус - уа )ЦУс - УБ ), (14)
/о = (/с + /А +А/)/2. (15)
Описанная методика определения частоты справедлива как для гауссовской, так и для негауссовской форм РИСС. Алгоритм определения частоты состоит в следующем:
• определяются два ФД, дающих наибольшие сигналы, но меньшие уровня насыщения;
• если ФД с наибольшим сигналом соответствует более низкой частоте, чем второй из определенных ФД, то с ним ассоциируется точка А, со вторым ФД - точка С (см. рис. 8, а), а расчет проводится по формулам (12) и (13);
• если ФД с наибольшим сигналом соответствует более высокой частоте, то с ним ассоциируется точка С, со вторым ФД - точка А (см. рис. 8, б), а расчет проводится по формулам (14) и (15);
• если сигналы, даваемые обоими выделенными ФД, равны, частота /о определяется по
Поскольку погрешность вычисления частоты описанным методом 5/ = РС - /о зависит от точности аппроксимации прямыми АВ (рис. 8, а) или СБ (рис. 8, б) соответствующих участков РИСС, то вычисления следует выполнять на участках РИСС с большой крутизной, поскольку этим участкам соответствуют минимальные погрешности аппроксимации. Знание формы РИСС, и в частности его коэффициента формы (для РИСС, описываемой гауссоидой, - коэффициент а), позволяет повысить точность определения /о . Так, в
этом случае для вычисления частотного интервала ОЕ можно воспользоваться формулами, описывающими участок РИСС между точками АВ.
На рис. 9 показан график зависимости методической погрешности 5/ от частотного
интервала А/, рассчитанный для коэффициента формы гауссоиды а = 0.51 МГц и частотного интервала между ФД АР = 518 кГц .
При асимметричной РИСС погрешность 5/ содержит систематическую составляющую, которую можно оценить и учесть в последующих вычислениях. Для получения систематической составляющей погрешности следует использовать калибровочный сигнал известной частоты.
Описанный алгоритм, так же, как и алгоритм для сигналов, не вызывающих насыщения ФД, экспериментально апробирован при обработке информации на выходе макета
Экспериментальные распределения сигналов на выходе макета АОИЧ отличались от форм, использованных при математическом моделировании (см. рис. 8). Сигналы, экспериментально полученные в режиме насыщения ФД, представлены на рис. 1о, на котором
ф°рмуле /0 = (fc + fA )/2.
АОИЧ [2].
Sf,
а (ÁF )2 = 0.624
>АЦП
2880 -
1920
960
0 130 260 390 Af, кГц
Рис. 9
0
1743
1746.5 1750 1753.5 f, МГц Рис. 10
цифрой 1 обозначен сигнал, соответствующий начальному участку насыщения ФД, цифрой 2 - сигнал при среднем уровне насыщения, а цифрой 3 - сигнал при глубоком насыщении ФД. Эти сигналы отличаются от теоретических главным образом тем, что в них наблюдается асимметрия относительно частоты радиосигнала Рс, причем она увеличивается с возрастанием степени насыщения. Кроме того, уровни сигналов насыщения от отдельных ФД различны, из-за чего этот участок не выражен явно.
С учетом указанных особенностей алгоритм вычисления частоты при наличии насыщения дополнен введением ограничения сигналов ФД, причем все сигналы, превышающие порог Р (см. рис. 10), полагались равными этому порогу. Порог Р выбирался в верхней части линейного участка амплитудной характеристики АОИЧ. После его установки вычисление частоты проводилось по описанному алгоритму.
Погрешности вычисления частоты для этого алгоритма, как и для предыдущего, оценивались по 50 реализациям РИСС в полосе АРА ~ 6ДР. Частота сигнала от реализации к реализации увеличивалась с дискретностью 65 кГц. Графики погрешностей для режимов ограничения (см. рис. 10, кривые 1, 2, 3) представлены на рис. 11. Все графики погрешностей включали в себя систематическую составляющую. В режимах 1 и 3 дополнительно наблюдалось влияние на погрешность дрейфа частотной шкалы.
Систематическая погрешность объясняется асимметрией формы сигнала и различной (из-за асимметрии) его крутизной на фронте и на спаде. Различие крутизны на этих участках привело к дополнительной погрешности вычисления частоты. После исключения систематической составляющей погрешность составила около ±110 кГц для кривой 1, ±70 кГц для кривой 2 и ±55 кГц для кривой 3.
При определении частоты сигнала как полусуммы частот крайних ФД АОИЧ, освещенных РИСС, погрешность единичного измерения не превышала половины частотной дискреты АР [4], [5]. При использовании более сложных алгоритмов (например, аппроксимации совокупности сигналов РИСС, снимаемых с к фотодиодов, некоторой известной кривой, координату экстремума которой отождествляют со значением измеряемой частоты) по теоретическим оценкам может быть реализована погрешность, примерно равная 0.004 АР .
Предложенные алгоритмы для определения положения РИСС требуют обработки информация с четырех ФД. Алгоритмы обеспечивают возможность измерения частоты и (потенциально) амплитуды сигнала при работе АОИЧ как в линейном, так и в нелинейном режимах, при этом максимальное теоретически оцененное значение погрешности измерения частоты для линейного режима, кото-
рую можно считать предельно достижимои, составляет величину порядка 0.01 АР . Однако как показывают экспериментальные исследования, искажения формы аппаратной функции АОИЧ не позволяют реализовать потенциальные возможности предложенных алгоритмов вычисления частоты.
Известия вузов России. Радиоэлектроника. 2008. Вып. 4======================================
Вместе с тем, даже при наличии искажающих факторов предложенные алгоритмы обеспечивают возможность измерения несущей частоты радиосигналов с погрешностью менее 10 % от величины частотной дискреты при работе АОИЧ в линейном и менее 20 % при работе в нелинейном режимах. Поэтому они могут быть рекомендованы для практического использования, причем как в АОИЧ, так и, по-видимому, в других измерителях.
Библиографический список
1. Роздобудько В. В. Широкополосные акустооптические измерители частотных и фазовых параметров радиосигналов // Радиотехника. 2001. № 1. С. 79-92.
2. Роздобудько В. В., Пелипенко М. И. Быстродействующий измеритель параметров СВЧ радиосигналов // Специальная техника. 2006. № 1. С. 28-36.
3. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. М.: Наука, 1978. 832 с.
4. Роздобудько В. В., Дикарев Б. Д. Высокоточный акустооптический приемник-частотомер комбинированного типа // Радиотехника. 2003. № 9. С. 31-36.
5. Grossman S. B., Emmons R. B. Performance analysis and size optimization of focal planes for point-source tracking algorithm applications // Opt. Engineering. 1984. Vol. 3, № 2. P. 167-176.
B. N. Wolfovsky, S. S. Shibaev, V. V. Rozdobudko Technology institute of the South federal university
Frequency computation algorithms in acousto-optic measurers of radiosignals parameters
Two heuristic algorithms offrequency computation, as abscissa of symmetry axis of video pulse, that is pickupped with photo receiver in acousto-optic measurers of radio signals parameters are offered and analyzed. Theoretical and experimental errors estimations are given. Experimentally evaluated error of frequency computation in enlarged dynamic range of input signals does not exceed 10 % of frequency grid step.
Algorithm, error, acousto-optic measurer, intensity distribution
Статья поступила в редакцию 14 апреля 2008 г.
