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

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

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

Представим выражение в круглых скобках в виде суммы про­ стых дробей

\/[р 2(1 Н“ ТоР)\ — Мр2T QIP -f- TQ!(\ -}- Top)-

Для всех простых дробей в табл. 8.1 имеются 2-преобразования. Поэтому, обращаясь к этой таблице, получим

 

Т У / „ \

K Q $ 2

z — 1

Z

 

А

 

 

А(2) =

 

z

 

14* Т^р I

 

_

0б2 2 - 1

 

Г

Тг________ TpZ

 

T0Z

1

 

 

 

61

г

 

L (z — 1 )\

{z— \

) ' z

— d \

 

 

 

_

б2/Со 1(Г -

Гр + 7VQ г +

(1 — й) Гр -

Td]

 

 

 

 

 

 

 

6 Х(z — 1 ) (г — d)

 

 

. где

d = е 'г/г" =

e-0,25 =

0,78.

 

 

 

 

Подставив

числовые значения входящих величин, получим

 

 

 

К (2) =

(Зг +

2,5)/(г2 -

1,782 +

0,78).

 

После

соответствующих

подстановок

z =

(1 +

ш)/(1 — w) и

w =

/7V 2

получим:

 

 

 

 

 

 

 

 

К и

=

12,5 (1

+ 0,09ny) (1 — w)/lw (1 + 8, lay) J;

 

 

 

Л(/Я) =

 

100 (1 + /0,011Ц (1 — ;0,13А.)

 

 

 

 

 

 

 

 

А (1 4“ Д.01Х)

 

 

 

§ 9.3. ПЕРЕДАТОЧНЫЕ ФУНКЦИИ

ЗАМКНУТЫХ ИМПУЛЬСНЫХ И ЦИФРОВЫХ СИСТЕМ

Знание дискретной передаточной функции разомкнутой им­ пульсной системы К (г) и системы с ЦВМ /С2 (z) = D (z) К (г) позволяет найти дискретные передаточные функции замкнутой системы. При этом, как и в непрерывных системах, здесь рассма­ тривают передаточные функции трех типов: по управляющему сигналу, по ошибке от управляющего сигнала и по ошибке от помехи.

В соответствии со схемой, приведенной на рис. 9.10, а, можно записать зависимость между 2-преобразованнями выходной вели­ чины у [п] и ошибкой х [п J

Y (2) = К (*) X (2).

(9.34)

Выходная и входная величины здесь рассматриваются в одни

и те же моменты времени t — пТ (п — 0,

1,2, ...), т. е. при отсут­

ствии смещения. Учитывая, что при этих условиях X (г) =

5 (г) —

Y (г), уравнение (9.34)

можно записать в виде

 

Y (г) -

К (z) [5 (2) -

Y (z) 1.

(9.35)

Отсюда следует, что

 

 

 

У (2) = IК (г)/И

+ К (г))| S (г) = Ф (г) S (г),

(9.36)

6 Лебедько Е. Г-

161

где выражение

 

Ф (г) — К (2)/S (2) = К (г)/[1 4- К (г))

(9.37)

называется дискретной передаточной функцией импульсной замк­ нутой системы по управляющему воздействию. Как следует из формулы, она равна отношению г-преобразований выходного сигнала у [п] к управляющему воздействию s [я].

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

^ S(z) ^X(z)

X(z)

Y(z)

 

Y{z)

 

Ks (z)

 

 

4

- ггНЕь

 

y[n!

Ф(2) =

 

 

 

 

\Ф ]

 

S (г)

 

1+ Ks (г)

 

 

 

 

 

 

 

D (z) К ( Z )

 

(9.38)

 

 

 

 

 

 

\+D(z)K (z)

 

 

 

HE(Z)

 

 

 

 

 

 

 

 

Учитывая, что в системах

 

 

 

 

 

о ЦВМ

всегда

имеется

за­

 

 

 

 

 

паздывание,

обусловленное

 

 

 

 

 

характером работы ЦВМ, сиг­

 

 

 

 

 

нал на выходе системы необ­

Рис. 9.10. Структурные схемы замкнутых

ходимо рассматривать в сме­

щенные

моменты

времени

ОЭС: а — импульсной; б — со встроен­

(п 4- е)

Т, в

то

время

 

как

 

 

 

ной ЦВМ

 

 

в

моменты

пТ (е — 0).

 

сигнал

на

входе

подается

Поэтому точная формула, согласно'ра­

боте [2], устанавливающая связь между 2-преобразованиями на

входе 5 (г, 0) и выходе У (г, е),

имеет вид

 

У(г,

D (z,

г) К (г, в)

? , п

(9.39)

1+D(z, 0) К (г,

0) °

 

 

Однако эта формула используется только в тех случаях, когда необходимо рассматривать поведение выходной величины у (£) между дискретными значениями отсчетов пТ.

Если в формулу (9.34) подставить Y (z) — S (2) — X (2), то получим

 

5(2) = Х (г)

[1 4-

K(Z)1

(9.40)

Отсюда следует,

что

 

 

 

 

X (2) =

5

(2)/П 4- К (2) ]

=

Фх (2) 5 (2),

(9.41)

где

 

X (г)/5 (2)

= 1/П

4* К (г)]

• ,

Фх (2) =

(9.42)

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

Для системы с ЦВМ дискретная передаточная функция по ошибке от управляющего воздействия имеет вид

Ф* (г) = 1 /П +

(г)] = 1 /1 1 + 0 (г) К. (г)).

