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

книги из ГПНТБ / Преснухин, Леонид Николаевич. Основы теории и проектирования приборов управления учебное пособие для машиностроительных и энергетических вузов

.pdf
Скачиваний:
12
Добавлен:
29.10.2023
Размер:
10.64 Mб
Скачать

§ 23. Выбор оптимального метода сглаживания в ЦСУ

139

Уравнение (226) в этом случае имеет вид

z2( —г2-)-а2-)-22)Р(т) = 2Д

при

или, что то же самое,

piv(T)_(a2_|-22)P"(T)=-2A (228)

Общее решение уравнения (228) имеет вид

Р (т) = - Дт2 + Дт + С +

Условие симметричности (225) дает

Ах = АТп- B^BeW'^2.

Остается вычислить три коэффициента: А, В и С функции веса:

Р(т) = Д-(7’н-т) + 5(^+ег(Гн-т)) + С

(229)

прит<7'н, где

r=]/a2-]-22.

 

Искомые постоянные будут найдены, если подставить

(227) и

(229) в исходное интегральное уравнение (224) и после

интегри­

рования приравнять нулю коэффициенты при функциях

времени.

В результате подстановки имеем

 

тн

 

B?q J

(Т'н - -) + В (е" 4- /(Гк-т)) + С] X

 

о

Используя свойства дельта-функции, получим

Гн

2аР (т) - г2 J [А-. (7'н - т) -\-В (еп + ег <гн~ ) +

О

+ С] е-‘V-I cos (2\t - ТI - arc tg —= -А- .

(230)

\

°2 — 22 /

 

Гн

t

 

Разбиваем интеграл J

на два: f + J, причем в первом ела.

о

о t

 

гаемом \t—1| заменяем на (t—т), а во втором \t — т| заменяем на — \t— т|. После этого интегрирование производится достаточно просто.

140

Глава III. Определение параметров движения цели

Произведя все необходимые упрощения, приводим равенство

(230)

к виду

 

4 — -J-£-aZ(Afcos Q^-H-^sin +

 

r2

 

+ ea(/“r«) [McosS(rH+0+WsinS(rH-0]=-^- , (231)

где

Л4=-ДГЙ-^-+ Ca 4-5 [a (/r« +1)(A-1)];

 

r2

 

L

J

 

AT=^+C2+^-[(a*-a2)(e'rH + 1) 4- ar (e^K - ]) ].

(232)

Г2

zS2

 

 

 

 

Уравнение (231)

обращается в тождество, если члены,

содержа­

щие t, обратятся в

нуль.

А это

возможно лишь при соблюдении

условий:

 

 

 

 

 

Тогда уравнение (231)

дает возможность найти дисперсию оши­

бок в сглаженных

параметрах движения, равную К.

 

 

 

лр"2 д~

 

 

(234)

 

к =—=£2. =£?.

 

 

г2

Ъц

аи

'

Условие (221) для функции

веса

(229) имеет вид

 

—н + С7'н4--(егГн-1) = 1.

(235)

6

 

г

 

 

 

Из уравнений (233) и (235) находим значения неизвестных ко­ эффициентов:

 

,

, -гГн1

 

 

 

 

— (г

—а)е

н]

 

 

 

(236)

Г /

4

\

/ 4

\

—тТ

 

+ Q2 (Г„Г2 4- 4a) П

------ТИ \-^-—+ тпу

«

 

 

2Q2 (7Hr2 + 4a) Ае

н .

 

 

 

 

 

Г1 [(г + a) — — а) е~гТн ] ’

 

 

 

 

 

С=^{[а7-нг2 + 2(а2-22)]

22(Гнг2 +

4а)(е

-гТ

\

 

H+U I

[(г + «)- (г-«)е-гГн] /'

£ 23. Выбор оптимального метода сглаживания в ДСУ

141

В случаях, когда гТИ^1, выражения для постоянных, .вхо­ дящих в функцию веса, можно несколько упростить, если в

суммах пренебречь членами, содержащими е~гГ* в качестве мно­ жителя. В результате получим

в(Гнг2 4-4д) Ле г7~н

 

Г4 (г+ а)

 

 

 

 

 

С = А [[а7нг2 + 2(а2 - 22) ] - 22(ГнГ2+-^Ц.

(237)

Г4

(

 

 

Г + а

J

 

Перепишем формулу (229) в виде

 

 

 

Р(,)=Л, ^=± + ct -L + S,

 

 

где

 

 

* И

* н

 

 

 

 

 

 

 

 

 

Л^з

С,==С7’Н;

В,=2ВЬе

т

1 — е"

 

Д, = —2-;

Тя\

Ь = !—^

 

1 .6

1

н-

1

 

г

Функция веса Р^)

уЗ

соответствует параболическо­

1 н

 

му методу сглаживания; функция веса /52(т) = — соответствует

методу конечных разностей; функция веса Рг(д) = е~” + е-г(Гъ-^

изображается цепной линией.

Отправляясь от уравнения (220), можно показать, что Pi(t) является оптимальной (в смысле минимума дисперсии) функцией веса сглаживания скорости, если ошибки в координате представ­ ляют собой чисто случайный процесс ’; Р2(т) —оптимальная функ­ ция веса метода сглаживания скорости, если ошибки в скорости — чисто случайный процесс, и, наконец, цепная линия Рз(т) — опти­ мальная функция веса для сглаживания параметров движения цели, искаженных ошибками с корреляционной функцией вида

R (т) = Е'2е-“ Iх cos 2т:.

Однако обычно такой вид имеют функции корреляции ошибок в координатах, а не в параметрах движения цели. При этом функ­

ции корреляции ошибок в скоростях

изменения координат

имеют

вид

(227).

 

 

1

Напомним, что функция корреляции

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

/?(т).—

= Х&(т).

142 Глава III. Определение параметров движения цели

Выкладки, приведенные в настоящем параграфе, позволяют сформулировать следующий вывод.

Если функция корреляции ошибок в координатах имеет вид

=cos St

исоответственно функции корреляции ошибок в параметрах движе­ ния цели изображаются формулой (227), то наивыгоднейшая функ­ ция веса сглаживания представляет собой линейную комбинацию трех рассмотренных выше функций Pi (т), Рг(т) и Рз(т). Это зна­ чит, что среднеквадратичная ошибка в сглаженной скорости будет минимальной, если ДСУ работает как схема из трех параллельно соединенных узлов, сглаживающих мгновенную скорость по мето­ дам, определяемым функциями веса Р1(т), /А (г) и Рз(т) соответ­ ственно. Искомая сглаженная скорость представляет собой линей­ ную комбинацию скоростей, определенных в каждом из этих узлов:

+ciДх2 4-ЛЧгЗ-

Коэффициенты Aj, В{ и Ci определяются постоянными а и Q функции корреляции и выбранным наблюдательным временем ТИ.

§ 24. Выбор наблюдательного времени ДСУ

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

а следовательно, и всего прибора в нормальный режим работы.

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

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

зависит от калибра артиллерийской установки, для которой проек­

тируется ПУАО. Если калибр большой, то обычно большими будут и времена полета снаряда до цели. А в этом случае малые ошибки в параметрах движения цели могут привести к большим погреш­ ностям в выходных величинах. При малом калибре времена полета снаряда будут также малыми, следовательно, в этом случае можно

допустить большие ошибки в определении параметров движения цели.

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

сошибками в измерении текущих координат цели ’. При этих усло-

1Строго говоря, допустимые ошибки прибора определяются ошибками всех других элементов артиллерийского комплекса. Однако обсуждение этого вопроса

выходит за рамки настоящего курса.

§ 24. Выбор наблюдательного времени ДСУ

143

виях ошибки ДСУ не будут являться доминирующими для ПУАО в целом. Они будут вызывать погрешности того же порядка, что и ошибки в измерении координат цели, неизбежные при слежении за быстро движущейся целью. Если это положение принять за основу при выборе допустимых ошибок, то можно сказать, что отношение ошибки в координате, поступающей на вход ДСУ, к ошибке в ско­

рости, получающейся на выходе, должно быть примерно равно' времени полета снаряда до цели.

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

отношение принять равным некоторому среднему времени

полета

снаряда, характерному для заданной артиллерийской

установ­

ки, т. е.

 

 

Е ~

 

(238)

8 <7

Тср-

Е,-

 

 

iq

 

 

В табл. 2 приведены предельные времена полета снаряда для зенитных орудий различных калибров. Эти времена определяют до­ сягаемость орудий соответствующих калибров. Поэтому top следует

брать меньше

величин,

приведенных в

табл.

2.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таблица 2

Калибр орудия

20

30

40

50

60

70

80

90

100

ПО

120

130

140

150

мм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ттах В сек.

3

5

8,5

12

15

19

24

30

36

42

48

55

61

68

 

таким

 

 

 

 

 

 

 

 

 

Е ~

 

 

Задавшись,

образом,

величиной

отношения —— и вос-

 

 

 

 

 

 

 

 

 

 

 

 

Е*ч

 

 

пользовавшись формулой (234), можно вычислить необходимое для ДСУ наблюдательное время. Из формулы (234) следует,

что

/£-\2

 

 

/ Ьд \

аЗ-ра2

(239)

 

 

 

Подставив значение А из формулы (236) и приняв отношение

равным его среднему значению

0,65 (см. гл. II),

а также

учитывая, что T„Q обычно

много

больше единицы,

получим

довольно простую формулу, из которой можно определить иско­ мое Тя.

(Е~ \ 2 о?’3

 

6 q | ,__ '"‘i н

(240)

 

144

Глава III. Определение параметров движения цели

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

(241)

На фиг. 91 приведены графики зависимости (241), рассчитанные

для

двух значений

Q.

 

 

 

 

 

Рассмотрим пример. Пусть необходимо спроектировать ДСУ

к ПУАЗО, обслуживающему артсистему, для которой характерно

 

 

 

тср =-^-4= 10 сек.,

а

координаты

 

 

 

цели поступают с автоматическо­

 

 

 

го радиолокатора, т. е.а=1,81/сек

 

 

 

и £2 = 2,6

1/сек.

 

 

по графику

 

 

 

Для этих условий

 

 

 

фиг. 91 получаем Та=7,5 сек., а по

 

 

 

формулам

(236)

подсчитываем зна­

 

 

 

чения коэффициентов функции веса:

 

 

 

А = 0,0137;

Д = 1,1-10-12;

 

 

 

 

С=0,00233.

Фиг. 91. График зависимости на­

Таким

образом,

 

оптимальная

блюдательного времени

от сред­

 

ней

продолжительности

полета

функция веса в данном случае имеет

 

снаряда.

 

вид

 

 

 

 

 

Р(ё) = 0,0137 т (7,5 —т)

1,1. Ю~12 (е3-17х-фе3-17 <7>5-") +

 

+ 0,00233

 

 

для

т 4 7,5;

 

Д(т) = 0

 

 

для

т>7,5.

Если принять, что координаты цели поступают

в

ПУАЗО от

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

вии с указаниями главы II полуавтоматические

приводы

и инди­

каторные устройства, то в этом

случае а=1,2 и £2 = 1,96. Для этих

данных согласно графику фиг.

91 рациональным наблюдательным

временем ДСУ будет 7^=8,3 сек., а по формулам (236)

получим

коэффициенты функции веса:

 

 

 

А = 0,0104; 5 = 2,24-Ю’10; С=-0,00186.

 

Следовательно, оптимальная функция веса для вторых условий

будет иметь вид:

 

 

 

Р (т) = 0,0104т (8,3 - т) + 2,24-10~10 (е2-3* -ф

 

— 0,00186

для т<ф8,3;

 

Р (т) = 0

для

8,3.

 

 

<£ 25. Аппроксимация оптимальной функции веса

145

На фиг. 92 сплошной кривой

 

изображен график функции веса,

 

рассчитанной для первых

усло­

 

вий,

а пунктиром —■ для

вторых

 

условий.

 

 

Осуществить сглаживание с функ­

 

цией

веса, точно совпадающей с ее

 

оптимальным видом, при помощи

счетно-решающих устройств непре­

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

ференцирующе-сглаживающих ус­ тройств.

Фиг. 92. Оптимальная функция веса для условий: А—рацио­ нального полуавтоматического

измерения координат; Б—ав­ томатического радиолокацион­ ного измерения координат.

§ 25. Аппроксимация оптимальной функции веса

ДСУ могут быть составлены из электрических контуров RC, автоматических фрикционов или электромеханических звеньев. Работа таких ДСУ может быть описана линейным дифферен­ циальным уравнением с постоянными коэффициентами.

Функция веса устройства определяется в этом случае соответ­

ствующим однородным дифференциальным уравнением. Решение

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

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

экспонент с действительными или комплексными показателями. По­ этому при аппроксимации оптимальной функции веса за функцию

веса реализуемого устройства следует принять

Р(т) = А1е-<^ + Д2<?-“»’+ . . .

 

(242)

Число п обычно ограничивается количеством звеньев,

из

кото­

рых предполагается составить ДСУ. Для того чтобы ДСУ

было

в работе устойчивым, а» должны быть положительными.

Если по­

казатели а* — комплексные числа, то должны быть положительны­

ми их действительные части. Экспоненты с действительными пока­

зателями соответствуют инерционным, а экспоненты с комплексны­ ми показателями — колебательным звеньям.

Задача аппроксимации оптимальной функции веса линейной комбинацией экспонент, так же как и задача усреднения парамет­ ров движения цели, не решается однозначно. Решение как той, так и другой задачи определяется тем, что подразумевают под «луч­

шим» приближением. Выбор «лучшего» метода сглаживания пара­ метров движения цели — один из решающих этапов при проекти­ ровании ДСУ. Различные же рациональные определения лучшего приближения при аппроксимации функции веса суммой экспонент,

10 604

146 Глава III. Определение параметров движения цели

как правило, не приводят к значительному разбросу характеристик ДСУ. Поэтому в каждом отдельном случае аппроксимация функции

веса производится по такому методу, который вызывает наимень­

шие расчетные трудности.

При этом всегда должны соблюдаться

некоторые определенные

условия.

Между коэффициентами

и показателями а* аппроксимирую­

щей функции веса (242) должна соблюдаться связь, определяемая коэффициентом усиления и наблюдательным временем ДСУ:

/(оо) = JP (т) ch = k;

О

 

Н!-/(О =

 

 

— k-q.

 

Применительно к функции веса

(242)

указанные условия приоб­

ретают вид

 

 

 

. 4-^ = ^;

(243)

 

Ctj

ot2

 

 

 

 

 

 

dl_ е-

_|_d2_ е~

4-

. . .

4. dm е-аптн = k-q.

(244)

СЦ

а2

 

 

 

“п

 

Как уже указывалось, обычно принимают ц = 0,05, а Л=1. Еще несколько условий, связывающих Д& и сц , можно получить исходя из того, что рациональное сглаживающее устройство должно сре­ зать колебания высоких частот, причем тем сильнее, чем больше частота. Иными словами, при со->оо амплитудно-частотная харак­

теристика рационального ДСУ должна стремиться к нулю не мед­

леннее, чем 1/со5, где $>1. Для получения связей между Ль и сц разложим функцию веса в ряд Маклорена.

P(t)=p(0)4-tP(0)+^p(0)4- . . .

Преобразование Лапласа от обеих частей равенства имеет вид

К (г)—Р(0)

1 Р(0)

г2

г3

Как известно,

К(/в>) = Д(<й)

Поэтому

. , .

... Р (0) ,

Р (0)

Р (0)

,

А (ш)е-^<ш) = -^2.-|----Y------- . . .

'

Уа)

С1)2

 

 

Отсюда видно, что только

при

Р(0)==0

Д(<в)—>0 при <в—»оо,

как — , а при

Р(0) = Р(0) = .

. .

= Р(~* 1>(0) — 0 А (ю) -» 0 при

Ы2

 

 

 

 

$ 25. Аппроксимация оптимальной функции веса

147

ш —> оо, как —г. Применительно к функции веса вида (242)

при-

0)

одно

веденные соображения позволяют получить соответственно

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

ми аппроксимирующей функции.

Л + Д2 + Лз+ . . . +Л = 0;

(245)

Aai+АдгЧ- ■ ■ • +^лал—0;

(246)

ДаГ2+Д2аГ2+ • • ■ +ДпаГ2=0..

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

Нередко удается получить удовлетворительное приближение реально осуществимых ДСУ к ДСУ с оптимальными характеристи­ ками, используя особенности амплитудно-фазовых характеристик,

отвечающих рациональным функциям веса. Такими дополнитель­

ными условиями могут

быть, например, совпадения

предельных

частот и динамических отставаний заданного метода

определения

параметров

движения

цели

и его

подходящего

приближе­

ния.

 

 

 

 

 

 

 

 

Функция веса

(242)

соответствует

случаю

различных

корней

характеристического уравнения ДСУ, реализующего

лучшее при­

ближение заданного

метода

выработки

сглаженных

скоро­

стей.

 

 

 

 

 

 

 

 

В случае кратных корней аппроксимирующая функция веса за­

писывается в

виде

 

 

 

 

 

 

Р(т)=(А + Л2г+ .

. . +Ллг«->)^.

 

(247)

Функция

веса

(242) имеет 2п параметров, которыми

можно

варьировать, а в функции (247)

можно распоряжаться выбором

только (п+1) величин.

Тем не менее опыт и

многочисленные ра­

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

интерес.

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

нения, определяющие искомые значения и а.

10*

148

Глава III. Определение параметров движения цели

Показатель

а удовлетворяет уравнению

1

1

2! .

.(га-1)!

У (а)

1

2!

3! .

.

га!

(~2а) У'(а)

2!

3!

4! .

(«+1)!

(-2а)2 У" (а) = 0, (248)

 

п! (а-Н)! (п+2)!

(2га- 1)!

(—2а)лК<") (а)

где

У (a), Y' (а), У" (а) —соответственно

 

значения заданной (оп­

тимальной) функции передачи и ее производных при z=а.

Коэффициенты Ай определяются из системы га-линейных уравне­

ний

по вычисленному из

(248)

значению а.

 

 

kl

л .

(* + D!

 

л

_

 

*(-2а) + 1

 

(-2а)&+2

2

’ ‘ •

• +(-Пл(-7^>4 = ^)М.

(249)

где А=0, 1,2,.

.

.

(.

^а)

 

га —1.

быть решено графически. Система

Уравнение

(248)

может

уравнений (249) просто решается аналитически.

Выше было показано, что амплитудно-частотная характеристи­

ка при увеличении со стремится к нулю, как 1/сог, если Р(0) =

=Р(0) = ... =Р(,-_9(0) =0. Для того чтобы определить коэффициен­ ты Ай и показатель а, обращающие в минимум среднеквадратич­

ное отклонение приближения от заданной функции веса и удовлетворяющие, кроме того, перечисленным дополнительным условиям, достаточно воспользоваться уравнениями (248) и (249).

В определителе (248) следует опустить (г—1) первых строк

и г первых столбцов, превратив его,

таким образом, в опреде­

литель (га—г+1) порядка. В системе

(249) следует

положить

А1 — А2 = А3^= . . . = А,_1 = 0 и найти

Аг; Аг+1, . . .

Ап из пер­

вых (га — г+1) уравнений.

 

 

Аппроксимируя, например, параболическую функцию веса выра­

жением вида (247) при дополнительном

требовании P(0)=Ai = 0,

получаем в зависимости от га следующие приближения:

а) при га = 2

 

 

ГЛ(~)=22,5^е

(250)

б) при га = 3

 

<251>

в) при га = 4

УЛ Ш=(740 й ~100 7? +15)

(252)

Соседние файлы в папке книги из ГПНТБ