Алгоритмы линейного вращения вектора численным интегрированием для реализации дискретного преобразования Фурье Текст научной статьи по специальности «Математика»

Зельманов Самуил Соломонович

Московский технический университет связи и информатики (Волго-Вятский филиал). E-mail: zelmanss@yandex.ru.

603011, г. Нижний Новгород, Менделеева, 15.

Тел.: 8312457505.

Zelmanov Samuil Solomonovich

Moscow Technical University of Communication and Information Sciens (Volgo-Vyatskiy Branch).

E-mail: zelmanss@yandex.ru.

15, Mendeleeva Street, Nizhny Novgorod, 603011, Russia.

Phone: 8312457505.

УДК 621.05.3

B.B. Владимиров, H.C. Звягинцев АЛГОРИТМЫ ЛИНЕЙНОГО ВРАЩЕНИЯ ВЕКТОРА ЧИСЛЕННЫМ ИНТЕГРИРОВАНИЕМ ДЛЯ РЕАЛИЗАЦИИ ДИСКРЕТНОГО ПРЕОБРАЗОВАНИЯ ФУРЬЕ

На основе вращения вектора численным интегрированием предлагается итерационный линейный алгоритм для рационального вычисления дискретного преобразования Фурье. Основное достоинство алгоритма - вычисление значений тригонометрического базиса за одну линейную итерацию в реальном времени. Как следствие, снижение вычислительного объема достигается за счет быстрого и параллельного вычисления элементов гармоник сигнала и соответствующих тригонометрический функций. Алгоритм дискретного квадратурного преобразования может представлять собой перспективный инструмент для реализации дискретного преобразования Фурье в задачах цифровой обработки сигналов.

Преобразования Фурье; тригонометрические функции.

V.V. Vladimirov, N.S. Zvyaginstev DISCRETE QUADRATURE ALGORITHMS FOR DISCRETE FURRIER TRANSFORMS COMPUTING

Herein introduced an iterative linear algorithm derived from vector's rotation by numerical integration that serves for rational computing of discrete Furrier transform. The main algorithm's benefit is the computing of trigonometric basis by single linear iteration performed in real time scale. As a result the computing process is simplified by fast and parallel work out of harmonic units together with relative trigonometric functions. The algorithm of discrete quadrature conversion can represent the perspectiv tool for implementation of discrete Furrier transform in problems of a digital processing of signals.

Furrier transform; trigonometric functions.

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

должны соответствовать условиям минимизации вычислительного объема (базовых операций сложения, умножения, сдвига) и простоты аппаратурной реализации при возможности организации параллельных вычислений [1]. Также отметим, что для ПЦОС характерным является наличие аппаратного умножителя, позволяющего выполнять умножение за один командный такт [2]. В универсальных же процессорах умножение обычно реализуется за несколько тактов, как последовательность операций сдвига и сложения. Другой особенностью ПЦОС является включение в систему специальных команд как, например, умножения с накоплением С;+1 = А;Б; + С; и с правилом изменения индексов, элементов массивов А и В [2].

Классическая запись ДПФ имеет вид [3]

Хк = I1 хш^тк , к = 0,1,2..^ -1, W = ехР- ]^

т=0 V N

где Хк - массив значений сигнала на выходе преобразования, Хк - массив координат сигнала. ДПФ обратимо, т.е. координатный массив определяется значениями сигнала по аналогичной зависимости.

Существует ряд алгоритмов, позволяющих значительно снизить вычислительный объем при реализации ДПФ. К ним относятся алгоритмы быстрого преобразования Фурье (БПФ), алгоритм Ре дера, алгоритм Виноградова и др. [3]. Не вдаваясь в подробности каждого способа, отметим, что снижение вычислительного объема достигается за счет рационального применения комбинаторных схем организации вычислений. В принятых оценках считается, что реализация ^точечного ДПФ классическим методом обладает вычислительной сложностью

N , БПФ по основанию два снижает сложность до N /о^2 N . Существуют также методы снижения сложности до N [3]. Однако указанная сложность выполнения БПФ и других преобразований рассчитывалась без учета вычисления значений показательных или тригонометрических функций, входящих в состав преобразования. Как правило, существует разграничение задач реализации преобразования как такового и определения значения функций, входящих в ДПФ. В доступной , -ции определяется разложением в показательный ряд (ряд Маклорена) или реализацией дискретных линейных преобразований методом цифра-за-цифрой - алгоритмом Волдера и его производными [4]. Для реализации этих преобразований необходимо выполнить вычисления слагаемых частных сумм рядов или фиксированное количество итераций алгоритма Волдера. Очевидно, что реализация этих вычислений несет с собой дополнительную нагрузку на вычислительные средства и существенно увеличивает вычислительную сложность ДПФ в целом.

