Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Теория и расчет импульсных и цифровых оптико-электронных систем

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
18.87 Mб
Скачать

Во втором способе оценки качества по точности работы рас­ сматривается установившаяся ошибка при гармоническом входном воздействии

s (0 = smsin

(10.15)

где sm— амплитуда входного непрерывного управляющего си­ гнала.

Если это воздействие представить решетчатой функцией

s [п] — smsin щпТ,

(10.16)

то верхнее граничное значение амплитуды ошибки в установив­ шемся режиме определится формулой

= SmIФ* (е'“ лТ, е)|,

(10.17)

где |ФЯ (е/(°ог , е)| — модуль дискретной частотной передаточной функции на частоте входного сигнала ю0.

Для замкнутых импульсных систем

Ф..с (е/®*г, е) =

1/(1 +

К

в)],

(10.18)

а для разомкнутых —

 

 

 

 

 

Ф*

е) = /См [1 -

К (е'“°г, е)],

(10.19)

где /См и К' (е'40»7, е)

имеют тот же смысл, что и в формуле

(10.14).

В некоторых случаях при проектировании импульсных систем

задаются лишь максимальные

скорость

sraax и ускорение

smax,

с которыми может

изменяться

входная

величина s (t).

В

этом

случае предварительная оценка точности работы системы может быть произведена по величине ошибки в установившемся эквива­ лентном гармоническом режиме.

Известно, что при гармоническом сигнале справедливы соот­ ношения :

s (t) — smsin со/;

s (t) = ds (t)/dt -f- smсо cos со/ = smcos to/;

s (t\= d2s (t)ldt2 = snlсо2 sin со/ = smsin cat.

Следовательно, если заданы максимальные значения smax, smax, то частота и амплитуда эквивалентного гармонического режима определяются выражениями:

 

(0Э= ^шзх/^щах»

(10.20)

5э ш

Siпах/^э — Smax/Sm a x *

( 10.21)

Подставив эти параметры в формулу (10.17), найдем

хэ m &эт

J

( 10.22)

е' э

Это и есть амплитуда ошибки в установившемся эквивалентном гармоническом режиме, при котором скорость и ускорение измене­ ния входного сигнала равны заданным максимальным значениям.

171

Глава 11

СЛУЧАЙНЫЕ ПРОЦЕССЫ И ИХ ПРОХОЖДЕНИЕ ЧЕРЕЗ ИМПУЛЬСНЫЕ ЛИНЕЙНЫЕ СИСТЕМЫ

§11.1. ОСНОВНЫЕ ХАРАКТЕРИСТИКИ СЛУЧАЙНЫХ ПРОЦЕССОВ

Если значение некоторой физической величины f (/) меняется непрерывно по случайному закону, то говорят, что она представ­ ляет собой случайный процесс. Такие процессы могут быть ста­ ционарными и нестационарными. В дальнейшем будут рассматри­ ваться только стационарные процессы, т. е. такие процессы, сред­ нестатистические характеристики и параметры которых не зависят от начала отсчета и не изменяются вдоль оси времени. Для непре­ рывных стационарных случайных процессов основными характе­ ристиками являются перечисленные ниже.

Одномерная плотность вероятности W (/), характеризующая вероятность распределения значений функции / (0 в любой точке на оси L

Двумерная плотность вероятности W (Д, /2, t2— А) — № (Д, /2> т), определяющая среднестатистическую связь между значе­ ниями fi — f (А) и /2 = / (У» разделенными на оси времени интер- 'валом т = /2 — h-

Корреляционная функция

/ 2) — В (т) — | h (0h it + Т)w Ifг (0, h (t + Т), т] dh dh,

ОО

(11.1)

описывающая изменение математического ожидания произведения двух значений случайной функции в зависимости от величины интервала Utx = т.

Благодаря эргодичности стационарного случайного процесса, корреляционная функция может быть найдена также осреднением одной достаточно протяженной его реализации

Bf (т) = lim

/ (0 f (t + Т) dt.

(11.2)

Т-*оо

 

 

Спектральная плотность мощности случайного процесса пред­ ставляет собой преобразование Фурье от корреляционной функции

+ 0 0

(?/(©)= j Bf (т) е~'шт dxt

(11.3)

-0 0

 

которому соответствует сопряженное преобразование

 

+оо.

 

В>М = 4т 1 0/(“)e/B,do)-

(П.4)

-00

 

172

Наиболее часто используемыми параметрами непрерывного случайного стационарного процесса являются математическое ожидание и дисперсия, которые определяются следующими выра­ жениями:

+оо

 

 

Т

!} (1) di =

/

В

; (т =

оо);

< Г ) =

f ( t J) W { f ( t ) ) d f = l i m ^ r

-оо

 

 

 

 

 

 

 

 

 

 

 

+ 0 0

 

 

 

 

 

 

 

о , =

а / (о - (пг-) = } if if) -

ini2 w if m df=

 

 

 

T

-oo

 

 

 

 

 

 

 

 

 

В,(т = 0) — Bf(i =

oo).

 

«= lim ф - f If(0 -

(/)]"-dl =

(11.6)

T-+oo

_JT

 

 

 

 

 

 

 

Если

непрерывный

случайный

процесс

 

f {t)

отсчитывать

в дискретные моменты времени с интервалом Т, то в результате будет получена случайная решетчатая функция f IпТ] = / [я]. При стационарном характере образующей случайной функции f (t) случайная решетчатая функция f In] также будет стационар­ ной и для оценки ее свойств могут быть использованы те же харак­ теристики и параметры, что и для непрерывных процессов:

одномерная плотность вероятности W (/);

двумерная плотность вероятности W (/х In], /2 In], пг] =

