
книги из ГПНТБ / Преснухин, Леонид Николаевич. Основы теории и проектирования приборов управления учебное пособие для машиностроительных и энергетических вузов
.pdf§ 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(т) = — соответствует
7н
методу конечных разностей; функция веса Рг(д) = е~” + е-г(Гъ-^
2Ь
изображается цепной линией.
Отправляясь от уравнения (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) |