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

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

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

§ 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(т) — импульсная функция.

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