Для снижения вычислительного объема и оптимизации процесса реализации ДПФ предлагается использовать дискретные линейные алгоритмы вращения вектора методом численного интегрирования - алгоритмы дискретного квадратурного преобразования (ДКП) [5]. Эти алгоритмы позволяют вычислять значения гармоник вместе с тригонометрическими функциями, не осложняя процесс дополни. -, . -мов состоит в следующем.

Приращение радиус-вектора Я(г) = {х( г),у( г ),г( г)} в декартовой системе координат с базисом {е1;,е2;,ез;}на текущем шаге аргумента Лг; = г;+1 -г; определяется углом вращения а(г) = {аХ( г ),ау( г ),а,7( г)}:

_ +ЛЛ

ЛЯ.(г) = |Щг)Са(г). (1)

г.

После выполнения преобразования (1) вектор

+1 = Щ г.+1) = {х( г.+1), у( г.+1), г( г.+1Л = Я. + Щ переходит в базис {вщ+1),е2(.+1 ),е3(.+1)} и таким образом сохраняется ортогональность преобразования (1). Это соответствует переиндексации . +1 = ..

При наличии начальных значений Я0 выражение (1) описывает вращение вектора вокруг произвольно направленной оси в системе координат, связанной с начальным значением вектора.

Такое вращение сводится к линейной комбинации двух, четырех или шести интегралов вида

О +Лг

1пт = | гп(г)Сат(г), где тФ п, и т,п е{х,у^} (2)

г.

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

а Д I ~1

Г(г) = гщ + £ —ТГ^П (а + в), (3)

I=1 1 ! /?=0

Л А] а .рг

ат (г) = атг + £ ~^Г~П (Я + ^, (4)

] =1 ] • >'=0

где А к/п~к = А к~1 /п~к+1~ ^~1 /п~к - нисходящие разности (/ ={а,г}),

г ~ г. и

а =------ , а и Ь - порядки разложения.

Лг

Ограничивая порядок разложения (3) и (4), запишем алгоритм вращения вектора Я0 в трехмерном пространстве методом дискретного квадратурного преобразования (ДКП):

х.+1 = х. ~ СгЛг( У,аг ) + £уЛг( г,ау ) ;

У.+1 = С7Лг( x, а7,). + У. ~ СхЛг( Ъ ах ).; (5)

Ъ+1 = ~СуЩ( х,ау ). + СхЛг( У,ах ). + ъ,

где /(г5) = Л, / е {х,у,г,ах,ау,аг }, 5 е {. +Ш ~1} , . е {oД,2...}, Сг е {1,0,~1} -

оператор направления вращения (введен для отработки углов в противоположен-

ном направлении и контроля текущих вычислений), а Лг(к,а1) , к, Iе {х, у,ъ}, кф1 определяются численной квадратурой (2):

2

Лг(к,а1). =Лаи £ gih , (6)

Н=0

где gi0 = ki, gil =AaaL (ki - ki-1 \ gi2 =Aak- (i - 2ki-1 + ki-2 ).

2 б

Основное качество алгоритма (5) в том, что при наличии координат R. (значения тригонометрических функций на момент tі) для вычисления вектора R.+1 (значения тригонометрических функций на момент t.+1 = t. + At) необходимо вы-

(5) (б).

'cos a:''

Ri =

sinai

, і = {0,1,2...n -1}значение пары

V

sinai

Ri+1 =

:

Ri+1 =

вычисляется алгоритмом ДКП за одну итерацию по одномерно-

і+1 у

/ \ хі+1 1

\Уі+1J A

ii 1 КУі

- ДКП по формуле прямоугольников;

i

Ri+1 =

;

/ \ хі+і 1

\Уі+1J A

ii JiJ

A a.

V Уі- Уі-1

= AiRi +Sf

ДКП по

Ri+1 =

/ \ хі+1

\Уі+ly

1 - Aa.

KAai 1 у

/ \ Xi

Ji;

Aa

- Х. + Х.

■i^ xi-1 Уі - У і-l

Aa

- хі + 2хі-і - хі-2 v Уі - 2 У.-і + Уі - 2 у

= А1Я1 + 8\г + 8?г - ДКП по формуле квадратичных парабол.

Пусть ак и Ьк - синусоидальная и косинусоидальная компоненты координа-

ты сигнала в функциональном базисе

V

, 2л

cosk-----1

T 2л sink-----1

T

. Тогда вектор Ri+1 определяет

значение базиса главной гармонической компоненты ДПФ, т.е. при k = 1 на мо-

мент t = t.+і. При начальных данных Rq =

1

0

алгоритму ДКП на текущем шаге

для определения указанного базиса надо выполнить одну итерацию, а для определения к-й гармоники необходимо выполнить к итераций с шагом А. При наличии значений гармонических базисов на текущем шаге вычислений определение самих гармоник и, как следствие, значения самой функции сводится к реализации 2N умножений и N слож ений.

Параллельные вычисления алгоритмом ДКП организуются следующим образом:

ak 0 bk О

О Ьк 0)t ;

