RU2436228C1 - Цифровой интеллектуальный многокаскадный фильтр - Google Patents

Цифровой интеллектуальный многокаскадный фильтр Download PDF

Info

Publication number
RU2436228C1
RU2436228C1 RU2010115924/08A RU2010115924A RU2436228C1 RU 2436228 C1 RU2436228 C1 RU 2436228C1 RU 2010115924/08 A RU2010115924/08 A RU 2010115924/08A RU 2010115924 A RU2010115924 A RU 2010115924A RU 2436228 C1 RU2436228 C1 RU 2436228C1
Authority
RU
Russia
Prior art keywords
block
output
information
information input
unit
Prior art date
Application number
RU2010115924/08A
Other languages
English (en)
Inventor
Андрей Александрович Костоглотов (RU)
Андрей Александрович Костоглотов
Антон Александрович Кузнецов (RU)
Антон Александрович Кузнецов
Сергей Валерьевич Лазаренко (RU)
Сергей Валерьевич Лазаренко
Дмитрий Сергеевич Андрашитов (RU)
Дмитрий Сергеевич Андрашитов
Игорь Владимирович Дерябкин (RU)
Игорь Владимирович Дерябкин
Анатолий Викторович Чеботарев (RU)
Анатолий Викторович Чеботарев
Original Assignee
Андрей Александрович Костоглотов
Антон Александрович Кузнецов
Сергей Валерьевич Лазаренко
Дмитрий Сергеевич Андрашитов
Игорь Владимирович Дерябкин
Анатолий Викторович Чеботарев
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Андрей Александрович Костоглотов, Антон Александрович Кузнецов, Сергей Валерьевич Лазаренко, Дмитрий Сергеевич Андрашитов, Игорь Владимирович Дерябкин, Анатолий Викторович Чеботарев filed Critical Андрей Александрович Костоглотов
Priority to RU2010115924/08A priority Critical patent/RU2436228C1/ru
Application granted granted Critical
Publication of RU2436228C1 publication Critical patent/RU2436228C1/ru

Links

Images

Landscapes

  • Complex Calculations (AREA)

Abstract

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

Description

