Цифровой многоитерациональный фильтр

Цифровой многоитерациональный фильтр

Иллюстрации

Показать все

Реферат

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

Известно устройство - расширенный калмановский фильтр [1, 2], недостатком которого является ограниченность функциональных возможностей, обусловленная линейной структурой обрабатываемых процессов и устройство [3], а так же эвристическое задание параметра регуляризации.

Наиболее близким по технической сущности к заявленному изобретению является цифровой интеллектуальный итеративный фильтр [3], содержащий первый, второй, третий блоки формирования разности; первый, второй, третий блоки коррекции; блок формирования и выдачи априорных данных; первый, второй, третий блоки формирования суммы; первый, второй, третий, четвертый, пятый и шестой блоки формирования матричных функций; первую, вторую и третью линии задержки; первый и второй блок расчета параметра регуляризации. Однако данное устройство в некоторых практически важных случаях не позволяет получить требуемой точности.

Улучшение точностных характеристик фильтрации параметров случайных процессов является актуальным направлением.

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

Одним из способов повышения точности оценки параметров динамических систем является использование методов решения некорректных задач на основе принципов регуляризации. Эффективность применения регуляризации для случая непрерывных систем доказана и получила свое выражение в виде метода А.Н.Тихонова [4] и его разновидности в виде метода итерационной регуляризации [5]. Получим уравнения фильтрации с использованием метода итерационной регуляризации для дискретной системы [6-9].

Пусть динамика оцениваемых параметров описывается системой разностных уравнений в дискретном времени

где - вектор состояния исследуемой системы; - вектор неизвестных внешних воздействий; переходная функция - непрерывная вместе с частными производными вектор-функция своих аргументов; G∈E^M×E^M - матрица интенсивности внешних воздействий; k, N, М - натуральные числа. Предполагается, что матрица имеет обратную матрицу.

Наблюдаемый сигнал, получаемый на выходе модели измерительной системы, описывается дискретным уравнением

где - вектор наблюдения,

- вектор дискретного белого гауссовского шума с известными локальными характеристиками

W - ковариационная матрица размерности L×L, δ(·) - векторная дельта - функция; - сигнальная вектор-функция, непрерывна вместе с частными производными; L, l - натуральные числа.

Поставим задачу синтеза рекуррентного фильтра оценки x^*(k), оптимального в смысле минимума функционала, характеризующего ошибку измерения

В силу непрерывности вектор - функции F(·), решение уравнения (1) непрерывно зависит от η(k), поэтому функционал ошибки (3) на каждом решении системы (1) непрерывно зависит от η(k). Таким образом, задача определения оценки x^*(k), доставляющей минимум (3), равносильна задаче определения

Задача (1), (2), (4) является некорректно поставленной обратной задачей [7]. Найти значения векторов х^*(k), η^*(k), путем решения совокупности уравнений (1), (2), (4) в условиях некорректности исходной задачи достаточно сложно, в связи с этим широкое распространение получили итерационные градиентные методы. Однако использование таких методов может привести к расходящейся последовательности приближений. Поэтому применение любого итерационного метода для решения задачи (1), (2), (4) требует определения регуляризирующего семейства операторов, в котором параметром регуляризации является номер итерации.

В соответствии с общим определением регуляризирующего семейства операторов по А.Н.Тихонову [7] будем говорить, что итерационный метод

в котором числовой параметр α_n удовлетворяет условиям:

где Δ(η_n) - невязка, порождает регуляризирующее семейство операторов, в котором параметром является номер итерации, если для любого начального приближения η₀ и для любого значения погрешности исходных данных σ, удовлетворяющего условию 0<|σ|<σ₀, σ₀=const, существует номер n(σ), такой, что

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

Согласно [8] выражение для определения градиента в точке η_n(k) имеет вид

где x_n(k) - решение задачи (1) при η_n(k), а вектор ψ_n(k) определяется из условий

Зная выражение для градиента (7) функционала (3), можно переписать выражение (6) для параметра регуляризации [9] в следующем виде

Выбор последовательности параметров α_n, удовлетворяющей условию (9), позволит реализовать фильтр цифровой обработки измерительной информации повышенной точности.

Для реализации итерационного метода (5) требуется определить градиент функционала (3), определяемый выражением (7). Приняв за нулевое приближение η₀(k)=0, запишем итерационную последовательность (5) в развернутой форме для

В результате, с учетом (7), имеем последовательность дискретных двухточечных краевых задач (ДДТКЗ) вида

Введем обозначение и помножим каждое из уравнений для сопряженных векторов ψ_i на величину α_i, тогда уравнения (10) принимают следующий вид

