Обобщенные функции и преобразования Хартли в системах счисления с постоянным основанием Текст научной статьи по специальности «Математика»

ИНФОРМАТИКА И ВЫЧИСЛИТЕЛЬНАЯ

ТЕХНИКА

УДК 519.216.1/2

ОБОБЩЕННЫЕ ФУНКЦИИ И ПРЕОБРАЗОВАНИЯ ХАРТЛИ В СИСТЕМАХ СЧИСЛЕНИЯ С ПОСТОЯННЫМ ОСНОВАНИЕМ

В.В. Сюзев

МГТУ им. Н.Э. Баумана, Москва, Российская Федерация e-mail: [email protected]

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

Ключевые слова: базисная функция, базисная система, преобразования Фурье, спектральный анализ, система счисления.

GENERALIZED FUNCTIONS AND HARTLEY TRANSFORMS IN NUMBER SYSTEMS WITH A PERMANENT BASE

V.V. Syuzev

Bauman Moscow State Technical University, Moscow, Russian Federation e-mail: [email protected]

To solve the problems of digital signal processing in the spectral domain, a new method for analytical synthesis ofdiscrete real parametric basis functions is proposed that uses the Hartmut procedure generalization to the case ofnumber systems with an arbitrary constant-value basis. The mathematical description of functions obtained is given; their basic properties and methods for construction of different orthogonalbasis systems and transformations on their basis are investigated. In terms of these functions, all the basic theorems of spectral analysis used in the theory and practice of digital processing are formulated and proved. The obtained results are original and form the basis of theory of the signal representation and conversion in the new orthogonal basis.

Keywords: basis function, a basis system, Fourier transform, spectral analysis, number system.

Вычислительная и функциональная эффективность решения многих задач цифровой обработки сигналов (ЦОС) спектральными методами существенно зависит от используемых систем базисных функций [1—4]. Поскольку ортогональных систем базисных функций существует неограниченное множество [5], то выбор рационального базиса является сложной теоретической и прикладной проблемой. В этих

условиях особенно полезными могут оказаться параметрические базисные функции, содержащие в своей структуре один или несколько изменяемых параметров, влияющих на их свойства. Известным и важным примером таких базисов служит класс комплексных экспоненциальных функций Виленкина-Крестенсона (ВКФ) [4, 6], управление свойствами которых осуществляется с помощью вариации основания используемой системы счисления. Изменяя величину основания и выбирая различные способы переупорядочения ВКФ, с их помощью можно получить широкое семейство полезных мультипликативных ортонормированных систем, для которых справедливы все теоремы спектрального анализа и существуют быстрые процедуры вычисления спектра [4, 6—9].

Однако, несмотря на указанные достоинства систем ВКФ, их комплексный характер требует использования в алгоритмах ЦОС трудоемкой комплексной арифметики, что может послужить весомым ограничением при практическом применении ВКФ, особенно при обработке высокочастотных многоразмерных сигналов в системах обработки жесткого реального масштаба времени. Поэтому целью настоящей работы поставлено решение теоретико-прикладной задачи синтеза и анализа нового действительного параметрического базиса со свойствами, близкими к свойствам базисов ВКФ, но оперирующего с вещественными числами и операциями. В основе математического подхода к разработке такого базиса положена процедура Хартли, использованная им при создании вещественной альтернативы комплексным экспоненциальным функциям Фурье, состоящим из обычных тригонометрических функций [10, 11], и распространенная здесь на обобщенные тригонометрические функции, образующие мнимые и действительные части ВКФ.

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

cos , используемых в ВКФ [4, 7],

и на интервале из N = рп точек имеют позиционные п-разрядные представления

где кт и {т являются т-ми разрядами этих представлений и лежат в диапазоне [0,р — 1]. Тогда из этих тригонометрических функций можно образовать следующие дискретные функции:

n

n

(1)

m= 1

m= 1

Cas(k, i) = cos

/2п A k . \ ./ 2п А \

I У kmim ) + sin ( / , kmim ) j (2)

VPi j v p i j

4 m=1 7 4 m=1 7

представляющие собой обобщение известных функций Хартли [10,11] на систему счисления с произвольным основанием. Для отражения этого факта в обозначении обобщенных функций Хартли (ОФХ) (2) использовано обозначение обычных функций Хартли, записанное с заглавной буквы. Из развернутой записи ОФХ в виде выражения (2) с помощью известных тригонометрических преобразований можно получить полезное более сжатое их представление:

/ 2П n \ Cas(k, i) = л/2 sin I — ^^ kmim + п/4 1

^ Р m=1 '

2п n

= ^соэ (— - п/4). (3)

^ Р т=1 '

Записанные таким образом ОФХ имеют ряд интересных свойств. Приведем основные из них.

