CC BY
87
УДК 621.396 + 681.5.015.52
М. А. Щербаков, В. В. Сазонов, Р. Н. Исянов
МАТРИЧНАЯ РЕАЛИЗАЦИЯ ДВУХМЕРНЫХ ДИСКРЕТНЫХ ФИЛЬТРОВ ВОЛЬТЕРРА
Аннотация.
Актуальность и цели. Объектом исследования являются многомерные дискретные фильтры Вольтерра (ФВ, полиномиальные фильтры), определяемые отрезком дискретного функционального ряда. Предметом исследования являются способы эффективной реализация ФВ, основанные на их представлении в матричном виде и выполнение процедуры декомпозиции. Целью работы является разработка способов эффективной реализации двухмерных дискретных ФВ за счет выполнения распараллеливания и унификации вычислительных операций.
Материалы и методы. Для описания процесса многомерной дискретной фильтрации использовано матричное представление ФВ. Введено понятие блочной свертки, позволяющей свести вычисление матричных нелинейных сверток к последовательности вычислений обычных двухмерных сверток, а также оценить объем необходимых вычислительных затрат.
Результаты. Предложен подход к реализации двухмерных дискретных ФВ, основанный на декомпозиции матричного представления нелинейных двухмерных сверток. Данный подход позволяет свести задачу реализации нелинейных ФВ к последовательно-параллельному выполнению обычных двумерных линейных сверток и может стать основой разработки параллельных вычислительных структур для реализации быстродействующих алгоритмов нелинейной фильтрации на основе использования матрицы параллельных процессоров.
Выводы. Метод реализации двухмерных дискретных ФВ, основанный на выполнении процедуры декомпозиции матричного представления нелинейных сверток, позволяет использовать для реализации ФВ известные методы и алгоритмы многомерной линейной фильтрации и дать им наглядную интерпретацию. Предложенный метод реализации ФВ за счет использования параллельно-последовательного выполнения процедуры вычислений обеспечивает возможность их реализации на основе использования структур процессоров систолического типа.
Ключевые слова: цифровая обработка сигналов и изображений, нелинейная фильтрация, фильтры Вольтерра, полиномиальная фильтрация, параллельные вычисления.
M. A. Shcherbakov, V. V. Sazonov, R. N. Isyanov
MATRIX REALIZATION OF TWO-DIMENSIONAL DISCRETE VOLTERRA FILTERS
Abstract.
Background. The object of the research is multidimensional discrete Volterra filters (VF, polynomial filters), determined by a section of a discrete functional series. The subject of the research is the methods of VF effective realization, based on representation thereof in matrix form and decomposition procedure perfomance. The article is aimed at development of ways of effective realization of two-dimensional discrete VF due to paralleling and unification of computing operations.
Materials and methods. To describe the process of multidimensional discrete filtration the authors used matrix representation of VF. The researchers introduced a notion of sectional convultion allowing to reduce calculation of matrix non-linear convultions to a sequence of calculations of regular two-dimensional convultions, and to estimate the volume of necessary calculations.
Results. The authors suggest an approach to realization of two-dimensional discrete VF, based on decomposition of matrix representation on non-linear twodimensional convultions. The given approach allows to reduce the problem of nonlinear VF realization to series-parallel execution of regular two-dimensional con-vultions and may serve as a foundation for development of parallel computational structures intedned for realization of high-speed algorithms of non-linear filtration through parallel process matrix application.
Conclusions. The method of realization of two-dimensional discrete VF, based on execution of decomposition of the matrix representation of non-linear con-vultions, allows to use the known methods and algorithms of multidimensional filtration for realization of VF and to give clear interpretation thereof. The suggested method of VF realization through using series-parallel execution of calculations provides an opportunity of realization thereof on the basis of application of structures of systolic processes.
Key words: digital signal and image processing, non-linear filtering, filters Volterra, polynomial filtering, parallel calculations.
Введение
С целью расширения спектра задач, решаемых средствами цифровой обработки сигналов и изображений, и преодоления ограничений, присущих методам линейной фильтрации, в настоящее время активно внедряются методы нелинейной фильтрации. Наиболее известными классами нелинейных фильтров являются: гомоморфные фильтры; морфологические фильтры; фильтры, основанные на порядковых статистиках, и их разновидности: L-, R-, M-фильтры, медианные фильтры; расширенные фильтры Калмана; нейронные фильтры; нелинейные фильтры Вольтерра [1].
Каждый из перечисленных классов имеет свои преимущества и область применения. Некоторые направления, такие, например, как фильтрация Кал-мана, гомоморфная фильтрация, имеют достаточно долгую историю. Другие направления, в частности фильтры Вольтерра, далее обозначаемые как ФВ [2, 3], появились сравнительно недавно и активно применяются в настоящее время в различных областях, таких как обработка сигналов, адаптивная эхо-компенсация, идентификация систем и управление и др.
Исследование данного класса фильтров, называемых также полиномиальными фильтрами [2], позволяет перенести известные методы многомерной линейной фильтрации на нелинейный случай. В то же время их реализация связана с существенным возрастание вычислительных затрат, что делает актуальной задачу разработки подходов, позволяющих использовать современные методы параллельных вычислений.
1. Математическое описание многомерных дискретных ФВ
В общем случае дискретный ФВ размерности r и порядка M определяется конечным дискретным рядом Вольтерра (функциональным полиномом) вида [3]
M M m
y(n) = h) + Z ym (n) = h0 + Z Z"Zhm (n1,...,nm)ПX(n -ni),
m=1
m=1 n
i=1
где hш(n1, пш) - многомерные импульсные характеристики (ядра) фильтра, зависящие от векторных аргументов п = \пі\...піг\.
Выходной сигнал у(п) таких фильтров представляет собой сумму составляющих, характеризующих нелинейности различного порядка: первая ^і(п) имеет вид линейной свертки, вторая у2(п) - квадратичной свертки и т.д. Составляющая уш(п) фильтра, определяемая нелинейной сверткой m-го порядка, является нелинейной относительно отсчетов входного сигнала, однако остается линейной по отношению к коэффициентам фильтра.
При m = 1 ядро hm(nl) представляет собой обычную импульсную характеристику многомерного линейного фильтра, в то время как при m = 2, ..., M ядра hm(n1, ..., пш) можно рассматривать как импульсные характеристики высших порядков, характеризующие нелинейные свойства многомерных полиномиальных фильтров \4\.
2. Матричное представление многомерных дискретных ФВ
Для реализации ФВ с конечной импульсной характеристикой удобным является эквивалентное матричное представление \5\. Для перехода к матричной форме записи введем в рассмотрение вектор входного сигнала:
xП =[хп(0) xn(1) ... xn^-1)].
Тогда векторное представление составляющей уш(п) m-го порядка может быть записано в виде
уш (п)=^ 4^
где вектор hш коэффициентов фильтра содержит лексикографически упорядоченные значения нелинейной импульсной характеристики hш(n1, ..., пш).
Для двухмерных полиномиальных фильтров вместо векторной формы можно использовать матричную. Воспользуемся простым отображением для преобразования двухмерной индексации точек опорной области в одномерную. Для маски 3 X 3, например, это преобразование будет выглядеть так:
"(0,0) (0,1) (0,2)' "0 1 2"
*2 = (1,0) (1,1) (1,2) 3 4 5
.(2,0) (2,1) (2,2)_ 6 7 8
Для такого представления опорной области входной сигнал может быть записан в виде N X ^-матрицы:
Xn =
Xn (0) Xn (1) Xn (2)
Xn (3) Xn (4) Xn (5)
Xn (6) Xn (7) Xn (8)
(1)
Здесь входной сигнал не упорядочивается в вектор, а представляется в более естественном для двухмерной фильтрации матричном виде.
Для формирования произведений отсчетов входного сигнала также воспользуемся кронекеровской степенью матрицы Хп. Образуем т-упорядо-ченную, по Кронекеру, матрицу Нт, содержащую элементы ядра т-го порядка, расположенные в соответствии с индексацией произведений отсчетов хп(/1), ..., хп(/т), содержащихся в матрице Хп(т). Будем называть матрицу Нт т-упорядоченной, по Кронекеру. В частности, для т = 2 и Хп в виде (1) данная матрица будет выглядеть следующим образом:
*2(0,0) *2(0,3) *2(0,6) *2(0,1) *2(0,4) *2(0,7) *2(0,2) *2(0,5) *2(0,8) ! *2(1,0) | *2(1,3) і *2(1,6) *2(1,1) *2(1,4) *2(1,7) *2(1,2) j *2(1,5) j *2(1,8) j *2(2,0) *2(2,3) *2(2,6) *2(2,1) *2(2,4) *2(2,7) *2(2,2) *2(2,5) *2(2,8)
Н2 = *2(3,0) *2(3,3) *2(3,6) *2(3,1) *2(3,4) *2(3,7) *2(3,2) *2(3,5) *2(3,8) [*2(4,0) ! *2(4,3) ! *2(4,6) *2(4,1) *2(4,4) *2(4,7) *2(4,2) j *2 (4,5) ! *2(4,8) і *2(5,0) *2(5,3) *2(5,6) *2(5,1) *2(5,4) *2(5,7) *2(5,2) *2(5,5) *2(5,8)
*2(6,0) *2(6,3) *2(6,6) *2(6,1) *2(6,4) *2(6,7) *2(6,2) *2(6,5) *2(6,8) ! *2(7,0) ! *2(7,3) *2(7,6) *2(7,1) *2(7,4) *2(7,7) *2(7,2) ! *2(7,5) j *2(7,8) *2(8,0) *2(8,3) *2(8,6) *2(8,1) *2(8,4) *2(8,7) *2(8,2) *2(8,5) *2(8,8)
?_2_! _H_2 і_H2
H3 1 H4 1 H5
H 2 I H 2 I H 2
h2 ! h2 ! h2
і
і
В общем случае составляющая _ут(п) двухмерного полиномиального фильтра может быть представлена в следующей матричной форме:
У г.
» = ^{нm о Xim)} ,
(2)
где ° означает произведение Адамара (поэлементное произведение матриц); Е(Л} - сумма всех элементов матрицы А.
3. Реализация двухмерных дискретных ФВ на основе матричных сверток
Взаимосвязь между полиномиальными и многомерными линейными фильтрами позволяет синтезировать эффективные алгоритмы на базе уже известных алгоритмов цифровой линейной фильтрации.
Декомпозицию двухмерного дискретного ФВ произвольной размерности т можно получить на основе его матричного представления. Для наглядности положим размер маски N = 3 и запишем (2) в блочном виде:
Ут (n) = 2
н0 ' н1 ' н2
_т_ і _ _ т_ і __т
НП H^ Н
пт _1 пт пт
н6 і н7 і Но
пт І пт I пт
"r„ (0)Xim-') 1 r'XS”-'» Г2 Xim-1)'
r3Xim-1) U Xir'* Г6 Xim-'» |r7 Xi"-'» Г5 X<”M) I Xl,m-11
o' I Xn (1) 'J y°m -1(n) ym-1(n) ym-1(n) 8
xn (3) Xn (4) Xn (5) О ym -1(n) 4 , \ ym-1(n) ym-1(n) > = Z Xiym -1(n)
_ Xn (6) Xn (7) Xn (8) _ У6т -1(n) УПт -1(n) ym -1(n) i=0
Данное выражение определяет выходной сигнал ут(п) ФВ т-го порядка через выходные сигналы у'т-1(п) полиномиальных фильтров (т - 1)-го порядка, равные
Шо) h|(2)" Xn (0) Xn (1) Xn (2)
h1(3) h1(4) hi(5) О Xn (3) Xn (4) Xn(5)
_h1(6) h1(7) ВД_ _Xn (6) Xn (7) Xn (8) _
ут-1(п)=х{нт ° хпт-1)}.
В частности, для линейного случая выражение (3) приобретает вид У1(п) = Х{Н1 о Хп }
и представляет собой матричную форму записи обычной двухмерной свертки.
Для описания матричного алгоритма реализации нелинейного двухмерного ФВ введем понятие блочной свертки.
Определение. Пусть имеются матрицы Н и X с размерностями соответственно N х N2 и Ы\ х М2 . Представим матрицу Н в виде блочном матрицы размерностью Ь\ х Ь2 вида
H =
H11 !
H Li !
H
1L2
H
L1L2
где каждый блок имеет размерность М\ х М2. Тогда под блочной сверткой будем понимать матрицу У = {уу} размерности Ь\ х Ь2, где каждый элемент представляет собой линейную двухмерную линейную свертку вида
Уч =^{Нч ох} •
Для обозначения блочной свертки будем использовать обозначение
У = Н □ Х.
Выполнение такой свертки, согласно определению, требует вычисления Ь\ х Ь2 обычных двухмерных сверток.
С учетом введенной операции блочной свертки и уравнения (3) декомпозиции процедуру вычисления выходного сигнала ут(п) двухмерного ФВ т-го порядка можно представить как последовательность вычисления т блочных сверток.
Для наглядности продемонстрируем данную процедуру на примере двухмерного ФВ третьего порядка с маской 3 х 3, выходной сигнал которого определяется выражением вида
Уз(п) = z{h оXi3)} ,
где Н3 - матрица коэффициентов фильтра размерностью 27 х 27; Хп - матрица отсчетов входного сигнала фильтра размерности 3 х 3.
В этом случае вычисление выходного сигнала ,у3(п) фильтра происходит путем последовательного вычисления следующих трех сверток:
1) блочной свертки матриц Н3 и Х, в результате которой формируется матрица У3[1] = Н3 п Х размерностью 9 х 9;
2) блочной свертки матриц У3[1] и X, в результате которой формируется матрица У3[2] = У3[1] п Х размерностью 3 х 3;
3) обычной свертки матриц У3[2] и X, в результате которой формируется скалярная величина ,у3[3] = .у3(п), представляющая собой выходной сигнал ФВ
Согласно приведенной процедуре для вычисления выходного сигнала ФВ третьего порядка в точке п требуется выполнить 9 х 9 + 3 х 3 + 1 = 91 операцию линейной двухмерной свертки с маской размерности 3 х 3. В общем случае реализация фильтра т-го порядка с маской N х N потребует выполнения
операций двухмерной свертки.
Таким образом, предложенная процедура матричной декомпозиции позволяет реализовать ФВ произвольного порядка путем последовательнопараллельного выполнения операций двухмерной линейной свертки, которые, в свою очередь, допускают высокоэффективную реализацию в виде структур систолического типа, матричных и волновых процессоров [6].
Предложенный метод реализации двухмерных дискретных ФВ, основанный на выполнении процедуры декомпозиции матричного представления нелинейных сверток, позволяет свести задачу реализации нелинейных ФВ к последовательно-параллельному выполнению обычных двумерных линейных сверток. Это делает возможным использовать при реализации ФВ известные методы и алгоритмы многомерной линейной фильтрации и дать наглядную интерпретацию процессу нелинейной фильтрации через введенное понятие блочной свертки. Предложенный подход может стать основой разработки параллельных вычислительных структур для реализации быстродействующих алгоритмов нелинейной фильтрации на основе использования методов параллельных вычислений на базе структур систолического типа.
Список литературы
1. Pit as, I. Nonlinear digital filters: principles and applications / I. Pitas, A. N. Venetsa-nopoulos. - Kluver Academic Publishers, 1990. - 391 p.
2. Mathews, V. J. Polynomial signal processing / V. J. Mathews, G. L. Sicuranza. -John Wiley & Sons, Inc., 2000. - 452 р.
3. Щерб аков, М. А. Цифровая полиномиальная фильтрация: теория и приложение / М. А. Щербаков. - Пенза : Изд-во Пенз. гос. техн. ун-та, 1997. - 246 с.
в точке n.
N2 -1
1=1
Заключение
4. Щерб аков, М. А. Метод синтеза цифровых полиномиальных фильтров с помощью базисных частотных функций / М. А. Щербаков, С. В. Сорокин // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2007. -№ 4. - С. 74-86.
5. Щерб аков, М. А. Матричное представление цифровых полиномиальных фильтров. / М. А. Щербаков // Надежность и качество : тр. Междунар. симпозиума. - Пенза : Изд-во Пенз. гос. ун-та, 2002. - С. 469-472.
6. СБИС для распознавания образов и обработки изображений / под. ред. К. Фу. -М. : Мир, 1988. - 248 с.
References
1. Pitas I., Venetsa-nopoulos A. N. Nonlinear digital filters: principles and applications. Kluver Academic Publishers, 1990, 391 p.
2. Mathews V. J., Sicuranza G. L. Polynomial signal processing. John Wiley & Sons, Inc., 2000, 452 р.
3. Shcherbakov M. A. Tsifrovaya polinomial’naya fil'tratsiya: teoriya i prilozhenie [Digital polynomial filtering: theory and application]. Penza: Izd-vo Penz. gos. tekhn. un-ta, 1997, 246 р.
4. Shcherbakov M. A., Sorokin S. V. University proceedings. Volga region. Engineering sciences. 2007, no. 4, pр. 74-86.
5. Shcherbakov M. A. Nadezhnost' i kachestvo: tr. Mezhdunar. simpoziuma [Reliability and Quality: Proceedings of the International Symposium]. Penza: Izd-vo Penz. gos. un-ta, 2002, pp. 469-472.
6. SBIS dlya raspoznavaniya obrazov i obrabotki izobrazheniy [VLSI for Pattern Recognition and Image Processing]. Ed. K. Fu. Moscow: Mir, 1988, 248 р.
Щербаков Михаил Александрович
доктор технических наук, профессор, заведующий кафедрой автоматики и телемеханики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: mashcherbakov@yandex.ru
Сазонов Владимир Васильевич
кандидат технических наук, доцент, кафедра автоматики и телемеханики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: mashcherbakov@yandex.ru
Исянов Руслан Николаевич
кандидат технических наук, доцент, кафедра автоматики и телемеханики, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: mashcherbakov@yandex.ru
Shcherbakov Mikhail Aleksandrovich Doctor of engineering sciences, professor, head of sub-department of automation and remote control, Penza State University (40 Krasnaya street, Penza, Russia)
Sazonov Vladimir Vasil'evich Candidate of engineering sciences, associate professor, sub-department of automation and remote control, Penza State University (40 Krasnaya street, Penza, Russia)
Isyanov Ruslan Nikolaevich Candidate of engineering sciences, associate professor, sub-department of automation and remote control, Penza State University (40 Krasnaya street, Penza, Russia)
УДК 621.396 + 681.5.015.52 Щербаков, М. А.
Матричная реализация двухмерных дискретных фильтров Воль-терра / М. А. Щербаков, В. В. Сазонов, Р. Н. Исянов // Известия высших учебных заведений. Поволжский регион. Технические науки. - 2014. -№ 1 (29). - С. 64-71.