Для получения алгоритма рекуррентного оценивания вектора состояния необходимо воспользоваться методом инвариантного погружения в дискретном варианте. Заметим, что уравнение для вектор - функций λ в ДДТКЗ (9) записано в обратном времени. Это требует его преобразования к виду, отражающему зависимость λ_n(k+1) от λ_n(k) и x_n(k). Производя соответствующие преобразования с учетом выражения для x_n(k+1) из (11) и используя разложение в ряд Тейлора в окрестности F(x_n(k),k), получим следующую последовательность ДДТКЗ

где функции β и γ вводятся для сокращения записи.

Заменим условие на конце λ_n(N)=0 более общим условием λ_n(N)=с и пусть N и с - переменные величины. Тогда значение вектора x_n(N) определяется как функция величин N и с

x_n(N)=r[с,N].

Изменение величины N Ha N+1 дает приращение Δс, тогда

Запишем выражение для r(с+Δс, N+1), используя аппарат конечных разностей

r(с+Δс, N+1)=r(c, N)+Δr_c(с, N)+Δr_N(с, N)+Δ²r_cN(с, N),

или, учитывая (13), получим

где

Согласно (12) выражения для Δx_n и Δс имеют вид

Разрешить разностное уравнение (14) относительно r(c,N), т.е. найти общее аналитическое решение не удается, и обычно обращаются к приближенным методам. Предположим, что r(c,N) линейна по с

где - оценка вектора состояния в момент N, P_n(N) - некоторая матрица размерности М×М.

Вычислим разности, входящие в выражение (14), используя выражение (16)

Подставив выражения (15), (16), (17) в (14), получим

Разлагая β и γ в ряд Тейлора в окрестности и пренебрегая членами порядка выше первого, можно записать уравнение (18) в виде

Соотношение (19) выполняется при с→0, поэтому, приравнивая коэффициенты при первой и нулевой степени с, получим разностные уравнения для и

Запишем ДДТКЗ (12) для случая, когда k=N при этом учтем, что это все ДДТКЗ для i=0,…,n-1 являются разрешенными и соответственно оценки x_i представляют собой известные функции параметра k. Имеем

Поскольку разницы в переобозначении матриц P_n и нет, запишем последовательность уравнений для оценки процесса (1), предполагая, что N постоянно изменяется и k=N, а также учитывая условие (7) налагаемое на параметр регуляризации, в виде

Последовательность уравнений (24) представляет собой цифровой интеллектуальный многокаскадный фильтр, который позволяет осуществить процесс цифровой обработки измерительной информации для дискретных динамических систем. Если сравнить полученные уравнения с уравнениями цифрового рекурсивного фильтра [3], то становится ясно, что они отличны друг от друга за счет введения дополнительной итерации, а также дополнительных связей последовательности параметров α_i, изменяющих общий коэффициент обратной связи в уравнении для оценки с оценками сигнальной вектор-функцией Н и матрицами 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 дискретной нелинейной системы третьего порядка

где параметр τ имеет смысл интервала времени, через который поступает измерительная информация в виде

Таким образом, как следует из соотношений (24), введение новых структурных элементов и связей позволяет в совокупности с общими признаками получить технический результат, состоящий в уменьшении дисперсии ошибок полученных на выходе фильтра оценок входных процессов.

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

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

Рассчитанный весовой коэффициент параметра регуляризации имеет вид [10]:

На основе прототипа ниже приводится описание изобретения расчитаного по формуле (27).

Сущность изобретения поясняется на рисунке 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 представлена структурная схема блока расчета точностных характеристик, входящего в первый, второй, третий и четвертый блоки коррекции, который содержит блок 20 формирования частных производных матричной функции, блок 21 транспонирования матриц, блок 22 формирования произведения, линию 23 задержки, блок 24 транспонирования матриц, блок 25 формирования произведения, блок 26 формирования суммы, блок 27 формирования частных производных матричной функции, блок 28 формирования произведения, блок 29 формирования разности, устройство 30 обращения матриц, блок 31 формирования произведения.

На рисунке 5 представлена структурная схема блока расчета параметра регуляризации, который содержит блок формирования матричной функции 26.1, линии задержки 26.2, 26.4, 26.8, 26.9, 26.13, 26.14, 26.18, 26.19, 26.20, 26.26, 26.27, 26.28, 26.32, 26.33, 26.34, 26.35, 26.43, 26.44, 26.45, 26.46, блоки 26.3, 26.10, 26.12, 26.21, 26.24, 26.25, 26.36, 26.40, 26.41, 26.42 формирования произведения числителя параметра регуляризации, блоки 26.5, 26.15, 26.29, 26.47 транспонирования, блоки 26.6, 26.16, 26.30, 26.48 формирования произведения, блоки 26.7, 26.17, 26.31, 26.49 формирования разности, блоки 26.11, 26.22, 26.23, 26.37, 26.38, 26.39 расчета прогноза, блоки 26.51, 26.50 формирования суммы, блок 26.52 формирования произведения знаменателя параметра регуляризации, блок 26.53 формирования отношений.