(9.43)

162

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

Рис. 9.11. Структурная схема замкнутой импульсной си­ стемы при воздействии помех

Для того чтобы можно было воспользоваться 2-преобразова­ нием, необходимо сначала помеху перенести на вход импульсного элемента. В рассматриваемой схеме для этого достаточно перенести помеху к началу системы. Если в точке приложения найти изобра­ жение помехи по Лапласу

ои

 

N (р) = J п (0 е-"' dt,

(9.44)

О

 

то на выходе системы ее изображение определится выражением

JV, (Р) = N (р) Kt (Р)■

Следовательно, если эту помеху отсчитывать в дискретные моменты замыкания ключа t — пТ, то можно найти ее z-преобра- вование

Nx (z)= Z \ N ip) Кг (Р)1

(9.45)

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

У„ (г) = - Z |N (р) Кг Ш П + К (г) I

(9.46)

Следует обратить внимание на то обстоятельство, что такой характер прохождения помехи имеет место в том случае, если в си­ стеме имеется реальное устройство, преобразующее все непрерыв­ ные сигналы в последовательность коротких импульсов (например, как в системах с ЦВМ, в которых ввод и вывод информации в про­ цессоре осуществляется в дискретные моменты времени с интерва­ лами 7\ и Ту. Многие ОЭС (в частности, системы слежения за

6*

a) {sfnTj
1

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

помехи п (t), приложенной так, как это изображено на рис. 9.11, определится по формуле

п и п

Т 2Т ЗТ 4 Т 5Т 6Т t

б) ! УМ

r^tf[nT]

1тI, Т т

ТZT ЗТ 47* 5Т 6Т t

4>1х[пТ]

Ya (p )= * - K 2(p) N {р)П\.+

 

+ К ( р ) J= Ф„ (р) N

(р), (9.47)

где

К (р) =

/Со (Р) =

Ki

(р)

Кг (р)

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

функция

непрерывной

части системы;

Yn (р) — изображение

по

Лапласу

ошибки

на

выходе си­

стемы от помехи

п (t);

 

 

 

 

Ф„ (Р) =

К (P)!N (р)

=

1 -

Т

ЗТ 4 Г ^ ^

 

-------Кг(р)Ш +

К(р)]

(9.48)

 

 

— передаточная функция непрерывной

 

 

 

5Т 6Т t

 

 

 

системы по ошибке от

помехи п (t).

 

 

 

 

 

 

 

 

 

 

В некоторых случаях ОЭС,

будучи

Рис.

9.12.

Графики

упра­

замкнутой по управляющему

воздей­

вляющего

сигнала и

реак­

ствию s (t), является разомкнутой по

ции

на

него в импульсной

помехе. Тогда сигнал уп(t) на

выходе

 

 

системе

 

системы от помехи п (t),

приложенной

 

 

 

 

 

так, как это изображено на рис. 9.11, определится формулой

^ (Р) = К2 (Р) N (р).

Для замкнутых по сигналу и по помехе систем с ЦВМ г-преобра- зование сигнала от помехи, приложенной в некоторой точке непре­ рывной системы, на выходе определится формулой