Изобретение относится к цифровой вычислительной технике и может быть использовано в системах цифровой обработки сигналов при решении задач оптимальной нелинейной фильтрации.
Известно устройство - расширенный калмановский фильтр [1, 2], недостатком которого является ограниченность функциональных возможностей, обусловленная линейной структурой обрабатываемых процессов, и устройство [3], а также эвристическое задание параметра регуляризации.
Наиболее близким по технической сущности к заявленному изобретению является цифровой интеллектуальный итеративный фильтр [3], содержащий первый, второй, третий блоки формирования разности; первый, второй, третий блоки коррекции; блок формирования и выдачи априорных данных; первый, второй, третий блоки формирования суммы; первый, второй, третий, четвертый, пятый и шестой блоки формирования матричных функций; первую, вторую и третью линии задержки; первый и второй блоки расчета параметра регуляризации. Однако данное устройство в некоторых практически важных случаях не позволяет получить требуемой точности.
Улучшение точностных характеристик фильтрации параметров случайных процессов является актуальным направлением.
Заявленное изобретение направлено на повышение точности при формировании оценки параметров информационного процесса в измерительных системах, что весьма важно при радиолокационном сопровождении целей, содержит блоки: первый, второй, третий, четвертый блоки формирования разности; первый, второй, третий, четвертый блоки коррекции; блок формирования и выдачи априорных данных; первый, второй, третий, четвертый блоки формирования суммы; первый, второй, третий, четвертый, пятый, шестой, седьмой и восьмой блоки формирования матричных функций; первую, вторую, третью и четвертую линии задержки; первый, второй и третий блоки расчета параметра регуляризации, при этом первый, второй, третий и четвертый выходы блока формирования и выдачи априорных данных соединены соответственно со вторым, третьим, четвертым и пятым информационными входами первого, второго, третьего и четвертого блоков коррекции, кроме того, первый, второй, пятый выходы блока формирования и выдачи априорных данных соединены соответственно с третьим, вторым и пятым информационными входами первого, второго и третьего блоков расчета параметра регуляризации, первый информационный выход первого блока коррекции соединен с первым информационным входом первого блока формирования суммы, выход которого соединен с информационным входом первого блока формирования матричной функции и четвертым информационным входом первого блока расчета параметра регуляризации, выход которого соединен с девятым информационным входом второго блока коррекции, выход первого блока формирования матричной функции соединен с информационным входом первой линии задержки, выход которой соединен со вторым информационным входом первого блока формирования суммы, с седьмым информационным входом первого блока коррекции и с информационным входом второго блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом первого блока коррекции, первый информационный выход второго блока формирования матричной функции соединен со вторым информационным входом блока формирования разности, выход которого соединен с первым информационным входом первого блока коррекции; второй информационный выход первого блока коррекции соединен с восьмым информационным входом второго блока коррекции, первый информационный выход которого соединен с первым информационным входом второго блока формирования суммы, выход которого соединен с информационным входом третьего блока формирования матричной функции и с четвертым информационным входом второго блока расчета параметра регуляризации, выход которого соединен с девятым информационным входом третьего блока коррекции, выход третьего блока формирования матричной функции соединен с информационным входом второй линии задержки, выход которой соединен со вторым информационным входом второго блока формирования суммы, седьмым информационным входом второго блока коррекции и с информационным входом четвертого блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом второго блока коррекции, первый информационный выход четвертого блока формирования матричной функции соединен со вторым информационным входом второго блока формирования разности, выход которого соединен с первым информационным входом второго блока коррекции; второй информационный выход второго блока коррекции соединен с восьмым информационным входом третьего блока коррекции, первый информационный выход которого соединен с первым информационным входом третьего блока формирования суммы, выход которого соединен с информационным входом пятого блока формирования матричной функции, выход которого соединен с информационным входом третьей линии задержки, выход которой соединен со вторым информационным входом третьего блока формирования суммы, седьмым информационным входом третьего блока коррекции и с информационным входом шестого блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом третьего блока коррекции, первый информационный выход шестого блока формирования матричной функции соединен со вторым информационным входом третьего блока формирования разности, выход которого соединен с первым информационным входом третьего блока коррекции; второй информационный выход третьего блока коррекции соединен с восьмым информационным входом четвертого блока коррекции, первый информационный выход которого соединен с первым информационным входом четвертого блока формирования суммы, выход которого является выходом устройства и соединен с информационным входом седьмого блока формирования матричной функции, выход которого соединен с информационным входом четвертой линии задержки, выход которой соединен со вторым информационным входом четвертого блока формирования суммы, седьмым информационным входом четвертого блока коррекции и с информационным входом восьмого блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом четвертого блока коррекции, первый информационный выход восьмого блока формирования матричной функции соединен со вторым информационным входом четвертого блока формирования разности, выход которого соединен с первым информационным входом четвертого блока коррекции.
Одним из способов повышения точности оценки параметров динамических систем является использование методов решения некорректных задач на основе принципов регуляризации. Эффективность применения регуляризации для случая непрерывных систем доказана и получила свое выражение в виде метода А.Н.Тихонова [4] и его разновидности в виде метода итерационной регуляризации [5]. Получим уравнения фильтрации с использованием метода итерационной регуляризации для дискретной системы [6, 9].
Пусть динамика оцениваемых параметров описывается системой разностных уравнений в дискретном времени
Figure 00000001
где x(k)=[x1(k), x2(k),…, xM(k)]T∈EM - вектор состояния исследуемой системы; η(k)=[η1(k), η2(k),…, ηM(k)]T∈EM - вектор неизвестных внешних воздействий; переходная функция F(x(k), k)=[F1(x(k), k), F2(x(k), k),…, FM(x(k), k)]T∈EM - непрерывная вместе с частными производными вектор-функция своих аргументов; G∈EM×EM - матрица интенсивности внешних воздействий; k, N, M - натуральные числа. Предполагается, что матрица
Figure 00000002
имеет обратную матрицу.
Наблюдаемый сигнал, получаемый на выходе модели измерительной системы, описывается дискретным уравнением
Figure 00000003
где y(k)=[y1(k), y2(k),…, yL(k)]T∈EL - вектор наблюдения, n(k)=[n1(k), n2(k),…, nL(k)]T∈EL - вектор дискретного белого гауссовского шума с известными локальными характеристиками
M[n(k)]=0,
M[n(k)nT(l)]=Wδ(k-l),
W - ковариационная матрица размерности L×L, δ(·) - векторная дельта-функция; H(x(k))=[H1(x(k)), H2(x(k)),…, HL(x(k))]T∈EL - сигнальная вектор-функция, непрерывна вместе с частными производными; L, l - натуральные числа.
Поставим задачу синтеза рекуррентного фильтра оценки х*(k), оптимального в смысле минимума функционала, характеризующего ошибку измерения
Figure 00000004
В силу непрерывности вектор-функции F(·), решение уравнения (1) непрерывно зависит от η(k), поэтому функционал ошибки (3) на каждом решении системы (1) непрерывно зависит от η(k). Таким образом, задача определения оценки х*(k), доставляющей минимум (3), равносильна задаче определения
Figure 00000005
Задача (1), (2), (4) является некорректно поставленной обратной задачей [7]. Найти значения векторов х*(k), η*(k),
Figure 00000006
путем решения совокупности уравнений (1), (2), (4) в условиях некорректности исходной задачи достаточно сложно, в связи с этим широкое распространение получили итерационные градиентные методы. Однако использование таких методов может привести к расходящейся последовательности приближений. Поэтому применение любого итерационного метода для решения задачи (1), (2), (4) требует определения регуляризирующего семейства операторов, в котором параметром регуляризации является номер итерации.
В соответствии с общим определением регуляризирующего семейства операторов по А.Н.Тихонову [7] будем говорить, что итерационный метод
Figure 00000007
в котором числовой параметр αn удовлетворяет условиям:
Figure 00000008
где Δ(ηn) - невязка, порождает регуляризирующее семейство операторов, в котором параметром является номер итерации, если для любого начального приближения η0 и для любого значения погрешности исходных данных σ, удовлетворяющего условию 0<|σ|<σ0, σ0=const, существует номер n(σ) такой, что
Figure 00000009
,
то есть полученные приближения сходятся к точному решению в норме пространства
Figure 00000010
при стремлении погрешности исходных данных к нулю.
Согласно [8] выражение для определения градиента в точке ηn(k) имеет вид
Figure 00000011
где xn(k) - решение задачи (1) при ηn(k), а вектор Ψn(k) определяется из условий
Figure 00000012
Зная выражение для градиента (7) функционала (3), можно переписать выражение (6) для параметра регуляризации [9] в следующем виде
α0=1,
Figure 00000013
.
Выбор последовательности параметров αn, удовлетворяющей условию (9), позволит реализовать фильтр цифровой обработки измерительной информации повышенной точности.
Для реализации итерационного метода (5) требуется определить градиент функционала (3), определяемый выражением (7). Приняв за нулевое приближение
Figure 00000014
, запишем итерационную последовательность (5) в развернутой форме для
Figure 00000015
,
Figure 00000016
В результате, с учетом (7), имеем последовательность дискретных двухточечных краевых задач (ДДТКЗ) вида
Figure 00000017
Figure 00000018
Figure 00000019
,
Figure 00000020
,
Figure 00000021
,
Figure 00000022
.
Введем обозначение
Figure 00000023
и помножим каждое из уравнений для сопряженных векторов Ψi на величину αi, тогда уравнения (10) принимают следующий вид
Figure 00000024
Для получения алгоритма рекуррентного оценивания вектора состояния необходимо воспользоваться методом инвариантного погружения в дискретном варианте. Заметим, что уравнение для вектор-функций λ в ДДТКЗ (9) записано в обратном времени. Это требует его преобразования к виду, отражающему зависимость λn(k+1) от λn(k) и xn(k). Производя соответствующие преобразования с учетом выражения для xn(k+1) из (11) и используя разложение в ряд Тейлора в окрестности F(xn(k), k), получим следующую последовательность ДДТКЗ
Figure 00000025
xi(k+1/k)=F(xi(k),k),
xi(0)=x*(0), λn(N)=0,
Figure 00000021
,
Figure 00000022
,
где функции β и γ вводятся для сокращения записи.
Заменим условие на конце λn(N)=0 более общим условием λn(N)=с и пусть N и с - переменные величины. Тогда значение вектора xn(N) определяется как функция величин N и с
xn(N)=r[с, N].
Изменение величины N на N+1 дает приращение Δс, тогда
Figure 00000026
Запишем выражение для r(c+Δc, N+1), используя аппарат конечных разностей
r(c+Δc, N+1)=r(c, N)+Δrc(c, N)+ΔrN(c, N)+Δ2rcN(c, N),
или, учитывая (13), получим
Figure 00000027
где
Δrc(c, N)=r(c+Δc, N)-r(c, N),
ΔrN(c, N)=r(c, N+1)-r(c, N),
Δ2rcN(c, N)=Δrc(c, N+1)-Δrc(c, N).
Согласно (12) выражения для Δxn и Δс имеют вид
Figure 00000028
Разрешить разностное уравнение (14) относительно r(с, N), т.е. найти общее аналитическое решение, не удается, и обычно обращаются к приближенным методам. Предположим, что r(c, N) линейна по с
Figure 00000029
где
Figure 00000030
- оценка вектора состояния в момент N, Pn(N) - некоторая матрица размерности М×М.
Вычислим разности, входящие в выражение (14), используя выражение (16)
Figure 00000031
Подставив выражения (15), (16), (17) в (14), получим
Figure 00000032
Разлагая β и γ в ряд Тейлора в окрестности
Figure 00000033
и пренебрегая членами порядка выше первого, можно записать уравнение (18) в виде
Figure 00000034
Соотношение (19) выполняется при с→0, поэтому, приравнивая коэффициенты при первой и нулевой степени с, получим разностные уравнения для
Figure 00000035
и Pn(N+1)
Figure 00000036
Запишем ДДТКЗ (12) для случая, когда k=N, при этом учтем, что это все ДДТКЗ для i=0,…, n-1 являются разрешенными и соответственно оценки xi представляют собой известные функции параметра k. Имеем
Figure 00000037
Тогда уравнение (20) преобразуется следующим образом
Figure 00000038
где
Figure 00000039
.
Введем обозначения
Figure 00000040
,
Figure 00000041
.
Тогда уравнение (22) запишем в виде
Figure 00000042
Поскольку разницы в переобозначении матриц Pn и
Figure 00000043
нет, запишем последовательность уравнений для оценки
Figure 00000044
процесса (1), предполагая, что N постоянно изменяется и k=N, а также учитывая условие (7), налагаемое на параметр регуляризации, в виде
Figure 00000045
Последовательность уравнений (24) представляет собой цифровой интеллектуальный многокаскадный фильтр, который позволяет осуществить процесс цифровой обработки измерительной информации для дискретных динамических систем. Если сравнить полученные уравнения с уравнениями цифрового рекурсивного фильтра [3], то становится ясно, что они отличны друг от друга за счет введения дополнительной итерации, а также дополнительных связей последовательности параметров αi, изменяющих общий коэффициент обратной связи в уравнении для оценки
Figure 00000046
, с оценками
Figure 00000047
, сигнальной вектор-функцией Н и матрицами G, W. То есть в фильтре (24), в отличие от фильтра [3], параметр регуляризации согласован в соответствии с принципом невязки с ошибками измерений, что позволяет получить более точную процедуру оценивания. Алгоритм (24) дает оптимальную для нелинейной системы оценку параметров в смысле минимума функционала, характеризующего среднеквадратическую ошибку измерительного канала.
Вычисление параметра регуляризации организовано следующим образом: пределы интегрирования [0, t] изменены на [t-4s, t], где t означает текущее время, s - шаг вычислений; для дискретного времени необходимо брать [N-1, N-5]. Исследования показали [6], что использование данных больше чем на четыре шага назад обеспечивает повышение точности получаемых оценок на доли процента при увеличении количества необходимых арифметических операций в несколько раз. Метод, использующий более четырех итераций, будем называть многокаскадным.
Следует отметить, что для реализации 2-й итерации разработанного алгоритма (24) необходимо, чтобы 1-я итерация была реализована, для реализации 3-й - 1-я и 2-я, для реализации 4-й - 1-я, 2-я и 3-я. То есть количество арифметических операций, необходимых для вычисления итерации алгоритма, включат в себя число арифметических операций предыдущих итераций. Анализ вычислительных затрат, требуемых на реализацию разработанного алгоритма оценки, позволяет сделать вывод о возможности его реализации в реальном масштабе времени на основе использования современной вычислительной техники.
Исследования показали, что эффективность каждой повторной обработки по отношению к предыдущей падает. Между тем объем вычислений резко увеличивается. На основании анализа результатов численного моделирования можно рекомендовать четырехкратную итеративную обработку данных.
Оценка эффективности функционирования разработанного фильтра произведена на основе численного моделирования задачи определения неизвестного постоянного параметра d дискретной нелинейной системы третьего порядка
Figure 00000048
где параметр τ имеет смысл интервала времени, через который поступает измерительная информация в виде
Figure 00000049
Таким образом, как следует из соотношений (24), введение новых структурных элементов и связей позволяет в совокупности с общими признаками получить технический результат, состоящий в уменьшении дисперсии ошибок полученных на выходе фильтра оценок входных процессов.
Заявленное устройство может быть применено в информационных системах, связанных со сбором и обработкой информации, например в информационных системах радиолокационных и радионавигационных комплексов.
Сущность изобретения поясняется фиг.1-7, где представлены структурные схемы интеллектуального многокаскадного цифрового фильтра, первого и второго блоков коррекции, блока расчета точностных характеристик, блока расчета параметра регуляризации, первого блока формирования произведения числителя параметра регуляризации, блока расчета прогноза.
На фиг.1 представлена структурная схема цифрового интеллектуального многокаскадного фильтра. Устройство содержит первый блок 1, второй блок 7, третий блок 14 и четвертый блок 20 формирования разности; первый блок 2, второй блок 8, третий блок 15 и четвертый блок 21 коррекции; первый блок 3, второй блок 9, третий блок 16 и четвертый блок 22 формирования суммы; первую линию задержки 5, вторую линию задержки 11, третью линию задержки 18 и четвертую линию задержки 24; первый блок 4, второй блок 6, третий блок 10, четвертый блок 12, пятый блок 17, шестой блок 19, седьмой блок 23 и восьмой блок 25 формирования матричных функций; блок 13 формирования и выдачи априорных данных; первый блок 26, второй блок 27 и третий блок 28 расчета параметра регуляризации.
На фиг.2 представлена структурная схема первого блока коррекции, который содержит блок 2.1 формирования частных производных, блок 2.2 транспонирования матричной функции, блок 2.3 формирования произведения, блок 2.4 вычисления точностных характеристик, блок 2.5 формирования произведения.
На фиг.3 представлена структурная схема второго блока коррекции, который содержит блок 8.1 формирования частных производных, блок 8.2 транспонирования матричной функции, блок 8.3 формирования произведения, блок 8.4 формирования суммы, блок 8.5 расчета точностных характеристик, блок 8.6 формирования произведения.
На фиг.4 представлена структурная схема блока расчета точностных характеристик, входящего в первый, второй, третий и четвертый блоки коррекции, который содержит блок 14 формирования частных производных матричной функции, блок 21 транспонирования матриц, блок 22 формирования произведения, линию 23 задержки, блок 24 транспонирования матриц, блок 25 формирования произведения, блок 26 формирования суммы, блок 27 формирования частных производных матричной функции, блок 28 формирования произведения, блок 29 формирования разности, устройство 30 обращения матриц, блок 31 формирования произведения.
На фиг.5 представлена структурная схема блока расчета параметра регуляризации, который содержит блок формирования матричной функции 26.1, линии задержки 26.2, 26.4, 26.7, 26.8, 26.12, 26.13, 26.16, 26.17, 26.18, 26.24, 26.25, 26.26, 26.29, 26.30, 26.31, 26.32, 26.40, 26.41, 26.42, 26.43, блоки 26.3, 26.9, 26.11, 26.19, 26.22, 26.23, 26.33, 26.37, 26.38, 26.39 формирования произведения числителя параметра регуляризации, блоки 26.5, 26.14, 26.27, 26.44 транспонирования, блоки 26.6, 26.15, 26.28, 26.45 формирования произведения, блоки 26.10, 26.20, 26.21, 26.34, 26.35, 26.36 расчета прогноза, блоки 26.46, 26.47 формирования суммы, блок 26.48 формирования отношения.
На фиг.6 представлена структурная схема первого блока формирования произведения числителя параметра регуляризации, входящего в первый, второй и третий блок формирования параметра регуляризации, который содержит блок 26.3.1 формирования матричной функции, блок 26.3.2 частных производных матричной функции, блок 26.3.3 транспонирования матриц, блок 26.3.4, блок 26.3.5 формирования произведения, блоки 26.3.6 и 26.3.7 транспонирования матриц и блок 26.3.8 формирования произведения.
На фиг.7 представлена структурная схема первого блока расчета прогноза, который входит в состав блока формирования произведения числителя параметра регуляризации, он включает блоки 26.10.1 и 26.10.5 формирования произведения, блоки 26.10.2 и 26.10.6 формирования суммы, блок 26.10.3 формирования частных производных матричной функции, блок 26.10.4 формирования матричной функции.
Информационные входы линий задержки 26.2, 26.7, 26.16, 26.29 (фиг.5) соединены с информационным выходом блока 3 формирования суммы, первые информационные входы линий задержки 26.4, 26.12, 26.24, 26.40 соединены с первым информационным выходом блока 26.1 формирования матричной функции
Figure 00000050
, информационный вход которого соединен с информационным выходом блока 3, информационный выход первой линии задержки 26.2 соединен со вторым информационным входом первого блока 26.3 формирования произведения числителя параметра регуляризации, первый и третий его информационные входы соединены со вторым и первым информационными выходами блока 13 формирования и выдачи априорных данных соответственно, а четвертый информационный вход со входом устройства, первый информационный выход первого блока 26.3 формирования произведения числителя параметра регуляризации соединен с шестым и седьмым информационными входами второго блока 26.46 формирования суммы, первый информационный выход которого соединен с первым информационным входом блока 26.48 формирования отношения, первый выход блока 26.48 формирования отношения является выходом блока 14; первый информационный выход второй линии задержки 26.4 соединен с первым информационным входом первого блока транспонирования 26.5 и первым информационным входом первого блока 26.6 формирования произведения, первый информационный выход блока 26.5 соединен со вторым информационным входом блока 26.6 формирования произведения, первый информационный выход которого соединен с четвертым информационным входом первого блока 26.47 формирования произведения, первый информационный выход которого соединен со вторым информационным входом блока 26.48 формирования отношения; первый информационный выход третей 26.7 линии задержки соединен с первым информационным входом четвертой 26.8 линии задержки, первый информационный выход которой соединен со вторым информационным входом второго блока 26.9 формирования произведения числителя параметра регуляризации и вторым информационным входом первого блока 26.10 расчета прогноза, второй информационный выход которого соединен со вторым информационным входом третьего блока 26.11 формирования произведения числителя параметра регуляризации, третьи и четвертые информационные входы блоков 26.9, 26.11 соединены с первым информационным выходом блока 13 и входом устройства, первый информационный выход блока 26.9 соединен с третьим и четвертым информационными входами блока 26.46, первый информационный выход блока 26.11 соединен с десятым информационным входом блока 26.46; первый информационный выход пятой 26.12 линии задержки соединен с первым информационным входом шестой 26.13 линии задержки, первый информационный выход которой соединен с первым информационным входом второго блока 26.14 транспонирования и первым информационным входом второго блока 26.15 формирования произведения, первый информационный выход второго 26.14 блока транспонирования соединен со вторым информационным входом блока 26.15, первый информационный выход которого соединен с третьим информационным входом блока 26.47 формирования суммы; первый информационный выход седьмой линии задержки 26.16 соединен с первым информационным входом восьмой линии задержки 26.17, первый информационный выход которой соединен с первым информационным входом девятой 26.18 линии задержки, первый информационный выход которой соединен со вторыми информационными входами четвертого 26.19 блока формирования произведения числителя параметра регуляризации и второго блока 26.20 расчета прогноза, второй информационный выход которого соединен со вторыми информационными входами третьего блока 26.21 расчета прогноза и 26.23 формирования произведения числителя параметра регуляризации, второй информационный вход блока 26.22 расчета прогноза соединен со вторым информационным выходом третьего блока 26.21 вычисления, третьи и четвертые информационные входы блоков 26.19, 26.22 и 26.23 соединены с первым информационным выходом блока 13 и входом устройства соответственно, первый информационный выход блока 26.19 соединен с одиннадцатым и двенадцатым информационными входами второго блока 26.46 формирования суммы, первые информационные выходы блоков 26.23, 26.22 соединены с тринадцатым и четырнадцатым информационными входами первого блока 26.46 соответственно; первый информационный выход десятой линии задержки 26.24 соединен с первым информационным входом одиннадцатой линии задержки 26.25, первый информационный выход которой соединен с первым информационным входом двенадцатой 26.26 линии задержки, первый информационный выход которой соединен с первым информационным входом третьего блока 26.27 транспонирования матричной функции и первым информационным входом третьего блока 26.28 формирования произведения, первый информационный выход третьего блока 26.27 транспонирования соединен со вторым информационным входом третьего блока 26.28 формирования произведения, первый информационный выход которого соединен с первым информационным входом первого блока 26.47 формирования суммы; первый информационный выход тринадцатой линии задержки 26.29 соединен с первым информационным входом четырнадцатой линии задержки 26.30, первый информационный выход которой соединен с первым информационным входом пятнадцатой 26.31 линии задержки, первый информационный выход которой в свою очередь соединен с первым информационным входом шестнадцатой линии задержки 26.32, выход которой соединен со вторым информационным входом седьмого блока 26.33 формирования произведения числителя параметра регуляризации и вторым информационным входом четвертого блока 26.34 расчета прогноза, второй информационный выход которого соединен со вторым информационным входом девятого блока 26.38 формирования произведения числителя параметра регуляризации и вторым информационным входом пятого блока 26.35 расчета прогноза, второй информационный выход последнего соединен со вторым информационным входом десятого блока 26.39 формирования произведения числителя параметра регуляризации и со вторым информационным входом шестого блока расчета прогноза 26.36, второй информационный выход которого соединен со вторым информационным входом восьмого блока 26.37 формирования произведения числителя параметра регуляризации, третьи и четвертые информационные входы блоков 26.33, 26.37, 26.38, 26.39 соединены соответственно с первым информационным выходом блока 13 и входом устройства, первый информационный выход блока 26.33 соединен с восьмым и девятым информационными входами блока 26.46 формирования суммы, первый информационный выход блока 26.37 формирования произведения числителя параметра регуляризации соединен с пятым информационным входом блока 26.46, первые информационные выходы блоков 26.38 и 26.39 формирования произведения числителя параметра регуляризации соединены со вторым и первым информационными входами блока 26.46 соответственно; первый информационный выход семнадцатой линии задержки 26.40 соединен с первым информационным входом восемнадцатой линии задержки 26.41, первый информационный выход которой соединен с первым информационным входом девятнадцатой 26.42 линии задержки, первый информационный выход которой в свою очередь соединен с первым информационным входом двадцатой линии задержки 26.43, информационный выход которой соединен с первыми информационными входами четвертого блока 26.44 транспонирования и четвертого блока 26.45 формирования произведения, первый информационный выход блока 26.44 соединен со вторым информационным входом блока 26.45, первый информационный выход которого соединен со вторым информационным входом блока 26.47; первые информационные входы второго 26.9, четвертого 26.19 и седьмого 26.33 блоков формирования произведения числителя параметра регуляризации соединены со входом устройства, а первые информационные входы третьего 26.11, пятого 26.22, шестого 26.23, восьмого 26.37, девятого 26.38 и десятого 26.39 блоков формирования произведения числителя параметра регуляризации соединены соответственно с первыми информационными выходами первого 26.10, третьего 26.21, второго 26.20, шестого 26.36, четвертого 26.34 и пятого 26.35 блоков расчета прогноза; первые информационные входы первого 26.10, второго 26.20 и четвертого 26.34 блоков расчета прогноза соединены со вторым информационным выходом блока 13, а первые информационные входы третьего 26.21, пятого 26.35, шестого 26.36 блоков расчета прогноза соединены с первыми информационными выходами соответственно второго 26.20, четвертого 26.34 и пятого 26.35 блоков расчета прогноза, третьи информационные входы блоков расчета прогноза соединены с пятым информационным выходом блока 13. Структура второго блока расчета параметра регуляризации аналогична.
Первый, второй и пятый информационные выходы блока 13 формирования и выдачи априорных данных (фиг.1) соединены соответственно с третьим, вторым и пятым информационными входами первого блока 26, второго блока 27 и третьего блока 28 расчета параметра регуляризации, выход которого соединен с девятым информационным входом четвертого блока 21 коррекции, а первый, второй, третий, четвертый информационные выходы блока 13 формирования и выдачи априорных данных соединены соответственно со вторым, третьим, четвертым, пятым информационными входами первого блока 2, второго блока 8, третьего блока 15 и четвертого блока 21 коррекции, второй информационный выход первого блока 2 коррекции соединен с восьмым информационным входом второго блока 8 коррекции, второй информационный выход второго блока 8 коррекции соединен с восьмым информационным входом третьего блока 15 коррекции, второй информационный выход третьего блока 15 коррекции соединен с восьмым информационным входом четвертого блока 21 коррекции, выход которого соединен с первым информационным входом четвертого блока 22 формирования суммы, выход которого является выходом устройства и соединен с информационным входом седьмого блока 23 формирования матричной функции, выход которого соединен с информационным входом четвертой линии задержки 24, выход которой соединен со вторым информационным входом четвертого блока 22 формирования суммы, с седьмым информационным входом четвертого блока 21 коррекции и с информационным входом восьмого блока 25 формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом четвертого блока 21 коррекции; первый информационный выход восьмого блока 25 формирования матричной функции соединен со вторым информационным входом четвертого блока 20 формирования разности, выход которого соединен с первым информационным входом четвертого блока 21 коррекции; выход первого блока 26 расчета параметра регуляризации соединен с девятым информационным входом второго блока 8 коррекции, выход которого соединен с первым информационным входом второго блока 9 формирования суммы, выход которого соединен с четвертым информационным входом второго блока 27 расчета параметра регуляризации и с информационным входом третьего блока 10 формирования матричной функции, выход которого соединен с информационным входом второй линии 11 задержки, выход которой соединен со вторым информационным входом второго блока 9 формирования суммы, с седьмым информационным входом второго блока 8 коррекции и с информационным входом четвертого блока 12 формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом блока 8 коррекции; первый информационный выход четвертого блока 12 формирования матричной функции соединен со вторым информационным входом второго блока 7 формирования разности, выход которого соединен с первым информационным входом второго блока 8 коррекции; первый информационный выход первого блока 2 коррекции соединен с первым информационным входом первого блока 3 формирования суммы, выход которого соединен с четвертым информационным входом первого блока 14 расчета параметра регуляризации и с информационным входом первого блока 4 формирования матричной функции, выход которого соединен с информационным входом первой линии 5 задержки, выход которой соединен со вторым информационным входом первого блока 3 формирования суммы, с седьмым информационным входом первого блока 2 коррекции и информационным входом второго блока 6 формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом первого блока 2 коррекции; первый информационный выход второго блока 6 формирования матричной функции соединен со вторым информационным входом блока 1 формирования разности, выход которого соединен с первым информационным входом первого блока 2 коррекции; первый информационный вход первого блока 1 формирования разности, первый информационный вход второго блока 7 формирования разности, первый информационный вход третьего блока 26 формирования разности и первый информационный вход четвертого блока 20 формирования разности, а также первый информационный вход первого блока 26 расчета параметра регуляризации, первый информационный вход второго блока 27 расчета параметра регуляризации и первый информационный вход третьего блока 28 расчета параметра регуляризации являются входами устройства.
Первый и четвертый информационные выходы блока 13 формирования и выдачи априорных данных соединены с третьим и четвертым информационными входами блока 2.3 формирования произведения (фиг.2); информационный выход первого блока 1 формирования разности соединен с первым информационным входом блока 2.3 формирования произведения; второй информационный выход второго блока 6 формирования матричной функции соединен с информационным входом блока 2.1 формирования частных производных, выход которого соединен с информационным входом блока 2.2 транспонирования матричной функции, выход которого соединен со вторым информационным входом блока 2.3, выход которого соединен с первым информационным входом блока 2.4 расчета точностных характеристик, выход которого соединен с первым информационным входом блока 2.5 формирования произведения, выход которого является выходом первого блока 2 коррекции; информационный выход линии задержки 5 (фиг.1) соединен со вторым информационным входом блока 2.4 расчета точностных характеристик; второй и третий выходы блока 13 формирования и выдачи априорных данных соединены с третьим и четвертым информационными входами блока 2.4 расчета точностных характеристик; выход блока 2.3 формирования произведения соединен со вторым информационным входом блока 2.5 формирования произведения и вторым информационным входом блока 8.4 формирования суммы (фиг.3).
Информационный выход второго блока 7 формирования разности соединен с первым информационным входом блока 8.3 формирования произведения (фиг.3). Первый информационный выход блока 13 формирования выдачи априорных данных соединен с третьим информационным входом блока 8.3 формирования произведения; выход блока 26.48 соединен со вторым информационным входом блока 8.3 формирования произведения; второй информационный выход четвертого блока 12 формирования матричной функции соединен с первым информационным входом блока 8.1 формирования частных производных, выход которого соединен с информационным входом блока 8.2 транспонирования матричной функции, выход которого соединен с четвертым информационным входом блока 8.3, выход которого соединен с первым информационным входом блока 8.4 формирования суммы, второй информационный вход которого соединен с информационным выходом блока 2.3 (фиг.2), выход блока 8.4 формирования суммы соединен с первым информационным выходом блока 8.5, а также со вторым информационным выходом блока 8.6 формирования произведения; третий и четвертый информационные выходы блока 13 формирования выдачи априорных данных соединены со вторым и третьим информационными входами блока 8.5 расчета точностных характеристик; информационный выход линии задержки 11 соединен с четвертым информационным входом блока 8.5 расчета точностных характеристик, выход которого соединен с первым информационным входом блока 8.6 формирования произведения, выход которого является выходом второго блока 8 коррекции (фиг.1).
Выход блока 2.3 формирования произведения (фиг.2) соединен с информационным входом блока 27 формирования частных производных (фиг.4), выход которого соединен с информационным входом блока 28 формирования произведения, выход которого соединен с первым информационным входом блока 29 формирования разности, выход которого соединен с информационным входом устройства 30 обращения матриц, выход которого соединен с первым информационным входом блока 31 формирования произведения, выход которого является выходом блока 2.4 (фиг.2) расчета точностных характеристик; четвертый выход блока 13 формирования и выдачи априорных данных соединяется с вторым информационным входом блока 29 формирования разности; выход блока 31 формирования произведения соединен с информационным входом линии 23 задержки, выход которой соединен с третьим информационным входом блока 22 формирования произведения, выход которого соединен со вторым информационным входом блока 26 формирования суммы, выход которого соединен со вторыми информационными входами блоков 28 и 31 формирования произведения; выход первой линии задержки 5 (фиг.1) соединен с информационным входом блока 14 формирования частных производных, выход которого соединен с информационным входом блока 21 транспонирования матриц и первым информационным входом блока 22 формирования произведения, второй информационный вход которого соединен с выходом блока 21; третий информационный выход блока 13 формирования и выдачи априорных данных соединен с информационным входом блока 24 транспонирования матриц и вторым информационным входом блока 25 формирования произведения, первый информационный вход которого соединен с информационным выходом блока 24; выход блока 25 соединен с информационным входом блока 26 формирования суммы.
Блок 8.5 расчета точностных характеристик (фиг.3) имеет структуру и связи, аналогичные блоку 2.4.
Первый и второй информационные выходы блока 13 формирования и выдачи априорных данных соединены с четвертым информационным входом блока 26.3.5 (фиг.6) и информационным входом блока 26.3.6 транспонирования матричной функции, информационный выход которого соединен с третьим информационным входом блока 26.3.5, информационный выход которого соединен, с информационным входом блока 26.3.7 транспонирования матричной функции и вторым информационным входом блока 26.3.8 формирования произведения, информационный выход которого является выходом первого блока формирования произведения числителя параметра регуляризации, информационный выход блока 26.3.7 соединен с первым информационным входом блока 26.3.8; первый информационный вход блока 26.3.5 соединен с информационным входом блока 26.3.4, второй информационный вход блока 26.3.5 соединен с информационным выходом блока 26.3.3 транспонирования матричной функции, информационный вход которого соединен с информационным выходом блока 26.3.2 формирования частных производных матричной функции, информационный выход которого в свою очередь соединен с информационным входом блока формирования матричной функции; структура всех блоков формирования произведения числителя параметра регуляризации идентична, они различаются только связями между информационным входом блока 26.3.1 (соединен либо с информационным выходом соответствующей линии задержки, либо первым информационным выходом блока расчета прогноза) и блока 26.3.6 (информационный вход которого соединен либо со вторым информационным выходом блока 13 формирования и выдачи априорных данных, либо вторым информационным выходом блока расчета прогноза).
Первый информационный вход блока формирования 26.10.1 произведения, который является первым информационным входом первого блока расчета прогноза (фиг.7), соединен с информационным выходом линии задержки, второй информационный вход блока 26.10.1 соединен с пятым информационным выходом блока 13 формирования и выдачи априорных данных, с этим же информационным выходом соединен первый информационный вход блока 26.10.2 формирования суммы и второй информационный выход блока 26.10.5, третий информационный вход блока 26.10.1 соединен с информационным выходом блока 26.10.3 формирования частных производных матричной функции, второй информационный вход которой соединен с пятым информационным входом блока 13 формирования и выдачи априорных данных, а первый информационный вход с информационным первым выходом блока 26.10.4 формирования матричной функции, второй информационный выход которой соединен с первым информационным входом блока 26.10.5 формирования произведения и вторым информационным входом блока 26.10.6 формирования суммы, информационный выход блока 26.10.5 соединен с первым информационным входом блока 26.10.6; информационный выход блока 26.10.1 соединен со вторым информационным входом блока 26.10.2, его информационный выход является вторым информационным выходом блока расчета прогноза и соединен с третьим информационным входом блока 26.3.6 второго блока формирования произведения числителя параметра регуляризации (фиг.1); структура последующих блоков расчета прогноза аналогична, различаются только связями между третьим информационным входом блока 26.3.5 (он соединен либо со вторым информационным входом предыдущего блока прогноза, либо с пятым информационным выходом блока 13 формирования и выдачи априорных данных) и первым информационным входом блока 26.3.1 (он соединен или с информационным выходом линии задержки, или первым информационным выходом предыдущего блока прогноза).
Устройство работает следующим образом (фиг.1). В исходном состоянии в блоке 13 формирования и выдачи априорных данных записаны значения матриц W-1, G, I, а также значения α0, Δt - значение шага дискретизации. Значение оценки информационного процесса в (k+1)-й момент времени
Figure 00000051
с выхода блока 3 формирования суммы поступает на вход блока 26 вычисления параметра регуляризации и вход блока 4 формирования матричной функции, на остальные входы блока 26 поступают значения y, G, W-1, с выхода блока 4 формирования матричной функции значение
Figure 00000052
поступает на вход линии задержки 5, с выхода которой значение
Figure 00000053
поступает на вход блока 3 формирования суммы, вход блока 2 коррекции и вход блока 6 формирования матричной функции
Figure 00000054
, значение которой с выхода блока 6 поступает на вход блока 2 и на вход блока 1, на выходе блока 1 формируется значение
Figure 00000055
невязки измерения, которое поступает на вход блока 2 коррекции, на остальные входы которого поступают значения α0, G, W-1, I; в блоке 2 формируется произведение матричного коэффициента усиления и невязки измерения, которое суммируется в блоке 3 со значением
Figure 00000056
; на выходе блока 26 начиная с четвертого момента времени формируется значение α1 параметра регуляризации (до этого имеем α1=0.5), которое поступает на вход второго блока 8 коррекции, на остальные входы блока 26 поступают значения α0, G, W-1, I; с одного из выходов блока 2 коррекции значение М0(k+1/k) поступает на вход блока 8 коррекции, в котором формируется значение
Figure 00000057
которое поступает на вход блока 9 формирования суммы; значение
Figure 00000058
на выходе блока 9 поступает на вход блока 27 вычисления параметра регуляризации и вход блока 10 формирования матричной функции, с выхода блока 10 значение
Figure 00000059
поступает на вход линии задержки 11, на выходе которой формируется значение
Figure 00000060
, которое суммируется со значением (27) в блоке 9, на выходе которого формируется значение
Figure 00000061
; с выхода блока 11 значение
Figure 00000060
поступает на вход блока 8 коррекции, на вход блока 9 и на вход блока 12 формирования матричной функции, на выходе которого формируется значение
Figure 00000062
которое подается на вход блока 8 и вход блока 7, на другой вход которого поступает входное колебание; невязка измерения
Figure 00000063
с выхода блока 7 поступает на вход блока 8, на остальные входы которого поступают значения α1, G, W-1, I; с одного из выходов блока 8 коррекции значение M1(k+1/k) поступает на вход блока 15 коррекции, в котором формируется значение
Figure 00000064
которое поступает на вход блока 16 формирования суммы; значение
Figure 00000065
на выходе блока 16 поступает на вход блока 17 формирования матричной функции, с выхода блока 17 значение
Figure 00000066
поступает на вход линии задержки 18, на выходе которой формируется значение
Figure 00000067
, которое суммируется со значением (28) в блоке 16, на выходе которого формируется значение
Figure 00000068
; с выхода блока 18 значение
Figure 00000069
поступает на вход блока 15 коррекции, на вход блока 16 и на вход блока 19 формирования матричной функции, на выходе которого формируется значение
Figure 00000070
которое подается на вход блока 15 и вход блока 14, на другой вход которого поступает входное колебание; невязка измерения
Figure 00000071
с выхода блока 14 поступает на вход блока 15, на остальные входы которого поступают значения α2, G, W-1, I, с одного из выходов блока 15 коррекции значение М2(k+1/k) поступает на вход блока 21 коррекции, в котором формируется значение
Figure 00000072
которое поступает на вход блока 22 формирования суммы; значение
Figure 00000073
на выходе блока 22 поступает на вход блока 23 формирования матричной функции, с выхода блока 23 значение
Figure 00000074
поступает на вход линии задержки 24, на выходе которой формируется значение
Figure 00000075
, которое суммируется со значением (29) в блоке 22, на выходе которого формируется значение
Figure 00000076
; с выхода блока 24 значение
Figure 00000077
поступает на вход блока 21 коррекции, на вход блока 22 и на вход блока 25 формирования матричной функции, на выходе которого формируется значение
Figure 00000078
, которое подается на вход блока 21 и вход блока 20, на другой вход которого поступает входное колебание; невязка измерения
Figure 00000079
с выхода блока 20 поступает на вход блока 21, на остальные входы которого поступают значения α2, G, W-1, I.
Первый блок 2 коррекции работает следующим образом (фиг.2). Значения матричной функции
Figure 00000080
поступают на вход блока 2.1 формирования частных производных, с выхода которого значения
Figure 00000081
поступают на вход блока транспонирования матричной функции 2.2, с выхода которого значения
Figure 00000082
, невязка измерения
Figure 00000083
и значения W-1, α0 поступают на входы блока 2.3 формирования произведения, с выхода которого значение М0(k+1/k) поступает на вход блока 2.4 расчета точностных характеристик, на другие входы которого поступают значения G, I,
Figure 00000084
и на выходе которого формируется значение P1(k+1), которое поступает на вход блока 2.5 формирования произведения, на другой вход которого поступает значение М0(k+1/k) с выхода блока 2.3. Выход блока 2.5 является выходом блока 2 коррекции.
Второй блок 8 коррекции работает следующим образом (фиг.3). Значение матричной функции
Figure 00000085
поступает на вход блока 8.1 формирования частных производных, с выхода которого значения
Figure 00000086
поступают на вход блока транспонирования матричной функции 8.2, с выхода которого значения
Figure 00000087
, а также значение невязки измерения
Figure 00000088
, α1, W-1 поступают на вход блока 8.3 формирования произведения, с выхода которого значение
Figure 00000089
поступает на вход блока 8.4 суммирования, на другой вход которого поступает значение М0(k+1/k); значение
Figure 00000090
с выхода блока 8.4 поступает на вход блока 8.5 расчета точностных характеристик, на другие входы которого поступают значения G, W-1,
Figure 00000060
, а на выходе формируется значение P1(k+1), которое умножается на (28) в блоке 8.6 формирования произведения; выход блока 8.6 является выходом блока 8 коррекции. Третий блок 15 коррекции и четвертый блок коррекции 14 работают аналогично.
Первый блок 2.4 расчета точностных характеристик работает следующим образом (фиг.5). Значение матричной функции
Figure 00000091
поступает на вход блока 14 формирования частных производных, с выхода которого значение
Figure 00000092
поступает на вход блока 21 транспонирования матричной функции и на вход блока 22 формирования произведения, на вход которого также поступает значение
Figure 00000093
с выхода блока 21 и значение P0(k) с выхода линии задержки 23, на вход которой с выхода блока 31 формирования произведения, являющегося выходом блока 2.3, поступает значение Р0(k+1); значение
Figure 00000094
с выхода блока 22 поступает на вход блока 26 формирования суммы, на другой вход которого поступает значение GGT, сформированное в блоке 25 формирования произведения, на вход которого поступает значение G и значение GT, сформированное в блоке 24 транспонирования матриц, на вход которого также поступает значение G; значение матричной функции М0(k+1/k) поступает на вход блока 27, формирования частных производных, с выхода которого значение
Figure 00000095
поступает на вход блока 28, на другой вход которого поступает значение P0(k+1/k), сформированное на выходе блока 26; значение
Figure 00000096
поступает на вход блока 29 формирования разности, на другой вход которого поступает значение I; с выхода блока 29 значение
Figure 00000097
поступает на вход устройства 30 обращения матриц, выход которого соединен с входом блока 31 формирования произведения, на другой вход которого поступает значение Р0(k+1/k) с выхода блока 26. Блок расчета точностных характеристик 8.5 второго, третьего и четвертого блоков коррекции работает аналогичным образом. На выходе блока 8.5 формируется значение Р1(k+1).
Блок расчета параметра регуляризации работает следующим образом. В исходном состоянии на вход блока расчета параметра регуляризации поступает
Figure 00000098
- значение оценки информационного процесса, где k - текущий момент времени, которое соответственно передается на линии задержки 26.2, 26.7, 26.16, 26.29 (фиг.5) и блок формирования матричной функции
Figure 00000099
, с выхода которого значение
Figure 00000099
поступает на вход 26.4, 26.12, 26.24 и 26.40 линий задержки; в следующей итерации фильтра на эти же блоки поступает текущее значение оценки информационного процесса, с выходов линий задержки 26.2, 26.7, 26.16, 26.29 и 26.4, 26.12, 26.24, 26.40 на линии задержки 26.8, 26.17, 26.30 и 26.13, 26.25, 26.41 поступают соответственно k-1 значение оценки и матричной функции
Figure 00000100
, эти же значения оценок и матричной функции поступают на входы блоков 26.3 формирования произведения числителя коэффициента регуляризации и блоков 26.5 транспонирования матричной функции, 26.6 формирования произведения; далее на входы линий задержки 26.2, 26.7, 26.16. 26.29 опять поступает текущее значение оценки информационного процесса, с этих линий задержки k-1 значение оценки поступает на 26.8, 26.17. 26.30 линий задержки, с выхода которых k-2 значение оценки информационного процесса поступает дальше на вход блока 26.9 и линий задержки 26.18, 26.31; аналогично и для 26.4, 26.12, 26.24, 26.40, а также 26.13, 26.25, 26.41 линий задержки, только для матричной функции Н; в результате на n+5 и последующих итерациях будем иметь на выходах линий задержки 26.2, 26.8, 26.18, 26.32 значения k-1, k-2, k-3, k-4 оценок информационного процесса соответственно и на выходах 26.4, 26.13, 26.26, 26.43 линий задержки значения матричной функции для тех же оценок. Эти значения оценок информационного процесса поступают на блоки 26.3, 26.9, 26.19, 26.33 формирования произведения числителя коэффициента регуляризации и 26.10, 26.20, 26.34 формирования прогноза; в блоки 26.3, 26.9, 26.19, 26.33 поступают из блока 13 G, W-1 и у со входа устройства, с выхода этих блоков значения
Figure 00000101
для
Figure 00000102
, хηn(j) для j=k, поступают на вход блока 26.46 формирования суммы; с выхода 26.10 блока на 26.11 блок поступает значение прогноза для k-1 момента времени, в блок 26.11 поступают значения G, W-1, y, с выхода 26.11 блока на вход 26.46 блока поступает выражение, аналогичное (29), с 26.20 блока в 26.23 блок и 26.21 блок поступает значение прогноза для k-2 момента времени, в блок 26.23 передаются значения G, W-1, y, с 26.23 блока на вход 26.45 блока поступает выражение, аналогичное (29), с выхода блока 26.21 поступает значение прогноза для k-1 момента времени, в блок 26.22 поступают значения G, W-1, y, с выхода 26.22 блока на 26.46 блок поступает выражение, аналогичное (29), на входах 26.34, 26.35, 26.36 блоков имеем
Figure 00000103
, xηn(k-3), хηn(k-2), с выходов этих блоков имеем значения для k-3, k-2, k-1 момента времени, которые поступают на блоки 26.38, 26.39, 26.37 соответственно, кроме того, блок 26.36 использует результаты работы 26.35 блока, а 26.35 блок 26.36 блока, и, наконец, в 26.34 блок передаются начальные значения, на входы блоков 26.37, 26.38, 26.39 поступают значения G, W-1, y, с выхода блоков 26.37, 26.38, 26.39 в блок 26.46 поступают выражения, аналогичные (29); с линий задержки 26.4, 26.13, 26.26, 26.43 на вход блоков 26.5, 26.14, 26.27, 26.44 транспонирования матричной функции и 26.6, 26.15, 26.28, 26.45 формирования произведения передается
Figure 00000104
, после чего с выхода блоков 26.6, 26.15, 26.28, 26.45 на вход блока 26.45 формирования сумм поступают значения
Figure 00000105
,
Figure 00000106
; с блока 26.46 в блок 26.47 формирования отношения передается значение
Figure 00000107
, на вход 26.47 блока с 26.45 блока поступает
Figure 00000108
на выходе блока 26.47 формирования отношения получаем значение коэффициента регуляризации αn+1. Второй блок 27 и третий блок 28 расчета параметра регуляризации работают аналогично.
Первый блок формирования произведения числителя параметра регуляризации 26.3 (фиг.6) работает следующим образом. На вход блока 26.3.1 формирования матричной функции поступает значение оценки информационного процесса
Figure 00000102
, с выхода блока на вход блока 26.3.2 формирования частных производных матричной функции и блока 26.3.4, далее с выхода блока 26.3.2 значения поступают на вход блока 26.3.3 транспонирования матричной функции, откуда на вход блока 26.3.5 формирования произведения поступает
Figure 00000109
, на другие входы поступают значения GT, W-1 и значение
Figure 00000110
с выхода блока 26.3.4, на один из входов которого поступает y, значение GT поступает на вход блока 26.3.5 с выхода блока 26.3.6 транспонирования матричной функции, на вход которого соответственно с выхода 13 блока поступает значение G; с выхода блока 26.3.5 имеем значение
Figure 00000111
, которое поступает соответственно на вход блоков 26.3.7 транспонирования матричной функции и 26.3.8 формирования произведения, в результате на выходе устройства будем иметь значения
Figure 00000112
.
Первый блок расчета прогноза работает следующим образом. На вход блоков 26.10.4 формирования матричной функции и 26.10.3 формирования частных производных поступает значение оценки информационного процесса
Figure 00000102
, с выхода блока 26.10.4 значение
Figure 00000113
поступает на другой вход блока 26.10.3 и блоков 26.10.5 формирования произведения, 26.10.6 формирования суммы, с выхода блока 26.10.3 значение
Figure 00000114
поступает на вход блока 26.10.1 формирования произведения, на другие входы которого поступают значения Gn(k) и Δt, с выхода этого блока на вход блока 26.10.2 формирования суммы передается значение
Figure 00000115
на другой вход блока 10.2 поступает значение Gn(k), в результате на выходе 26.10.2 блока, выход которого является первым выходом блока расчета прогноза, имеем прогноз Gn(j), j=k+1; на второй вход блока 26.10.5 поступает значение Δt, с выхода которого значение
Figure 00000116
поступает на второй вход блока 26.10.6, выход которого является вторым выходом блока расчета прогноза, с выхода которого имеем хηn(j). Остальные блоки расчета прогноза работают аналогично.
Источники информации
1. А.А.Венгеров, В.А.Щаренский. Прикладные вопросы оптимальной линейной фильтрации. - М.: Энергоиздат, 1982. - 96 с.
2. А.Фарина, Ф.Студер. Цифровая обработка радиолокационной информации. Сопровождение целей. - М.: Радио и связь, 1993. - 118 с.
3. Патент №2362255. Россия. 2009. Цифровой интеллектуальный рекурсивный фильтр // Костоглотов А.А., Кузнецов А.А., Пономарев А.С., Лазаренко С.В.
4. Костоглотов А.А. Синтез интеллектуальных измерительных процедур на основе принципа регуляризации А.Н.Тихонова // Измерительная техника, №1, 2001. С.8-12.
5. Костоглотов А.А. Метод последовательных приближений в теории фильтрации // Автоматика и вычислительная техника, №3, 2000. С.53-63.
6. Костоглотов А.А., Кузнецов А.А. Синтез интеллектуальной измерительной процедуры на основе метода минимальных ошибок // Измерительная техника, №7, 2005. С.8-13.
7. Тихонов А.Н., Арсенин В.Я. Методы решения некорректных задач. - М.: Наука, 1986. - 288 с.
8. Васильев Ф.П. Методы решения экстремальных задач. - М.: Наука, 1981. - 106 с.
9. Костоглотов А.А. Цифровая интеллектуальная измерительная процедура // Измерительная техника, №7, 2002. С.16-21.