1. ОФХ являются действительными функциями, принимающими только р различных значений. Справедливость этого свойства следует из самого аналитического описания ОФХ в виде (2) или (3).

2. В ОФХ переменные к и г являются равноправными, поэтому, если их поменять местами, функция не изменится, т.е.

Сав(к,г) = Сав(г,к).

В этом проявляется свойство двойственности ОФХ относительно своих аргументов, которое приводит к симметричности матрицы значений ОФХ.

3. ОФХ являются периодическими функциями с периодом N = рп,

п = 1, 2,... .

Это свойство следует из того, что при смещении числа г на N единиц младшие п разрядов в р-ичном представлении числа остаются без изменения.

4. Среднее значение любой ОФХ, кроме нулевой, равно нулю, т.е.

1 N-1

Оав(к, г) = 0, к = 0. (4)

N

i=0

Действительно

N Х1 Сав(к, г) = Х1 ("" Е ^т + п/4

¿=0 ¿=0 ^ Р т=1

^2 р-1 р-1 р-1 2П п 2П

= ^Е Е ктгт + п/4) + —к1г1]. (5)

г„ =0 г„-1=0 ¿1=0 т=2

Рассмотрим в этом выражении внутреннюю сумму по индексу г1. Она является табличной и представляется в виде произведения трех

