CC BY
33
УДК 519.673
ПРИМЕНЕНИЕ ЯЧЕЕК ОДНОРОДНОЙ СТРУКТУРЫ ДЛЯ ВЫЧИСЛЕНИЯ НЕПРЕРЫВНОГО ВЕЙВЛЕТ-ПРЕОБРАЗОВАНИЯ
А.А. Хамухин
Томский политехнический университет E-mail: [email protected]
Рассмотрено применение ячеек однородной структуры для параллельного синхронного вычисления непрерывного вейвлет-преобразования. Показана возможность реализации функционального наполнения ячейки независимо от вида используемого вейвлета. Представлены результаты моделирования матрицы ячеек размером 16х64для тестового входного сигнала, состоящего из двух гармонических колебаний разной частоты.
Ключевые слова:
Ячейка однородной вычислительной структуры, параллельные вычисления, непрерывное вейвлет-преобразование, математическая модель.
Key words:
Cell of homogeneous computing structure, parallel computation, continuous wavelet transform, mathematical model.
Введение
Методология однородных вычислительных структур (ОВС), разработанная в нашей стране Э.В. Евреиновым с 60-х по 80-е гг. прошлого века [1], при существенном увеличении количества транзисторов на кристалле становится вновь предметом внимания разработчиков [2, 3]. Если в качестве базовой вычислительной ячейки используется микропроцессор, то это позволяет легко перенастраивать ОВС на решение разных задач, но не дает такого увеличения производительности, которое можно ожидать от аппаратной реализации ячейки.
Можно выделить достаточно широкий класс задач, для которых затраты на аппаратную реализацию будут оправданы, если требуемое увеличение производительности будет достигнуто. Так, для решения дифференциальных уравнений в частных производных разработана схема вычислительной ячейки [4], показана возможность реконфигурирования ОВС, состоящей из таких ячеек [5], с помощью математического моделирования продемонстрировано увеличение производительности на порядок [6].
Объект и цель исследования
Непрерывное вейвлет-преобразование (НВП) отличается от решения дифференциальных уравнений в частных производных, поэтому прямое использование разработанных ячеек однородной структуры невозможно. Но, используя методику разработки ячейки для решения дифференциальных уравнений в частных производных, можно аналогично разработать и ячейку однородной структуры для НВП, что и является объектом нашего исследования.
Целью исследования является найти математически обоснованный переход от известного математического описания непрерывного вейвлет-преобразования к системе уравнений, содержащих элементы, известные как математические описания функциональных блоков ячеек однородной струк-
туры (интегратор, сумматор, масштабный интегратор и др.).
Пример построения объекта
Непрерывное вейвлет-преобразование при вычислении на цифровой технике требует дискретизации и потому не является строго непрерывным. Но этот термин сохраняется, чтобы отличать его от дискретного вейвлет-преобразования, которое основано на децимации сигнала после каждой пары фильтров высокой и низкой частоты [7]. Иногда для большей точности НВП называют дискретным вейвлет-преобразованием с произвольным шагом дискретизации [8].
Непрерывное вейвлет-преобразование произвольного входного сигнала £(/) описывается уравнением:
W (a, b) = ^= J 5 (t)¥
(1)
где Ж(а,Ь) - искомая функция; / - время; Ь - сдвиг по времени; а - масштаб по времени (обратно проч ^ - Ь
порциональный частоте сигнала);
- не-
которая функция с определенными свойствами (вейвлет).
Для оцифрованного с равномерным шагом сигнала, имеющего конечные пределы по времени, интеграл в уравнении (1) аппроксимируется конечной суммой:
П Л г \
w a, ьк )^ (t. )y ■\jaj ‘=0
At,
(2)
j У
где I, ], к - индексы по времени I, по масштабу а, по сдвигу по времени Ь; п - количество шагов по времени; А/ - шаг по времени; $(/,.) - оцифрованный входной сигнал $(/) (данные).
Если принять начальный момент времени нулевым, а сдвиг по времени равным шагу по времени (что соответствует наибольшей разрешающей спо-
а
собности НВП), то уравнение (2) можно записать в виде системы уравнений:
W(й =
At
&
a,
V i
S (0)T +S (2At)T
+ S(At)T
a
V i /
a
V i
+... + S( nAt) T( —)
W (a b^) =
At
4a,
S (0)T
a
V i
+ S(At)T
a
V i
+S (2At )T
+...
w a, Ьз)=
At
S (0)T
a,
\ i J
+ S(At)T
a,
\ i J
+S (2At )T
+...
W (</ bk) =
S(0)T (k -1) At + S(At)T k At +
At l_ a, J l_ a, J
+S (2At)T
(k + 1)At
(3)
в которой количество уравнений равно количеству сдвигов по времени. Очевидно, что все выражения, на которые должны умножаться данные 5(/А/), могут быть вычислены заранее и занесены в виде коэффициентов для блока умножения создаваемой ячейки. Систему уравнений (3) можно записать в виде итерационной формулы по индексу г.
W (a,, bk ) =
= W‘ (a,, bk ) + S (iAt) <
At
T
(i + k -1)At
или
где
Wi+1 (a bk) = Wi (a , bk) + S(iAt)F>,
IVa,
(i + k -1)At
(4)
(5)
(6)
Формула (5) легко реализуется с помощью цифрового и масштабного интеграторов, применяемых в ячейках однородной структуры. Коэффициенты (6) определяются только видом выбранного материнского вейвлета, не зависят от результатов
текущих вычислений, поэтому могут быть заранее подготовлены и записаны в регистры ячейки.
Следовательно, объединяя упомянутые функциональные блоки в одну вычислительную ячейку, можно вычислить одну точку двумерного массива дискретизированного вейвлет-преобразования, соответствующего точному преобразованию с погрешностью, зависящей от шага дискретизации At.
Объединяя далее такие ячейки в однородную вычислительную структуру (матрицу) размером JxK, где J - максимальное количество масштабов по времени, K- максимальное количество сдвигов по времени, можно получать дискретизированное непрерывное вейвлет-преобразования данных в полном объеме.
Еще одна составляющая погрешности вычисления НВП определяется разрядностью представления всех указанных переменных в ячейках однородной структуры, которая легко может масштабироваться в зависимости от поставленной задачи и требуемой точности расчета.
Таким образом, вычисление непрерывного вейвлет-преобразования дискретизированного входного сигнала распределяется между синхронно параллельно работающими ячейками однородной структуры. На вход каждой ячейки синхронно подаются оцифрованные данные, подлежащие НВП, которые умножаются на заранее подготовленные коэффициенты (6), хранящиеся в каждой ячейке. Результат накапливается в цифровом интеграторе с учетом (4). Каждая ячейка матрицы дает одно значение дискретизированного непрерывного вейвлет-преобразования W(aj,bk), поэтому вся однородная структура ячеек размером JxK даст полный спектр значений НВП с погрешностью, зависящей от шага дискретизации и разрядности регистров ячеек. Параллельность и синхронность работы всех ячеек однородной структуры являются залогом достижения максимальной производительности при прочих равных условиях.
Результаты моделирования
Математическая модель (1)-(6) была реализована с помощью MS Excel в таблице размером 16x64, что соответствовало однородной вычислительной структуре (матрице) того же размера (JxK). В каждой клетке таблицы MS Excel производилось независимое от соседних клеток вычисление одного значение НВП по формуле (5). Таким образом, моделировалась параллельность и синхронность работы ячеек ОВС. Для коэффициентов (6) использовался материнский вейвлет, известный под названием «мексиканская шляпа» (рис. 1), который является второй производной функции Гаусса и описывается уравнением
T( х) =
1
•Jiña''
f х2^ 1 -О-
a
(7)
где х=(—Ь)/а; а - параметр, позволяющий изменять масштаб вейвлета в соответствии с видом входного сигнала.
Рис. 2. Вид дискретизированного тестового входного сигнала (8): Д=А=/; со=0,627; М=1
В качестве тестового входного сигнала использовалась часто применяемая для тестирования вейвлет-преобразований функция
Б (?) = Л^тШ?) + Л2$т(5ш), (8)
представляющая собой сумму двух гармонических колебаний с частотами ш и 5ш Вид дискретизированного входного сигнала, для которого выполнялось НВП с помощью описанной выше модели, представлен на рис. 2.
Параметры вейвлета при моделировании НВП изменялись в пределах а от 1 до 16, Ь от 1 до 57. Количество отсчетов тестового входного сигнала п=64. Полученное в результате компьютерного моделирования однородной вычислительной структуры непрерывное вейвлет-преобразование Ж(а,Ь) входного сигнала £(?) представлено на рис. 3 в виде поверхности в трехмерном пространстве и на рис. 4 в виде проекции линий равного уровня на плоскость (а, Ь).
На рис. 5 представлены сечения полученного НВП на разных временных масштабах, хорошо показывающие составляющие частоты тестового входного сигнала.
Исследования математической модели ячейки однородной структуры для вычисления НВП проводились для различных наборов параметров а, Ь, п, А1, А2, /, К, ш, а. От выбора соотношения этих параметров работоспособность созданной модели не зависела, но наглядность полученных результатов была разной. На рисунках, представленных выше, приведены результаты моделирования для одного набора параметров. По графикам видно полное совпадение расчетных данных с результатами, представленными в первоисточнике [7].
Выводы
Показано, что для параллельного вычисления непрерывного вейвлет-преобразования после соответствующей доработки возможно применение ячеек, аналогичных ячейкам однородной структуры для решения дифференциальных уравнений в частных производных.
Функциональное наполнение ячейки однородной структуры для непрерывного вейвлет-преобразования можно реализовать независимым от вида используемого вейвлета, который рассчитывается заранее для всего диапазона масштабов и сдвигов
Рис. 3. Результат моделирования непрерывного вейвлет-преобразования тестового входного сигнала (8) с помощью формул ОВС (4)-(6) на вейвлете (7)
&^:§^!£58£!»ВЕ2ЙИИ7///ЛНШ||М1|1]|||1111|,11|,1ШПШ\\УК«К£№ЕевЗаг22г^;г!£:г^
«?5?5ййй5йкй3==я55йи0ии»//лшш11111итттумм8к«£в=55йяйвйй:й%:й;2!»
й&«ка»й5;Щ!г=:й==йЯ5йИ000«шш1111т11илт\И!'гкйй5а==ййг:й'ййаййг?й;й'й
я^8й«&йй^^=:==-==г:=га«й^ии,ллш111к\л:|Х«^«^^ь=;==-=5=й==ййа8;йг«я^
в«¥®еее«5;5=зй====£зеззавщ!*и1111,и*ей'К5е^эй====55ййеаеззайия
ка$а5£г£58&ЗЁ»ггггг8!!:11!1!1,и1111тгккЕ5г^г^гз:1:Еагггггг99аа1
■Е*=^;!»якяввг㧧^=^зй1яя1,н«иипбЁЕ§§§§5г№!!шга9яВ11
ЯЖ1Вё=^ШШаН№^:^^за1;8&;8М1ПК№88Я1!£^§§^НП1!Ш11Г1!1Ш^ё€а!11ЕЯ
!6»№Иг!^ШЖ4Ва№Ш^13^ШМбаКШШВ»П№>г!£€^:Ш1115ВП1Н№!Шё:а«П17/<вга
»Ш№111111П№№№ЮаИ№7Л|М!1ИВВМ)»1Ш1иЖПГ1111111111И111И,ГЛ»№№7>Я0№№ЮЮ!вВ1И!№(М№ГЛаИ1И[М№«!а
за№£=аа«££^У1Б^гаи1!£г=г>!П^=гзш1£=зв!1К£г=^!аг1г>£=гг1^=
Рис. 4. Проекция непрерывного вейвлет-преобразования Ща,Ь), полученного с помощью формул ОВС (4)-(6) на плоскость (а,Ь)
по времени и в виде коэффициентов хранится и используется в регистрах соответствующих ячеек.
Моделирование матрицы ячеек размером 16х64 показало соответствие известным ре-
зультатам непрерывного вейвлет-преобразования тестового входного сигнала, состоящего из двух гармонических колебаний разной частоты.
Рис. 5. Сечения полученного непрерывного вейвлет-преобразования тестового входного сигнала на разных временных масштабах
СПИСОК ЛИТЕРАТУРЫ
1. Евреинов Э.В. Однородные вычислительные системы, структуры и среды. - М.: Радио и связь, 1981. - 208 с.
2. Каляев И.А., Левин И.И., Семерников Е.А., Шмойлов В.И. Реконфигурируемые мультиконвейерные вычислительные структуры. - Ростов на Дону: ЮНЦ РАН, 2008. - 393 с. 2010. URL: http://parallel.ru/FPGA/papers/rmvs.pdf (дата обращения: 21.08.2010).
3. Giefers H., Platzner M. A Many-Core Implementation Based on the Reconfigurable Mesh Model // IEEE Xplore DIGITAL LIBRARY. 2010. URL: http://ieeexplore.ieee.org/Xplore/defdeny.jsp? url=http://ieeexplore.ieee.org/stamp/stamp.jsp%3Ftp%3D%26ar-number%3D4380623&denyReason = -134&arnum-ber=4380623&productsMatched=null (дата обращения: 21.08.2010).
4. Ячейка однородной структуры для решения дифференциальных уравнений в частных производных: пат. 2359322 Рос. Федерация. № 2007141832/09; заявл. 12.11.2007; опубл. 20.06.2009, Бюл. № 17. - 6 с.
5. Хамухин А.А. Реконфигурирование однородной вычислительной структуры с непрограммируемыми ячейками для решения дифференциальных уравнений в частных производных // Известия Томского политехнического университета. - 2010. -Т. 316. - № 5. - С. 68-72.
6. Хамухин А.А. Ячеечная модель устройства для решения дифференциальных уравнений в частных производных // Известия Томского политехнического университета. - 2010. -Т 316. - № 5. - С. 62-67.
7. Яковлев А.Н. Введение в вейвлет-преобразования. - Новосибирск: Изд-во НГТУ, 2003. - 104 с.
8. Способ и устройство быстрого вычисления дискретного вейвлет-преобразования сигнала с произвольным шагом дискретизации масштабных коэффициентов: пат. 2246132 Рос. Федерация. № 2003100794/09; заявл. 09.01.2003; опубл. 10.02.2005, Бюл. № 4. - 20 с.
Поступила 01.09.2010 г.