Claims (1)

  1. Цифровой интеллектуальный многокаскадный фильтр, содержащий первый, второй, третий и четвертый блоки формирования разности, первый, второй, третий и четвертый блоки коррекции, блок формирования и выдачи априорных данных, первый, второй и третий блоки формирования суммы, первый, второй, третий, четвертый, пятый, шестой, седьмой и восьмой блоки формирования матричных функций, первую, вторую, третью и четвертую линии задержки, первый, второй и третий блоки расчета параметра регуляризации, при этом первый, второй, третий и четвертый выходы блока формирования и выдачи априорных данных соединены соответственно со вторым, третьим, четвертым и пятым информационным входами первого, второго, третьего и четвертого блока коррекции, первый информационный выход первого блока коррекции соединен с первым информационным входом первого блока формирования суммы, выход которого соединен с информационным входом первого блока формирования матричной функции, и с четвертым информационным входом первого блока расчета параметра регуляризации, первый информационный вход которого является входом устройства, третий, второй и пятый информационные входы первого блока расчета параметра регуляризации соединены соответственно с первым, вторым и пятым выходами блока формирования и выдачи априорных данных, выход первого блока расчета параметра регуляризации соединен с девятым входом второго блока коррекции, выход первого блока формирования матричной функции соединен с информационным входом первой линии задержки, выход которой соединен со вторым информационным входом первого блока формирования суммы, с седьмым информационным входом первого блока коррекции и с информационным входом второго блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом первого блока коррекции, первый информационный выход второго блока формирования матричной функции соединен со вторым информационным входом первого блока формирования разности, выход которого соединен с первым информационным входом первого блока коррекции, второй информационный выход первого блока коррекции соединен с восьмым информационным входом второго блока коррекции, первый информационный выход которого соединен с первым информационным входом второго блока формирования суммы, выход которого соединен с информационным входом третьего блока формирования матричной функции и с четвертым входом второго блока расчета параметра регуляризации, первый информационный вход которого является входом устройства, второй, третий и пятый информационные входы второго блока расчета параметра регуляризации соединены соответственно с первым, вторым и пятым выходами блока формирования и выдачи априорных данных, выход второго блока расчета параметра регуляризации соединен с девятым входом третьего блока коррекции, выход третьего блока формирования матричной функции соединен с информационным входом второй линии задержки, выход которой соединен со вторым информационным входом второго блока формирования суммы, седьмым информационным входом второго блока коррекции и с информационным входом четвертого блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом второго блока коррекции, первый информационный выход четвертого блока формирования матричной функции соединен со вторым информационным входом второго блока формирования разности, выход которого соединен с первым информационным входом второго блока коррекции, второй информационный выход второго блока коррекции соединен с восьмым информационным входом третьего блока коррекции, первый информационный выход которого соединен с первым информационным входом третьего блока формирования суммы, выход которого соединен с информационным входом пятого блока формирования матричной функции, выход которого соединен с информационным входом третьей линии задержки, выход которой соединен со вторым информационным входом третьего блока формирования суммы, седьмым информационным входом третьего блока коррекции и с информационным входом шестого блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом третьего блока коррекции, первый информационный выход шестого блока формирования матричной функции соединен со вторым информационным входом третьего блока формирования разности, выход которого соединен с первым информационным входом третьего блока коррекции, второй информационный выход третьего блока коррекции соединен с восьмым информационным входом четвертого блока коррекции, первый информационный выход которого соединен с первым информационным входом четвертого блока формирования суммы, выход которого является выходом устройства и соединен с информационным входом седьмого блока формирования матричной функции, выход которого соединен с информационным входом четвертой линии задержки, выход которой соединен со вторым информационным входом четвертого блока формирования суммы, седьмым информационным входом четвертого блока коррекции и с информационным входом восьмого блока формирования матричной функции, второй информационный выход которого соединен с шестым информационным входом четвертого блока коррекции, первый информационный выход восьмого блока формирования матричной функции соединен со вторым информационным входом четвертого блока формирования разности, выход которого соединен с первым информационным входом четвертого блока коррекции, первый информационный вход первого блока формирования разности, первый информационный вход второго блока формирования разности, первый информационный вход третьего блока формирования разности и первый информационный вход четвертого блока формирования разности являются входами устройства, в устройство введен третий блок расчета параметра регуляризации, первый информационный вход которого является входом устройства, второй, третий и пятый информационные входы третьего блока расчета параметра регуляризации соединены соответственно с первым, вторым и пятым выходами блока формирования и выдачи априорных данных, четвертый информационный вход третьего блока расчета параметра регуляризации соединен с выходом третьего блока формирования суммы, выход третьего блока расчета параметра регуляризации соединен с девятым входом четвертого блока коррекции.
