Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

SYNTHESIS OF OPTIMAL FILTERS ANALOG ELECTRONICS USING A DEVICE MATRIX OPERATORS

Melnikov D.V. 1
1 Bauman Moscow State Technical University (Kaluga Branch)
In the article the method of synthesis of non-stationary filters for the analog electronics. For the synthesis filters is proposed to use the unit matrix operators, which is well established for the solution of control problems for a class of linear stationary and non-stationary systems. Algorithm method involves the implementation of two quite time-consuming stages: the calculation of the matrix operator of the filter and restore the differential equation of the filter is found in the matrix operator. In the first phase of the algorithm uses the same ideology, as in the authors ‘ projection matrix method of synthesis of nonlinear control systems – replacement of nonlinear elements of the equivalent matrix operators are using methods of mathematical programming. In the second stage, we solve the following problem. Known Mat-hex filter operator in the class of stationary or non-stationary linear systems, it is necessary to build a differential equation equivalent to this matrix operator. The structure of the differential equation of the regulator is chosen for reasons of hardware implementation, so the problem reduces to determining the coefficients of the equation. The unknown coefficients are calculated in the form of expansions in an orthonormal basis. It is shown that the algorithm of calculation of these coefficients is reduced to solving an overdetermined system of linear algebraic equations for the desired Fourier coefficients. The algorithm of the method is demonstrated on the example of constructing an optimal filter of the electronic device. The proposed approach to the construction of filters based on the matrix operators allows us to find the differential equation of optimal filters for a class of stationary and nonstationary random processes, with no restriction on the stationarity of the object.
matrix operator
orthonormal basis
synthesis
filter
stochastic process
variance
correlation function
1. Aung Ch.S., Makarenkov A.M., Mo P.S. Identifikacija sluchajnyh parametrov matema-ticheskoj modeli jelektrogidravlicheskogo sledjashhego privoda // Fundamentalnye is-sledovanija. 2016. no. 2–2. рр. 231–235.
2. Martjanov P.S. Sintez analogovyh filtrov na jelementarnyh zvenjah dlja radioteh-nicheskih sistem: dissertacija ... kandidata tehnicheskih nauk: M., 2010. 156 р.
3. Melnikov D.V., Egupov N.D., Min Ch.Tu. Metod matrichnyh operatorov sinteza nesta-cionarnyh filtrov // Voprosy radiojelektroniki. 2013. T. 1, no. 4. рр. 65–75.
4. Melnikov D.V., Kornjushin P.Ju., Min Ch.T., Chzho T.A., Okar M. Proekcionno–matrichnyj podhod k analizu i sintezu sistem upravlenija jelektrojenergeticheskih si-stem // Nauchnoe obozrenie. 2015. no. 2. рр. 88–97.
5. Melnikov D.V., Kornjushin Ju.P., Min Ch.Tu. Spektralnyj metod sinteza nelinej-nyh radiojelektronnyh sistem upravlenija // Voprosy radiojelektroniki. 2013. T. 1, no. 4. рр. 116–124.
6. Melnikov D.V., Min Ch.Tu. Issledovanie radiotehnicheskih ustrojstv v stohasticheskih rezhimah // Voprosy radiojelektroniki. 2013. T. 1, no. 4. рр. 76–88.
7. Melnikov D.V., Chzho Tu.A., Okar M. Algoritm statisticheskogo analiza linejnyh ne-stacionarnyh sistem so sluchajnymi parametrami s ispolzovaniem apparata matrich-nyh operatorov i teorii chuvstvitelnosti // Nauchnoe obozrenie. 2015. no. 20. рр. 185–190.
8. Min Ch.Tu., Melnikov D.V., Algoritm stohasticheskogo issledovanija nelinejnyh si-stem upravlenija s ispolzovaniem proekcionno–matrichnogo apparata // Nauchnoe obo-zrenie. 2014. no. 10–1. рр. 87–93.
9. Mishin G.T. Sovremennaja analogovaja mikrojelektronika: teorija i praktika. M.: Iz-datelstvo «Radiotehnika», 2007. 208 р.
10. Okar Min, Melnikov D.V. Algoritm rascheta nelinejnyh sistem upravlenija proek-cionno-matrichnym metodom // Inzhenernyj zhurnal: nauka i innovacii. 2014 no. 12 (36). рр. 17; http://engjournal.ru/catalog/it/asu/1268.html.
11. Pupkov K.A., Egupov N.D., Kolesnikov L.V., Lukashenko Ju.L., Melnikov D.V., V.M. Rybin, Trofimov A.I. Matrichnye metody rascheta i proektirovanija slozhnyh sistem avtomaticheskogo upravlenija dlja inzhenerov / Pod red. K.A. Pupkova i N.D. Egupova. M.: Izdatelstvo MGTU im. N.Je. Baumana, 2007. 664р.
12. Pupkov K.A., Egupov N.D., Kolesnikov L.V., Melnikov D.V., Trofimov A.I. Vysoko-tochnye sistemy samonavedenija. Raschet i proektirovanie. Vychislitelnyj jeksperi-ment / Pod. red. Pupkova K.A., Egupova N.D. M.: Fizmatlit, 2011. 512 р.
13. Chzho Tu Aung, Melnikov D.V. Algoritm issledovanija nelinejnyh sistem avtomati-cheskogo upravlenija v stohasticheskih rezhimah // Inzhenernyj zhurnal: nauka i innovacii. 2014. no. 4 (28). рр. 3; http://engjournal.ru/catalog/it/asu/1270.html.

