книги / Теория оптико-электронных следящих систем
..pdfпрозрачной части модулятора, когда выполняется условие |
|
*м/ > 0. |
(3.2.47) |
В этом случае вычисляется сигнал от данной точки g(xif уу ) в соответствии с принятой аппроксимацией ФРТ и засылается в сумматор, моделирующий выражение (3.2.37). В противном случае, при хМ1- < 0, точка находится на непрозрачной части модулятора и сигнал от нее равен нулю (не вычис ляется) .
Таким образом проводятся вычисления для всех точек изображения в каждый момент времени t. Аналогичные вычисления повторяются для следующего момента времени с шагом At. Индексы /, / в (3.2.45) опреде ляются размерами изображения. Для объектов заданной геометрической формы удобнее производить перебор точек, изменяя координаты х, у . от
*min до |
*тах |
и от y min до у тах |
с шагом А. Например, для круглого |
||||
изображения |
|
|
|
________ |
|
||
Утах |
~ |
'О, |
*тах |
“ |
\ / г 0 |
—у 2, |
(3.2.48) |
Уmin |
= ~ гО, |
xmin |
= —V Г% ~ |
У2. |
|
Подробнее рассмотренный алгоритм представлен на рис. 3.9. Для опреде ленности распределение освещенности внутри кружка рассеяния аппрокси мируется гауссоидой (3.2.16). Как видно из блок-схемы, программа состоит из ряда проверок попадания точек изображения за границу поля зрения, границу изображения, границу модулятора. Аналогичные алгорит мы получаются для моделирования анализаторов изображения других ти пов, поэтому в дальнейшем не будем приводить полную блок-схему, а огра ничимся рассмотреннием основных моментов при моделировании анализа торов других типов.
В р а щ а ю щ и й с я р а д и а л ь н о - щ е л е в о й м о д у л я т о р . Анализатор изображения данного типа приведен на рис. 3.10. Здесь произ водится модуляция потока с двойной частотой (двойная амплитудная модуляция), что позволяет произвести пространственную селекцию точеч ного объекта [67].
Алгоритм проверки попадания точки изображения а/у на прозрачную или непрозрачную часть модулятора распадается на две задачи. Во-первых,
Рис. 3.10. Радиально-щелевой модулятор
191
необходимо определить, на какую половину модулятора (непрозрачную или полупрозрачную) попадает выбранная точка. Во-вторых, в случае попадания на полупрозрачную половину нужно определить, в каком секто ре (прозрачном или непрозрачном) находится точка.
Решение первой задачи производится по алгоритму, рассмотренному выше. Вторая задача решается путем сравнения угла otMij в подвижной системе модулятора с величиной углового размера сектора. Величина угла равна
о-м// = arctg ( y Mj / x Mi). |
(3.2.49) |
Пронумеруем секторы от 0 до и (на рис. 3.9 п = 3) в каждой из четвертей полупрозрачной половины.
Операция выделения целой части |
|
|
|
||
|
|
а |
|
|
|
лг |
= int Iе ) |
|
(3.2.50) |
||
дает |
номер |
сектора. Дополнительная |
проверка |
на знак y Mj дает |
номер |
квадранта |
(1 или IV), в котором |
находится |
рассматриваемая |
точка. |
В итоге точка оказывается на прозрачной части модулятора в нечетных секторах первого квадранта и в четных секторах IV квадранта, т.е. при выполнении условий
у М1> 0 , |
п |
нечетное, |
(3 2 51) |
y Mi < 0, |
п |
четное. |
|
Для этой точки вычисляется величина сигнала £(*/,.У/), засылается в сум матор и производится переход к анализу следующей точки изображения. После просмотра всех точек изображения вычисляется сигнал на выходе анализатора изображения
и = S A xA y HgiXi.yj). |
|
|
(3.2.52) |
А н а л и з а т о р с |
п е р е н о с о м и з о б р а ж е н и я . В этом АИ |
||
(рис. 3 .1 1 ) изображение |
с помощью вращающегося наклонного зеркала |
||
|
переносится по плоскости |
||
|
неподвижного модулято |
||
|
ра [73]. Если объект нахо |
||
|
дится |
на оптической оси, |
|
|
центр окружности перено |
||
|
са радиуса Rn совпадает с |
||
|
центром |
модулятора, а |
|
|
смещение |
объекта вызы |
|
|
вает смещение центра Ом |
||
|
Окружности переноса. |
||
|
При моделировании та |
||
|
кого |
API |
целесообразно |
Рис. 3.11. Анализатор с пере носом изображения
192
перейти от подвижной системы координат Омх м у м к неподвижной Ок хк у к . В произвольный момент времени координаты центра изображе ния в подвижной системе равны
Хц |
-^п cos comfit, |
Ущ |
-^п sin |
(3.2.53) |
|
Координаты центра изображения в системе координатора Окх^ук |
|
||||
Х К. |
~ * Ц |
"*■^ М ) |
УК ” |
Уц + Ум. 5 |
(3.2.54) |
где хы, Ум |
— координаты |
центра окружности переноса при |
смещении |
||
объекта от оптической оси. |
|
|
|||
Координаты точки изображения ац в системе Окхку к |
|
||||
xKi |
~ хъс ■*" |
Ук.] ~ Ук “*■/А. |
(3.2.55) |
Как видно из рис. 3.11, плоскость АИ можно разбить на внешнюю, коль цевую, зону и центральную, границей между которыми является окруж ность радиуса / 3. Поэтому анализ положения точки ау состоит из трех про верок: попадания в плоскость АИ, попадания в кольцевую зону и попада ния в центральную зону АИ. Попадание в плоскость АИ задается условием
I = -у/х2 + у 2 < /4 , |
(3.2.56) |
где для сокращения записи обозначено xKi =х>ук/- - у . Попадание в коль цевую зону происходит при выполнении неравенств
/з < / < /4)
а попадание в центральную зону происходит при / < / 3.
Если точка изображения ац попадает в кольцевую зону, то производится
сравнение ее углового положения |
|
|
| a if I = arctg У |
|
(3.2.57) |
х |
|
|
последовательно с углами <pl9 $ 2 , |
$ 4° Выполнение условия | |
| < <рп |
определяет номер кольцевого сектора, который имеет прозрачность 0 или 1. При попадании точки на прозрачный сектор вычисляется сигнал g(x, у) (здесь х, у —координаты в системе Оху) , который засылается в сумма тор. Далее производится переход к следующей точке изображения.
При попадании точки в центральную часть АИ производится проверка
координат x Ki, |
xKj |
на величину отклонений l t |
и / 2. Точка попадает на |
прозрачную часть, если одновременно выполняются условия |
|||
h < U K I I < |
/2, |
h < Ij'K/ I < /2 . |
(3.2.58) |
Таким образом, все проверки можно разбить на две группы: проверка по углу и проверка по величине отклонения. В связи с симметричностью рисунка модулятора все проверки можно производить для одной четверти, взяв все углы и координаты по абсолютной величине.
М а т р и ч н ы й а н а л и з а т о р и з о б р а ж е н и я . Анализатор изо бражения, выполненный в виде матрицы фотоприемника, приведен на
13. Ю.М. Астапов |
193 |
Рис. 3.12. Матричный анализатор изображения
рис. 3.12. Для определенности изобра жена матрица размера 6 X 6 элементов. Частным случаем такого фотоприемни ка является четырехплощадочный фо топриемник. Изображение, как и ра нее, рассматривается в виде кружка рассеяния; координаты элементарной точки суть х { = хц + /А, у; = у п + /Д .
Можно построить алгоритм после довательной проверки попадания точ ки изображения в элементы матрицы,
однако такой алгоритм содержит большое число вычислений, особенно при большой размерности матрицы. Поэтому здесь целесообразно построить ал горитм, позволяющий сразу указать номер элемента, в который попадает точка изображения.
Назовем один из элементов матрицы, примыкающий к началу коорди нат 0 К, центральным. Тогда номер элемента, в который попадает точка изображения ау с координатами xif уу , определяется выражением
N = АГц.э |
+ Ny + NNX. |
(3.2.59) |
Здесь |
— номер центрального элемента; |
N —количество элементов |
в строке; Nx, Ny —количество целых элементов, укладывающихся на со ответствующих отрезках:
Nx = int |
(3.2.60) |
где Дм —размер элемента матрицы.
При использовании этого алгоритма необходимо учитывать, что цен тральный элемент примыкает к началу координат справа и сверху.
Поэтому при положительных координатах |
(х* > 0, yj > |
0) |
в (3.2.59) |
||
подставляются Nx и Ny непосредственно из |
(3.2.60), а при отрицательных |
||||
координатах. |
|
|
|
|
|
Nx = |
int |
1, |
|
|
(3.2.61) |
Такой |
алгоритм |
дает большую экономию вычислений |
по |
сравнению |
с последовательным просмотром. На каждый элемент матрицы отводится свой сумматор Xg(xf, y j ), в который заносится значение яркости попав шей на этот элемент точки.
Т е л е в и з и о н н ы й а н а л и з а т о р и з о б р а ж е н и я . В данном анализаторе производится последовательный просмотр всего поля зрения по определенной траектории сканирования. При использовании обычного телевизионного стандарта это линейная развертка. В связи с большим
194
объемом информации, содержащейся в телевизионном кадре, моделиро вание полного телевизионного сигнала весьма затруднительно в связи с необходимостью больших вычислительных затрат. Поэтому одной из основных задач при математическом моделировании является воспроизве дение телевизионного сигнала с разумным ограничением избыточной информации.
Свяжем количество строк телевизионной развертки при моделировании с верхней воспроизводимой частотой телевизионного сигнала. Длительность одной строки равна
Гс |
(3.2.62) |
где п - количество строк разложения в модели; Гк —период кадровой развертки.
При одинаковом количестве элементов разложения по строкам и столбцам шаг по времени между двумя элементами разложения
ДгР |
Т± |
Тк |
(3.2.63) |
|
п |
|
|||
|
|
|
||
По теореме Котельникова |
[135] верхняя частота сигнала связана с ша |
|||
гом счета выражением |
|
|||
|
1 |
|
(3.2.64) |
|
/в = 2 A t p |
2 Тк |
|||
|
Отсюда получаем необходимое количество элементов разложения в строке при заданной верхней частоте телевизионного сигнала
п |
(3.2.65) |
где F K = 1/Гк —кадровая частота.
Обычно в отико-электронных системах не используется .полная ширина спектра телевизионного сигнала, поэтому размер п X п матрицы разложе ния кадра получается в разумных пределах. Например, при / в = 300 кГц, FK = 25 Гц получим п = 155. Тогда шаг дискретизации по координатам равен
Ахр = А ур = Atp Vp ~ — , |
(3.2.66) |
где Vp —скорость развертки; L — размер мишени телевизионного анализа тора. Фильтрующее свойство развертывающего элемента характеризуется его импульсной характеристикой gp( х , у ) , а траектория движения задается координатами хр(Г), y p(t). Если форма развертывающего элемента не из меняется в процессе развертки (пространственно-инвариантная импульсная характеристика), то сигнал на выходе телевизионного анализатора описы-
1 3 * |
195 |
вается выражением |
|
||||
«О ) = u[x(t),y(t)] = |
|
||||
= SS Ви(х, y)gp [х —Xp(t), у —У р ( 0 ] dxdy, |
(3.2.67) |
||||
где Ви(х, у) |
—изображение, сформированное в фокальной плоскости опта- |
||||
ческой системы. |
|
|
|
||
Для перехода к дискретной форме положим |
|
||||
х = iАх, |
у |
= jAy, хр = аАхр = kVpxA t , |
(3.2.68) |
||
Ур = РДУр |
= ^ VpyAt. |
||||
|
|||||
Тогда можно записать в матричной форме |
|
||||
u(kAt) —2 |
2 |
gptoi' jpbnijу |
(3.2.69) |
||
|
а. |
(3 |
|
|
|
гДе &pia, / р, |
b<nij |
—элементы матриц развертки изображения. |
|
Здесь необходимо отметить два обстоятельства. Во-первых, полезно сравнить ширину импульсных характеристик объектива и развертывающего элемента. Если разрешение одного из рассмотренных элементов преобразо вания сигнала значительно выше, чем другого, то сигнал на выходе будет определяться в основном тем элементом, разрешающая способность кото рого ниже. Поэтому при моделировании можно учитывать только этот эле мент. Например, при хорошей оптике и грубой развертке (например, мозаичный фотоприемник) объектив можно считать идеальным и не учиты вать его фильтрующие свойства. В этом случае математическая модель становится более простой.
Во-вторых, шаги дискретизации Ах, Ау, Ах р , Аур изображения и развертки в общем случае не равны друг другу (обычно Ар < А). Поэтому при моделировании телевизионного анализатора приходится при менять интерполяцию функции двух переменных для нахождения проме жуточных значений. Если изображение задано в дискретных точках
В„(хо + г А х , у 0 +кАу) |
= f ik |
0\ к = 0, ±1, . . . ) , |
|
|
а изменения аргумента равны |
|
|
||
а = |
*о |
У ~ Уо |
|
(3.2.70) |
0 = |
Ау |
|
||
|
Ах |
|
|
|
то, например, формула билинейной интерполяции имеет вид [148] |
|
|||
f i x , у) |
= [/о о (1 -а ) + a / i o ] ( l - « + [/o i(l-< * ) + <*/ii]0. |
(3.2.71) |
Таким* образом, в результате рассмотренных преобразований получаем модель телевизионного сигнала (видеосигнала) в дискретные моменты времени.
Моделирование электрического сигнала на выходе фотоприемного устройства. В рассмотренные анализаторы изображения входит как состав ная часть приемник излучения, а в некоторых случаях, например в телеви зионных системах, сам фотоприемник (телевизионная трубка) является
196
анализатором изображения. Поэтому при моделировании необходимо учи тывать характеристики фотоприемников. На выходе приемника излучения получается электрический сигнал, а в математической модели —одномер ная функция времени, которая с учетом шумов фотоприемника является случайным процессом.
Х а р а к т е р и с т и к и п р и е м н и к о в и з л у ч е н и я . Свойства приемника излучения наиболее полно могут быть описаны системой харак теристик, выражающих зависимость сигнала и шума, вырабатываемых приемником, от различных* факторов: величины, спектрального состава и частоты модуляции излучения, падающего на приемник, температуры окружающей среды, полосы пропускания усилителя и т.д.
Наиболее распространенными характеристиками приемника излучения являются: амплитудная, спектральная, частотная, шумовая [81]. При про ведении математического моделирования не всегда имеется полный набор характеристик, поэтому во многих случаях характеристики заменяются числовыми параметрами, выражающими свойства приемника для опреде ленных условий. В соответствии с принципом функционального моделиро вания будем рассматривать только те характеристики и параметры фото приемников, которые позволяют описать его как преобразователь потока излучения в электрический сигнал. С этой точки зрения приемник можно рассматривать как четырехполюсник, имеющий крутизну преобразования S , которая определяет величину сигнала, вырабатываемого приемником и приходящегося на единицу падающего на него потока излучения. Крутизна преобразования обычно называется чувствительностью и в зависимости от типа фотоприемника выражается в различных единицах (вольт на ватт, ом на ватт, ампер на люкс и т.д.). Рассмотрим подробнее важнейшие характеристики.
Амплитудная характеристика (энергетическая или световая) дает зави симость сигнала от величины потока излучения, падающего на приемник, т.е. и - / ( Ф). Эта зависимость является нелинейной, поэтому описать амплитудную характеристику одним значением крутизны преобразования невозможно. Различают три значения крутизны (чувствительности).
Статическая крутизна или крутизна для смодулированного сигнала
SCT = X . |
(3.2.72) |
Ф |
|
Цифференциальная крутизна или чувствительность для модулированного сигнала
*^д |
d u |
(3.2.73) |
|
|
аФ |
Средняя крутизна S cp определяется отношением амплитуды первой гар моники сигнала с фотоприемника к амплитуде синусоидально модулиро ваниюго потока излучения
Аи |
(3.2.74) |
*^ср ” . - • |
197
Значение крутизны преобразования зависит от уровня фоновой засветки (смодулированного потока излучения). При низких уровнях постоянной составляющей, т.е. при работе на линейном участке амплитудной характе ристики,все три типа крутизны преобразования совпадают. Крутизна преоб разования измеряется при интегральном облучении от абсолютно черного тела или другого эталонного источника излучения, поэтому является интегральным параметром.
Для описания зависимости сигнала от длины волны падающего излуче ния вводится спектральная крутизна
5( Х) = |
du |
(3.2.75) |
---------- , |
ФэxdX
где X - длина волны; Фэ\ —спектральная плотность потока излучения. Интегральная и спектральная крутизна связаны между собой со
отношением
оо
/ Ф Э*5(Х)</Х |
|
S = —-------------------. |
(3.2.76) |
оо |
|
о
Формула (3.2.76) дает значение крутизны преобразования приемника при облучении его потоком сложного спектрального состава.
Частотная характеристика показывает зависимость дифференциальной крутизны от частоты модуляции потока излучения 5 Д = ( / м ). Эта харак теристика отражает динамические свойства приемника —его способность реагировать на быстрые изменения потока излучения. В наиболее простом случае динамическим эквивалентом приемника является апериодическое звено с частотной характеристикой
ИЧ/со) = т А ~ ^ ’ |
(3‘2*77) |
1 + /соТ |
|
где Т —постоянная времени фотоприемника.
Вообще говоря, процессы, происходящие в реальных фотоприемниках, более сложны, и дальнейшее развитие модели (3.2.77) при функциональ ном моделировании заключается в введении двух постоянных времени Тх и Т2 для возрастания и убывания потока излучения.
Ш у м ы п р и е м н и к о в и з л у ч е н и я . Одной из важнейших ха рактеристик шумов фотоприемников является энергетический спектр шума. Однако в большинстве случаев математического моделирования эта характеристика отсутствует, имеются сведения лишь о среднеквадратиче ской величине шума в определенной полосе частот. Поэтому о спектре шу ма судят зачастую на основании общих сведений о его природе, позволяю щих указать распределение мощности шума по частотам для приемника того или иного типа с учетом условий эксплуатации.
1 9 8
Основными видами шумов приемников излучения являются: тепловой, дробовой, токовый, генерационно-рекомбинационный. Рассмотрим мате матическое описание перечисленных шумов [39, 65, 81, 119].
Тепловой шум обусловлен хаотическим движением носителей заряда и имеет равномерный спектр до частот порядка 1011 Гц. Спектральная плот
ность теплового шума описывается формулой |
|
Nr ( f ) = 4kTR В2/Гц, |
(3.2.78) |
где R — сопротивление, ом; Т — абсолютная температура, |
К; к = |
= 1,38-1(Г2* Дж-К*"1 —постоянная Больцмана. |
|
ЭДС эквивалентного генератора шума в полосе частот А/ равна
ет = V NT( f ) A f . |
(3.2.79) |
Дробовой шум определяется тем, что электрический ток представляет собой поток дискретных частиц, заряд которых равен или кратен заряду электрона. Этот шум проявляется в основном в электровакуумных прибо рах. Спектральная плотность тока постоянна и определяется формулой
/ ДР( П |
= 2 eJ0 а2 / Гц, |
(3.2.80) |
где е —заряд электрона; J 0 —уровень постоянной составляющей. |
|
|
Величина тока эквивалентного генератора шума равна |
|
|
/др |
yj 2 еJQAf . |
(3.2.81) |
Токовый шум зависит от протекающего через полупроводник тока (природа возникновения этого шума не совсем ясна). Спектральная
плотность выражается формулой |
|
|
const |
|
|
N r A f ) = |
0 |
(3.2.82) |
/ |
|
где параметр 0 = 0,8-М,5 и определяется экспериментально для каждого типа фотоприемника. Токовый шум сказывается в основном на низких частотах. ЭДС эквивалентного генератора шума равна
^тк |
[ / ^ т к ( / ) ^ / ] |
(3.2.83) |
|
Л |
|
Генерационно-рекомбинационный шум связан с флуктуациями процессов генерации и рекомбинации носителей зарядов в полупроводниках. Он аналогичен дробовому шуму в вакуумных приборах. Спектральная плотность равномерна в широкой полосе частот и определяется опытным путем для каждого типа полупроводника.
Из рассмотренных шумов для полупроводниковых фотоприемников на низких частотах преобладает токовый шум, на средних —генерационно-ре комбинационный, на высших частотах —тепловой шум. Суммарная спек тральная плотность приведена на рис. 3.13.
199
Обычно спектральная плотность шумов фотоприемника неизвестна или задана довольно приблизительно. Поэтому при моделировании поступают следующим образом. Для реального фотоприемника измеряется средне квадратическое значение шумов на выходе усилителя с полосой А/ при подключенном на входе фотоприемнике. Полученное значение пересчиты вается ко входу усилителя по формуле
^ ш . в х |
мш. вых |
(3.2.84) |
|
к |
|||
|
|
Таким образом учитываются также шумы первых каскадов усилителя. Величина пропорциональна площади под кривой спектра в полосе про пускания усилителя (см. рис. 3.13). В математической модели ОЭС шумы фотоприемника вводятся аддитивно, т.е. прибавляются к сигналу, полу ченному после анализатора изображения с учетом крутизны преобразова ния S. Шумы задаются с датчика случайных чисел с нормальным законом распределения, равномерным спектром в полосе частот А / (дискретный белый шум) и дисперсией о2 = Мщ ВХ. Методика генерации на ЭВМ слу чайных процессов рассмотрена в § 3.5.
В ряде случаев приходится учитывать неравномерность чувствительности по полю фотоприемника и задавать ее как функцию координат S(x, у). Обычно чувствительность является случайной функцией координат, т.е. случайным полем. Методы генерации случайных полей на ЭВМ также рас смотрены в § 3.5. Здесь отметим только, что неравномерность чувстви тельности по полю эквивалентна мультипликативному шуму, так как сигнал умножается на случайную крутизну преобразования.
Рис. 3.13. Шумы приемников излучения: 1 - генерационно-рекомбинационный; 2 - токо вый; 3 —тепловой: Еш - спектральная плот ность шума
Статистические характеристики сигнала на выходе сканирующей оптиче ской системы. Как отмечалось выше, теоретически вычислить сигнал на вы ходе анализатора изображения, как правило, не удается и приходится применять численные методы. Однако в ряде случаев можно определить аналитически статистические характеристики сигнала на выходе анализато ра изображения и, таким образом, избежать прямого моделирования АИ. При необходимости можно затем произвести генерацию случайного про цесса по известным статистическим характеристикам.
200