=W (/ь /2, т ); корреляционная функция

 

 

 

 

+ 0 0

 

 

 

Bf [пТ, (n +

т ) Т] = Bf [пг] = JJ /1/ 2^

(/ъ /а>

nz) rf/i df2 —

 

 

 

 

N

 

 

 

 

~

 

'2W + 1 '

S n « ] / t « + mi;

(и .7)

 

 

 

 

n = —N

 

 

 

математическое ожидание

 

 

 

 

 

+oo

 

 

 

+N

 

 

(/> =

J/ [ft] ^ (/ M ) df = jirn

-jNl+ t

^

f [n] =

1/ Bf [m = oo];

 

-0 0

 

 

 

n=-?-N

 

(11.8)

дисперсия

 

 

 

 

 

 

 

 

 

 

 

+00

 

 

 

 

УУ

(f[n] - (f)f =

0/= J(fin)- if))*Wif W)df =

 

2

-00

 

 

 

 

Г2=—iV

 

 

=

Bf [m =

0] — Bf [tn — 00].

(11.9)

Спектральная

плотность случайного стационарного •процесса,

согласно работе

[2],

может

быть

представлена

в нескольких

формах

записи:

 

 

 

 

 

 

m

в виде двустороннего 2-преобразования корреляционной функ­ ции

 

 

т=+оо

 

 

< ?(*)-

Е Bf [т] z~m=

F (2) - f F (2-1) — Bf [0],

(11.10)

 

 

т=—оо

 

 

где F (2) — 2-преобразование от корреляционной функции;

в виде функции от круговой частоты

 

 

4 -0 0

 

 

 