RU2010115924/08A 2010-04-21 2010-04-21 Цифровой интеллектуальный многокаскадный фильтр RU2436228C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2010115924/08A RU2436228C1 (ru) 2010-04-21 2010-04-21 Цифровой интеллектуальный многокаскадный фильтр

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2010115924/08A RU2436228C1 (ru) 2010-04-21 2010-04-21 Цифровой интеллектуальный многокаскадный фильтр

Publications (1)

Publication Number Publication Date
RU2436228C1 true RU2436228C1 (ru) 2011-12-10

Family

ID=45405771

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2010115924/08A RU2436228C1 (ru) 2010-04-21 2010-04-21 Цифровой интеллектуальный многокаскадный фильтр

Country Status (1)

Country Link
RU (1) RU2436228C1 (ru)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2599970C2 (ru) * 2012-05-10 2016-10-20 Долби Лабораторис Лайсэнзин Корпорейшн Многокаскадный iir-фильтр и распараллеленная фильтрация данных таковым
RU2653485C1 (ru) * 2017-06-19 2018-05-08 Публичное акционерное общество "Российский институт мощного радиостроения" (ПАО "РИМР") Способ адаптивного выбора оптимального параметра алгоритма коррекции сигналов

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2599970C2 (ru) * 2012-05-10 2016-10-20 Долби Лабораторис Лайсэнзин Корпорейшн Многокаскадный iir-фильтр и распараллеленная фильтрация данных таковым
RU2653485C1 (ru) * 2017-06-19 2018-05-08 Публичное акционерное общество "Российский институт мощного радиостроения" (ПАО "РИМР") Способ адаптивного выбора оптимального параметра алгоритма коррекции сигналов