На рисунке 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.11.1 и 26.11.5 формирования произведения, блоки 26.11.2 и 26.11.6 формирования суммы, блок 26.11.3 формирования частных производных матричной функции, блок 26.11.4 формирования матричной функции.

Для первого блока расчета параметра регуляризации (рисунок 5) информационные входы линий задержки 26.2, 26.8, 26.18, 26.32, соединены с информационным выходом блока 3 формирования суммы, (рисунок 1) первые информационные входы линий задержки 26.4, 26.13, 26.26, 26.43, соединены с первым информационным выходом первого блока 26.1 формирования матричной функции , информационный вход которого соединен с информационным выходом блока 3, информационный выход первой линии задержки 26.2 соединен со вторым информационным входом первого блока 26.3 формирования произведения числителя параметра регуляризации, первый и третий его информационные входы соединены со вторым и первым информационными выходами блока 13 формирования и выдачи априорных данных соответственно, а четвертый информационный вход - со входом устройства, первый информационный выход первого блока 26.3 формирования произведения числителя параметра регуляризации соединен с шестым и седьмым информационными входами второго блока 26.50 формирования суммы, первый информационный выход которого соединен с первым информационным входом блока 26.53 формирования отношения, первый выход блока 26.53 формирования отношения является выходом блока 26; первый информационный выход второй линии задержки 26.4 соединен с первым информационным входом первого блока транспонирования 26.5 и первым информационным входом первого блока 26.6 формирования произведения, первый информационный выход блока 26.5 соединен со вторым информационным входом блока 26.6 формирования произведения, первый информационный выход которого соединен со вторым информационным входом первого блока 26.7 формирования разности, первый информационный вход которого соединен с четвертым информационным входом блока 26.3, первый информационный выход блока 26.7 соединен с четвертым информационным входом первого блока 26.51 формирования суммы, первый информационный выход которого соединен с первым информационным входом блока 26.52 формирования произведения знаменателя параметра регуляризации, а второй вход блока 26.52 соединен с пятым выходом блока 13 формирования и выдачи априорных данных; первый информационный выход блока 26.52 формирования произведения знаменателя параметра регуляризации соединен со вторым входом блока 26.53 формирования отношения; первый информационный выход третей 26.8 линии задержки соединен с первым информационным входом четвертой 26.9 линии задержки, первый информационный выход которой соединен со вторым информационным входом второго блока 26.10 формирования произведения числителя параметра регуляризации и вторым информационным входом первого блока 26.11 вычисления прогноза, второй информационный выход которого соединен со вторым информационным входом третьего блока 26.12 формирования произведения числителя параметра регуляризации, третьи и четвертые информационные входы блоков 26.10, 26.12, соединены с первым информационным выходом блока 13 и входом устройства, первый информационный выход блока 26.10 соединен с четвертым и третьим информационными входами блока 26.50, первый информационный выход блока 26.12 соединен с десятым информационным входом блока 26.50 формирования суммы; первый информационный выход пятой 26.13 линии задержки соединен с первым информационным входом шестой 26.14 линии задержки, первый информационный выход которой соединен с первым информационным входом второго блока 26.15 транспонирования и первым информационным входом второго блока 26.16 формирования произведения, первый информационный выход второго 26.15 блока транспонирования соединен со вторым информационным входом блока 26.16, первый информационный выход которого соединен со вторым информационным входом второго блока формирования разности 26.17, первый информационный вход которого соединен с четвертым информационным входом блока 26.10 и первый информационный выход блока 26.17 соединен с третьим информационным входом блока 26.51 формирования суммы; первый информационный выход седьмой линии задержки 26.18 соединен с первым информационным входом восьмой линии задержки 26.19, первый информационный выход которой соединен с первым информационным входом девятой 26.20 линии задержки, первый информационный выход которой соединен со вторым информационным входом четвертого 26.21 блока формирования произведения числителя параметра регуляризации и второго блока 26.22 вычисления прогноза, первый информационный выход которого соединен с первыми информационными входами третьего блока 26.23 вычисления прогноза и 26.25 формирования произведения числителя параметра регуляризации, третьи и четвертые информационные входы блоков 26.21, 26.24 и 26.25 соединены с первым информационным выходом блока 13 и входом устройства, первый информационный выход блока 26.21 соединен с двенадцатым и одиннадцатыми информационными входами второго блока 26.50 формирования суммы, первые информационные выходы блоков 26.25, 26.24 соединены с тринадцатым и четырнадцатым информационными входами первого блока 26.50 соответственно; первый информационный выход десятой линии задержки 26.26 соединен с первым информационным входом одиннадцатой линии задержки 26.27, первый информационный выход которой соединен с первым информационным входом двенадцатой 26.28 линии задержки, первый информационный выход которой соединен с первыми информационным входом третьего блока 26.29 транспонирования матричной функции и первым информационным входом третьего блока 26.30 формирования произведения, первый информационный выход третьего блока 26.29 транспонирования соединен со вторым информационным входом третьего блока 26.30 формирования произведения, первый информационный выход которого соединен с первым информационным входом третьего блока формирования разности 26.31, первый информационный выход которого соединен с первым информационным входом блока 26.51 формирования суммы, и первый информационный вход третьего блока 26.31 соединен с четвертым информационным входом блока 26.21; первый информационный выход тринадцатой линии задержки 26.32 соединен с первым информационным входом четырнадцатой линии задержки 26.33, первый информационный выход которой соединен с первым информационным входом пятнадцатой 26.34 линии задержки, первый информационный выход которой в свою очередь соединен с первым информационным входом шестнадцатой линии задержки 26.35, выход которой соединен со вторым информационным входом седьмого блока 26.36 формирования произведения числителя параметра регуляризации и вторым информационным входом четвертого блока 26.37 вычисления прогноза, второй информационный выход которого соединен со вторым информационным входом девятого блока 26.41 формирования произведения числителя параметра регуляризации и вторым информационным входом пятого блока 26.38 вычисления прогноза, второй информационный выход последнего соединен со вторым информационным входом десятого блока 26.42 формирования произведения числителя параметра регуляризации и со вторым информационным входом шестого блока вычисления прогноза 26.39, второй информационный выход которого соединен со вторым информационным входом восьмого блока 26.40 формирования произведения числителя параметра регуляризации, третьи и четвертые информационные входы блоков 26.36, 26.40, 26.41, 26.42 соединены соответственно с первым информационным выходом блока 13 и входом устройства, первый информационный выход блока 26.36 соединен с девятым и восьмым информационными входами блока 26.50 формирования суммы, четвертый вход этого блока соединен с первым входом четвертого блока формирования разности 26.49; первый информационный выход блока 26.40 формирования произведения числителя параметра регуляризации соединен с пятым информационным входом блока 26.50, первые информационные выходы блоков 26.41 и 26.42 формирования произведения числителя параметра регуляризации соединены с первым и вторым информационными входами блока 26.50 соответственно; первый информационный выход семнадцатой линии задержки 26.43 соединен с первым информационным входом восемнадцатой линии задержки 26.44, первый информационный выход которой соединен с первым информационным входом девятнадцатой 26.45 линии задержки, первый информационный выход которой в свою очередь соединен с первым информационным входом двадцатой линии задержки 26.46, информационный выход которой соединен с первыми информационными входами четвертого блока 26.47 транспонирования и четвертого блока 26.48 формирования произведения, первый информационный выход блока 26.47 соединен со вторым информационным входом блока 26.48, первый информационный выход которого соединен со вторым информационным входом блока формирования разносности 26.49, который соединен со вторым информационным входом блока 26.51 формирования суммы; а также первый информационный вход блока 26.49 соединен с четвертым информационным входом блока 26.36; первые информационные входы второго 26.10, четвертого 26.21 и седьмого 26.36 блоков формирования произведения числителя параметра регуляризации соединены со вторым информационным выходом блока 13 формирования и выдачи априорных данных, а первые информационные входы третьего 26.12, пятого 26.24, шестого 26.25, восьмого 26.40, девятого 26.41 и десятого 26.42 блоков формирования произведения числителя параметра регуляризации соединены, соответственно, с первыми информационными выходами первого 26.21, третьего 26.23, второго 26.22, шестого 26.39, четвертого 26.37 и пятого 26.38 блоков вычисления прогноза; первые информационные входы первого 26.11, второго 26.22 и четвертого 26.37 блоков вычисления прогноза соединены со вторым информационным выходом блока 13, а первые информационные входы третьего 26.23, пятого 26.38, шестого 26.39 блоков вычисления прогноза соединены с первыми информационными выходами соответственно второго 26.22, четвертого 26.37 и пятого 26.38 блоков вычисления прогноза, третьи информационные входы блоков вычисления прогноза соединены с пятым информационным выходом блока 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 формирования суммы, выход которого соединен с четвертым информационным входом в