G(е'“ г ) =

S

[т] ег!ытТ= F (ЫшТ) + F (е-/шГ) Bf [0].

(11.11)

m = — 00

 

 

 

Если учесть, что корреляционная решетчатая функция явля­

ется четной, то можно записать

ОО

 

 

 

 

 

 

G(е'шГ) = В, [0] + 2

S В, [ml cos amT.

(11.12)

 

 

m=1

 

Обратное преобразование имеет вид

 

 

л/Т

 

л/Т

 

Bf [m] = - ~

J G(e/w7’)e/Mm7’ rfco =

-^- j G(&<*T)cos ®mTda.

(11.13)

 

-л/т

 

о

 

При этом средний квадрат решетчатой функции, рассматривае­ мой в дискретные моменты времени, определяется интегралом

 

л/Т

 

(Р[п, е]) =

j G(&oT, e)dco.

(11.14)

 

о

 

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

Для получения такой зависимости сначала необходимо получить

спектральную плотность С/ (г) согласно формуле (11.10),

а затем

сделать подстановку

г =

(1 4* а>)/(1 — ш)|и,-дг/2.

В результате

будет

получена функция

 

 

 

 

 

G^ = G i - ^

m

) = 2 R e F ( ^

w k ) ~

B>l0]-

(11Л5)

Здесь

F обозначено

2-преобразование

корреляционной

функции

с последующей подстановкой.

 

 

 

Обратный переход от спектральной плотности к корреляцион­ ной функции может быть осуществлен согласно выражению

+00

G (X) dX

(11.16)

11 + (ХТ/2 р

 

00

 

При m — 0 из этой формулы можно получить выражение для расчета среднего квадрата случайного решетчатого процесса

174

в виде

 

 

 

 

 

 

 

</■ И ) = (772л)

х

 

 

+0 0

 

 

 

 

 

 

X

I G(X) dX/11 + jXT/212.

 

 

 

 

 

 

 

(11.17)

 

В

табл.

11.1

приведены

корреляционные

 

функции

и

спектральные

плотности

некоторых

наиболее

часто

встречающихся

 

стационар­

ных

случайных

решетчатых

процессов.

 

 

 

 

 

Величина среднего ква­

драта

случайного

решетча­

того

процесса,

представлен­

ная формулой (11.17),

обыч­

но

может

быть

записана

в следующей форме:

 

 

 

(1Чп])=Т1п.

(11.18)

 

Подынтегральная

функ­

ция в интеграле

 

 

 

 

 

 

+ 0 0

 

 

 

 

у

_

1

г

G ( X ) c a

_

 

л

J

|1-ЬД 7У 2| 2

 

 

 

- 0 0

 

 

 

 

 

_

1

+ 0 0

В (%)d%

 

 

Г

 

 

 

J

А (А ) А (—/Я)

-о о

(11.19)

может быть представлена от­ ношением полиномов, кото­ рые определяются согласно выражениям:

В(Х)=Й0(А)2(П-■>+

+ &1 (А)2 <"-2>+ • • • + ЬЛа

A (jX) = а0 ЦХУ +

+fli (jX)”-' Н---------Ь ап.

(11.20)

Если процесс является устойчивым, то интеграл (11.19) находится по значе-

VO

- ^

К

 

eg

 

(0772)

 

Ь

 

£

 

N

 

S EN

 

 

С1

 

sina—

 

 

0

а

II

•-

С5

 

^

-2.

s ^

 

<N

со

сч

 

Cl

“Ь ^

-

 

со

 

 

 

 

сГ

 

 

+

+ +

 

 

 

а(772)

С4 Е- £

 

 

 

 

CI

 

 

 

 

С1

SZ

■S

£

 

1

+ sh

 

+ 3=

 

 

1

оо

 

 

 

 

 

1

н

о

о

н

к,

иО

3

Р

^ I

О

Е.3

о

4) N

3

2 а«=:

нс .

с25

U

+W|

03

2 g

о

2оо Ю

к 2 4 3

о £

0Q

к

з

со Cl^3 О

CN

04

тз

+

к

8 IО

N 43N

сГ

X*

N

8

Q

S

0Q

8^ S

CCL

оо

I

Is—- 3

'—к

Q

17S

ниям коэффициентов а, и b,t входящих в полиномы

(11.20). В ча­

стности, при

1 < п с

4 эти интегралы определяются

по формулам:

т _

bp .

j _

а0&1 — Дг^о .

