Главная -> Книги

(0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (74) (75) ( 76 ) (77) (78) (79) (80) (81) (82) (83) (84) (85) (86) (87) (88) (89) (90) (91) (92) (93) (94) (95) (96) (97) (98) (99) (100) (101) (102) (103) (104) (105) (106) (107) (108) (109) (110) (76)

Матрица корреляционных моментов ошибки вектора

E(t/t) = G(t)-G(t/t) в выражении (11.59) удовлетворяет векторному дифференциальному уравнению

Rh(OC(OP- it)C (OX

R£(0 = А(0 R£(0 + R£(0A(0

X Rh(0 + B(OQ(OB(0

(11.60)

с начальными условиями Re(0) =Reo-

Выражение (11.60) в теории оптимальных систем называют уравнением Риккати. Это нелинейное уравнение,

G(t)

fJ>

Рис. 11.5. Структурная схема оптимального непрерывного фильтра

которое в общем виде не решается. Поэтому вычисление элементов матрицы корреляционных моментов ошибки фильтрации сводится к решению системы из дифференциальных уравнений первого порядка, которая получается путем приравнивания в уравнении (11.60) элементов матрицы слева соответствующим элементам матрицы справа. Так как матрица Re (О симметричная, то число уравнений равно 0,5п(п+1), где п - порядок формирующего фильтра. Из найденных решений следует отобрать только то, при котором матрица ReH) положительно определена (см. приложение П. 4), такое решение является единственным.

В установившемся режиме А(0=А, 8(0 = В, С(/)==

= С, Р(0 = Р, Q(0=Q и R.e(0 = Re, поэтому Re=0 и уравнения (11.58), (11.59) и (11.60) принимают такой вид:

G (О = AG (О + К [ F (О - О (01; (и .61) К = Re CP-; (11.62)

О = ARf -f RePJ - Re СР-» С" Re -f BQBr . <11.63)



уравнение (11.63) в общем виде также не решается, поэтому элементы матрицы Re определяются из решения системы алгебраических уравнений, получаемых из 1(11.63) после выполнения умножения и сложения матриц в правой части уравнения (11.63). Из найденных решений необходимо отобрать только то решение, при котором матрица Re положительно определена.

Отметим, что передаточная функция оптимального фильтра, соответствующая уравнениям (11.61) и (11.62), совпадает с передаточной функцией, получаемой при синтезе оптимального фильтра по методу Винера.

Найденный оптимальный фильтр может быть использован для проектирования оптимальной системы РА. Если проектируемая система является счетно-решающим устройством, то при технической реализации могут быть применены однотипные интеграторы и сумматоры, используемые в вычислительной технике. Если проектируемая система предназначена для управления динамическим объектом, то реализация оптимальной системы сводится к определению структуры и параметров корректируемого устройства, подключение которого к объекту управления позволяет получить оптимальную систему. При стационарных сигнале и помехе реализация оптимальной системы упрощается. В этом случае матрицы в уравнениях (11.58) и (11.59) не зависят от времени, поэтому можно найти передаточные функции оптимальной системы и затем, используя методику гл. 7, определить передаточные функции корректирующих устройств, которые в данном случае будут стационарными.

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

Пример 11.4. Определить из условия минимума среднеквадратической ошибки в установившемся режиме структурную схему и параметры оптимальной системы, на которую действует смесь сигна-

* См.; Казаков И. Е. Статистическая теория систем управления в пространстве состояний.- М.: Наука, 1975.



ла И помехи со спектральными плотностями S;t(cu)=---I

Решение. Передаточная функция формирующего фильтра сигнала-это инерционное звеио с постоянной времени Тх и коэффициентом усиления, равным единице. Такой передаточной функции формирующего фильтра соответствует следующее дифференциальное уравнение:

где ао=Ь(,= 11Тх.

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

gi (0=-flogi(0 + boHO; gi (0 = (0-

Таким образом, в рассматриваемой задаче А=-оо; В=1; С=1; P=iV„; Q = blN.

В установившемся режиме из уравнения (11.61) следует, что Л л л

оценка сигнала g{t)=-aog\U)+4f{t)-eiit)].

Для вычисления коэффициента усиления оптимальной системы необходимо найти дисперсию ошибки фильтрации, для чего нужно решить уравнение (11.63), которое в данном примере следующее:

Rl+ 2aoNnRE-blNxNn = 0.

Согласно последнему квадратному уравнению, дисперсия ошибки фильтрации

R£ = Val Nl + bl JV„ -%N= {УТП - 1) NJT, где p=NxfNn.

Коэффициент усиления рассчитаем по формуле (11.62): fe=

На рис. 11.6 показана структурная схема оптимальной системы. Передаточная функция системы

1 + рТо

где йо=1-1 Т+р; 7-0=7-j/TTp.

Полученная передаточная функция оптимальной системы совпадает с передаточной функцией, найденной по методу Винера в примере 9.2.

Пример 11.5. Спектральные плотности сигнала и помехи характеризуются выражениями 5;t((o)=giVJ(o; 5n(cu) = iVn- Определить структурную схему и параметры оптимальной системы.

Решение. Формирующий фильтр состоит из двух интегрирующих звеньев, соединенных последовательно. Такому фильтру соответствует дифференциальное уравнение

где ц (О -белый шум интенсивностью Nx<



(0) (1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) (14) (15) (16) (17) (18) (19) (20) (21) (22) (23) (24) (25) (26) (27) (28) (29) (30) (31) (32) (33) (34) (35) (36) (37) (38) (39) (40) (41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) (65) (66) (67) (68) (69) (70) (71) (72) (73) (74) (75) ( 76 ) (77) (78) (79) (80) (81) (82) (83) (84) (85) (86) (87) (88) (89) (90) (91) (92) (93) (94) (95) (96) (97) (98) (99) (100) (101) (102) (103) (104) (105) (106) (107) (108) (109) (110)