Similar Documents

Publication Publication Date Title
Rego et al. Guaranteed methods based on constrained zonotopes for set-valued state estimation of nonlinear discrete-time systems
Ahmad et al. On a coupled system of fractional differential equations with coupled nonlocal and integral boundary conditions
Schoukens et al. Identification of Wiener–Hammerstein systems by a nonparametric separation of the best linear approximation
Nestler et al. Fast Ewald summation based on NFFT with mixed periodicity
Daid et al. On the convergence of the unscented Kalman filter
Lozi et al. A new reliable numerical method for computing chaotic solutions of dynamical systems: the Chen attractor case
Ebrahimi et al. Design of a robust central difference Kalman filter in the presence of uncertainties and unknown measurement errors
Ponomareva et al. Higher order sigma point filter: A new heuristic for nonlinear time series filtering
CN107766293B (zh) 部分采样数据规则性缺失时的信号频谱分析方法及系统
RU2436228C1 (ru) Цифровой интеллектуальный многокаскадный фильтр
RU2357357C2 (ru) Цифровой интеллектуальный рекурсивный фильтр
Heydari et al. Numerical solution of distributed-order time fractional Klein–Gordon–Zakharov system
RU2452080C1 (ru) Цифровой многоитерациональный фильтр
Wu et al. Reduced-order H2 filtering for discrete linear repetitive processes
Madureira et al. Numerical analysis for a thermoelastic diffusion problem in moving boundary
RU2362265C1 (ru) Цифровой интеллектуальный итеративный фильтр
Yong et al. An asynchronous sensor bias estimation algorithm utilizing targets’ positions only
RU2257667C2 (ru) Цифровой рекурсивный фильтр
RU2393522C1 (ru) Устройство управления подвижным объектом
Maier et al. Distribution of a cotangent sum related to the Nyman–Beurling criterion for the Riemann Hypothesis
Nagar et al. A bivariate distribution whose marginal laws are gamma and Macdonald
Carvajal-Rubio et al. Discretization of the robust exact filtering differentiator based on the matching approach
Straka et al. Design of discrete second order filters for continuous-discrete models
Borisov Analysis and estimation of the states of special jump markov processes. ii. optimal filtration in wiener noise
CN100520442C (zh) 采用总体最小二乘和均衡算法的到达时差定位方法

Legal Events

Date Code Title Description
MM4A The patent is invalid due to non-payment of fees

Effective date: 20120422