j _

Дз (QQ^I ~ а2^о) — ДрД^г »

1

2а0а, ’

2

2a{sa1a2

'

3

2а0Дз (ДрДз — Д1Д2) *

 

Д40 (Д2Д3 —Д1Д4) —Дц(Дз^1

Д^г)! ~НДр^з (ДрДз —Д1Д2)

 

 

2Д|)Д. (Д0Д5 +

Ф л

— а 1а2ал)

 

( 11.21)

 

 

 

 

Общее выражение для вычисления интегралов (11.19) по вели­ чине коэффициентов at и bt при любом п можно найти в работе [2 ].

§ 11.2. ПРОХОЖДЕНИЕ СЛУЧАЙНОГО ИМПУЛЬСНОГО СИГНАЛА ЧЕРЕЗ ЛИНЕЙНУЮ ДИСКРЕТНУЮ СИСТЕМУ

Если на вход разомкнутой импульсной системы, имеющей частотные передаточные функции К (е'и)Г) и К (IXT/2), поступает стационарный случайный решетчатый сигнал х [п] со спектраль­ ной плотностью Gx (е/“ г) или Gx (X), то спектральная плотность случайной выходной величины у [/г] определяется формулами:

Gy(У<*т) =

|К (&аТ) |2 Gx (е/®7');

(11.22)

Gy (X) =

|К (/Ш2>|2 Gx (X).

(11.23)

Воспользовавшись формулами (11ЛЗ) и (11.16), можно найти

корреляционную функцию

Ви [m] выходного

случайного

про­

цесса.

значение выходной

величины

у Iп]

Среднее квадратическое

определяется интегралами:

 

п/т

 

 

{У2М) =

j

|К (е'аТ) Р Gx (е'“П А»;

(11.24)

 

- Я

 

 

 

 

+ 0 0

1к a m 2) |« GX(X)

 

 

{У2[п}) = ~

j

dx.

(11.25)

 

 

11 + j%T/2 I*

 

 

Если импульсная система является замкнутой и к ее входу приложено случайное управляющее воздействие s (/) со спектраль­ ной плотностью мощности Ga (X), то спектральную плотность мощ­ ности выходного сигнала Gy (К) можно определить по формуле

Gy(X) = \<S(jXTI2)\*Gs (X) =

К (А772)

2G8(%).

(11.26)

1 + к QXT/2)

Спектральная плотность динамической ошибки от управляю­ щего воздействия имеет вид

G,j

шш

|Фх

Ц'ХТ'

G s (X) =

1 + К (jXT/2) * Gs

(11.27)

 

 

/2) |2

 

 

Если на вход импульсной системы поступает импульсная помеха со спектральной плотностью мощности Gn (X), то спектраль­

176

ная плотность мощности ошибки на выходе от помехи равна произведению

Gl/n(X) = \G>(jkT12)\*GM =

К ( /М 7 2 )

<?„(*)• (11.28)

1+ К (АГ/2)

В том случае, когда на вход замкнутой импульсной системы одновременно поступают случайные импульсные некоррелирован­ ные управляющий сигнал и помеха со спектральными плотностями Gs (А,) и Ga (Л), то суммарная средняя квадратическая ошибка представляет собой сумму интегралов

+оо

 

 

 

1

Gs (^) т-

 

 

1

+ К

 

 

 

( /М 7 2 )

 

К ФТ/2)

|2„

1

dX

(11.29)

+ 14- К ( j l T l 2)

I

Un *

'J

114- iXT/2 |2 *

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

б|щ(&) = |/Сп(/Ш 2 )/[1 + /С ( /Ш 2 ) ] |* О п(А,),

(Н .30)

где Кп (jXTl2) — частотная передаточная функция от

места ее

приложения до выхода.

 

