Адаптивное устройство для обработки данных

.Бурыи Бонджиова ание фильт ичии неопре мы, - Ракет м 11, М 5,1 Ап полито оектиров па при нал ров систе автика, то на973 54) АДАПТИ АБОТКИ ДА ДЛЯ ОБНОЕ УСТРОИНЫХ СУДАРСТВЕННЫЙ КОМИТЕТИЗОБРЕТЕНИЯМ И ОТКРЫТИИ ГКНТ СССР СПИСАНИЕ ИЗО К АВТОРСКОМУ СВИДЕТЕЛЬС(56) Хачинсон, ДМинимаксное прКалмановского тиленности параметтехника и космонс. 150-157. Ж, 1741 О 98 А 1(57) Сущность изобретения: устройство содержит блок 1 элементов И, минимаксный фильтр 2 Калмана, блок 3 расчета коэффициента коррекции, блок 4 расчета апостериорной матрицы ковариации, два блока 5 и 25 элементов задержки, девять умножителей 6, 9, 10, 15, 17, 19, 22. 23 и 24, четыре сумматора 11, 12, 13 и 18, блок 14 памяти матрицы коэффициентов системы, блок 16 памяти измерительной матрицы и матрицы шумов измерений, блок 20 памяти границы неопределенности переходной матрицы, блок 21 памяти предполагаемой переходной матрицы, блок 26 управления, блок 27 памяти переходной матрицы системы.4 ил,Изобретение относится к линейным устройствам калмановского типа для обработки данных в реальном масштабе времени в условиях неопределенности задания априорной информации.Известен фильтр Калмана для обработки данных, содержащий сумматоры, блоки умножения, деления и блоки памяти.Недостатком этого фильтра является низкая точность обработки данных вследствие расходимости .оценок вектора состояния при наличии неопределенности в задании априорной информации о моделях систем и измерений,Наиболее близким по технической сущности к предлагаемому устройству является устройство, содержащее дискретный минимаксный фильтр Калмана, блок расчета коэффициента коррекции, блок расчета апостериорной матрицы ковариации, первый блок оперативной памяти, первый, второй и третий блоки памяти, первый, второй,. третий и четвертый блоки умножения и первый сумматор, первый вход которого подключен к выходу четвертого блока умножения, первый вход которого подключен к ыходу второго блока памяти, а второй вход - к выходу третьего блока умножения, первый и второй входы третьего блока умножения подключены соответственно к выходам третьего и второго блоков памяти, второи вход первого сумматора подключен к выходу второго блока умножения, первый вход которого поключен к выходу первого блока памяти, а второй вход - к выходу первого блока умножения, первый вход блока умножения подключен к выходу первого блока памяти, а второй - к выходу первого блока оперативной памяти, вход первого блока оперативной памяти подключен к выходу блока расчета апостериорной матрицы ковариации, первый вход которого подключен к выходу блока расчета коэффициента коррекции, выход первого сумматора подкл ючен к входу блока расчета коэффициента коррекции, выход которого подключен к первому входу дискретного минимаксного фильтра Калмана.Известное устройство проводит обработку данных в соответствии с формулойХ(1 с/1 с) = Ф(1 с, 1 с - 1) Х (1 с - 1/1 с - 1)++ Км(1 с/1 с) (м Н Ф(1 с, 1 с - 1) Хм (1 с - 1/1 с - 1, (1)где Хм(1 с/1 с) - минимаксная оценка векторасостояния в 1 с-й момент времени;Н - измерительная матрица; Ф(1 с,1 с) - переходная матрица системы;Е 1 с - вектор измерений в 1 с-й момент времени;5 Км(1 с/1 с) - минимаксный матричный коэффициент усиления (коррекции), определяемый соотношением Км(1 с/1) =(Мк/к)м Н ( Н (Мк/с)м Н+ В) . (2) 10где В - матрица интенсивности шумов измерения;(Мкуп) - минимаксная априорная матрица ковариации ошибки оценки;15 знак Т - транспонирование матриц,Его недостатком является то, что обработка данных осуществляется с учетом неопределенности задания характеристик шумов системы и измерений и не учитыва ются неопределенности в задании переходной матрицы системы Ф(1 с,1 с). Поэтому при наличии неопределенности в задании переходной матрицы системы точность обработки данных начинает ухудшаться, а дискретный 25 минимаксный фильтр Калмана расходится,Цель изобретения - повышение точности устройства.Поставленная цель достигается тем,что в адаптивное устройство для обработ ки данных, содержащее последовательносоединенные блок памяти матрицы интенсивности шумов системы, первый и второй умножители, первый сумматор, блок памяти матрицы коэффициентов системы, подклю ченный выходом к вторым входам первого ивторого умножителей, блок памяти и редполагаемой переходной матрицы, соединенный выходом через третий умнокитель с первым входом четвертого умножителя и не посредственно с вторым его входом, выходчетвертого умнокителя подключен к второму входу первого сумматора, минимаксный фильтр Калмана, соединенный первым входом с выходом блока расчета коэффициента 45 коррекции и первым входом блока расчетаапостериорной матрицы ковариации, второй вход которого подключен к первому входу блока расчета коэффициента коррекции и первому выходу блока памяти измерительной 50 матрицы и матрицы шумов измерений, третий вход - к второму выходу блока памяти измерительной матрицы и матрицы шумов измерений и вторым входом блока расчета коэффициента коррекции и минимаксного 55 фильтра Калмана, а выход через первый блокэлементов задержки - к второму входу третьего умножителя, введены блок элементов И, квадратор, три сумматора, с пятого по девятый умножители, блок памяти границы неопределенности нормы переходной матрицы, блок памяти переходной матрицы системы, второй блок элементов задержки, первый и второй входы второго сумматора соединены соответственно с выходами первого и третьего сумматора, а выход - с входом второго блока элементов задержки, третьим входом блока расчета коэффициента коррекции и четвертым входом блока расчета апостериорной матрицы ковариации, выход блока памяти предполагаемой переходной матрицы через пятый умножитель подключен к первому входу четвертого сумматора, к второму входу которого через шестой умножитель подключен выход блока памяти границы неопределенности нормы переходной матрицы, соединенный также через последовательно соединенные квадратор и седьмой умножитель с первым входом третьего сумматора, подключенного вторым входом к выходу четвертого сумматора, выход второго блока элементов задержки соединен с вторым входом седьмого умножителя, выход блока памяти предполагаемой переходной матрицы подключен также к первому входу девятого умножителя, соединенного вторым входом с выходом первого блока элементов задержки и первым входом восьмого умножителя, а выходом - с вторым входом шестого умножителя, второй вход пятого умножителя подключен к выходу восьмого умножителя, соединенного вторым входом с выходом блока памяти границы неопределенности нормы переходной матрицы, первый выход блока управления подключен к управляющему входу блока расчета коэффициента коррекции, первому входу блока элементов И и входу разрешения считывания блока памяти переходной матрицы системы, соединенному адресным входом с вторым выходом блока управления, а выходом - с третьим входом минимаксного фильтра Калмана, подключенного четвертым входом к выходу блока элементов И, второй вход блока элементов И и выход минимаксного фильтра Калмана являются соответственно информационным входом и выходом устройства,Такая совокупность признаков предлагаемого устройства позволяет, по сравнению с прототипом, повысить точность обработки данных в условиях неопределенности задания переходной матрицы системы за счет учета неопределенности в задании переходной матрицы системы Ф(К Е - 1) и дополнительного введения в устройство второго блока задержки, двух блоков памяти, квадратора, блока элементов И,. пяти умножителей, блока управления, второго, третьего, четвертого сумматоров, соединенных особым образом, указанным в формуле изобретения.Указанные существенные отличия позволяют на каждом шаге оценивания рассчиты вать минимаксную априорную матрицуковариации ошибки оценки с учетом неточности задания переходной матрицы системы, При этом повышается вес каждого текущего измерения, в силу чего устраняется расходи масть дискретного минимаксного фильтраКалмана.Устройство осуществляет расчет минимаксной априорной матрицы ковариации по формуле15+ Ь Ф(М-а-)м (4)25Г - матрица системы;Л - положительное число, определяющее верхнюю границу неопределенностинормы матрицы перехода, а именно 30аФ Ь(Мк/к)ч - минимаксная апостериорная матрица ковариации ошибки оценки, опре 35 деляемая выражением(Мо)м =- (1 - (КОВм НЗ (Мкй-)м (1 -- (К(й, Н+ (К(В),.В(К(В) , (5)40где 1 - единичная матрица;К - матрица интенсивности шумов измерений;0 - матрица интенсивности шумов системы.На фиг,1 изображена функциональная. схема предлагаемого устройства; на фиг,2 -схема блока расчета коэффициента коррекции; на фиг.3 - схема блока расчета апостериорной матрицы ковариации; на фиг,4 - схема минимаксного фильтра Калмана,Адаптивное устройство для обработкиданных содержит блок 1 элементов И, минимаксный фильтр 2 Калмана, блок 3 расчета коэффициента коррекции, блок 4 расчета апостериорной матрицы ковариации, второй блок 5 элементов задержки, седьмой умножитель б, квадратор 7. блок 8 памяти матрицы интенсивности шумов системы, первый 9 и второй 10 умножители, первый -третий сумматоры 11 - 13, блок 14 памятиматрицы коэффициентов системы, четвертый умножитель 15, блок 16 памяти измерительной матрицы и матрицы шумовизмерений, пятый умножитель 17, четвертый сумматор 18, шестой умножитель 19,блок 20 памяти границы неопределенностипереходной матрицы, блок 21 памяти предполагаемой переходной матрицы, третий22, восьмой 23, девятый 24 умножители,первый блок 25 элементов задержки, блок26 управления, блок 27 памяти переходнойматрицы системы.Блок 3 содержит первый - третий умножители 28-30, блок 31 регистров, сумматор32, делитель ЗЗ (фиг,2).Блок 4 содержит пероый - третий умножители 34-36, вычитагель 37, сумматор 38(фиг,З),В минимаксном фильтре 2 Калмана вкачестве усилителя используется матричный умножитель (цифровой аналог усилителя), Он содержит (фиг.4) вычитатель 39,первый 40, второй 41, третий 42 умножители, сумматор 43, блок 44 элементов задержки. Первая группа входов фильтра 2подк,. очена к умножителю 40, вторая - кумножителю 41, третья - к умножителю 42,четвертая - к оычитателю 39,Блок 26 управления реализован в видепоследовательно соединенных генераторатактовых импульсов и счетчика, Генераторосуществляет управление поступлением обрабатываемой информации в устройство,записью данных в блок 31, а счетчик обеспечивает управление адресными входамиблока 27.Блок 27,может быть выполнен на МКСЕ 155 РУ 2, Запись данных в блок 27 производится до начала работы устройства, а сигналом разрешения считывания данных можетслужить сигнал с выхода генератора, блока26; в прямом виде на вводе 3 указанноймикросхемы и в инвертированном - на вводе 2.Блок 16 памяти представляет собой объединение двух блоков памяти (регистров памяти), в каждом из которых хранитсяинформация, записанная перед началомцикла обработки,Второй вход блока 1 является информационным входом всего устройства, первыйвход блока 1 соединен с первым выходомблока 26, который соединен также с входомсчитывания блока 27, адресный вход которого соединен с вторым выходом блока 26,выход блока 27 подключен к третьему входуфильтра 2, четвертый вход которого соединен с выходом блока 1, а второй вход соединен с вторым выходом блока 16, которыйтакже соединяется с вторым входом блока 3и с третьим входом блока 4, управляющий вход блока 3 подключен к первому выходу блока 26, а первый вход подключен к первому выходу блока 16, который также соединяется с вторым входом блока 4, четвертый вход которого, как и третий вход блока 3,подключен к выходу сумматора 12, с которым также соединяется вход блока 5, выход10 которого подключен к второму входу умно-.жителя 6, первый вход его подключен к вы 15 20 25 30 35 40 45 50 55 ходу квадратора 7, а выход соединен через первый вход сумматора 13 с вторым входом сумматора 12, первый вход которого подключен к выходу сумматора 11; на входы которого подключены выходы умножителей 10 и 15, на первьй вход первого из них подключается через умножитель 9 выход блока 8, на другие входы умножителей 9 и 10 подключен выход блока 14, на первый вход умножителя 15 подключен выход умно- жителя 22, а на оставшиеся входы умножителей 15 и 22 подключен выход блока 21, который соединяется также через первый вход умножителя 17 с первым входом сумматора 18 и через первый вход умножителя 24 - с вторым входом умножителя 19, на первый вход которого подключен выход блока 20, который соединяется также с входом квадратора 7 и вторым входсм умножителя 23, на первый вход которого, как и на соответствующие входы умножителей 24 и 22, подключен выход блока 25, вход которого соединен с выходом блока 4, выходы умножителей 23 и 19 соединены с вторыми входами блоков 17 и 18 соответственно, а выход последнего подключен к второму входу сумматора 13, первые входы фильтра 2 и блока 4 соединены с выходом блока 3, а выход фильтра 2 является выходом всего устройства.Устройстоо работает следующим образом.В блоки 21, 14, 8, 16, 20 и 27 записываются пакеты сигналов, пропорциональные соответственно элементам предполагаемой матрицы системы Ф матрицы системы Г, матрицы интенсивностей шумов системы О, измерительная матрица Н и матрица интенсивностей шумов измерений Н, верхняя граница неопределенности нормы матрицы перехода Л, переходная матрица системы Ф,В устройстве организован синхронный принцип накопления, вычисления и хранения информации и использован конвейерный режим работц, При этом режим "накачки" входными сигналами определяется временем выполнения самого процессамирует пакет сигналов, пропорциональных 55 произведению матриц Г Ы, который поступает в умножитель 10, формирующий пакет сигналов, пропорциональных произведению матриц Г О Г, который поступает в вычисления. С такой же скоростью выдаются результаты вычислений. В приведеннойструктуре производится обработка сигналов в каждом блоке по мере их накопления во входных регистрах, Таким образом, различные элементы устройства осуществляютодновременно обработки сигналов по своим алгоритмам с различным временем окончания работы, и затем передают результаты оешений на следующие блоки, работа которых начинается после получения сигналов, пропорциональных предыдущим значениям вычисленных данных.В момент включения в блоках памяти записаны соответствующие матрицы, в регистрах блока 3 - "0", счетчик в блоке 26 обнулен. С появлением первого импульса на первом выходе блока 26 открываются элементы И блока 1, после чего входной сигнал 2 к поступает на фильтр 2, на третий вход которого поступают составляющие матрицы с блока 27 под управлением на его адресные входы с выходов (2 - и + 1) блока 26, На второй и первый входы фильтра 2 поступают сигналы с блоков 16 (матрица Н) и 3 - Км(кй), Пакет сигналов, пропорциональных элементам матрицы М(ка), сформированной в блоке 4, поступает в блок 25 элементов задержки. Причем блок 4 осуществляет преобразования в соответствии с уравнением (5). На выходе вычитателя 37 вырабатывается сигнал, пропорциональный выражению (1- К(1/1 Д, Операция транспортирования в умножителях 35 и 36осуществляется за счет соответствующей коммутации второй руппы входов умножителей относительно первой. На выходе сумматора 38 вырабатывается матрица (Мй)м, которая поступает в блок 25, Одновременнос выхода блока 25 в умножители 22 - 24 поступает кодовая комбинация, пропорциональная элементам матрицы (М 1-1 Й)м, соответствующая предыдущему такту работы. Элементы матрицы Фс блока 21 действуют на соответствующие входы умножителей 22 и 24. После умножения в умножителе 15 фоРмиРУетсЯ матРица Ф(Мка)л ч, сигнал которой поступает на сумматор 11, Пакет сигналов, пропорциональных элементам матрицы системы Г, с блока 14 поступает одновременно в умнокители 9 и 10. Сигналы, пропорциональные элементам матрицы интенсивностей шумов системы О, с блока 8 поступают в умножитель 9, который фор 10 15 20 2530354050 сумматор 11, Сумматор 11 формирует пакетсигналов, пропорциональных сумме произведений матриц Ф(М-1)мФ + ГО Г,который поступает в сумматор 12,Сигнал, пропорциональный величинеЛ, с блока 20 памяти одновременно поступает в умножители 23 и 19. Умножитель 23формирует пакет сигналов, пропорциональных элементам матрицы Л(ма)л, который поступает на умножитель 17. Пакетсигналов, пропорциональных элементамматрицы Ф, поступает в умножитель 17,формирующий пакет сигналов, пропорциональных элементам произведения матрицд(А(к)м Ф, который поступает на сумматор 18,Умножитель 24 формирует пакет сигналов, пропорциональных элементам пролзведениЯ матРиц Ф(ма)л, котоРыйпоступает в умножитель 19, формирующийпакет сигналов, пропорциональных элементам произведения матриц ЛФ(Мй)м,который поступает в сумматор 18. Сумматор18 формирует пакет сигналов, пропорциональных элементам суммы произведений матрицЛ (Мк/к-фФ+ Л Ф(Мк)л, котоРыипоступает на сумматор .13. С выхода блока 5пакет сигналов, пропорциональных элементам матрицы (Мкк)м, рассчитанныйна предыдущем такте, поступает в умножитель б. Сигнал, пропорциональный величине Л, с квадратора 7 поступает в2умножитель б, формирующий пакет сигналов, пропорциональных элементам матри 2цы Л (Мк/-2)м, который поступаетна сумматор 13, Сумматор 13 формирует пакет сигналов, пропорциональныхэлементам суммы произведений матриц+ ФЛ (М(к)м, который поступает в сумматор 12, Сумматор 12 формирует пакет сигналов, пропорциональных элементамминимаксной матрицы (Му)м, которыйпоступает в блок 4 расчета минимакснойапостериорной матрицы ковариации, блок 3расчета коэффициента коррекции и в блок 5для задержки на один такт, В блоке 3 осуществляется преобразование в соответствии свыражением (2). На выходе умножителя 29вырабатывается сигнал, пропорциональныйвыражению (Мкл)м Н, а на выходе делителя ЗЗ - сигнал, пропорциональный выражению в фигурных скобках уравнения (2).Таким образом, в умножителе 30 формируется сигнал, соответствующий Кл(М/1), который записывается в блок 31 подуправлением тактового сигнала с первоговыхода блока 26 управления. С блока 31 ве 1741098 12личина КчЕЬ) поступает на первый вход блока 2 для обеспечения обработки в соответствии с алгоритмом (1) и на первую группу входов блока 4,Предлагаемое устройство по сравнению с прототипом обладает более высокой точностью обработки входных сигналов, обусловленной дополнительным изменением на каждом шаге обработки матричного коэффициента коррекции фильтра Калмана за счет учета неточности задания переходной матрицы системы. При этом процессе формирования оценки сигнала придается больший вес априорным данным, в частности модели системы, неверное задание переходной матрицы которой приводит к расхождению фильтра Калмана.Можно показать, что ограничение снизу коэффициента коррекции К(В) величиной, например, 1/Е приводит к следующей зависимости погрешности оценкиЕ(ХМ =: (1 в 1 /3 ) С(1+ )/2,где С,- некоторые постоянные величины,Тогда при увеличении количества измерений погрешность оценки стремится к нулю и фильтр не расходится.Величина добавки к коэффициенту коррекции зависит от верхней границы неопределенности задания нормы переходной матрицы системы и выбирается исходя из минимизации следа матрицы ковариации ошибки оценки при максимально возмокной неопределенности задания нормы переходной матрицы системы,Формула изобретения Адаптивное устройство для обработки данных, содержащее последовательно соединенные блок памяти матрицы интенсивности шумов системы, первый и второй умножители, первый сумматор, блок памяти матрицы коэффициентов системы, подключенный выходом к вторым входам первого и второго умножителей, блок памяти предполагаемой переходной матрицы, соединенный выходом через третий умножитель с первым входом четвертого умнокителя и непосредственно с вторым его входом, выход четвертого умножителя подключен к второму входу первого сумматора, минимаксный Фильтр Калмана, соединенный . первым входом с выходом блока расчета коэффициента коррекции и первым входом блока расчета апостериорной матрицы ковариации, второй вход которого подключен к первому входу блока расчета коэффициента коррекции и первому выходу блока памяти измерительной матрицы и матрицы шумов измерений, третий вход - к второму выходу блока памяти измерительной матрицы и матрицы шумов измерений и вторым входом блока расчета коэффициентов кор рекции и минимаксного фильтра Калмана, авыход через первый блок элементов задержки - к второму входу третьего умножителя, о т л и ч а ю щ е е с я тем, что, с целью повышения точности устройства, введены 10 блок элементов И, квадратор, три сумматора, с пятого гю девятый умножители, блок памяти границы неопределенности нормы переходной матрицы, блок памяти переходной матрицы системы, второй блок элемен тов задеркки, первый и второй входывторого сумматора соединены соответственно с выходами первого и третьего сумматоров, а выход - с входом второго блока элементов задеркки, третьим входом блока 20 расчета коэффициента коррекции и четвертым входом блока расчета апостериорной матрицы ковариации, выход блока памяти предполагаемой переходной матрицы через пятый умножитель подключен к первому 25 входу четвертого сумматора, к второму входу которого через шестой умножитель подключен выход блока памяти границы неопределенносги нормы переходной матрицы, соединенный также через последова тельно подключенные квадратор и седьмойумножитель с первым входом третьего сумматора, подключенного вторым входом к выходу четвертого сумматора, выход второго блока элементов задеркки соединен с 35 вторым входом седьмого умножителя, выход блока памяти предполагаемой переходной матрицы подключен также к первому входу девятого умножителя, соединенного вторым входом с выходом первого блока 40 элементов задержки и первым входом восьмого умножителя, а выходом - с вторым входом шестого умножителя, второй вход пятого умнокителя подключен к выходу восьмого умножителя, соединенного вто рым входом с выходом блока памяти границы неопределенности нормы переходной матрицы, первый выход блока управления подключен к управляющему входу блока расчета коэффициента коррекции, первому 50 входу блока элементов И и входу разрешения считывания блока памяти переходной матрицы системы, соединенного адресным входом с вторым выходом блока управления, а выходом - с третьим входом минимак сного фильтра Калмана, подключенногочетвертым входом к выходу блока элементов И, второй вход блока элементов И и выход минимаксного фильтра Калмана являются соответственно информационным входом и выходом устройства.Составитель В,Бета новедактор М, Кобылянская Техред М.Моргентал Корректор Т.Ланская атент", г. Ужгоро арина, 1 ельский комбин роизводствен аз 2084ВКИИПИ Госуд Тираженного комитет13035, Москва,Подписноео изобретениям и открытиям при ГКНТ С5, Раушскэя наб., 4/5