В настоящее время в области анализа и проектирования устройств аналоговой электроники (устройства отображения информации, системы автоматического управления, аналоговые фильтры) наблюдаются попытки перехода на уровень выпуска системных микросхем, однако отсутствие единого общепризнанного математического аппарата решения задач анализа и синтеза сдерживают развитие данного направления. Выход из создавшегося положения наметился в работе [2], где за естественнонаучную базу аналоговой микроэлектроники принята математическая теория систем дифференциальных уравнений (ДУ) в обобщенной форме. Переход к естественнонаучным представлениям позволяет сформировать новый элементный базис. Необходимым шагом перехода на новую элементную базу является замена модели радиоэлектронного устройства в виде принципиальной электрической схемы на равноценную модель в виде структурной схемы, которая должна состоять только из элементарных звеньев предлагаемой элементной базы [2, 9]. Важным шагом перехода на предлагаемую элементную базу является обоснование процедуры синтеза электронного устройства: выбор структуры дифференциального уравнения; расчет его коэффициентов и сборка синтезируемого устройства. Всё это, вместе взятое, определяет новый подход к анализу и синтезу электронных устройств на элементарных звеньях. На основе вновь введенной элементной базы возможно сформировать как МАБИС (матричные аналоговые большие интегральные схемы), так и ПАИС (программируемые аналоговые интегральные схемы). Такая возможность обусловлена малочисленностью элементов в базе и позволяет спроектировать уникальный чип для данной МАБИС, или единый перепрограммируемый элемент для ПАИС [2].

Однако те большие возможности, которые открываются введением предлагаемой элементной базы, требуют обширной предварительной работы, связанной с исследованием особенностей проектирования различных типов аналоговых устройств. В настоящей статье предлагается алгоритм синтеза таких широко применяющихся устройств, как фильтры. Для синтеза фильтров предлагается использовать аппарат матричных операторов, который хорошо зарекомендовал себя для решения задач управления для класса линейных стационарных и нестационарных систем [1, 3, 7], а в последнее время нашел развитие и на класс нелинейных систем [4, 5, 6, 10, 13]. Рассмотрим алгоритм решения задачи синтеза оптимального фильтра на конкретном примере. Предположим, что система задана следующей структурной схемой (рис. 1).

melnik1.wmf

Рис. 1. Структурная схема системы

Динамика объекта описывается уравнением

mel01.wmf (1)

где

mel02.wmf

mel03.wmf mel04.wmf

mel05.wmf

На вход системы поступает аддитивная смесь центрированных не коррелированных сигналов m(t) (полезный сигнал) и n(t) (помеха) с известными корреляционными функциями

mel06.wmf (2)

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

Этап 1. Выберем в качестве ортонормированного базиса (ОНБ) ортонормированные на промежутке [0, 6] полиномы Лежандра mel08.wmf

Этап 2. Проведем статистический анализ системы без фильтра. Структурная схема системы в операторной форме представлена на рис. 2 (mel09.wmf – спектральные характеристики входного сигнала, сигнала ошибки, управления и выходного сигнала соответственно; Аф, А0 – матричные операторы фильтра и объекта).

melnik2.wmf

Рис. 2. Структурная схема системы в операторном виде

Задав размер базиса l = 12, одним из способов, изложенных в [11, 12], найдем матричный оператор объекта А0 в выбранном базисе. Далее определяем матричный оператор всей системы без регулятора: mel10.wmf Применяя двумерное преобразование Фурье, рассчитываем спектральные характеристики корреляционных функций полезного сигнала и помехи:

mel11.wmf

mel12.wmf

По следующим зависимостям определяем спектральные характеристики корреляционной функции выходного сигнала, обусловленные действием: полезного сигнала и помехи – mel13.wmf полезного сигнала – mel14.wmf помехи – mel15.wmf.