В рассмотренных случаях прохождения случайного сигнала через импульсные системы предполагалось, что непрерывные слу­ чайные управляющий сигнал s (I) и помеха п (t) прерываются ключом одновременно и в непрерывную часть они поступают в виде одинаковых по форме импульсов со случайно изменяющейся амплитудой. При этом, когда ключ разомкнут, в систему не проходит ни сигнал, ни помеха.

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

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

177

При этом для замкнутых по помехе систем спектральная плот­ ность мощности ошибки от каждой составляющей непрерывной помехи определится выражением

 

(СО) = I Кп (/»)/(! + К (/со)]I2 G" (<о),

(11.31)

где

(©) — спектральная плотность мощности непрерывной по­

мехи

в точке ее приложения; /Сп (/со) — частотная

передаточная

функция непрерывной замкнутой системы от точки приложения помехи до выхода системы; К (/со) — частотная передаточная функция непрерывной системы в разомкнутом состоянии.

Для разомкнутых по помехе систем спектральная плотность мощности ошибки на выходе от непрерывной помехи определится

формулой

 

G"a(<») = \K„ (/со)|2Gn(со),

(11.32)

где Ка (/<я) и СЙ С®) имеют тот же смысл, что и в формуле (11.31). Таким образом, если по управляющему сигналу замкнутая система является импульсной, а по помехе — замкнутой непрерыв­ ной, причем управляющий сигнал имеет случайный характер со спектральной плотностью G* (X), то суммарная средняя квадрати­

ческая ошибка системы представляет собой сумму интегралов

 

4-00

р

 

 

1

Kni (/М72)

X

 

1 + К (А772)

i=\ i + к a m 2)

 

 

X

dX

Knha®)

dco, (11.33)

11 + Д 7721

1+ К (/CD)

где p — число импульсных помех, действующих на систему (обычно это флуктуационные помехи, обусловленные статистическим ха­ рактером числа фотонов в импульсах принимаемого излучения и преобразования их фотоприемниками в импульсы электрического сигнала);.? — число непрерывных помех, по отношению к которым система является непрерывной.

При указанных условиях средняя квадратическая ошибка разомкнутой по сигналу и по помехе импульсной системы будет определяться суммой интегралов

+00 Г

 

<*2> =

-ЕГ 1

* 1 11-

к

