Дискретная система оценивания

 

Изобретение относится к радиотехнике и автоматике и может быть использовано в радиолокации для оптимальной линейной оценки параметров движения летательных аппаратов. Цель изобретения - повышение динамической точности системы. Дискретная система оценивания содержит первый вычитатель 1, блок 2 коррекции, сумматор 3, первый блок 4 задержки, блок 5 экстраполяции , второй 6 и третий 7 вычитатели, второй блок 8 задержки, усилитель 9, счетчик 10, двухканальный мультиплексор 11, компаратор 12; умножитель 13, сумматор 14, элементы 15-17 памяти коэффициентов. Цель изобретения достигается за счет введения элементов и блоков 6-17. 3 ил.

СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК (я)5 G 05 В 17/00, G 05 F 1/00

ГОСУДАРСТВЕННЫЙ КОМИТЕТ

ПО ИЗОБРЕТЕНИЯМ И ОТКРЫТИЯМ

ПРИ ГКНТ СССР, ",)О il i! g

ОПИСАНИЕ ИЗОБРЕТЕНИЯ л

X(n) X(n) К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ (21) 4756238/24 (22) 03.11.89 (46) 15.04.92. 6юл. М 14 (72) Ю.А.Пушкарев, А.А.Шалаевский и

В.И.Ищенко (53) 62-50(088.8) (56) Гришин Ю.П., Казаринов Ю,М. Динамические системы, устойчивые к отказам. М.:

Радио, 1985, с. 19. (54) ДИСКРЕТНАЯ СИСТЕМА ОЦЕНИВАНИЯ (57) Изобретение относится к радиотехнике и автоматике и может быть использовано в

° ЯЛ 1727110 А1 радиолокации для оптимальной линейной оценки параметров движения летательных аппаратов. Цель изобретения — повышение динамической точности системы. Дискретная система оценивания содержит первый вычитатель 1, блок 2 коррекции, сумматор 3, первый блок 4 задержки, блок 5 экстраполяции, второй 6 и третий 7 вычитатели, второй блок 8 задержки, усилитель 9, счетчик 10, двухканальный мультиплексор 11, компаратор 12; умножитель 13, сумматор 14, элементы 15 — 17 памяти коэффициентов. Цель изобретения достигается за счет введения элементов и блоков 6-17. 3 ил.

1727110

Изобретение относится к радиотехнике и автоматике, в частности может использоваться в радиолокации для оптимальной линейной оценки параметров движения л етател ьн ых объектов.

Известны дискретные системы оценивания, предназначенные для оценивания параметров движения объектов.

Недостатком систем является отсутствие возможности качественно управлять динамической точностью в переходном режиме.

Известна также дискретная система оценивания, содержащая последовательно соединенные элемент сравнения, блок коррекции, сумматор, блок задержки, блок экстраполяции, блок обратной связи, Данная система имеет возможность оптимально оценивать параметры движения объектов в установившемся режиме. Недостатком данной системы является невысокая динамическая точность в переходном режиме.

Наиболее близкой к предлагаемой является дискретная система оценивания, которая имеет воэможность оптимальным образом оценивать параметры движения обьектов в установившемся режиме. Она содержит последовательно соединенные элемент сравнения, блок коррекции, сумматор, блок задержки, блок экстраполяции, выход которого соединен с вторыми входами элемента сравнения и сумматора.

Недостатком данной системы является низкая динамическая точность в переходном режиме.

Цель изобретения — повышение динамической точности системы в переходном режиме.

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

На фиг. 1 изображена структурная схема системы оценивания; на фиг. 2 — структурная схема блока экстраполяции.

Дискретная система оценивания (фиг. 1) содержит вычитатели 1, 6 и 7, блок 2 коррекции, сумматоры 3 и 14, блоки 4 и 8 задержки с операторами задержки z, блок 5 экстраполяции, усилитель 9, счетчик 10, двухканальный мультиплексор 11, компаратор 12, умножитель 13, элементы 15 — 17 памяти коэффициентов, В статическом состоянии система оценивания характеризуется следующими связями: содержит последовательно соединенные вычитатель 1, блок 2 коррек5

