
книги из ГПНТБ / Свешников А.А. Вероятностные методы в прикладной теории гироскопов
.pdf§ 9 . 2 ] |
ОП РЕДЕЛЕН И Е СИСТЕМЫ ПРИ БЕСКОНЕЧНОМ ВРЕМЕНИ |
477 |
§ 9.2. Определение оптимальной динамической системы при бесконечном времени наблюдения
Как было указано в § 9.1, в задаче Винера на вход динамической системы поступает сумма управляющего воздействия (полезного сигнала) U (t) и возмущающего воздействия (помехи) V (t), кото рые являются стационарными случайными функциями времени (будем считать й=г;=0). Рассматриваемая система должна на своем выходе воспроизводить некоторый требуемый сигнал Z (t), который связан с U (t) заданным функциональным соотношением
Z(t) = NU(t), |
(9.5) |
где N — некоторый преобразующий оператор. |
осуществимой, |
Оптимальная система должна быть физически |
т. е. ее реакция на единичную импульсную функцию должна обра титься в нуль для if 0. Следовательно, импульсная переходная функция I (t) (весовая функция) должна удовлетворять условию
l(t) — 0 при t < 0. |
(9.6) |
Погрешность е (t) системы определяется разностью между действительным выходным сигналом Y (t) и требуемым выходным сигналом Z (t), т. е.
s(t) = Y (t) — Z (t). |
(9.7) |
В связи с этим условие минимума средней квадратической ошибки может быть записано в виде
М [е2 (f)] = М {[F (і) — Z (f)]2} = min. |
(9.8) |
Задача синтеза оптимальной системы считается решенной, если определена импульсная переходная функция 1(х) или связан ная с ней взаимно однозначными соотношениями передаточная функция системы L (іш):
СО
L (іш) = J I (t) e~imidt,
|
|
° |
„ |
|
|
(9-9) |
|
|
|
eimiL (іш) d<n. |
|
||
|
|
|
—CO |
|
|
|
Спектральная |
плотность |
входной |
случайной функции |
X (t) |
||
в соответствии с |
(1 ) |
определяется |
соотношением |
|
||
s x И = |
S U («>) + |
S V (®) + |
S |
« , (№) + Sna(а»), |
(9.10) |
где Su (ш) и (о>) — спектральные плотности случайных процес сов U (f )n F (t) соответственно, а Suv ( ш) — их взаимная спектраль ная плотность,
§ 9 . 2 1 ОП РЕДЕЛЕНИ Е СИСТЕМЫ ПРИ БЕСКОНЕЧНОМ ВРЕМЕНИ |
479 |
Если Z (<) является результатом применения оператора N к текущим или будущим значениям ординат функции U (t), то I65]
/у (/«>) : |
1 |
(?„(<■>) |
1 И . |
(9.19) |
а'1 |
Рт Н |
|||
|
а |
Іг |
|
(9.20) |
|
|
|
|
г=1к=1
d lr~k |
(9.21) |
|
(lr-k)'- dJr~ |
||
|
||
где \ r — полюс кратности lr |
Ql (<•>) |
|
выражения p*'(-;)) SX!! (со), лежащий |
в верхней полуплоскости.
Если функция Z (t) является результатом применения опера тора N к ординатам функции U (t) в прошлые моменты времени, то
|
|
|
|
|
’ |
|
а2 |
^ п Н |
|
|
|
(9.22) |
|
|
|
|
|
Іг |
|
|
|
|
|||
|
|
|
|
|
аг |
|
|
|
|
|
|
|
|
|
|
*F (« |
|
_ |
ск |
_ |
|
|
|
(9.22) |
|
|
|
|
|
|
г—Л/.—] ѵ |
|
|
|
|
|
|
|
|
5*г — (/; _ к) ! |
Д*--* г, |
|
:«(<• |
|
|
|
(9.24) |
||||
|
уг_к[(“ |
КУГp ^ (J) s » Н ] |
о = Х - |
|
||||||||
где |
., |
|
кратности |
,, |
|
|
QZM с |
/ \ |
лежащий |
|||
Аг — полюс |
Zr выражения |
" |
' охг (ш), |
|||||||||
в нижней полуплоскости. |
|
|
|
|
являющаяся мини |
|||||||
|
Дисперсия ошибки оптимальной системы, |
|||||||||||
мальной для |
принятых исходных данных, |
определяется соотноше |
||||||||||
нием |
|
|
|
|
|
|
|
|
|
|
|
|
M[s2 (7)lmin = |
D[Z(0 ] - D ( F ( 0 ] = |
|
|
|
|
|
|
|||||
|
|
|
|
|
с о |
|
|
с о |
|
|
|
|
|
|
|
|
= |
J 5,(«o)Äo— |
J I /у (іш)(2 ^ (ш) |
(9.25) |
|||||
или |
|
|
|
|
|
|
|
|
|
|
(9.26) |
|
|
|
М |
[® 2 |
(7)]щіп= |
0 |
|
|
2 ( 0 |
dt, |
|
||
|
|
|
Кг ( ) |
|
2 л |
ß |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
где |
ß (t) определяется формулой |
(16). |
|
|
когда управ |
|||||||
|
Обобщение приведенных результатов на случай, |
ляющее и возмущающее воздействия приложены к различным точ кам системы, дано в [47], [49]. При этом рассматривается случай, когда кроме полезного сигнала U (t) и помехи V (t) на систему действует возмущение, точка приложения которого не совпадает с точками приложения U (t) ж V (t).
480 |
ВЫ БОР |
ОПТИМАЛЬНЫХ |
СХКМ И ПАРАМЕТРОВ |
ГУ |
[ГЛ. 9 |
|
§ 9.3. Определение оптимальной динамической системы |
||||
|
|
при конечном времени наблюдения |
|
|
|
|
Излагаемый |
ниже метод оптимизации динамических систем |
|||
отличается от |
рассмотренного |
в предыдущем параграфе |
двумя |
||
существенными |
особенностями: |
во-первых, время |
наблюдения |
конечное, т. е. система должна обладать, как говорят, конечной «памятью», реагируя только на входные данные, поступающие в си стему не ранее чем за время Т до настоящего момента времени. Во-вторых, входной полезный сигнал, помимо случайной состав
ляющей, содержит также составляющую в |
виде полинома g (t) |
|
от времени, т. е. математическое |
ожидание |
полезного сигнала |
g (t) может быть задано в виде многочлена от времени |
||
г (0 = 2 |
Ѵ г- |
(9-27) |
#= 0 |
|
|
В зависимости от имеющихся сведений о составляющей g (t) следует различать три случая: а) о коэффициентах кд ничего не известно; б) значения коэффициентов кд заданы; в) коэффициенты kq являются случайными и известны их вероятностные харак теристики.
Таким образом, в отличие от рассмотренного к § 9.2 случая Винера входная функция системы X (t) содержит неслучайное слагаемое g (t), т. е. выражается формулой
X ( t ) = g ( t ) + U ( t ) + V ( t ) , |
(9.28) |
где [g (t) -\-U ([)] и F (t) по-прежнему обозначают полезный сигнал и помеху, а функции U (t) и V (t) будем считать стационарными, несвязанными случайными функциями с нулевыми математиче скими ожиданиями.
На выходе динамической системы получаем некоторую вели чину Y (t). Для импульсной переходной функции системы I (<) вместо одного условия (6 ) для рассматриваемой системы с конеч ной «памятью» должны выполняться следующие два условия:
I (t) — 0 |
при |
t О 0, |
(9.29) |
I (t) = 0 |
при |
t^>T. |
(9.30) |
Первое из них, аналогичное (6 ), является условием физической осуществимости системы; второе — отражает сделанное предполо жение о невосприимчивости системы к сигналам, поступающим в систему ранее чем за время Т до настоящего момента времени (конечная «память» системы).
Входной сигнал X (t) должен быть преобразован в некоторый выходной сигнал Y (t), который должен, в определенном смысле
482 |
ВЫ БОР ОПТИМАЛЬНЫХ СХЕМ И ПАРАМЕТРОВ ГУ |
Ггл. 9 |
Рассмотрим первый случай, когда коэффициенты кд в выра жении (27) являются неизвестными.
Требование отсутствия систематической ошибки записывается в виде
‘с (0 = 0. |
(9.40) |
В соответствии с (38) есис (t) является полиномом от t степени г, обращение в нуль которого возможно только в том случае, когда все коэффициенты полинома обращаются в нуль. Приравнивая эти коэффициенты нулю, получим (г + 1 ) дополнительное условие, которому должна удовлетворять импульсная переходная функ ция I (т) оптимальной системы
Тсо
[т)г = ^ z4 (т) dz == |
j тѵп(т)йт |
(ѵ = 0, 1 , . .. ,r) . |
(9.41) |
0 |
— CD |
|
|
Следовательно, в рассматриваемом случае необходимо найти импульсную переходную функцию I (t) системы, обращающую в минимум величину М [е2 (£)] и одновременно удовлетворяющую (r-j-1) условию (41), т. е. нужно решить изопериметрическую за дачу вариационного исчисления (найти условный экстремум не которого функционала от искомой функции, на которую наложен ряд ограничений). В вариационном исчислении доказывается, что эту задачу на условный экстремум можно свести к задаче опре деления безусловного экстремума функционала I, если положить
I = |
М [е2 (t)] — 2у0т]0 — ... |
— 2угт)г, |
(9.42) |
где постоянные у0, |
у15 . . . , уг — так |
называемые |
множители |
Лагранжа. |
|
|
|
В случае, рассмотренном Заде и Рагаззини [90], предполагается, что спектральная плотность Sx( u>) случайной функции X (t) — ■—g (t) представляет собой дробно-рациональную функцию от со, т. е.
|
|
С л ,Л |
|
Л , І ^ И 1 2 _ |
д а Р „ И ^ Н |
(9.43) |
|||
|
|
|
|
||||||
где |
полиномы |
Pm(о>), |
Qu(a>), |
Р*т (<о) и (?*(«>) |
по-прежнему опреде |
||||
ляются формулами (18), а |<2 , » |
| и |
о)Р |
полиномы степени п |
||||||
и соответственно т от |
|
с постоянными коэффициентами. |
|||||||
|
При сделанных допущениях оптимальная импульсная переход |
||||||||
ная |
функция |
I (і) имеет |
вид |
|
|
|
|
||
|
|
2 т |
|
|
|
|
|
|
|
l(t) = H D/ j + 2 С.еѴ + 2 |
[А 8 (м) (0 + |
Вд§<«’ (t - |
Т)] + |
||||||
|
j-o ' г |
' Р а ~ г ‘ |
' Р L‘ г |
|
|
|
|||
|
+ |
1 <?„М12 |
N (im)SM (іо)etatdw |
( 0 < і < Г ) , |
(9.44) |
||||
|
ЛпИІS |
2 |
484 |
ВЫ БОР ОПТИМАЛЬНЫХ СХЕМ И ПАРАМЕТРОВ ГУ |
[ГЛ. 9 |
Расчеты по определению оптимальных импульсных переходных функций динамических систем с конечной «памятью» существенно упрощаются при использовании приведенных в [ 70 ] таблиц, составленных для различных выражений спектральных плотно стей случайной составляющей полезного сигнала и помехи.
§9.4. Определение оптимальных параметров системы
сзаданной структурой
Применительно к гироскопическим устройствам широкое при менение может получить следующая задача синтеза: известна структурная схема ГУ и его динамические характеристики, опре деляющиеся дифференциальными уравнениями или соответствую щими передаточными функциями; требуется определить оптималь ные значения параметров ГУ. Решение подобной задачи может потребоваться, например, при выборе оптимальных значений пара метров гиротахометра, акселерометра, поплавкового интегрирую щего гироскопа и других ГУ, структурная схема которых задана.
При определении оптимальных параметров системы с заданной структурой речь идет не о выборе всех параметров системы, а лишь тех из них, которые могут быть легко изменены в широких преде лах и которые существенно влияют на точность системы.
Обозначим через р15 р2, . . ., уя параметры системы. Поскольку структура системы (передаточная функция линейной системы) и вероятностные характеристики входного сигнала и помехи из вестны, то критерий качества системы, например дисперсия ее ошибки е (t), может быть выражен в виде функции параметров системы
D [> (01 = F (рі, р2, .. .,рл). |
(9.49) |
Оптимальные значения параметров системы определяются из системы уравнений, которая получается приравниванием нулю частных производных критерия качества по неизвестным пара метрам системы
(9.50)
Решения ^ = |а°лт, ц2 = р.®лт, . . ., |ія = р.°пт системы уравне ний (50) определяют оптимальные параметры рассматриваемой динамической системы.
Указанный подход был применен Р. Филиппсом [19] для определения оптимальной динамической системы, передаточная функция которой выбрана в виде дробно-рациональной функции с неизвестными коэффициентами, значения которых определя лись при оптимизации качества рассматриваемой системы.
4 8 6 |
Е Е В Ы БО Р |
ОПТИМАЛЬНЫХ СХЕМ И ПАРАМЕТРОВ ГУ |
[ГЛ . 9 |
|
Как известно, |
корректирующее устройство может быть после |
довательным или параллельным. На рис. 9.1, а показано коррек тирующее устройство с передаточной функцией LK(s), включен ное последовательно с объектом (с неизменяемой частью системы),
имеющим передаточную функцию L0 (s); на рис. |
9.1, |
б |
показано |
||||||
параллельное |
включение корректирующего устройства. |
Предпо |
|||||||
|
|
лагается, |
что оптимальная |
||||||
LH(s) |
LJs) |
передаточная функция L (s) |
|||||||
замкнутой системы опреде |
|||||||||
|
|
||||||||
|
|
лена, например, с помощью |
|||||||
|
|
методов, |
§ |
излагающихся |
|||||
|
а) |
в § 9.2 |
и |
9.3. |
Зная пе |
||||
|
редаточную функцию L0 (s) |
||||||||
|
|
неизменяемой части систе |
|||||||
|
|
мы, требуется |
определить |
||||||
|
|
передаточною |
|
функцию |
|||||
|
|
Ls(s) |
корректирующего |
||||||
|
|
устройства, |
его структур |
||||||
|
|
ную схему и параметры. |
|||||||
|
|
Эта |
задача |
|
решается |
||||
|
6) |
сравнительно |
просто |
для |
|||||
|
стационарных |
линейных |
|||||||
Рис. 9.1. Схемы |
последовательного (а) и |
||||||||
систем. Принимая во внима |
|||||||||
параллельного (б) включения корректирую |
ние введенные выше обоз |
||||||||
щих устройств. |
|||||||||
|
|
начения, можно найти сле |
|||||||
|
|
дующие |
|
выражения |
для |
передаточных функций корректирующих устройств при последо вательном и параллельном их включении
т |
L (s) |
(9.52) |
|
г |
_Ln (s) L (s) |
(9.53) |
|
> |
L ( s ) L 0(s) |
||
’ |
Обычно оптимальная передаточная функция L (s) системы является дробно-рациональной или трансцендентной. В последнем случае и передаточная функция корректирующего устройства также будет трансцендентной, и ее не представляется возможным реализовать с помощью элементарных звеньев. Для приближенной реализации передаточной функции корректирующего устройства ее необходимо аппроксимировать дробно-рациональной функцией. Для этого обычно пользуются методом логарифмических частотных характеристик [70]. Существуют также различные приближенные аналитические методы аппроксимации трансцендентных функ ций (см., например, [36]).