Кп (г) =

- Z \N (р) Кг {р)\![\ + D (z) К (г) J.

(9.49)

По дискретным

передаточным функциям

замкнутых систем

посредством подстановок г = efaT, z — (1

w)/(1 — до),

до =

= jXT/2 можно затем получить соответствующие частотные ха­ рактеристики Ф (е'407) и Ф (/X). Зная дискретные передаточные функции Ф (г) и входное воздействие 5 (г), можно найти также решетчатые функции у Iп] и х [п ].

Рассмотрим, например, какими будут решетчатые функции у [п] и х [п ] для импульсной системы, структурная схема которой изображена на рис. 9.4. Как было показано ранее, дискретная передаточная функция этой системы в разомкнутом состоянии имеет вид

* «

=

1Г -л7(Г--дГ - ? -

(Ж + 5 .7Г •

<9-50>

Пусть на

вход

импульсной

системы

прикладывается

сигнал

в виде единичной

ступенчатой

функции

s (t) =* 1 (t). Если этот

164

сигнал отсчитывать в дискретные моменты времени с интерва­ лом 7\ то будем иметь единичную ступенчатую функцию s Iп] (рис. 9.12, а), 2-преобразование которой, согласно табл. 8.1, опре­ делится формулой

 

5

(2) =

2/(2 — 1).

 

(9.51)

Подставив значения К (г) и

5 (г) в выражение (9.37), найдем

Y (г) =

К (г) 5 (г)/[1 +

К (г)] = 0,1122/(г2

— l,67z +

0,78) X

X

(z — 1) = 0,1123/(г3 — 2,67г2 + 2,45z — 0,78).

Разделив числитель на знаменатель, получим ряд Лорана

функции Y (z) в виде

 

 

 

 

 

Y (г) = 0, Иг"1 -|- 0,29г-2 + О.Бг'3 +

 

 

+ 0,81г‘ 4 +

1,16г-5 + 1,512-°+ ...

 

Следовательно, у 10) =

0,

у [1 ] = 0,11,

у [2 1 =

0,29, ...»

...,# [6] =

1,51, ...

 

 

 

 

Если рассмотреть более продолжительную последовательность

решетчатой функции у [п],

то можно убедиться, что ее изменение

носит затухающий колебательный характер (рис. 9.12,6).

Подставив значения К (г) из выражения (9.50) и 5 (г) из выра­

жения (9.51) в формулу (9.41), найдем г-преобразование ошибки следящей системы

 

X (г) = 5

(г)/[1

+ /<(г)]

=

 

=

(г2 — 1,782 + 0,78) г/(г2

— 1,67г +

0,78) (г — 1).

Однако для нахождения

решетчатой

функции

х Iп 1 более

просто

воспользоваться

разностью

х

[п] = s

[/г] — у [п]

(рис. 9.12, в).

 

 

 

 

 

Глава 10

УСТОЙЧИВОСТЬ И ТОЧНОСТЬ РАБОТЫ ИМПУЛЬСНЫХ И ЦИФРОВЫХ СИСТЕМ

§ 10.1. УСТОЙЧИВОСТЬ ИМПУЛЬСНЫХ СИСТЕМ

Для обеспечения устойчивости импульсной системы необходимо и достаточно, чтобы корни характеристического уравнения, полу­ чаемого путем приравнивания нулю знаменателя дискретной передаточной функции замкнутой системы (9.37):

\ + К (*рТ) = 0 ,

(10.1)

располагались в левой полуплоскости комплексной плоскости корней р = с -Ь /со (рис. 10.1, а). Для систем с ЦВМ характери­ стическое уравнение определяется из формулы

1 + D (&т) К (tpT) = 0.

(10.2)

165

Корни

этого

уравнения

 

 

 

 

Н, а = 0,835 ± /(0,835)* -

0,78 =

0,835 ±

/0,292.

Модуль

корней

 

 

 

 

 

 

 

 

I г,,„| = V (0.835)2 +

(0.292)2

=

0,88 <

1.

Следовательно,

система

устойчива.

условиям (10.6). Здесь

Тот же

результат дает

проверка по

А = —1,67,

В — 0,78.

Поэтому получаем 1 — 1,67 0,78 =

= 0,И >

0;

1 +

1,67 +

0,78 =

3,45 >

0;

0,78 <

1.

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

зование. Переход к таким

уравнениям осуществляется

подста­

новкой z =

(1 + w)/(l — до) в передаточную функцию К (z) с по­

следующим

составлением

уравнения

 

 

1 +

К (до) = 0.

(10.7)

Для устойчивости линейной импульсной системы необходимо и достаточно, чтобы все корни характеристического уравнения, которое получается из выражения (10.7), находились в левой полу­ плоскости на комплексной плоскости корней (рис. 10.1, в).

Однако, располагая характеристическим уравнением, выра­ женным через до-преобразование, для решения вопроса об устой­ чивости системы корни находить не обязательно, так как для этой цели можно воспользоваться критерием Гурвица, широко приме­ няемым при анализе непрерывных систем. Для исследования воп­ роса устойчивости системы по критерию Гурвица необходимо составить матрицу из коэффициентов характеристического урав­ нения

a 0Wm -(- Ц1дот "1 + • • • - { - а т-1Ш + а т = 0

в виде

Oi I а3j

а0 CI4j а4

0 «1 а3

• •

• •

• •

Для устойчивости системы необходимо и достаточно, чтобы при

Оо >

0 все определители частных матриц, очерченных пункти­

ром,

были больше нуля.

Пусть, например, имеется характеристическое уравнение вто­

рого

порядка 22 + Аг -+• В = 0.

Врезультате замены переменной z — (1 4- до)/(1—до) получаем (1 — А 4- В) до2 -Ь 2 (1 — В) до + (1 + А + В) = 0. (10.8)

167

Таким образом,

здесь

а0 =

1

Л + В, й\ — 2 (1 — В),

At = 1 4- А

В.

имеет

вид

 

 

Матрица

Гурвица

0

 

 

 

 

а1 1

(10.9)

 

 

 

а0

#2

 

 

 

 

Отсюда получаем условие устойчивости в виде неравенств (10.6), Из этих условий можно определить, при каких значениях неко­ торых параметров система может находиться на границе устой­ чивости. Пусть, например, дискретная передаточная функция

разомкнутой

системы

определяется выражением

 

К (z) =

уТК (1 — d) г!(г - 1) (г — d),

где у = 0,1;

Т — 0,05 с; d = 0,78.

Определим коэффициент передачи К непрерывной части си­ стемы, при котором система будет находиться на границе устой­

чивости. Из формулы 1 +

К (z) = 0 находим характеристичес­

кое уравнение

 

 

 

(г — 1) (г — d) + уТК (1 — d) z =

г2 — [1 +

d

— уТК (1

— d ) ] z - M

= 0.

 

Подстановкой z = (1 +

w)l{\ w) переходим к

до-преобра-

8 0 В З Н И Ю

 

2 (1 — d) w +

[2 (1 -f d) — уТК (Г — d)) w2 +

+ уТК (1 — d) = 0.

Из матрицы Гурвица (10.9) следует, что все коэффициенты

характеристического

уравнения

должны быть

больше нуля,

 

т.

е.

2 (1 + 'd ) — уТК (1 — d) >

0; (1 —

 

— d) > 0;

уТК (1

— d)

> 0 .

 

 

 

 

В

результате

получаем

 

 

 

 

_

2(1 + d )

 

3,56

 

=

3240 с-1

 

 

 

yT(\-d) ~ 0,1-0,05-0,22,

 

 

 

 

Переход к амплитудно-фазовой частот­

 

ной

характеристике

К (/А.)

позволяет

 

оценить устойчивость

системы путем ис­

 

пользования критерия Найквиста. В со­

 

ответствии с этим

критерием при

измене­

Рис. 10.2. Амплитудно­

нии частоты А от 0 до оо годограф ампли­

тудно-фазовой частотной

характеристики

фазовая характеристика

на

комплексной

плоскости

не

должен

системы, устойчивой по

критерию Найквиста

охватывать точку

с координатами Jm =0,

 

Re =

—1.

 

 

 

 

 

 

Пусть, например, дискретная передаточная функция имеет вид

К (z) = КТЦг — 1).

Подставив г — (1+ ш)/(1 — w),

найдем

К (w)—KT (1 -\-w)l2w.

В результате

подстановки w =

/А772

получим

К (А) =

К (1 — jXTl2)/jX = —у {КЩ — КТ12.

168

На комплексной плоскости эта функция имеет вид, изображен­ ный на рис. 10.2. Следовательно, критическое значение величины kT, при котором система находится на границе устойчивости, определяется величиной КТ — 2.

§10.2. ОШИБКИ ЗАМКНУТЫХ ИМПУЛЬСНЫХ СИСТЕМ

ВУСТАНОВИВШИХСЯ РЕЖИМАХ РАБОТЫ

Вконечном счете основным требованием к автоматическим системам (в том числе и к импульсным) является заданная точность воспроизведения или оценки управляющего сигнала. Поэтому, найдя ошибки таких систем при некоторых типовых входных сигна­ лах, можно судить о динамических качествах системы. Наиболее полное представление о них можно получить по кривой переход­ ного процесса, получаемого, например, при входном сигнале в виде единичной ступенчатой функции. Такой переходный процесс позволяет найти и максимальное значение ошибки, и ее изменение во времени, позволяющее судить о склонности системы к колеба­ тельности. Однако такой путь оценки качества импульсных авто­ матических систем связан с необходимостью выполнения большого объема вычислений и последующей обработки их результатов. Поэтому для оценки качества автоматических систем пользуются различными критериями. В настоящее время разработано не­ сколько групп таких критериев, приведенных в работах [2, 9, 101. Здесь мы рассмотрим лишь два способа оценки качества, базирую­ щихся на оценке точности работы в установившихся типовых режимах.

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

изводных s, s, ...» s<m> = dms (t)/dtm. Тогда дискретное изображе­ ние ошибки для импульсной замкнутой системы определится выра­ жением

X (ерТ) = фя (ерТ) S (ерГ) = S (ерТ)/[1 + К (ерТ)], (10.10)

где Ф* (еРт) — дискретная передаточная функция замкнутой системы по ошибке от управляющего воздействия; 5 (ерГ) — дискретное преобразование входного сигнала s (t). Разложив эту передаточную функцию в ряд Тейлора по возрастающим степеням р , получим

X (е"г) = [с, + clP + -if- р * + £

р> + • • • ] S(е»>0.

(Ю.11)

Переходя от изображения к оригиналу, найдем ошибку в уста­

новившемся режиме работы в виде ряда

 

* [я] = c0s [п] + c,s [я] +

If-Si [л] + •••

(10.12)

Коэффициенты этого ряда определяются по формуле

 

сп= [dnOx (epT)/(dp'%=0-

(10.13)

169

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

К (г) = 0,llz/(z — 1) (z d),

если на вход системы подано управляющее воздействие с постоян­ ной скоростью изменения входного сигнала s (t) — s0t. Решетчатая

функция этого сигнала s [п] =

s0tiT, где s0 =

ds (t)/dt — постоян­

ная

скорость

изменения

входного

сигнала,

и, следовательно,

S [fi]

= s0.

что г =

ерТ, найдем

 

 

Учитывая,

(ерГ — \)(ерТ 0,78)

Ф* (е*т) = ! / [ ! + К (z)] |г==е,г =

 

(ерт— 1) (ерГ — 0,78) + 0 ,11ерТ

Из предыдущего уравнения находим:

 

 

 

 

со =

IФх (ерТ) |р=о = 0 ;

 

 

йФ*(ерГ)

_

0,1\Терт(е2рт— 0,78)

 

 

dp

 

[(еРГ - 1) (ерГ -

0,78) +

0,11ерТ]2

Ci = № х (epT)/(dp)]p=о = 2Т.

Если полученный результат применить к примеру системы слежения за уголковым отражателем, рассмотренному ранее, где 0,11 = уКТ (1 — d), d — 0,78 = е~7’/7\ то будем иметь сг =

=1/уК = 0,1 с.

Подставив выражения s [л], s [л], с0 и сг в формулу (10.12),

найдем значение установившейся ошибки слежения

*уст [п\ — Cis Iп1 = 2 Т SQ-

Д ля примера системы слежения за уголковым отражателем имеем

■^уст [^] — sJyK — 0,150-

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

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

Фх (ерг) = -Км [1 — К' (ерГ)]»

(Ю-14)

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

Коэффициенты ошибок с0, clf ... в этом случае получаются дифференцированием выражения (10.14), согласно формуле (10.13).

170