задается начальный вектор размерности 4 К, алгоритм ДКП выполняет вращение с матрицей А в 2К независимых плоскостях, и каждая координата преобразуется в своей плоскости.

+

2

+

2

б

Гармоники вычисляются суммированием первой и четвертой компоненты каждой пары плоскостей, начиная с первой, т.е. к-я гармоника определяется как хи + Ун , а сам0 значение массива ДПФ равно

К

Хк = а0 + X( х11 + у21 ) .

I =1

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

определяются с множителем, соответствующим координате сигнала.

Структура выполнения обратного дискретного преобразования Фурье полно. . 1, н.

Таблица 1

Состав преобразований для реализации ДПФ

К

Вид ДКФ Кол-во итераций ДКП (n) Кол-во операций сложения Кол-во операций умножения

Последовательный n = K K2 2K

Параллельный 1 < n < K K2 G

2

Состав ДКП на простой итерации

ДКФ Сложение Умножение

Прямоугольник 2 2

Трапеция б 2

Квадратная парабола 12 2

Сравнительный анализ таблиц показывает, что вычислительная сложность алгоритма ДКП при реализации классического ДПФ сравнима с известными алгоритмами БПФ и др. При этом в алгоритме ДКП нет необходимости усложнять вычисления определением тригонометрических или показательных функций. Также отметим ориентацию структуры алгоритма ДКП на базовые операции ПЦОС и возможность применения последовательно алгоритма ДКП для вычисления тригонометрических функций в классических схемах БПФ.

, -струмент для реализации ДПФ в задачах ЦОС, отличающийся простотой и широкими возможностями организации параллельных вычислений.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Букашкин С.А., Лисицын Г.Ф., Миронов ВТ. Применение цифровых процессоров обработки сигналов - магистральный путь развития современных средств телекоммуникаций // Сборник докладов 3-й Межд. конф. «Цифровая обработка сигналов и ее применение». Т. 1. - М., 2000. - C. 3-4.

2. Корнеев В.В., Киселев А.В. Современные микропроцессоры. - М.: Нолидж, 1998. - 240 с.

3. Нуссбаумэр Г. Быстрое преобразование Фурье и алгоритмы вычисления сверток. - M.: Радиосвязь, 1985. - 248 с.

4. Voider J.E. The CORDIC Trigonometric Computing Technique // IRE Trans. On Electronic Computers. - 1959. - Vol. EC-8 (3). - P. 330-334.

5. .. ..

// . - -

. , 2004.

Статью рекомендовал к опубликованию д.т.н., профессор Л.П. Фельдман.

Владимиров Виктор Владимирович

Государственная морская академия им. Ф.Ф. Ушакова.

E-mail: vev-71@mail.ru.

353918, . , . , 93.

Тел.: 88617232322.

Звягинцев Николай Сергеевич

Новороссийское морское пароходство.

E-mail: zns1@yandex.ru.

353918, . , . 93.

.: 88617601763.

Vladimirov Victor Vladimirovich

Novorossiysk State Marine Academy.

E-mail:vev-71@mail.ru.

93, Lenina Pr., Novorossiysk, 353918, Russia.

Phone: +78617232322.

Zvyaginstev Nikolay Sergeevich

Quality division, JSC Novoship.

E-mail: zns1@yandex.ru.

93, Lenina Pr., Novorossiysk, 353918, Russia.

Phone: +78617601763.

УДК 681.518.54

C^. Байдаров

-

РАБОТЫ ДАТЧИКОВОЙ АППАРАТУРЫ

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

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

Вейвлет-преобразования; диагностика; неразрушающий контроль; датчик давления.

S.U. Baydarov

USING WAVELET TRANSFORM FOR DIAGNOSIS OF SENSOR

EQUIPMENT

The article is devoted to the nondestructive testing of piezoelectric smart sensors rapidly and acoustic pressures at the example of research vibrochuvstvitelnogo element with separate electrodes.As a diagnostic sensor vibrochuvstvitelnym element was used direct and inverse piezoelectric effect, and for signal analysis - wavelet transform, followed by research skeylolgram-my.This method allows to increase the information content of the diagnosis and monitoring of various units and units of special equipment.

Wavelet transform; diagnosis; nondestructive testing; pressure sensor.

Для неразрушающей диагностики датчиковой аппаратуры (ДА) широкое распространение получил пьезоэффект, благодаря которому для виброчувствитель-