(А772) I*G. ( Ц +

 

 

 

- 0 0 L

 

 

 

 

 

 

Р

 

1

 

+°° Q

 

+

2 I ( А

7 7 2 ) I2 б

ш

W

А

J +2 Т[5 Г | (]ш ) |2 ° « » (C0) da>’

 

1 = 1

 

- I

 

- 0 0 k — \

 

 

 

 

 

 

 

(11.34)

где все обозначения соответствуют приведенным ранее.

 

по

Следует иметь в виду, что многие ОЭС, являющиеся замкнутыми

управляющему

сигналу

s (/),

могут быть разомкнутыми

по

помехе. Например, в рассмотренной

ранее системе слежения

за

178

уголковым отражателем фоны на выходе фотоприемника создают непрерывные флуктуации. Такие же флуктуации генерируются в электронном тракте. Как те, так и другие проходят только до исполнительного устройства, вызывая случайное движение опти­ ческой оси системы. При этом отрицательная обратная связь по помехам при отсутствии управляющего сигнала отсутствует. Аналогичная ситуация возникает и в телевизионной системе, предназначенной для определения положения точечного источника излучения. В этих случаях при оценке погрешностей работы ОЭС, обусловленных помехами, систему необходимо рассматривать разомкнутой, а обусловленных управляющим воздействием s (/), — замкнутой.

§ 11.3. ОПТИМИЗАЦИЯ ПАРАМЕТРОВ ИМПУЛЬСНОЙ СИСТЕМЫ ПО МИНИМУМУ

СРЕДНЕЙ КВАДРАТИЧЕСКОЙ ОШИБКИ

В общем случае работа любой автоматической системы (в том числе и импульсной) происходит под воздействием как полезного управляющего сигнала, так и помехи. При этом, как было показано выше, возникают динамическая ошибка воспроизведения полез­ ного сигнала хдин и ошибка от помехи хп. В линейных системах результирующая ошибка определяется суммой этих двух состав­ ляющих. По условиям минимизации их величины динамическая и помеховая ошибки вызывают появление противоположных требо­ ваний к динамическим характеристикам и параметрам системы. Для уменьшения динамической ошибки необходимо увеличивать быстродействие и уменьшать инерционность системы. Однако при этом система будет хорошо реагировать не только на полезный управляющий сигнал, но и на помеху и, следовательно, будет увеличиваться ошибка от помехи. Поэтому можно ожидать, что существуют оптимальные значения параметров линейной системы, при которых суммарная ошибка будет минимальной. Чаще всего при проектировании минимизации подвергается средний квадрат суммарной ошибки <дс|>. Если и управляющий сигнал, и помехи являются стационарными случайными процессами, то в зависи­ мости от конкретной системы для расчета (х^> можно воспользо­ ваться формулами (11.29), (11.33) и (11.34). Если же управляющий сигнал является детерминированным, то суммарную среднюю квадратическую ошибку можно представить в виде

( 4 ) = + <**>,

(п.35)

где Хуст— квадрат установившейся динамической ошибки от управляющего сигнала, ьеличина которой может быть определена, например, по формуле (10.12), (10.17) или (10.22); (я?,) — ошибка от помех, которая в общем случае может быть рассчитана с по­ мощью интеграла (11.33) или (11.34), в которых полагается

G3 (К) = 0.

179

В любом случае в результате расчета определяют средний квадрат суммарной ошибки как функцию <х|) = (х2> (аь а2, ...» ah) от некоторых параметров системы (общего коэффициента пере­ дачи КУпостоянных времени Тъ 7\>, ..., Тт , периода размыкания ключа Ту относительной длительности импульса у и др.), которые обусловливают динамические свойства системы. Для минимизации суммарной ошибки составляется система уравнений

d ( * s )

__ д ( * 1 ) ( а 1> а2» • • м a k) _ Q

даг

дал

з(4)

д(4) (ai' а2........ ак)

= а.

 

(11.36)

да2

да2

1

 

д ( * | ) _

d ( * | ) ( f l i> а 2» • • •» ак)

___Q

Q С

k

daq

daq

 

4 ^

Решение этой системы уравнений и позволяет найти те значения

параметров аъ аъ

..., aq (q <* &), при которых суммарная ошибка

будет минимальна.

 

 

 

 

Глава 12

ВЫБОР ОСНОВНЫХ ПАРАМЕТРОВ ЦВМ, ВХОДЯЩИХ В СОСТАВ О ЭС

К числу основных параметров ЦВМ, которые необходимо выбрать разработчику при проектировании ОЭС со встроенными микропроцессорами и микроЭВМ, можно отнести: разрядность процессора, а также входного и выходного преобразователей; быстродействие; объем оперативного и долговременного запоми­ нающих устройств; шаг квантования по времени во входном и выходном преобразователе. Выбор этих параметров является достаточно сложной проблемой, представляющей самостоятельный интерес, поэтому ниже даются некоторые рекомендации, позво­ ляющие выбрать указанные параметры лишь в первом приближе­ нии. Более точно эти параметры можно определить эксперимен­ тально в процессе отладки системы.

§ 12.1. ВЫБОР РАЗРЯДНОСТИ ПРОЦЕССОРА. ВХОДНОГО И ВЫХОДНОГО ПРЕОБРАЗОВАТЕЛЕЙ

При рассмотрении этого вопроса будем полагать, что ЦВМ, входящая в состав ОЭС, работает в двоичном коде.

Как было показано ранее, между шагом квантования по уровню и разрядностью двоичного кода тхвходного преобразователя Н—К справедлива" зависимость

— *im/(2mi — 1),

(12.1)

где х1т— максимальное значение сигнала на входе преобразова­ теля Н—к.

180