55 ции, сумматор 3, блок 4 задержки, блок 5 экстраполяции. Первый вход вычитателя 1 и выход сумматора 3 являются соответственно информационным входом и выходом дискретной системы оценивания. Кроме того, система оценивания содержит вычитатели 6 и 7, сумматор 14, элементы 15 — 17 памяти коэффициентов, счетчик 10, компаратор 12, усилитель 9, умножитель 13, блок 8 задержки и двухканальный мультиплексор 11. Выход вычитателя 1 через последовательно соединенные блок 8 задержки и усилитель 9 подключен к первой группе информационныхх входов мультиплексора 11, соединенного второй группой информационных входов с выходами. элементов 15 — 17 памяти коэффициентов, адресным входом — с выходом счетчика 10 и первым входом компаратора

12. Первый выход блока 5 экстраполяции подключен к первому входу умножителя 13, соединенного вторым входом с первым выходом мультиплексора 11, а выходом — с первыми входами сумматора 14 и вычитателя 6. Второй выход блока 5 экстраполяции подключен к первому входу вычитателя 7 и второму входу вычитателя 6, соединенного выходом с вторым входом вычитателя 1.

Второй выход мультиплексора 11 подключен к второму входу сумматора 14, соединенного выходом с вторым входом вычитателя 7, выход которого подключен к второму входу сумматора 3. Второй вход компаратора 12 соединен с шиной задания числа шагов, а выход — с информационным входом счетчика 10, тактовый вход которого является тактовым входом дискретной системы оценивания.

При работе устройства возможны два режима работы; переходный и установившийся.

В установившемся режиме устройство работает следующим образом.

Входной сигнал системы Х(п), представляющий собой сумму детерминированной составляющей вх(п), центрированной случайной составляющей X(n) и ошибок измерения f(n) (помехи), поступает на первый вход вычитателя 1, на его второй вход подается экстраполированное значение входного сигнала Хэ(п). В установившемся режиме блоки 8-17 системы не работают, поэтому информация на первом входе вычитателя 6 и втором входе вычитателя 7 отсутствует. В этом случае экстраполированное значение входного сигнала X>(n), сформированное блоком 5, проходит через вычитатели 6 и 7 без изменений. На выходе вычитателя 1 имеем сигнал невязки 0(п), который поступает на блок 2 коррекции. На выходе сумматора 3 имеем сигнал оценки Х(п):

1727110

30 (4) +U(Z)(1-bo(n)) = Xe(Z) +

55 (7) +U(Z)(1 — bp(n)), -2Ь1(п)) Z + Ь1(п)2, (9), л л л

X(n) = 2Х(п-1)-X(n-2)+ K(n) U(n). (1)

Первые два члена в правой части выражения (1) являются сигналом Хэ(п), сформированным блоком 5 экстраполяции. Блок 5 экстраполяции синтезирован на основании 5 известного метода "трех полиномов" из условия: входное воздействие на систему имеет вид

X(n) = mx(n) + X(n) + f(n), (2)

10 где mx(n) — детерминированная составляющая задающего воздействия X(n);

Х(п) — его центрированная случайная составляющая;

f(n) — центрированная возмущающая со- 15 ставляющая (ошибки наблюдения).

Для исключения динамической ошибки . в дискретной системе при оценивании коор- . динат объекта, описываемого выражением (2), необходимо иметь порядок астатизма 20

v = 2, а следовательно из метода "трех полиномов" определим полиномы

A(Z,n) = ao(nX1 — 2 ) и

B(Z,п) = bp(n)(1 Z ). (3) Характеристический полином выбираем из условия устойчивости системы.

C(Z,n) = ао(п) — 2bp(n)Z + bp(n)Z

