
книги из ГПНТБ / Преснухин, Леонид Николаевич. Основы теории и проектирования приборов управления учебное пособие для машиностроительных и энергетических вузов
.pdf§ 21. Устройства для определения параметров движения цели |
129 |
Частотные характеристики, отвечающие параболическому сгла живанию, определяются из соотношения
А (а))е-^(“) = С 6т <zh —Ч
' 7 |
1 |
уЗ |
|
V |
1 Н |
О
Несложные преобразования приводят к формуле
А (о) |
]/4ф-sin |
\ 2 |
со6/ |
г |
Фи1г. 86. Ступенчатая функция fiifco). |
Фиг. 87. Амплитудно-частот |
|||||
|
|
|
|
ная характеристика |
пара |
|
|
|
|
|
болического |
сглаживания. |
|
В силу того, |
что |
Л(ш) — существенно положительная функ- |
||||
А (ш) = |
“ 1 Н |
К4+“2/н Isin |
-агс te |
|; |
|
|
|
|
I |
' |
2 / \ |
|
|
где E.t (ш) — ступенчатая функция. |
|
|
|
|||
£1(w) = ^ |
при |
Ь<—— arctg^<(A-H)1c; |
|
|||
|
|
|
£=0, 1,2,. . . |
|
|
|
График £,1(ш) |
изображен на фиг. |
86. |
|
|
||
Графики амплитудно-частотной и фазово-частотной характери |
||||||
стик представлены на фиг. 87 и 88. |
|
|
не ско |
|||
Если синусоидальное |
колебание 3x=sin at искажает |
рость, а координату, то амплитуда ошибки в сглаженной скорости В(и) равна
В(и))==0)Д(<в) = ^1-|Л4+^1 sin&-arctg^)|.
Ш2/ „ г | \ 2 2 / I
9 604
130 |
Глава III. Определение параметров движения цели |
|
|
График В(®) |
изображен на фиг. 89. По характеристике -В (со) |
легко определить |
предельную частоту соо. Как уже отмечалось, |
в том случае, когда амплитудная характеристика из-за наличия тригонометрических функций состоит из ряда полуволн, в качест ве соо целесообразно принимать предельную частоту «огибающей», или первый максимум. Для параболического сглаживания наимень-
|
dQ |
= 0 является также наименьшим |
||
ший положительный корень--- |
- |
|||
1 |
d<i> |
|
|
|
корнем уравнения |
|
|
__ |
4<д0Гн |
to- |
|
|
||
Ь |
о |
|
Л |
,2т2 ’ |
|
2 |
|
4 |
— ш07 н |
Фиг. 88. Фазово-частотная ха |
Фиг. 89. Характеристика |
рактеристика параболического |
В (со) для параболического |
сглаживания. |
сглаживания. |
Показатель затухания амплитудно-частотной характеристики
параболического сглаживания s=2. Действительно, при о)->со и $=2
1/4+ <027-2
Д (<о)=—-------z------------
ограничено и отлично от нуля.
Время реакции на маневр при параболическом сглаживании,
как и для метода конечных разностей, равно наблюдательному
времени. Действительно, из выражения для фазово-частотной ха рактеристики параболического метода сглаживания получаем
Д, = 27^=2'/(0) = Гн.
Функция передачи параболического сглаживания определяется
непосредственно как лапласово изображение параболической функции веса.
Y(z) = -J- {(THz - 2) + |
(1\.z + 2)}. |
$ 22. Характеристики ошибок в параметрах движения цели |
131 |
|
На фиг. 90 изображены амплитудно-частотные |
характеристики |
|
, отвечающие параболическому сглаживанию |
(ПС) |
и методу |
конечных разностей (КР), в функции от отношения наблюдатель
ного |
времени |
(7\) |
к |
периоду |
|
|
|
|||
(■&) |
ошибок, искажающих |
ко |
|
|
|
|||||
ординату цели. Сравнивая эти |
|
|
|
|||||||
характеристики, можно сделать |
|
|
|
|||||||
следующее |
заключение. |
Если |
|
|
|
|||||
основная часть спектра ошибок |
|
|
|
|||||||
слежения расположена в райо |
|
|
|
|||||||
не малых частот (когда перио |
|
|
|
|||||||
ды ошибок не малы по сравне |
|
|
|
|||||||
нию с наблюдательным |
време |
|
|
|
||||||
нем ДСУ), лучшим-для опреде-. |
|
|
|
|||||||
ления параметров движения це |
|
|
|
|||||||
ли следует признать метод ко |
|
|
|
|||||||
нечных |
разностей. |
В |
другом |
|
|
|
||||
случае, |
когда периоды |
ошибок |
Фиг. 90. Сравнение амплитудно-частот |
|||||||
слежения |
гораздо |
меньше |
на |
ных |
характеристик |
параболического |
||||
блюдательного |
времени, |
луч |
сглаживания (ПС) и метода конечных |
|||||||
шим методом является |
парабо |
|
разностей |
(КР). |
||||||
лическое сглаживание. |
в |
приборах |
непрерывного |
действия нет |
||||||
В настоящее время |
устройств, осуществляющих принципиально точное параболическое сглаживание. Схема такого устройства, как и ДСУ на конечных
разностях, требует наличия запаздывающего звена с большим вре
менем задержки. Однако имеется много устройств, реализующих
с той или иной степенью приближения параболическое сглажива ние. Некоторые из них рассмотрены ниже.
§ 22. Статистические характеристики ошибок в постоянных параметрах движения цели
Основной задачей проектирования ДСУ является выбор такой схемы ДСУ и таких значений ее параметров, при которых посту пающие на вход узла ошибки вызывали бы наименьшие при задан ном наблюдательном времени ошибки на выходе.
В предыдущих параграфах рассмотрены основные методы опре деления параметров движения цели и их характеристики. В гла ве II приведены статистические.характеристики ошибок слежения.
Весь накопленный материал позволяет теперь получить корреля ционные зависимости между ошибками в параметрах движения, необходимые для выбора рациональных схем и расчета конструк тивных постоянных ДСУ.
Выпишем вначале общие формулы, определяющие статистиче ские ошибки в постоянных параметрах движения цели, и восполь-
*9
132 Глава III. Определение параметров движения цели
зуемся ими для сравнения и оценки рассмотренных в § 19 парамет
ров движения цели.
Начнем с выражения функции корреляции и спектральной плот ности ошибок в прямоугольных составляющих вектора скорости через характеристики ошибок слежения, получающихся при изме рении сферических координат цели. Имеем
х== Дсоэ & cos р,
у= Д cos е sin р,
Н—Д sins.
Отсюда следуют зависимости между ошибками в прямоуголь ных и сферических координатах цели:
8х = 8Дсоз ecos Р —Д8езй1 ecos р— Д8рсоэ s sin Р;
о_у = 8Дсоз е sin р—Д8е sin ssinp + Дор cos в cos В;
оН= 8Д sin е Д8е cos в.
Учитывая формулы, связывающие прямоугольные и сферические координаты, получим
8х = X — X tg е ■ 8е —у8р;
*y = -jy-y tge.8s4-x8p;
8Д = 8-^-/У—Л/ctg е-8г.
На большей части курса цели ошибки изменяются с гораздо большей скоростью, чем координаты. Поэтому при определении ста тистических характеристик ошибок можно пренебречь скоростью изменения координат, полагая их неизменными на отдельных участ ках курса.
Учитывая, что ошибки слежения в сферических координатах — практически не связанные случайные функции, можно записать сле дующие выражения для функции корреляции ошибок в прямо
угольных координатах:
&х W ~ Ri Д (х) + tg2 sA’se Д} + У2/?s? (т);
А £
д (х) +-У2 ^2 &Rb (х)
А*
Rih (т) = RtД(9 + Я2 cig2 е/?6б(т)
§ 22. Характеристики ошибок в параметрах движения цели |
133 |
и соответственно для |
спектральных плотностей |
О8л(«>) = ^ 084(<o) + x2tg2e Gse(w)+j/2GE (и); |
|
ДЛ |
|
Gt у (“) = |
(ш) -j-y2 tg2 е GJe(a>) + x2GEp (a>); |
*Д |
|
GS//(w) = :^-Дгд (<u) + 7/2ctg2s GSe(w).
Д1
Из полученных формул видно, что при
/?8 (") = ДД") = /?8Д (т)
функции корреляции ошибок на входе ДСУ имели бы тот же вид,
что и функции корреляции в отдельных сферических координатах,
и отличались бы от них только меняющейся вдоль курса цели дис персией.
В § 12 и 14 показано, что функция корреляции ошибок в сфе рических координатах цели может быть представлена в виде
/?(т) = £2е_“ Iх'cos Gt.
Из выводов указанных параграфов следует также, что при про ектировании систем слежения по любой координате целесообразно отдавать предпочтение полуавтоматическим приводам с постоянной времени 7=3 сек. и /^0,04. В этом случае всегда а= 1,2= 1,3 1/сек и, следовательно, Й = 1,96=2,12 (1/сек). При автоматическом сле жении за целью с помощью радиолокатора ошибки в угловьдх коор динатах имеют функции корреляции с 1,8 и Q=2,6 (см. табл. 1).
Таким образом, при рациональном проектировании систем слежения (и при пренебрежимо малых инструментальных ошиб ках приборов управления) пределы изменения характеристик
функций корреляции для различных |
координат довольно узкие: |
||
a =1,2-4-1,8 (1/сек) и 2 = 1,96-4-2,6 |
(1/сек). |
|
|
Отсюда ясно, что функции корреляции ошибок, поступающих на |
|||
вход ДСУ, обычно имеют вид |
|
|
|
Rtq |
e-“lTl cos 2т. |
(209) |
|
Здесь а и й имеют те же значения, что и для ошибок в сфери |
|||
ческих координатах. Величину |
следует вычислить для отдель |
ных участков курса в соответствии со значениями координат цели в середине каждого участка:
Функции корреляции |
(209) соответствует спектральная |
плот |
ность |
|
|
Ог (ш) = |
----------а2 + 82 + 0)2----------. |
(210) |
4 V 7 |
я («2 Q2)2 + 2 (<х2 _ Q2) ш2 ш4 |
134 Глава III. Определение параметров движения цели
Если — ошибки в прямоугольных координатах цели, то спек
тральная плотность ошибок в прямоугольных составляющих векто
ра скорости будет |
иметь вид |
|
|
Gi „ (ш) =--------- |
а2 -L <22 _|_ ш2 |
(211) |
|
------------;--------. |
|||
4 ' ' |
л |
(а2 + Q2)2 2 (а2 _ Q2) ш2 ш4 |
|
Аналогичные выражения получаются для статистических харак
теристик ошибок |
в |
других параметрах движения. В частности, |
|
имеем |
|
|
|
8т»г=соз (Q — ^){cossMZ— Ног-Г-d tg (Q — p) op}, |
|
||
^8(2 = sin (Q—Р){созе8Д —/У8е —dtg (Q — р)Вр}. |
(212) |
||
Следовательно, |
при |
принятых допущениях о характеристиках |
|
ошибок слежения в сферических координатах спектральные |
плот |
ности несглаженных ошибок в vr и Q будут также иметь вид (211). Согласно известным положениям теории случайных функций спек
тральная плотность ошибок в |
сглаженных значениях параметров |
|||||
движения цели Gsu (со) может |
быть |
получена |
по формуле |
|||
Gs« (<о) = о>2А2 (ш) G (<о), |
|
(213) |
||||
где <o2G («) — спектральная плотность ошибок в |
мгновенных зна |
|||||
чениях сглаживаемого параметра; |
|
сглаживаю |
||||
А (со)— амплитудно-частотная |
характеристика |
|||||
щего устройства. |
|
известно также, |
|
|||
Из теории случайных функций |
что средне |
|||||
квадратичная ошибка Eiu определяется из соотношения |
||||||
ос |
|
|
|
|
|
|
= jGSu(<o)c/<°. |
|
|
(214) |
|||
о |
|
|
|
|
|
|
Воспользуемся выражениями (213) |
и |
(214) |
для вычисления |
|||
среднеквадратичных значений |
ошибок |
в |
параметрах движения |
цели, усредненных в соответствии с рассмотренными в предыдущих параграфах методами. Опуская промежуточные громоздкие вы кладки, приведем лишь окончательные результаты. Среднеквадра тичная ошибка в скорости, полученной по методу конечных разно
стей, равна
— |
2£. |
_________________ |
(215) |
yi—e-“rHcos 2ГН. |
|||
Гн |
|
|
|
Здесь Еьч — среднеквадратичное значение ошибки в |
координа |
||
те. Из формулы (215) видно, |
что при аТ'нЭ 1 среднеквадратичная |
ошибка в скорости, вычисленной по методу конечных разностей, за
§ 23. Выбор оптимального метода сглаживания в ДСУ |
135 |
висит только от амплитуд ошибок в координате и не зависит от их
частотного содержания. |
ошибка |
в параметрах |
движения |
цели, |
|||
Среднеквадратичная |
|||||||
сглаженных с |
экспоненциальной |
|
1 |
~ ~т |
, |
||
функцией веса Р(т) = — е |
|||||||
определяется по формуле |
|
|
|
|
|||
£2 |
=» |
тИ2 (°2 + Р-2)2 + (Q2~а2)1~аГ2(°2в- Д2)-1! |
/216) |
|
|||
S“ |
Т |
|
Л (а2 -|_ 22)2 + 2Z2 (22 „ а2) + 1 |
‘ |
1 |
' |
|
Аналогичные, но несколько более громоздкие |
формулы |
полу |
|||||
чаются для n-кратного |
экспоненциального и для параболического |
сглаживания.
Как видим, при одних и тех же статистических характеристиках ошибок в координатах различные методы определения параметров движения цели приводят к существенно различным ошибкам. От сюда вытекает естественная задача: какой метод определения па раметров движения цели при заданных характеристиках ошибок во
входных координатах определяет наименьшие среднеквадратичные величины ошибок в параметрах движения цели? Решению этого вопроса посвящается следующий параграф.
§ 23. Выбор оптимального метода сглаживания в ДСУ
При решении этой задачи наивыгоднейшим будем считать метод сглаживания, при котором среднеквадратичное значение ошибки в сглаженной скорости (при заданных статистических характери стиках ошибок слежения) будет наименьшим.
Дисперсия ошибок в скорости (или другом параметре движения цели), сглаженной устройством с функцией передачи У (г), может
быть найдена, как указывалось ранее, из |
формулы |
^2o=Jl^(j“>)l2>)^G~(<* . |
(217) |
о |
|
где G~((u)—спектральная плотность скорости изменения ошибок
в величине q, поступающей на вход ДСУ:
Q -(ио) = <i)2Gs9 (ш).
Выбор наивыгоднейшего метода определения параметров дви жения цели сводится, таким образом, к нахождению функции Y(z), обращающей в минимум интеграл (217). Иногда удобнее опериро вать не со спектральной плотностью, а с функцией корреляции оши бок слежения или скорости их изменения. В таких случаях для ре
136 Глава 111. Определение параметров движения цели
шения поставленной задачи формуле (217) можно придать более удобный вид.
Ошибка в параметре, сглаженном ДСУ с функцией веса Р(т), выражается через ошибки несглаженного параметра посредством соотношения
|
Гн ~ |
|
= |
т)Р(т)(/с, |
(218) |
и
при этом Р(т)=0, если т>7’н.
Согласно теореме Пугачева функция корреляции /?8И(^ —£2) будет получена из /?8~(^ —^), если произвести над ней дважды
операцию (218) |
применительно к каждому аргументу, т. |
е. |
Г |
г |
(219) |
|
р?~ [(^-12)~ (Ъ-т2)] P^)P(r2)d^ dx2. |
|
о |
о |
|
Чтобы получить |
дисперсию 8и, следует положить ti = t2. В этом |
|
случае получим |
г |
|
Г |
|
|
|
нн |
(220) |
Л2а=я«(0)=у |
||
о |
о |
|
Условие обращения в минимум Еъа определяет |
функцию веса |
Р(т), отвечающую наиболее выгодному методу выработки парамет ров движения цели.
Искомая функция Р(т), как и всякая |
функция веса, должна, |
кроме того, удовлетворять соотношению |
|
т |
(221) |
р>(т)<Д=1. |
|
о |
|
Если сглаживаемый параметр в соответствии с гипотезой не ме няется во времени, то (221) —единственное условие, ограничиваю щее выбор Р(т). Определение относительного минимума функцио нала вида (220) при условии (221) сводится к отысканию абсолют ного минимума выражения
н |
Тн |
= |
[P^)dx-1 , (222) |
ОС |
1-0 |
где Xi — постоянное число.
Пусть Р(т) — искомая функция веса. Тогда всякая другая функ ция Q (т) =Р(т) +аДР, как бы мало по абсолютной величине ни было постоянное число а, увеличивает значение функционала /.
/(С) = /(Р+адР)>/(Р).
§ 23. Выбор оптимального метода сглаживания в ДСУ |
137 |
||
Это значит, что /(Р+'аДР) |
достигает минимума при а=0. |
||
Для этого необходимо, чтобы |
|
||
р/(Р+аЛР)-1 =0 |
(223) |
||
[ |
да |
]а=О |
|
Запишем подробнее |
условие |
(223). |
|
тН тН
Z(P+ аДР) = J J [Р(г,) + аДР^)] [Р (т2) + аДР (т2)] R-(т, -t2) X
О о
/н
X dxx Йт2-|-^11J [Р(т)4-аДР(т)] rfx—1
д
или
т т
нн
I(Р-фаДР)==а2 J *j ДР(т2) ДР^)/?-^ —т2) dtx dxi-\- О о
г~Тнгн
-|-2а j J Р (Tj) ДР (т2) р ~ (?! — т2) dt! а1т2 + -о о
of-c + /(Р).
Дифференцируя это выражение по а и приравнивая [в соответ
ствии с условием (223)] полученный результат нулю при а = 0, по лучим
тН |
* н |
= 0. |
J ДР (т2) JР (хх) р~ - т2) d-4 |
оLo
Это равенство должно соблюдаться при любой функции ДР(т). Для этого необходимо, чтобы второй множитель подынтегрального выражения был тождественно равен нулю.
Следовательно, оптимальная функция веса устройства, сглажи вающего ошибки постоянного по гипотезе параметра, должна при 0<т<Тн удовлетворять уравнению
|
f Р (l)P~ (t — т) rf-c = k, |
(224) |
|
0 |
|
где |
—Х =Величина К определяется из условия (221). Из фор |
|
мул |
(220) и (224) определяем физический |
смысл постоянной X. |
Имеем
тнгт1
т2 Т1) *^(^2)^2
ОО
138 Глава III. Определение параметров движения цели
Согласно (224) |
j52u= J кР(т2) |
|
|
|
|
|
о |
|
|
В силу условия |
(221) Е?и = 1, т. |
е. к |
в условии (224) равна |
|
дисперсии |
ошибок |
в скорости, сглаженной наивыгоднейшим |
||
образом. |
Заменяя в |
формуле (224) |
t на |
Тя — t и т на Тя — х, за |
мечаем, что функция веса наивыгоднейшего метода сглаживания параметров движения цели всегда симметрична, т. е.
Р^ = Р(ТЯ-^. |
|
(225) |
В случае, когда спектральная плотность |
ошибок на входе СУ |
|
представляется рациональной дробью, а |
этот |
случай именно и |
интересен для нас, для решения интегрального |
уравнения (224) |
можно применить следующий искусственный метод. Пусть спект ральная плотность ошибок в несглаженном параметре представ ляется дробью
0^) = ^,
'S(u2)
где Q(co2) и S (со2) —многочлены от ®2. Тогда решение интеграль
ного уравнения (224) на интервале (0; 7Н) совпадает |
с решением |
линейного дифференциального уравнения |
(226) |
Q (— z1) Р (т) = const, |
|
где z — оператор дифференцирования, т. е. Z—d/dt. |
Значения про |
извольных постоянных могут быть вычислены при подстановке об щего решения уравнения (226) в интегральное уравнение (224) и условие (221).
Применим приведенное правило к определению оптимальной функции веса устройства, сглаживающего неизменные по гипотезе параметры движения цели. Пусть ошибки, искажающие координа ту, имеют функцию корреляции вида (214), т. е.
/?г 9 (т) = е-а И cos 2т.
При этом спектральная плотность ошибок в скорости изменения координаты или в параметрах движения цели определяется соотно
шением (216).
„ . , . _ |
q>2 (щ2 а2 + Q2) |
SH03!— |
(u4 + 2(a2_22)<o2 + (a2 4-Q2)2 • |
Этой спектральной плоскости соответствует функция корре
ляции
/?^(т) = 2аЕ2?8(т)-
- (a2 + 22) Efq е~Л N cos [2 hl - arc tg -^^-1. (227)
Здесь 3(т) — импульсная функция.