УДК 517, 621.391
С.Я. Новиков, М.Е. Федина1
ПРИНЦИПЫ НЕОПРЕДЕЛЕННОСТИ НА ГРУППАХ И ВОССТАНОВЛЕНИЕ СИГНАЛОВ2
Показано, как принципы неопределенности гармонического анализа переносятся на конечные абелевы группы. Выделены недавние результаты Т. Тао и его соавторов о циклических группах простого порядка. Найдены аналоги гауссовых функций на конечных абелевых группах, индикаторные функции подгрупп. Доказан конечномерный вариант формулы суммирования Пуассона. Намечены возможности применения полученных результатов для восстановления дискретных сигналов по неполному набору коэффициентов. Сформулирован принцип частичной изометрии, в соответствии с которым можно определить минимальное количество измерений для устойчивого восстановления сигнала.
Ключевые слова: принципы неопределенности, циклические конечные группы, восстановление, разреженный сигнал, индикаторные функции, формула Пуассона.
Принцип неопределенности в гармоническом анализе утверждает: если функция / : О ^ С на абелевой группе О сконцентрирована на маленьком множестве, то ее преобразование Фурье / : О ^ С имеет "достаточно большой" носитель. Известно много результатов, точно выражающих этот принцип.
Например, для вещественной прямой О = М и стандартного преобразования Фурье /(£) = /(х)в-2п1Х^3х справедлив [1]
Принцип неопределенности Гейзенберга. Если \\/\\ь2(к) = \\ /\\ь2(К) = 1 и хо, е М, то
\\(х - хо)/\\ь2(к)\\(£ - £о)/\\ь2(ж) >
Лаконичная запись этого принципа имеет следующий вид:(Д^x)(Af£) ^ • Равенство достигается на центрированных гауссовых функциях
пАх2. ?= с в-п£2/А
/(х) = ; / = -=е
при хо = Со = 0.
х© Новиков С.Я., Федина М.Е., 2015
Новиков Сергей Яковлевич (nvks@samsu.ru), кафедра теории вероятностей и математической статистики, Самарский государственный университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
Федина Мария Ефимовна (phedina75@gmail.com), кафедра безопасности информационных систем, Самарский государственный университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
2Работа выполнена при финансовой поддержке Минобрнауки России в рамках базовой части государственного задания, проект № 204.
Для конечной (абелевой) аддитивной группы (О, +) пусть #О обозначает мощность О. Нормированная мера Хаара на этой группе — это нормированная считающая мера
/ (х^х :=#оТ,/(х).
# хео
Получаем конечномерное гильбертово пространство Ь2(О) функций / : О ^ С со скалярным или внутренним произведением
{/,9)ь2(о) := /(х)д(х)Лх.
о
(х)
о
Так как группа О конечна, все нормы эквивалентны, так что Ь2(О) = Ь1(О) = = Ь(О). Каждый элемент группы у € О порождает оператор сдвига Тгапву : Ь2(О) ^ Ь2(О), определенный соотношением
ТгапВу/(х) := /(х - у).
Каждый такой оператор является унитарным на Ь2(О), кроме того, отображение у ^ Тгапву является гомоморфизмом:
Тгапву Тгапв2 = Тгапву+2. (1)
Другими словами, отображение у ^ Тгапву является унитарным представлением группы О, действующей на гильбертовом пространстве, это Ь2(О), т. н. регулярное представление. Одно из следствий (1) состоит в том, что все сдвиги коммутируют между собой.
Свертка функций /,д € Ь2(О) обозначается / * д, / * д € Ь2(О) и определяется формулой
/ * д := /(у)д(х - у) Лу.
о
Свертка определяет билинейную, ассоциативную и коммутативную операцию (последнее только для абелевой группы). Существует единичный элемент 6 € Ь2(О), определяемый в виде 6(х) := (#О)1щ, таким образом, / * 6 = 6 * / = 1.
Для фиксированной функции д € Ь2(О) можно определить оператор свертки Тд : Ь2(О) ^ Ь2(О) следующим образом: Тд/ := д * / или в интегральной форме
Тд = д(у)Ттапву Лу.
о
Таким образом, операторы свертки оказываются линейными комбинациями операторов сдвига. Операторы свертки коммутируют между собой и инвариантны относительно сдвига. И обратно каждый инвариантный относительно сдвига оператор является оператором свертки (отображение д ^ Тд можно рассматривать как представление сверточной алгебры Ь2(О) на себя).
Определим подпространство V С Ь2(О), инвариантное относительно сдвига, как подпространство V, инвариантное относительно всех сдвигов Тгапву, т. е. ТгапвуV = V для всех у. Каждое такое подпространство есть компонента регулярного представления О. Эквивалентно V сохраняется при действии всех операторов свертки Тд. Кроме тривиальных примеров {0}, Ь2(О), простыми примерами являются подпространства констант {с : с € О} и подпространство {/ € Ь2(О) : О / = 0} функций с нулевым средним. Другая важная пара примеров: ядро {/ € Ь2(О) : Тд/ = 0} и образ {Тд/ : / € Ь2(О)} оператора свертки Тд (основано на том, что операторы свертки коммутируют со сдвигами). Векторная сумма и пересечение двух инвариантных относительно сдвига подпространств
инвариантно относительно сдвига. Так как операторы сдвига унитарны, ортогональное дополнение подпространства, инвариантного относительно сдвига, также инвариантно относительно сдвига.
Пример 1. Циклическая группа G = Z/NZ. Для данного £ G Z/NZ можно построить одномерное инвариантное относительно сдвига подпространство V^, порожденное характером e^ : x ^ e2nlxt/N. Подпространство инвариантно относительно сдвига, тогда оно является прямой суммой некоторых из этих подпространств Vt, или имеет вид {f G L2(G) : /\e = 0} для некоторого фиксированного множества частот E с Z/NZ.
Так как Tg5 = g для всех g G L2(G), получаем, что инвариантное относительно сдвига подпространство, содержащее сверточную единицу 5, будет совпадать со всем пространством L2(G), т. е. инвариантное относительно сдвига подпространство является собственным тогда, когда оно не содержит 5.
Из перечисленных выше фактов следует: (а) каждое собственное инвариантное относительно сдвига подпространство содержится в максимальном инвариантном относительно сдвига подпространстве; (б) каждое инвариантное относительно сдвига подпространство может быть представлено (возможно, разными способами) в виде прямой суммы неприводимых инвариантных относительно сдвига подпространств, т. е. ненулевых подпространств, которые не могут быть нетривиальным образом разбиты на сумму двух меньших инвариантных относительно сдвига подпространств. Можно показать, что инвариантное относительно сдвига подпространство неприводимо тогда, когда его ортогональное дополнение максимально.
Замечание 1. Ортогональная проекция на инвариантное относительно сдвига подпространство является оператором, инвариантным относительно сдвига, и поэтому является сверткой с некоторой функцией в частности, ¡л * ¡л = Такие функции называются идемпотентными мерами. Они играли значительную роль в развитии гармонического анализа.
Предложение 1 (частный случай теоремы Гельфанда — Мазура [2]). Все максимальные инвариантные относительно сдвига подпространства являются гиперплоскостями (т. е. имеют коразмерность один).
Предложение 2. Все неприводимые инвариантные относительно сдвига подпространства одномерны.
Доказательство. Приведенные выше два предложения эквивалентны друг другу.
Пусть V — неприводимое инвариантное относительно сдвига подпространство, т. е. одномерное. Это означает, что каждый оператор сдвига Transy действует на V как умножение на комплексную константу xv(у). Так как Transy унитарный, получаем, что \xv(y)\ = 1. Также из (1) видим, что xv(у) гомоморфизм: xv(у + +z) = XV (y)xv (z). Другими словами, xv : G ^ S1 является мультипликативным характером G; обратно каждый мультипликативный характер х : G ^ S1 порождает одномерное подпространство, неприводимое, инвариантное относительно сдвига. Таким образом, неприводимые инвариантные относительно сдвига подпространства находятся во взаимно-однозначном соответствии с мультипликативными характерами. Используя экспоненциальную функцию e : R/Z ^ S1, e(x) := := e2mx, можно записать каждый мультипликативный характер х как х = e(£),
где £ : G ^ R/Z — аддитивный характер, т. е. аддитивный гомоморфизм из G в R/Z. Таким образом, неприводимые инвариантные относительно сдвига подпространства находятся также во взаимно-однозначном соответствии с аддитивными характерами. □
Замечание 2. Как следствие, имеем, что каждое максимальное инвариантное относительно сдвига подпространство является ортогональным дополнением мультипликативного характера х.
Определим группу G, дуальную к G (по Понтрягину), как пространство всех аддитивных характеров G; G образует аддитивную группу. Будем пользоваться обозначением £■x для £(x) € R/Z, x € G,£ € G. Элементы будем называть частотами. Для фиксированной частоты £ соответствующее неприводимое инвариантное относительно сдвига подпространство V, является линейной оболочкой мультипликативного характера e, : x ^ e(£ ■ x).
Лемма (ортогональность). Если £,п — две различные частоты, то соотвеству-ющие инвариантные относительно сдвига подпространства V, и Vn ортогональны.
Доказательство. Достаточно показать, что e, и en ортогональны, другими словами, надо показать, что выражение
I := e(£ ■ x)e(-tf ■ x) dx
Ja
равно нулю. Сдвигая x на y, видно, что
1 = e(£ ■ y)e(-n ■ y)I
для всех y € G. Но так как £, п различны, то существует y такой, что £ ■ y = п ■ У, и, следовательно, I = 0. □
Разбивая регулярные представления на неприводимые (а также ортогональные) компоненты, получаем
Следствие 1 (теорема Петера — Вейля [3], случай конечной абелевой группы). Имеет место равенство L2(G) = 0,е(з V,с, из которого следует, что #G = #G. Пространство {e, : £ € G} мультипликативных характеров образует ортонормиро-ванный базис в L2(G).
Заметим, если f € L2(G) и £ € G, проекция f на V, задается как f (£)e,, где f (£) = {f,eс)l2(g) = Jq f (x)e(-£ ■ x)dx. Таким образом, получаем
Следствие 2 (формула обращения Фурье). Для каждого f € L2(G) имеем
f = £f (£)eс.
Как еще одно следствие получается т. н. тождество Планшереля
\\f h4G) = иЬф)
и более общее тождество Парсеваля
if,g)L2(G) = {f,g )i2{G).
Имеем тождество свертки
/ * 9 = ¡9
и дуальное тождество
¡9 = I * 9,
где * справа означает дискретную свертку (использующую считающую меру на О вместо нормированной считающей меры на О).
Заметим, что каждый х € О может рассматриваться как характер х ^ £ • х
на О, таким образом получается каноническое отображение из О в О. Это отображение инъективно. Действительно, предположим, что х принадлежит ядру этого отображения, тогда £ • х = 0 для всех £ € О, или эквивалентно Тгап8; оставляет неподвижным каждый из характеров в^. В силу формулы обращения Фурье Тгап8; оставляет неподвижными все функции. Из теоремы Петера — Вейля известно, что мощности О и О совпадают. Таким образом, отображение биективно. Другими словами дуальная по Понтрягину к О канонически отождествляется с самой О.
Пример 2. Если О = Ъ/ЫZ, каждый £ € Ъ/ЫZ порождает характер по формуле £ • х := £х/Ы. Они образуют N различных характеров и, следовательно, по теореме Петера — Вейля, других характеров не существует. Таким образом, абстрактное преобразование Фурье совпадает с обычным конечным преобразованием Фурье на Ъ/ЫЪ.
Имеет место стандартное соотношение между сдвигом и модуляцией: для каждого f € L2(G) и y е G, £ € G
Transy f = Mod_y f; Mod f = Trans5 f ,
где Mod_yF(£) := e(-£ ■ y)F(£) и Mod5f (x) := e(£ ■ x)f (x).
Аналогами гауссовых функций на конечных абелевых группах являются индикаторные функции подгрупп.
Если H С G — подгруппа G, определим ортогональное дополнение HЧ С G как
HЧ := {£ е G : £ ■ x = 0 для всех x € H}. Имеет место формула суммирования Пуассона
Г H11
1н = щ н±
(в частности, преобразованием Фурье единицы является функция Дирака и наоборот). Отсюда и из формулы Планшереля получаем основное тождество
\H \x\H Ч = \G\.
Для конечной абелевой группы G принцип неопределенности Донохо —
Старка [4] принимает вид: для любой ненулевой функции f : G ^ C имеем \supp(f)\\supp(f)\ > \G\.
Доказательство. Объединить теорему Планшереля с неравенством Гельдера
\\f\\li(g) < \supp(f)\1/2\G\-1/2\\f\\L4G); \\f\\nG) < \supp(/)\1/2\\fy^)
и неравенство Римана — Лебега
IIf II,-(G) < f МО-
□
Можно показать, что равенство достигается на индикаторах 1н подгрупп с точностью до сдвига, модуляции и умножения на константы.
Функциональные аналоги: если f G L2(R), suppf С T и f С Q то \T\\Q\ ^ 1.
Если f G L1(R), suppf С T, suppf С Q и \T\\Q\ < то, то f = 0.
Для произвольных групп G и произвольных функций f невозможно получить лучшие оценки, чем приведенные выше принципы неопределенности, об этом говорят примеры с f = 1н, если H — подгруппа G.
С другой стороны, для групп и функций специального вида можно ожидать улучшенных оценок.
Например, для циклической группы G = Z/pZ простого порядка, которая не имеет нетривиальных подгрупп, имеем
Принцип неопределенности для Z/pZ [5]. Пусть p — простое число. Если f : Z/pZ ^ C нетривиальна, то
\supp(f )\ + \supp(f)\ > p +1.
Обратно, если A и B — непустые подмножества Z/pZ такие, что A + \B \ ^ p+ + 1, то существует функция f с supp(f) = A и supp(f) = B.
В [5] найдены связи этого принципа с результатом Чеботарева о том, что все миноры матрицы Фурье (e(x£/p))i^x^p невырождены.
Этот принцип неопределенности имеет приложения к арифметической комбинаторике, например, он влечет неравенство Коши — Давенпорта
\A + B\ > min(\A\ + \B\- 1,p)
для подмножеств A, B множества Z/pZ.
Доказательство. Применить принцип неопределенности к функциям вида f * g, где f расположена на A, g расположена на B, и supp(f ), supp(g ) выбираются так, чтобы они имели наименьшее пересечение.
□
Принцип неопределенности для Z/pZ можно проинтерпретировать следующим образом:
Если сигнал f — S-разреженный (т. е. \supp(f)\ < S ), а \Q\ > S, то f имеет на Q, по крайней мере, один ненулевой коэффициент.
С точки зрения цифровой обработки сигналов, произвольные S коэффициентов Фурье позволяют обнаружить присутствие S-разреженного сигнала.
Для составных p это неверно. Если, например, N является точным квадратом, то Z/NZ содержит подгруппу из VN элементов, и индикатор этой подгруппы (т. н. Дирак-комб сигнал), являясь %/N-разреженным, имеет N — %/N нулевых коэффициентов Фурье.
Еще одно следствие принципа неопределенности: если f — неизвестный разреженный сигнал, и удалось измерить 2S коэффициентов Фурье, то сигнал может быть точно восстановлен. В противном случае, если два S-разреженных сигнала f и g имеют совпадающие коэффициенты на Q, то 2S-разреженная разность f — g имеет нулевые коэффициенты на Q, что невозможно.
Утверждения такого типа характерны для сжатого зондирования (compressed sensing): возможность восстановления разреженного или сжатого сигнала, используя небольшое количество измерений, не имея сведений о расположении носителя сигнала (в стандартной ситуации для восстановления сигнала требуются все p
коэффициентов Фурье, разреженные сигналы имеют меньшую энтропию и, как следствие, восстанавливаются значительно меньшим количеством измерений, чем общий сигнал).
Практически сформулированные выше результаты не могут считаться удовлетворительными по двум причинам:
1) теоретически сигнал восстанавливается по 2S коэффициентам Фурье, однако если сигнал имеет S + 1 ненулевых компонент, то его уже не удается, вообще говоря, восстановить по 2S измерениям;
2) процедура восстановления не является робастной (устойчивой) к возмущениям: важна простота числа p, нет устойчивости и по отношению к малым возмущениям f.
Обе эти проблемы решаются, если множество частот Q удовлетворяет т. н. "принципу частичной изометрии" с параметрами S и S:
(1 - S) H||f\\14Z/NZ < WfWhn) < (1 + S) H||f \\14Z/NZ
для всех S-разреженных сигналов f.
Заметим, что множитель N согласован с теоремой Планшереля. Он показывает, что Q всегда схватывает "порядочную долю" энергии разреженного сигнала. Таким образом, Q не только обнаруживает присутствие S-разреженного сигнала, но "вылавливает" большую его часть.
"Принцип частичной изометрии" полезен в сжатом зондировании.
Теорема [6]. Если Q С Z/NZ удовлетворяет "принципу частичной изометрии" c параметрами 4S и S = 1/4, то любой S-разреженный сигнал f является единственным решением g : Z/NZ ^ C задачи g |q= f |q с минимальной 1\(С)-нормой. В частности, сигнал f может быть восстановлен по коэффициентам Фурье f |п с помощью выпуклой оптимизации.
Литература
[1] Grochenig K. Foundations of Time-Frequency Analysis. Boston; Basel; Berlin: Birkhauser. 2000. 360 p.
[2] Рудин У. Функциональный анализ. М.: Мир, 1975. 443 с.
[3] Понтрягин Л.С. Непрерывные группы. М.: Наука: Физматлит. 1973. 527 с.
[4] Donoho D.L., Stark P.B., Edidin D. Uncertainty principles and signal recovery // Journal Applied Mathematics (SIAM). 1989. V. 49. I. 3. P. 906-931.
[5] Tao T. An uncertainty principle for cyclic groups of prime order // Mathematical Research Letters. 2005. V. 12. P. 121-127.
[6] Candes E.J., Romberg J., Tao T. Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information // IEEE Trans. Inform. Theory. 2004. V. 52. P. 489-509.
References
[1] Grochenig K. Foundations of Time-Frequency Analysis. Boston-Basel-Berlin, Birkhauser, 2000, 360 pp. [in English].
[2] Rudin U. Functional analysis. M., Mir, 1975, 443 p. [in Russian].
[3] Pontryagin L.S. Continuous groups. M., Nauka. Fizmatlit, 1973, 527 p. [in Russian].
[4] Donoho D.L., Stark P.B., Edidin D. Uncertainty principles and signal recovery. SIAM Journal Applied Mathematics, 1989, V. 49, I. 3, pp. 906-931 [in English].
[5] Tao T. An uncertainty principle for cyclic groups of prime order. Mathematical Research Letters, 2005, Vol. 12, pp. 121-127 [in Russian].
[6] Candes E.J., Romberg J., Tao T. Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information. IEEE Trans. Inform. Theory, 2004, Vol. 52, pp. 489-509 [in Russian].
S.Y. Novikov, M.E. Fedina3
UNCERTAINTY PRINCIPLES FOR GROUPS AND RECONSTRUCTION OF SIGNALS4
Uncertainty principles of harmonic analysis and their analogues for finite abelian groups are considered in the paper. Special attention is paid to the recent results of T. Tao and coauthors about cyclic groups of prime order. It is shown, that indicator functions of subgroups of finite Abelian groups are analogues of Gaussian functions. Finite-dimensional version of Poisson summation formula is proved. Opportunities of application of these results for reconstruction of discrete signals with incomplete number of coefficients are suggested. The principle of partial isometric whereby we can determine the minimum number of measurements for stable recovery of the signal are formulated.
Key words: uncertainty principles, cyclic finite groups, reconstruction, sparse signal, indicator functions, Poisson formula.
Статья поступила в редакцию 28/V/2015. The article received 28/V/2015.
3Novikov Sergey Yakovlevich (nvks@samsu.ru), Department of Probability Theory and Mathematical Statistics, Samara State University, 1, Acad. Pavlov Street, Samara, 443011, Russian Federation.
Fedina Maria Efimovna (phedina75@gmail.com), Department of Security of Information Systems, Samara State University, 1, Acad. Pavlov Street, Samara, 443011, Russian Federation.
4The work was supported by the Russian Ministry of Education as part of the basic part of a state task, project № 204.