Используя выражение для оценки метода "трех полиномов" (.(Z,и1- S(Z,n)

Х(2 = u(z) ) (5)

А(2,и) получим ар(и)-2bo(h)Z +Ьр(Ы2 -Ъ (n)+2 (n)Z - Ьр(и)2

)((z)

ao(n)-Qa,(n)Z +а,()Z

ao(nl — ЬО(ь) . (--Ос - (6), U(Z) = tJ(z) 6

aр(й) 2aр()Z +0 (г ) {-2Z +Z

4М ао и где U(n) — невязка, а— о и

После преобразования (6) получим выражение для оценки в операторном виде 50

X(Z) = 2X(Z)Z — X(Z)Z +

/ где Xs(2) = 2X(Z)Z X(Z)Z — экстраполированное значение координаты;

0(2)(1-Ь (и)) = (.)(2)К(п) — корректировочное значение.

Блок 5 экстраполяции, структурная схема которого изображена на фиг. 2, содержит блоки 18 и 19 задержки с оператором задержки Z усилители 20 и 22 с коэффициентами усиления "два" и вычитатели 21 и 23.

На основании проведенного синтеза уравнения(7) построена часть блока 5, включающая в себя блок 18 задержки, усилитель

20, вычитатель 21. Эта часть блока 5 используется устройством при работе как в установившемся, так и в переходном режиме.

Блоки 19, 22 и 23 используются устройством при работе только в переходном режиме.

Входным сигналом для блока 5 экстраполяции является сигнал с выхода блока задержки Х(п-1), Сигнал с первого выхода блока 5 поступает на первый вход умножителя 13, а сигнал с второго выхода в виде

X (n) — на первый вход вычитателя 7 и на второй вход вычитателя 6, через которые поступает соответственно на вторые входы сумматора 3 и вычитателя 1. Сигнал Х(п) формируется в виде (1), сигнал невязки U(n) — в следующем виде: ю/ и л

U(n) = Х(п) — 2X(n — 1) + X(n-2). (8) Прежде, чем описывать работу устройства в переходном режиме, необходимо отметить, что в замкнутых системах оценивания, использующих для сглаживания (фильтрации) и управления только одни и те же весовые коэффициенты, существует противоречие между условиями уменьшения флюктуационных ошибок (условиями фильтрации) и условиями повышения динамической точности в переходном режиме.

Это противоречие может быть разрешено с помощью введения дополнительных весовых коэффициентов по каналу управления.

Тогда можно добиваться высокого качества переходного процесса как по каналу управления, так и по каналу оценивания, не увеличивая при этом флюктуационные ошибки.

В связи с изложенным полиномы A(Z,ï), B(Z,n), C(Z,n) представим в виде

A(Z,n) =(1-2 .1 (1+ а)(п)2 ), В(2,п) = (1 — 2 ) (Bo(n) + Ь1(п)2 1, C(Z, п) = 1 + (31(n) — 2 Ьp(n))Z 1 + (bp(n) 1727110 (10) (мин = Ч1мин + Чгмин + ° ° °

U(n) = X(n) — Хэ(п), (15) 10 имеем

A(Z,n1 (е1) X (n).

Из(11) при условии, что (12) имеем

C(-1)(Z,n) U (n) = U,),(n), A(Z,n) Хо(п) = (.4ын(п), (13) Так как Г1+ F1= 1, то (16) Поэтому

При оценивании процесса(2) обеспечим требуемое качество переходного процесса по критерию где V1M» = мин U „(1), Чгм» = мин U „(2) и г " г т.д„

U() — переходная составляющая ошибки наблюдения.

Для отыскания оптимальных коэффициентов а1 р (и) обратимся к методу "трех полиномов" для невязки

Хо, ЕСЛИ и -0

Х() = Х (п) „,й О

Up(n) = С(-1)(Z,n) Up(n) + A(Z,n) Хо(п), где полином C(-1)(Z,n) — это полином C(Z,n) без свободного слагаемого.

Следовательно, есть свободная и вынужденная составляющие динамической ошибки в переходном режиме. И та и другая составляющие зависят от коэффициентов а1(п), причем противоречивым образом. Поэтому существуют оптимальные значения коэффициентов а1оР (n), удОВЛЕтВОряЮщИЕ КрИтЕрИЮ (10).

Для отыскания минимальных значений

VM» критерия качества I вводятся весовые коэффициенты Г и F по свободной и вынужденной составляющей соответственно.

Из (11) с учетом (12) и (9) получим следующее разностное уравнение:

° U(n)< = Ь U(n 1)п AU(n 2)п f3U(n 3)п +

Xp(n) + (а 1(п) - 2) Xp(n — 1) + (1 — 2а1(п)) Xp(n — 2) + а1(п) Хо(п — 3), (13)

ГдЕ Ь = a1(n)-2Ьо(П), (г = bp(n) — 2Ь1(п), (14) (з = b1(n).

Принимая во внимание (12), найдем переходную составляющую ошибки в начальный момент времени.

Учитывая, что канал наблюдения описы5 вается выражением по методу "трех полиномов", U(0). = Х(О)-Х.(О) = Х.(О)-О = Х..

Так как в первую часть уравнения (13)

15 входит множитель Хо(п) Ф О при n = 1, 2, 3, то из данного выражения получим следующие три уравнения: и()„=-Е,й(а)„.(а,щ-2) Х.(o), 2О ц(р „=-l,Ц(1„-6,О(О „+(.(-йа (91(,(О), ()()„=-Я,ц(Я1„-(дЯ(<)„-0 0(0)„+а,(ъ) х,(0, Здесь слагаемые, содержание Un в правой части, образуют свободную составляю25 щую ошибки, а Хо(0) — вынужденную.

Для квадратов ошибок имеем ц (()„=jL,(11ц(01„) ц,+(с(,( (((2l„.((е,«(),ß2,é(o),) )(+

+ t(— а,(65 х, (о1, (U((ú),-((е,«д,) (е "(i).)

35 (Е, U(o) „) I(..; ; a,(з) х,(o) Fâ.

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

40 получается относительно предыдущего состояния, полученного оптимальным образом, На первом шаге оптимизации имеем

U (il„=U, мак((а,(1)-2b.(<)) X.(O)r,+ а,(1)

+(а,(Ц-2) X (o) Г,)

OO а.,(„=X (O)(2(o (×-2>.«))«2(а,(П-2)фа. а1 Р (1) = 2 (Ьо(1) Г1+ F1).

U п(1) р = (а1 р (1) — 2Ьо(1)) Х p(0)l + (а1ор (1) — 2)г Х p(0)F1.

Второй шаг

1727110

10 (17) Третий шаг (20) ll ф„-Н -мин ((а,(21-2Ь.(й)) Й И„ +

a,(i)а,(2)

+(b,(21-2Ь,(2)) Х, ÎI)r (i-2а,(2)) х,(O1F„ о . а 2

--R(a,(21-2ь,(2))Й (цд Г

+a((-2a,®3x оp(0)(-aF,) =D Р а1ь,(6Ю () „P Ã,. 4,3

° ь(0 (Ь)„=Чз- мин ((а,(3) -2Ь (М) (у (2) Р«

«(Ь,(Н-2Ь,(э)) О () р +Ь,(э)х, (а) а,())((p)F

3у

Ba,(ц= Р ьИ ("-). (ь а () ((о)г,=О ор1 2.bo()u (2) „ Гз ° 7, (1 8) При выбранных полиномах А(2,п), B(Z,n), C(Z,ï) коэффициент а1(п) не входит в управление. Это приводит к ухудшению качества фильтрации. Чтобы снять этот недостаток, необходимо в разнице полиномов

C(Z,n)-B(Z,ï) исключить слагаемые, содержащие коэффициенты а1(п). Это можно сделать лишь в том случае, если полином C(z,ï) и B(Z, N ) включают в себя коэффициенты

a<(n).

Если полиномы A(Z,n) и B(Z,ï) задать в виде

A(z,n)=((-z ) (1 а,()z- j, 5(z,è) =(i- z ) (b,(n)+b|(n)z )(ha<(n)z ), (19) то из выражения по методу "трех полиномов" для C(z,ï) Ь(2,в — Ь(2, ) э(7)

С(2,n)— э(7) найдем выражения для алгоритмов оценивания и управления в переходном режиме.

Согласно (20) C(Z,п) — BZ, и) = 1 — bp(n) — Ь )(n)z, (21) Поэтому выражение для оценки примет

BRA h и и и

Х(п) =2X(n — 1) — X(n — 2) — а1(п)Х(п — 1) +

+ 2а11(п)Х(п — 2)-а1(п)Х(п-3) +

+K(n)U(n) — Ь1(п)(.)(п-1).. (22) 5 Для того, чтобы система была устойчивой, необходимо выполнить условие (1 — 2 ) (1 + a >(n)Z ) = 1-)=,(2).

10 Здесь выражение для функции F (z) (программы управления) будет иметь вид

F,(2) = (2 — a)(n)2 1+ (2а1(п)-1)Z -a>(n)2

15 и, следовательно, алгоритм управления примет вид п h h

Хэ(п) = 2X(n-1)-X(n-2)-а 1(п)Х(п-1) + л A

20 + 2а1(п) Х (n — 2)-а1(п)Х(п-3). (23) Таким образом, выбор полиномов А и В в соответствии с выражением (19) создает условия для оптимальных переходных про25 цессов без нарушения условий фильтрации

В переходном режиме устройство работает следующим образом, В момент включения устройства на вход системы поступает ступенчатое воздейст30 вие вида (12), обусловленное априорной неопределенностью состояния обьекта. Для того, чтобы улучшить качество переходного процесса заранее по формулам (16) — (1 8), вычисляются коэффициенты а1(1), а1(2), а1(3) и

35 записываются в блоки 15-17. Количество шагов, равное трем, определяется порядком характеристического полинома при ступенчатом воздействии.

Процесс управления в переходном ре40 жиме реализуется устройством при помощи блока 8 задержки, усилителя 9 с коэффициентом усиления Ь1(п), счетчика 10, двухканального мультиплексора 11, компаратора

12, умножителя 13, сумматора 14, элемента

45 15 хранения коэффициента а1(1), элемента

16 хранения коэффициента а1(2), элемента

17 хранения коэффициента а1(3).

Перед началом работы устройства на второй вход компаратора 12 подается код N

50 числа шагов работы системы в переходном режиме, Для рассматриваемого случая N =

3. На вход счетчика 10 поступает последовательность импульсов пТ, синхронная моментам выдачи информации на обработку.

56 Выходные сигналы счетчика 10 поступают на адресный вход мультиплексора 11 (в качестве которого может быть использована, например, микросхема К155КП2). В первом такте коэффициент а1(1) считывается из элемента 15 памяти и проходит через мульти1727110

35

50

55 плексор 11 на второй вход умножителя 13, в котором умножается на сигнал X(n-I) — 2X(n2)+ Х(п — 3), сформированный блоками 18, 19 и 20 и вычитателем 23 блока 5. На выходе умножителя 1 формируется сигнал

a1(1)X(n-1) — 2а1(1)Х(п — 2) + а1(1)Х)п — 3). B тоже время по второму каналу мультиплексора 11 поступает сигнал b<(1)U(n-1), сформированный блоками 8 и 9, на второй вход сумматора 14. Информация с выхода умножителя 13 поступает на первый вход сумматора 14 и первый вход вычитателя 6. На выходе сумматора 14 имеем сигнал л и л а (1) Х (n — 1) — 2а (1)Х(п — 2) + a1(1)X(n — 3) +

B результате на выходе вычитателя 6 имеем сигнал управления (23), а на выходе вычитателя 7 сигнал

21(— 1) — X(n — 2)-а <(1)X(n-1) + 2a1(1)X(n2)-à <(1)X(n — 3) — Ь (1)0(п — 1), который принимает участие в формировании сигнала оценки (22). В следующих двух тактах формирование рассмотренных сигналов осуществляется аналогично для коэффициентов а1(2), а (3), Ь1(2), Ь (3).

Коэффициент Ь1(п) определяется из условий устойчивости системы и в рассматриваемом случае является постоянным, т.е. Ь (1) = Ь (2) = Ь (3). При достижении счетчиком 10 числа

N компаратор 12 вырабатывает сигнал, блокирующий работу счетчика 10. Работа системы в переходном режиме на этом заканчивается. Количество шагов оптимизации N рассчитывается на основе принципов оптимальности Белмана, уравнения (16)(18).

С целью подтверждения эффективности предлагаемой системы было проведено цифровое моделирование на ЭВМ. Для расчета коэффициентов a>(n) принимались Х =

30 м, Т = 0,1 с. Результаты моделирования приведены на фиг. 3, где пунктиром показаны ошибки оценки в переходном режиме для системы с а (п) = О, а сплошной — с использованием а (п) О, где коэффициенты а (п) принимали значение; 2; 0,615; 2 соответственно, Данным значениям соответствуют весовые коэффициенты Г = F=

= 0,5.

Анализ результатов моделирования показал, что применение блоков 8 — 17 позволяет повысить динамическую точность системы в переходном режиме. Время регулирования переходного процесса уменьши5

30 лосьсср=3.6сдотр=2 5с,а перерегулирование — с о =. 80 до 40 . Ошибка оценки еуменьшилась с 24 м до 12 м. Видно, что для предлагаемого устройства динамическая точность в переходном режиме повысилась в два раза (при заданных условиях), что особенно важно в системах обработки информации, в частности, при выдаче целеуказанной системе оценивания при захвате цели.

Формула изобретения

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

1727110

:-IO

Составитель N3.nóøêàðåâ

Техред М.Моргентал Корректор С.Лыжова

Редактор Э.Слиган

Заказ 1278 Тираж Подписное

ВНИИПИ Государственного комитета по изобретениям и открытиям при ГКНТ СССР

113035, Москва, Ж-35, Раушская наб., 4/5

Производственно-издательский комбинат "Патент", г. Ужгород, ул.Гагарина, 101