При синтезе фильтра спектральную характеристику mel16.wmf (обусловленную действием только полезного сигнала) будем использовать в качестве желаемой (эталонной). По формуле

mel17.wmf

mel18.wmf

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

Найдем автокорреляционные функции, дисперсии и среднеквадратическое отклонение (СКО) выходного сигнала, обусловленные действием:

– полезного сигнала и помехи

mel19.wmf

mel20.wmf

– полезного сигнала

mel21.wmf

mel22.wmf

– помехи mel23.wmf

mel24.wmf

Графики СКО выходного сигнала, обусловленные совместным действием полезного сигнала и помехи, а также действием отдельно помехи и полезного сигнала, представлены на рис. 3.

melnik3.wmf

Рис. 3. Графики СКО выходного сигнала системы без фильтра

Этап 3. Расчет оптимального матричного оператора фильтра. Вычислим матричный оператор фильтра относительно двух критериев: критерия близости корреляционных функций и критерия близости дисперсий. Матричный оператор системы с фильтром имеет вид

mel25.wmf. (3)

Оптимальные матричные операторы найдем из следующих зависимостей:

- критерий близости корреляционных функций:

mel26.wmf (4)

где

mel27.wmf

- критерий близости дисперсий:

mel28.wmf (5)

где

mel29a.wmf

mel29b.wmf

mel30.wmf

Выбрав евклидову норму матриц (4), (5), получим следующие целевые функции, подлежащие минимизации относительно элементов матрицы Аф:

mel31.wmf

mel32.wmf

Положив в качестве начального приближения Аф = I, (I – единичная матрица) методом Гаусса – Ньютона вычисляем оптимальные матричные операторы фильтра mel33.wmf, mel34.wmf, при этом mel35.wmf mel36.wmf

Этап 4. Построение математической модели оптимального фильтра в форме дифференциального уравнения по известным матричным операторам mel37.wmf mel38.wmf Имеется несколько подходов к решению этой задачи, один их них основан на использовании понятия порождающих функций [11]. В случае если порядок ДУ невысокий, можно предложить следующий способ. Зададим следующую структуру ДУ, описывающее динамику фильтра:

mel39.wmf (6)

Определим коэффициенты mel40.wmf mel41.wmf в виде разложения по ОНБ так, чтобы матричный оператор уравнения (6) Аф как можно точнее совпадал с mel42.wmf в одном случае и с mel43.wmf в другом. Матричные операторы правой mel44.wmf и левой mel45.wmf частей уравнения (6) можно представить следующим образом (см. переход о ДУ к соответствующему интегральному уравнению в [11]):

mel46.wmf

где mel47.wmf – спектральная характеристика mel48.wmf

mel49.wmf – спектральная характеристика mel50.wmf элементы матрицы mel51.wmf согласно вычисляются следующим образом:

mel52.wmf,

а ее структура имеет размер (12*122):

mel54.wmf

Тогда матричный оператор фильтра определяется следующим образом:

mel55.wmf (7)

Оптимальные параметры mel56.wmf найдем исходя из следующих зависимостей:

– критерий близости корреляционных функций:

mel57.wmf (8)

– критерий близости дисперсий:

mel58.wmf (9)

Задавшись евклидовой нормой для матриц (8), (9), получим следующие целевые функции, подлежащие минимизации:

mel59.wmf

mel60.wmf

где mel61.wmf – совокупность искомых параметров.

В результате минимизации целевых функций методом Гаусса – Ньютона были найдены: mel62.wmf – спектральные характеристики соответствующих коэффициентов фильтра, найденные по критерию (8); mel63.wmf – спектральные характеристики, найденные по критерию (9).

Целевые функции при оптимальных параметрах принимают следующие значения:

mel64.wmf mel65.wmf

mel66.wmf mel67.wmf

Этап 5. Анализ системы с найденными фильтрами. Матричные операторы всей системы:

mel68.wmf

mel69.wmf

Спектральные характеристики выходного сигнала, обусловленные воздействием полезного сигнала и помехи:

mel70.wmf

mel71.wmf

Автокорреляционные функции, дисперсии и СКО выходного сигнала, обусловленные действием полезного сигнала и помехи:

mel72.wmf

mel73.wmf

mel74.wmf

mel75.wmf

На рис. 4 представлены графики СКО выходного сигнала с фильтром, построенным на основе равенства корреляционных функций выходного и эталонного процессов (фильтр 1) и на основе равенства дисперсий выходного и эталонного процессов (фильтр 2). В случае получения результата с недостаточной точностью, необходимо на этапе 4 задать более сложную структуру фильтра.

melnik4.wmf

Рис. 4. Графики СКО выходного сигнала системы

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

Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований и Правительства Калужской области (гранты № 14-48-03013, № 16-41-400701).