сомножителей [12]: ( 2п

p-i

S sin

ii=0

V p

y^ kmim + n/4 1 +--kiii

m=2

2n

P

(2n

= sin —

v p

V kmim + П + -(p - 1)ki ) sin(nki) cos ec( nki )

m=2 4 p / VP /

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

Среднее значение нулевой ОФХ равно единице, так как Сав(0, г) = 1

1

N-1

и

N

Ei = i.

i=0

5. Мощность Pk любой k-й ОФХ равна единице:

1

n- i

Pk = ^Е | Cas(k, i)

=1.

(6)

i=0

Для доказательства этого свойства представим квадрат ОФК с учетом соотношений (3) в виде

_.2 / 2П n Cas(k,i) = 2 sin i — ^^ kmim +

n 2n n

~т I cos ( ^ ^ kmir

m=i

4

П

P

m=i

4

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

Оав(к, г) = 1 + ( —^ ^ ктгг

V р 1

4 т=1

С учетом этого результата мощность Р^ можно записать так:

1 N-1 (—п ' Рк = 1 + NN Е ^п ( — У ктг

i=0

P

mbm

m=i

1 P-i P-i P-i

i+^E E ...Esi

sin | — S kmi^ . (7) P m=i

г„=0 г„_1=0 гх=0

В выражении (7) внутренняя сумма по индексу г1 также табличная и тоже равна произведению

81П

4п P

У^ kmim +---1— (р — 1)ki sin(2nki) cos ec( —k\ \,

m=2 р р J v р у

р р

которое из-за сомножителя в1п(2пк1) равно нулю при любом значении

2

к1. Поэтому вся многомерная сумма в уравнении (7) также будет равна нулю и Рк = 1. Справедливость соотношения (6) доказана.

6. Обобщенные функции Хартли являются ортогональными функциями, т.е.

1

N

N-1

^Cas(k,i)Cas(\,i) = 0, k = Л.

(8)

¿=0

Для доказательства этого свойства вновь используем описание ОФХ выражениями (3). Тогда произведение ОФК будет равно

Cas(k,i)Cas^,i) = "2п

= Э1П

Р

^ ^ (km + Лm)im

m=1

+ соэ

2п Р

^ ^ (km Лm)im

m=1

а сумма таких произведений по индексу г принимает следующий

вид:

N-1

N- 1

¿=0

¿=0

Cas(k, i)Cas(Л, i) = ^^ sin — ^^(km + Лт)1

Р N1

^m ~ /vm I hm

m=1

+

¿=0

2п

+ ^ ^ COs ^ ^ (km Лm)im

Р

m=1

. (9)

Рассмотрим каждое слагаемое этой суммы по отдельности. Запишем первое слагаемое в многомерном виде

N1

¿=0

2п

^ ^ Sin ^ ^ (km + Лm)im

Р

m=1

p— 1 p— 1 p—1

£

in=0 in-1=0 ¿1 =0

Й1П

- y~](km + Лm)im +--(k1 + Лх)г1

Р

m=2

Р

Внутренняя сумма по индексу г1 здесь будет табличной и представляется в виде произведения, один из сомножителей которого равен э1п[п(к1 + Л1)] и принимает нулевые значения при любых значениях к1 и Л1. Вследствие этого и внутренняя сумма, и вся многомерная сумма также будут равны нулю.

Для второго слагаемого в выражении (9) по аналогии имеем

N1

¿=0

2п

^ ^ COs ^ ^ (km Лm)im

Р

m=1

p— 1 p— 1 p— 1 Е Е ...Еcos

in=0 in-1=0 ¿1 =0

2^ n 2п

— Лт )im +--(k1 — Л1)i1

Р

m=2

Р

При = Л1 внутренняя сумма представляется произведением, один из сомножителей которого имеет вид з1п[п(к1 — Л1)] и равен нулю при любых не равных друг другу значениях к1 и Л1. В случае равенства значений младших разрядов кодов чисел к и Л (т.е. при к1 = Л1) внутренняя сумма во втором слагаемом выражения (9) равна

p-i

Е cos

n

Р ^

1 m=2

= p cos

n

Р ^

1 m=2

и второе слагаемое представляется уже (n — 1)-мерной суммой

p—1 p—1 p—1

cos

in=0 in-1=0 ¿2=0

Р^ У cos — - )im + — (&2 - ^¿2

РР

m=3

Внутренняя сумма этого выражения в зависимости от значений к2 и Л2 либо будет равна нулю, либо р, что снова приведет к уменьшению размерности представления второго слагаемого в общем выражении (9). Поскольку к = Л, то обязательно хотя бы одни значения кт и Лт не будут равны между собой. Поэтому и эта многомерная сумма в конечном счете также станет равной нулю. Таким образом, условие ортогональности доказано.

Используя свойства 5 и 6, можно записать свойство ортонормиро-ванности ОФХ:

1 ^ ( 1, к = Л,

— > Саз(к, г)Саз(Л, г) = < * ¿0 ( , ) ( , ) I 0, к = Л.

7. Объединение N первых ОФХ приводит к полной ортонорми-рованной базисной системе, пригодной для представления любых решетчатых сигналов конечной мощности, определенных на дискретном интервале [0,*). Полнота системы обеспечивается тем, что к ней невозможно добавить ни одной новой функции, которая была бы ортогональна ко всем остальным.

Системы дискретных ОФХ удобно записывать в виде матриц С их значений. Эти матрицы будут симметрическими и ортогональными. Обратные к ним матрицы будут совпадать с прямыми с точностью до постоянного множителя 1/* : С= С. Свойства матриц ОФХ одинаковы для строк и столбцов (в силу их симметричности). Матрицы содержат ровно р различных действительных элементов. Элементы нулевых строк и столбцов равны единице.

Для фиксированных значений р возможны матрицы ОФХ, отличающиеся порядком следования строк и столбцов, т.е. другими словами, возможны различные способы упорядочения функций Саз(к,г)в системе. Это свойство матриц ОФХ подобно аналогичному свойству матриц значений ВКФ [4, 7] и позволяет значительно расширить ассортимент действительных тригонометрических базисов.

Изменение порядка следования функций в базисной системе достигается путем применения различных замкнутых операций переупорядочения к номерам базисных функций. Наибольшее распространение получили операции инвертирования р-ичных кодов чисел к и их обобщенное кодирование Грея [4,6]. Применим их к ОФХ.

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

Оав(к, г) = соэ ( — ^ кп+1_тгт ) + эт ( — ^ кп+1_тгт ), (10) ^ Р т=1 ' \ Р т=1 /

Cas(k, i) = а/2 sin ^ k,n+i-mim + - J =

(2ъ

m=1 n

= v^cos ^kn+i-mim - -j, (11)

m = 1

а заменяя прямой код чисел k их кодом Грея, — обобщенную систему Хартли - Хармута:

Cas(k, i) =cos i— < km > imj +sin (— < km > imj, (12)

^ P m=1 ' \ P m=i /

Cas(k, i) = V2 sin f — < km > im + П j =

^ P m=1 '

= V-coS (- 2 < km > im - (13)

^ P m=1 '

где разряды кода Грея < km > вычисляются по правилу: < km >= km+ + km+1 ( mod p), kn+1 = 0. Фамилии Пэли и Хармута включены в название систем ОФХ (10)—(13) по аналогии с системами ВКФ.

Следует отметить, что все ранее приведенные свойства ОФХ будут справедливы для всех перечисленных базисных систем. Однако свойства спектров конкретных сигналов в системах с различным порядком следования ОФХ могут сильно отличаться.

Системы ОФХ носят обобщенный характер. Из них за счет выбора основания p системы счисления можно получить множество известных и неизученных систем. Так при p = N и n =1 все системы ОФХ переходят в систему обычных функций Хартли, так как в этом случае

Cas(k, i) = cos (-^kij + sin (-^kij = cas(k, i).

При p = 2 и n = 1 из систем ОФХ получаются различные системы Уолша [4, 13]: из ОФХ (2) — системы Уолша-Адамара

Cas(k,i) = cos (к kmim ) = (-1)Em=1 kmim = had(k,i);

^ m— 1 '

m=1

из ОФХ (10) — системы Уолша-Пэли

Cas(k, i) = cos kn+1-mim ) = (-1)^m=1 kn+i-mim = p a/(k, i);

^ m=1 '

из ОФХ (12) — системы Уолша-Хармута

Cas(k, i) = cos ln^2<km > im ) = (-1)^= <km>im = har(k, i).

^ m=1 '

Обобщенные преобразования Хартли и свойства обобщенных спектров Хартли. Обобщенные преобразования Хартли (ОПХ) представляются в виде следующей пары прямого и обратного дискретных преобразований Фурье

1 N—1

Xx(k) = x(i)Cas(M), (14)

i=0

N — 1

x(i) = £ Xx(k)Cas(k, i), (15)

k=0

где x(i) являются отсчетами дискретного входного сигнала, а XX(k) — составляющими его обобщенного спектра Хартли. Обе решетчатые функции x(i) и XX (k) в ОПХ являются действительными и определены на целочисленном интервале [0, N). Энергетическая взаимосвязь сигнала и его спектра в базисе ОФХ устанавливается равенством Пар-севаля

1 N—1 N — 1

i £ *2(i) = £ xx (k)

г=0 k=0

справедливость которого подтверждает к тому же полноту базисной системы ОФХ.

Системы ОФХ не обладают свойством мультипликативности, поскольку произведение двух любых функций Хартли не дает функцию той же системы. По этой причине в базисе ОФХ в прямом виде не выполняются теоремы преобразования спектров, записанные для мультипликативных базисов и имеющие важное значение в теории и практике дискретного спектрального анализа (теоремы об обобщенном сдвиге сигналов, о модуляции сигнала и спектра, о свертке и корреляции, об умножении сигналов и т.п.) [4,5]. Однако эти теоремы можно сформулировать и в терминах спектров ОФХ, если использовать взаимосвязь спектров ОФХ со спектрами мультипликативных

базисов Виленкина-Крестенсона, для которых теоремы спектрального анализа справедливы [4].

Обобщенные функции Хартли и ВКФ используют в своей струк-

туре одинаковые обобщенные тригонометрические функции: четные

системы этих функций являются родственными и их отличие состоит

только в том, что в действительных ОФХ эти функции используются

в качестве слагаемых, а в комплексных ВКФ — в качестве их действи-

тельной и мнимой частей соответственно.

соэ

m = 1

и нечетные эт

. В этом смысле обе

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

1 N-1 (2П п \

Хч(к) = — Е ж(«) СОй ( — Е кт«т ) ,

г=0 ^ Р т=1 '

1 N-1 ^ . (2п

Xя(k) = ^^ Е Х(^) Sin ( - E kmim

г=0 ^ P m=1

Тогда спектры Хх (k) Хартли и ХВк(k) Виленкина-Крестенсона можно записать в виде

Xx (k) = Хч (k) + Хя (k), (16)

Xbk(k) = Хч(k) - jXH(k). (17)

В последнем уравнении j означает мнимую единицу (j = л,/—1), а наличие знака минус перед мнимой частью спектра Виленкина-Крестенсона связано с комплексно-сопряженным характером прямого преобразования Фурье в базисе ВКФ [4]. Поскольку ХЧ(k) является четной функцией переменной k, а Хя(k) — нечетной функцией той же переменной, полезно ввести спектр Хартли для отрицательных значений его номера

Хх (—k) = Хч (k) — Хя (k), (18)

который в базисе ОФХ будет играть роль, аналогичную роли комплексно-сопряженного спектра ХВк (k) в базисе ВКФ. Последовательно суммируя и вычитая уравнения (16) и (18), после преобразования получаем, что

Хч(k) = [Хх(k) + Хх(—k)]/2, Хя(k) = [Хх(k) — Хх(—k)]/2, и спектр Виленкина-Крестенсона принимает следующий вид:

Хвк (k) = [Хх (k) + Хх (—k)]/2 — j [Хх (k) — Хх (—k)]/2. (19)

Это уравнение и является уравнением связи спектров в родственных базисах ОФХ и ВКФ. Его использование позволяет свойства спектров Виленкина-Крестенсона трансформировать в свойства спектра Хартли. Выполним такую трансформацию, записав основные свойства спектров в виде соответствующих теорем спектрального анализа, как это принято в ЦОС [1-5].

Теорема 1. О спектре сигналов с обобщенным сдвигом во времени. Спектр Хартли сигнала, сдвинутого по оси времени i на величину т по закону операций прямого i ф т и обратного i 0 т обобщенного сдвига, выполняемого с помощью поразрядного суммирования или вычитания по модулю p p-ичных кодов чисел i и т, равен спектру Хартли несдвинутого сигнала, модулированному обобщенными тригонометрическими функциями в момент времени т.

Доказательство. Если сигнал у (г) получается путем прямого обобщенного сдвига г ф т по оси времени г на величину т сигнала х(г) (т.е. у (г) = х(гфт)), то его спектр в базисе ВКФ имеет следующий вид [4, 6]:

г 2П п \

Убк(к) = Хвк(к) ехр ( з — Е кттт) , (20)

^ Р т=1 '

где все тт являются разрядами р-ичного кода числа т,

п

Т = ^ ТтРт-1.

т=1

Используя уравнения (17) и (19) и развернутую запись комплексной обобщенной экспоненты, выражение (20) можно преобразовать к виду

Ybk (k) =

XX(k) + Xx(-k)coJ^ V k T ) +

coS I / J kmTm +

2 \p Xx (k) - Xx (—k) . f 2n

+ XX(k) " XX(-k) sin ( - ^ kmTm

p

m=1

n

^ i

то ' m

Xx(k) - Xx(-k)co^2n , k T

I 7 lbrn I rn. I

m=l

n

mTm

m 1 m

km

2P

m=l

n

kmTm

Xx(k)+ Xx(-k)si/2% k t

Sil! 7 kmTm

2 V Р

т=1

Учитывая теперь, что спектр Хартли Ух (к) получается путем суммирования действительной и мнимой частей спектра Виленкина -Крестенсона Убк(к), в окончательном виде получаем

f 2П n \

Yx (k) = Xx(k) cos f — kmTm )

^ P m=1 '

р т=1

г 2П п \ - Хх (-к)эт(- У^ ктТт) . (21)

^ Р т=1 '

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

f 2П n \

Yx (k) = Xx(k) cos i — kmim ) +

^ P m=1 '

{2-K n \

+ Xx (-k)sini— У~] kmTm). (22)

^ P m=1 '

Выражения (21) и (22) и представляют собой аналитическую запись теоремы об обобщенном сдвиге сигнала в терминах спектров Хартли. Теорема 1 доказана. Из нее следует, что обобщенный сдвиг сигнала по оси времени приводит к модуляции спектра несдвинутого сигнала обобщенными тригонометрическими функциями (составляющие спектра с положительными номерами — четными, а составляющие спектра с отрицательными номерами — нечетными тригонометрическими функциями). Сами спектры сдвинутого сигнала в этом случае выражаются при прямом сдвиге через разность, а при обратном сдвиге — через сумму модулированных составляющих.

Теорема 2. О модуляции сигнала базисной функцией. Умножение сигнала на базисную функцию приводит к изменению порядка следования его спектральных составляющих.

Доказательство. Пусть сигнал у (г) получается путем умножения

сигнала х(г) на базисную функцию Саз(А, г). Тогда

1 ^

(к) = — Е[х(г)Сав(А,г)]Са8(М).

г=0

Используя развернутую запись ОФХ в виде (2) и применяя известные тригонометрические теоремы, после преобразования получаем

(к) = 1[Хх(к + А) + Хх(к - А) + ХХ(А - к) - ХХ(-к - А)]. (23)

Теорема 2 доказана. Из нее следует, что при модулировании сигнала базисной функцией сам спектр сигнала не меняется, но меняется порядок его следования.

Теорема 3. О свертке. Спектр Хартли сигнала, являющегося результатом обобщенной свертки двух других сигналов, равен с точностью до постоянного множителя произведению спектров этих сигналов.

Доказательство. Пусть сигнал у(г) есть обобщенная свертка двух сигналов х(г) и и(г):

N -1

у(г) = Е х(А)и(А 0 г).

А=0

В базисе ВКФ спектр такой свертки равен [4]

Увк (к) = (к)ивк (к).

Используя связь спектров Виленкина-Крестенсона и Хартли в форме уравнения (19), из этой формулы после преобразования получаем

Хек (к) = "4|[Хх (к) + Хх (-ЩЦх (к) + (-&)]-- [Хх(к) - Хх(-Щ^х(к) - их(-&)]}-

3^{[Xx (k)+ Xx (-k)][Ux (k)-

N 1

- их (-к)] + [Хх (к) - Хх (к)][Цх (к) + их (-к)] откуда следует, что

N

Yx (k) = ~2

Xx (k)Ux (-k) + Xx (-k)Ux (k) + Xx (k)Ux (k)- Xx(-k)Ux(-k)

= "21 [Хх(к) + Хх(-к)]их(к) + [Хх(к) - Хх(-к)]их(-к)

Эта зависимость представляет собой математическое описание теоремы об обобщенной свертке на языке спектров Хартли.

Теорема 3 доказана. Видно, что ее запись в базисе ОФХ сложнее, чем в базисе ВКФ. В частном случае, когда один или оба свертываемых сигнала являются либо четными, либо нечетными относительно середины интервала определения, запись теоремы 3 о свертке в базисах ОФК и ВКФ совпадает. Это связано с тем, что в этом случае спектры Хартли также являются четными или нечетными функциями номера к. Так, например, для четного сигнала х(г) = х(— - г) и Хх (к) = Хх (-к), поэтому Ух (к) = -Хх (к)Цх (к).

Теорема 4. О корреляции. Спектр Хартли обобщенной корреляционной функции сигнала равен спектру мощности этого сигнала.

Доказательство. Если обобщенная корреляционная функция двух сигналов имеет вид

! *-1

У(г) = "уЕ х(А)х(А ф г)

А=0

то ее спектр в базисе ВКФ совпадает со спектром мощности этого сигнала [4]:

Убк (к) = ХВ к (к)Хвк (к).

Используя уравнение (19), это энергетическое соотношение можно записать и в базисе ОФК:

Ух(к) = 2[ХХ (к) + ХХ (-к)].

Теорема 4 доказана.

Теорема 5. О независимости спектра мощности и энергетического спектра сигнала от его обобщенного сдвига по оси времени.

Спектр мощности и энергетический спектр сигнала в базисе Хартли не изменяются при его обобщенном сдвиге по оси времени.

Доказательство. Пусть сигнал у (г) получают из сигнала х (г) путем его обратного обобщенного сдвига на время т, т.е. у(г) = х(г 0 т).

Энергетический спектр и спектр мощности этого сигнала в базисе ВКФ имеют следующий вид записи [4]:

(к) = 1 Увк (ВД* (к) = 1 Хвк (к)ХВк (к),

Ру (к) = (к) = Ув* (ВД* (к) = Хв* (к)Х£ * (к).

Подставляя в эти формулы вместо спектров Виленкина-Крестенсона их представления в виде спектров ОФХ (19), после преобразования получаем

^ (к) = ^ УХ (к) + УХ (-к)] = ^ [XX (к) + XX (-к)] = зд, Ру (к) = 2 [УХ (к) + УХ (-к)] = 2[ХХ (к) + XX (-к)] = Рх(к).

Теорема 5 доказана. Она свидетельствует об инвариантности энергетических характеристик сигнала в базисе ОФХ относительно его обобщенного сдвига по оси дискретного времени.

Теорема 6. Об умножении сигналов. В базисе ОФХ спектр сигнала — произведения двух других сигналов представляется сумой обобщенных сверток спектров этих сигналов.

Доказательство. Пусть сигналы ж(г), у(г) и м(г) имеют одинаковый интервал определения [0,Х) и у(г)=ж(г)м(г).

В базисе ВКФ спектры этих сигналов связаны между собой соотношением обобщенной свертки [4], т.е.

N -1

Ув*(к) = £ Хв*(А)Цвк(А 0 к). (24)

А=0

Но в соответствии с выражением (19)

„ ... Хх (А)+ ХХ (-А) ,ХХ (А) - ХХ (-А)

(А) --^--•-2-,

ив* (А 0 к) - (их (А 0 к) + их (-(А 0 к)))/2-

- • (их (А 0 к) - их (-(А 0 к)))/2.

Поэтому с учетом этих соотношений из общей формулы (24) после преобразований получаем

1 *-1

Ybk(k) = ^ £ {[Xx(А) + Xx(-A)][Ux(А ö k) + Ux(-(А ö k))]-

Л=0

- [XxU) - Xx(-A)][Ux(A ö k) - Ux(-(A ö k))]}-

1 N-1

- j1 Е {[Xx(A) + Xx(-A)][Ux(A ö k) - Ux(-(A ö k))]+

Л=0

+ [XX (Л) - Хх (-Л)][их (л е к))] + их (-(Л е к))]}.

Выделяя в комплексном спектре Увк(к) действительную и мнимую части и суммируя их, получаем спектр сигнала у (г) в базисе ОФХ:

Yx (k) = 2

N — 1 N— 1

^ Xx(A)Ux(-(Л 0 k)) + ^ Xx(-A)Ux(Л ö k)+

A=0 Л=0

N 1 N 1

+ ^ Хх (Л)их (л е к) - Хх (-Л)их (-(ле))

А=0 А=0

Этот спектр выражается сумой четырех обобщенных сверток спектров Хартли сигналов-сомножителей.

Выражение (25) можно упростить и представить спектр Ух в виде суммы двух сверток, если учесть, что

[их (Л е к) + их (-(Л е к))]/2 = ич (Л е к), [их (Л е к) - их (-(Л е к))]/2 = ин (Л е к),

и поэтому

N— 1 N— 1

Ух (к) = ^ Хх (Л)ич (Л е к) + ^ Хх (-Л)ин (Л е к).

А=0 А=0

Теорема 6 доказана. Выбор конкретного варианта ее записи зависит от вида сигналов ж (г) и и(г).

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

кт = к*т = (р - кт)р = (р - кт) (шодр), (25)

где (р - кт)р означает вычет числа р - кт по модулю р. Соотношение (25) легко подтверждается, если учесть, что вычет разности двух чисел р и кт равен разности вычетов и вычет числа р по модулю р равен нулю, а вычет числа кт, меньшего р, равен самому числу кт.

Все приведенные теоремы обобщенного спектрального анализа Хартли записаны и выведены для систем обобщенных функций Хартли - Адамара. Однако они будут выполняться и для систем Хартли с другими способами упорядочения базисных функций (Пэли и Хар-мута). Эти теоремы носят обобщенный характер и при конкретных значениях параметра р могут приводить как к известным, так и к новым результатам. Так, например, при р = N и п = 1, когда ОФХ

переходят в обычные функции Хартли, приведенные теоремы совпадают с теоремами обычного спектрального анализа Хартли [10], а при p = 2 и N = 2n, когда ОФХ становятся функциями Уолша, они совпадают с теоремами спектрального анализа Уолша [4, 13].

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

Схожесть свойств полученных базисных систем со свойствами комплексных систем Виленкина-Крестенсона и однозначная взаимосвязь спектров этих систем позволяет рассматривать базис ОФК в качестве вещественной альтернативы комплексному базису ВКФ. Можно сказать, что вряд ли существуют такие задачи, для которых справедливо использование комплексных преобразований Виленки-на - Крестенсона и одновременно не может быть применено вещественное обобщенное преобразование Хартли.

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

ЛИТЕРАТУРА

1. Айфичер Э., Джервис Б. Цифровая обработка сигналов: практический подход, 2-е издание / пер. с англ. М.: Издательский дом "Вильямс", 2004. 992 с.

2. Залманзон Л.А. Преобразования Фурье, Уолша, Хаара и их применение в управлении, связи и других областях. М.: Наука, 1989. 496 с.

3. Цифровая обработка сигналов и изображений в радиофизических приложениях / под ред. В.Ф. Кравченко. М.: Физматлит, 2007. 554 с.

4. Трахтман А.М., Трахтман В.А. Основы теории дискретных сигналов на конечных интервалах. М.: Сов. радио, 1975. 208 с.

5. Трахтман А.М.Введение в обобщенную спектральную теорию сигналов. М.: Сов. радио, 1972. 352 с.

6. Сюзев В.В. Энергетические спектры сигналов в базисе функций Виленкина-Крестенсона, инвариантные к циклическому сдвигу // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. 2007. № 2. С. 64-77.

7. Сюзев В.В. Операторы взаимосвязи спектров в базисах комплексных экспоненциальных функций и функций Виленкина-Крестенсона // Инженерный журнал: наука и инновации. 2013. Вып. 11. URL:http://engjournal. ry/catalog /it/ hidden /1051.html.

8. Сюзев В.В., Савельев А.Я., Гудзенко Д.Ю. Методы представления и преобразования сигналов в базисе обобщенных функций Крестенсона // Наука и образование. 2012. № 3. http://technomag.edu.ru/doc/372760.html

9. Власенко В.А., Лаппа Ю.М., Ярославский Л.П. Методы синтеза быстрых алгоритмов свертки и спектрального анализа сигналов. М.: Наука, 1990. 180 с.

10. Брайсуэлл Р. Преобразование Хартли / пер. с англ. М.: Мир, 1990. 175 с.

11. Сюзев В.В. Теоретические основы спектрального анализа в базисе Хартли // Наука и образование. 2011. № 10. URL: http://technomag.edu.ru/doc/230816.html.

12. Двайт Г.Б. Таблицы интегралов и другие математические формулы / пер. с англ. М.: Наука, 1969. 227 с.

13. Сюзев В.В. Скалярный метод синтеза быстрых преобразований Уолша-Адамара // Вестник МГТУ им. Н.Э. Баумана. Сер. Приборостроение. Спец. вып. "Информационные технологии и компьютерные системы". 2011. С. 128-137.

REFERENCES

[1] Ifeachor E.C, Jervis B.W. Digital Signal Processing: A Practical Approach. 2nd Ed. 2002. 933 p. (Russ. Ed.: Ayficher E., Dzhervis B. Tsifrovaya obrabotka signalov: prakticheskiy podkhod. 2-e izd. Moscow, St. Petersburg., Kiev, Izd. dom "Vil'yams" Publ., 2004. 992 p.)

[2] Zalmanzon L.A. Preobrazovaniya Fur'e, Uolsha, Khaara i ikh primenenie v upravlenii, svyazi i drugikh oblastyakh [Fourier, Walsh and Haar transforms, and their application in control circuit, communication and other fields]. Moscow, Nauka Publ., 1989. 496 p.

[3] Kravchenko V.F. Tsifrovaya obrabotka signalov i izobrazheniy v radiofizicheskikh prilozheniyakh [Digital signal and image in radiophysical applications]. Moscow, Fizmatlit Publ., 2007. 554 p.

[4] Trakhtman A.M., Trakhtman V.A. Osnovy teorii diskretnykh signalov na konechnykh intervalakh [The heart of the theory of discrete signals on finite intervals]. Moscow, Sovetskoe Radio Publ., 1975. 208 p.

[5] Trakhtman A.M. Vvedenie v obobshchennuyu spektral'nuyu teoriyu signalov [Introduction to the generalized spectrum theory of signals]. Moscow, Sovetskoe Radio Publ., 1972. 352 p.

[6] Syuzev V.V. The energy spectra of signals in the Vilenkin-Christenson functions basis invariant to the cyclic shift. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Estestv. Nauki [Herald of the Bauman Moscow State Tech. Univ., Nat. Sci.], 2007, no. 2, pp. 64-77 (in Russ.).

[7] Syuzev V.V. Association functional of spectra in the complex exponential and Vilenkin-Christenson functions bases. Inzhenernyy zhurnal: nauka i innovatsii [Eng. J.: Science and Innovation], 2013, iss. 11. pp. 1-7 (in. Russ.). Sc. Educ: El. Sc. Techn. Publ. Available at: http://engjournal.ru/catalog/it/hidden/1051.html (accessed 22.02.2014).

[8] Syuzev V.V. Presentation methods and signal conversion in the generalized Christenson functions basis. Nauka i obrazovanie: elektronnoe nauchno-tekhnicheskoe izdanie [SCIENCE & EDUCATION: El. Sc. Techn. Periodical], 2012, no. 3, pp. 1-28 (in Russ.). Available at: http://technomag.edu.ru/doc/372760.html (accessed 22.02.2014).

[9] Vlasenko V.A., Lappa Yu.M., Yaroslavskiy L.P. Metody sinteza bystrykh algoritmov svertki i spektral'nogo analiza signalov [Synthesis technique of fast algorithms of fold and spectral analysis of signals]. Moscow, Nauka Publ., 1990. 180 p.

[10] Bracewell R. Hartley Transform: theory and applications. New York: Oxford University Press, 1986. (Russ. Ed.: Braysuell R. Preobrazovanie Khartli Moscow, Mir Publ., 1990. 175 p.)

[11] Syuzev V.V. Theoretical foundations of spectral analysis in the Hartley basis. Nauka i obrazovanie: elektronnoe nauchno-tekhnicheskoe izdanie [SCIENCE & EDUCATION: El. Sc. Techn. Periodical], 2011, no. 10, pp. 1-47. Available at: http://technomag.edu.ru/doc/230816.html (accessed 22.02.2014) (in Russ.).

[12] Dwight H.B. Tables of Integrals and Other Mathematical Data. 3th Ed. The Macmillan Co., N.Y., 1947. 198 p. (Russ. Ed.: Dvayt G.B. Tablitsy integralov i drugie matematicheskie formuly. Moscow, Nauka Publ., 1969. 227 p.)

[13] Syuzev V.V. Scalar method for the synthesis of the fast Walsh-Hadamard transform. Vestn. Mosk. Gos. Tekh. Univ. im. N.E. Baumana, Priborostr., Spetsvyp. "Informatsionnye tekhnologii i komp'yuternye sistemy" [Herald of the Bauman Moscow State Tech. Univ., Instrum. Eng., Spec. Issue "Information technologies and computer systems"], 2011, pp. 128-137 (in Russ.).

Статья поступила в редакцию 16.12.2013

Владимир Васильевич Сюзев — д-р техн. наук, профессор, заведующий кафедрой "Компьютерные системы и сети" МГТУ им. Н.Э. Баумана. Автор более 200 научных работ в области цифровой обработки сигналов и компьютерных информационно-управляющих систем.

МГТУ им. Н.Э. Баумана, Российская Федерация, 105005, Москва, 2-я Бауманская ул., д. 5.

V.V. Syuzev — Dr. Sci. (Eng.), professor, head of "Computer Systems and Networks" department of the Bauman Moscow State Technical University. Author of more than 200 publications in the field of digital processing of signals and computer information and control systems.

Bauman Moscow State Technical University, Vtoraya Baumanskaya ul. 5, Moscow, 105005 Russian Federation.