CC BY
24Физико-математические науки
УДК: 621.391.1
Кулешова А. А.
магистрант,
Самарский Национальный Исследовательский Университет имени Академика С. П. Королева» (Самарский университет)
ПОИСК АЛГОРИТМОВ ДЛЯ ВОССТАНОВЛЕНИЯ СИГНАЛА БЕЗ ФАЗЫ
Поиск алгоритмов для восстановления сигнала без фазы актуален в настоящее время. Алгоритмы восстановления важны в обработке разнообразных сигналов, в особенности в технологии распознавания речи, в томографии. Главное свойство фреймов, которое делает их настолько полезными в прикладных задачах — их избыточность. Хорошо выбранный фрейм может обеспечить численную устойчивость для восстановления сигнала и получение важных характеристик сигнала.
Ключевые слова: фрейм, восстановление сигнала без фаз, равномерные фреймы.
Дискретизация и квантование аналогового сигнала приводят к рассмотрению сигнала как элемента некоторого конечномерного пространства V. В таком пространстве, вообще говоря, комплексном, вводится комплексное скалярное произведение и соответствующая
эрмитова норма. По ортонормированному базису (ОНБ) {щ «сигнал» X е V единственным образом представляется суммой:
х = X (х, щ) щ.
I =1
Комплексные коэффициенты описываемого представления [X, и^ дают возможность полного восстановления сигнала и часто понимаются как «измерения» сигнала. Реальные измерения получаются вещественными, и зазор между (х, и1 ) и амплитудами измерений
X, щ
оказывается непреодолимым при восстановлении сигнала.
Представляя сигнал в различных базисах, можно получить о нем разностороннюю информацию. Так, переход от представления по ортам к представлению в базисе Фурье, позволяет получить частотные характеристики сигнала, дающие широкие возможности для его цифровой обработки.
Последние годы значительное количество работ посвящено решению следующей
задаче: построить такие системы «измерительных» векторов Е = , которые позволяют
восстановить произвольный сигнал X е V по набору вещественных чисел (х,
В классе ОНБ такая задача не имеет решения.
Основная проблема, поставленная в [13], до сих пор далека от окончательного решения:
Найти необходимые и достаточные условия на систему векторов представления Е = (т.н. «измерительных векторов»), которые обеспечивают инъективность и
устойчивость отображения «измерения амплитуды» сигнала х:
о
(А( X))(/ ):= (х,
Семейство векторов {}М^ называется фреймом гильбертового пространства Н^
если существуют такие константы 0 < A < B < ю, такие, что для всех x £ H выполняются следующие неравенства:
2
М
2
2
АЩ < X (х,/ ) < В||х|
г =1
A и B называются границами фрейма. Наибольшая из нижних границ называется оптимальной нижней границей, а наименьшая из верхних границ — оптимальной верхней границей. Если A = B, то фрейм называется А-жестким, а если A = B = 1, то фреймом Парсеваля-Отеклова.
Числа X,/г )}г_1 называются фреймовыми коэффициентами.
Пусть {}М=1 — фрейм, линейное отображение:
Т: НЩ ^ НМ = 12(I), Т(х) = {(х,/ )}М=1
называется оператором анализа. Линейное отображение:
Т * : НМ = 12(!) ^ НЩ, Т *(с) =
г =1
называется оператором синтеза.
Композиция Т и Т определяет фреймовый оператор — положительный, самосопряженный обратимый оператор:
5 = Т*Т : НЩ ^ НЩ : 5х = Т*Тх = £(х,/ )fi .
г =1
Он обеспечивает точную формулу для восстановления:
м
X = М(х, /> "/ .
г =1
Семейство векторов Е = {/ }М=Х является равномерным равноугольным жёстким фреймом, если
1) 3Р> 0: || / || =р V/ = 1,М ;
2) 3с > 0 : для любой пары векторов фрейма / ■ и / , ^к, фрейма мы имеем:
/], Ь) = с
Известно, что есть верхняя граница для числа векторов в равномерном равноугольном жёстком фрейме Е = {/ }М=Х на N-мерном Гильбертовом пространстве Н. В вещественном
^ +1) 2 случае это М <-, в комплексном случае — М < N2([10], [3]). Построение
максимального числа векторов для равномерного равноугольного жесткого фрейма очень сложная и нерешенная задача в теории фреймов.
Рассмотрим Р - нелинейное отображение, переводящее вектор в набор модулей фреймовых коэффициентов:
с л М
Р: Н ^ 12(1 ),Р(х) = \ X/
г =1
Множество От (N, М; Я) — это множество К-мерных линейных подпространств в
ЯМ , которое имеет структуру N(М — N) — мерного множества.
Множество От (N, М; Я) называется многообразием Грассмана.
Фрейм {}г=1 называется фреймом общего положения, если {}М=^ С Ь еП, где и — открытое по Зарисскому множество и П С От (N, М) . Теорема 1: (Вещественный случай) [13] Если М > 2N — 1, то для фрейма общего положения Е = {/1 }М, нелинейное отображение Р инъективно.
Множество {/^ }М С ЯN называется набором с полным спарком, если каждое его
О N
подмножество из N векторов полно в Я .
Для анализа инъективности важным является свойство альтернативной полноты.
Набор векторов Е = {fi в ЯN (С^) назовем альтернативно полным, если для
любого р с {1, . . . ,М} , либо {fi } ¡ер, либо {fi } . с полно в Яп (С* ).
Пусть Н / ~ — фактор-пространство, полученное отождествлением двух векторов, если они отличаются постоянным фазовым коэффициентом. Таким образом, х~у означает, что существует константа с:|с|=1, такая что у = сх.
Для вещественных Гильбертовых пространств с = + 1 , тогда Нт = Н /{+1}.
Для комплексных гильбертовых пространств с = е , тогда Нг = Н/{Т,
где Т1 - окружность единичного радиуса на комплексной плоскости. В квантовой механике эти проективные лучи определяют квантовые состояния [12].
Теорема 2: Пусть Е = {fi }М=1 — набор векторов в ЯN. Отображение А: Яп /{+1}
^ Я определено равенствами: (А(х))(I) := (х, fЛ , I = 1,..,М. Тогда отображение А инъективно тогда и только тогда, когда Е - альтернативно полно.
Теорема 3: Всякий набор с полным спарком Е = {fi }М1 в ЯN с М > 2N — 1 удовлетворяет свойству альтернативной полноты.
Теорема 4: В вещественном случае, если Е = {fi }М1 в и М < 2N — 2, то отображение А : RN /{+1} ^ ЯМ не является инъективным. Если М = 2N — 1, то
отображение А инъективно тогда и только тогда, когда Е = {fi }М1 — полный спарк.
Доказательство: Если М < 2N — 2, то множество {1,..,М} можно разбить на
с
множества р и р так, чтобы мощность каждого не превосходила N — 1. Ни одно из множеств {fi }iер, {fi } . с не может быть полным.
Если М = 2N — 1 и Е = {fi }М^ — полный спарк, то инъективность А следует из леммы 2 и теоремы 3.
Обратно, если А — инъективно, то Е = является альтернативно полным
семейством. Возьмем произвольное подмножество (С {1,.., М} с ( = N. Тогда
(рС = N — 1 и {fi }. с не может быть полным. Следовательно, {полно, и
Е = Iм — полный спарк.
Неизвестна точная минимальная граница для комплексного случая. Также, в вещественном случае существует простой прямой метод для проверки инъективности отображения А для соответствующего фрейма [4].
Следствие 1: Если Б это М-элементный фрейм в Я1^ с М > 2N — 1 и у Е = {fi
каждые К-элементов фрейма линейно независимы, то оператор А : Н^г ^ НМ = 12 (I) — инъективен.
Теорема 5: (Комплексный случай) Если М > 4К-2, то для фрейма общего положения Б, отображение А инъективно.
(4M - 4)-гипотеза: (Комплексный случай) Если М > 4К-2, то для фрейма общего положения Б, отображение А инъективно.
Пусть Е = {fi }М с С N и отображение А : CN /{Т1} ^ ЯМ определено соотношениями (А(л))(I) := (л,^ \ , I = 1,..,М. Если М > 2, то:
(а) Если M < 4 N — 4 , то отображение A не является инъективным (опровергнуто для N=4)
(б) Если M > 4N — 4, то отображение A может быть инъективным для некоторых фреймов.
В [1] показывается, что фрейм общего положения будет восстанавливать сигнал без фазы за полиномиальное число шагов.
Пусть H — фиксированное N-мерное векторное пространство (вещественное или комплексное), пусть } - ортонормированный базис для H.
Теорема 6:
(а) Если H е RN, M > —-) и F = If }MJ- фрейм общего положения,
нелинейное отображение P — инъективно. Тогда вектор x е H может быть восстановлен (с точностью до знака) из множества {| (x, f^ |} M1 модулей фреймовых коэффициентов за
полиномиальное число шагов (O(N6 )) .
(б) Если H е CN, M > N2 и F = Ifi }Mj — фрейм общего положения, нелинейное отображение P — инъективно. Тогда вектор x е H может быть восстановлен (с точностью умножения на корень из единицы) из множества {| (x, f^j |} M1 модулей фреймовых
коэффициентов за полиномиальное число шагов (O(N6 )) .
Литература
1. R. Balan, B. G. Bodmann, P. G. Casazza, D. Edidin, Fast algorithems for signal reconstruction without phase, Proceedings of SPIE-Wavelets XII, San Diego 67G1 (2GG7) 67G11192G-67G111932.
2. R. Balan. Equivalence relations and distances between Hilbert frames. Proc. Amer. Math. Soc., 127(8) : 2353-23бб, 1999.
3. R. Balan, B. G. Bodman, P. G. Casazza and D. Edidin, Painless reconstruction from magnitudes of frame coefficients, preprint.
4. R. Balan, P. Casazza, D. Edidin, On signal reconstruction without phase, Appl.Comput.Harmon.Anal. 20 (2006), 345-356.
5. A. R. Calderbank, P. J. Cameron, W. M. Kantor, and J. J. Seidel, Z4-Kerdock codes, orthogonal spreads, and extremal Euclidean line-sets, Proc. London Math. Soc. (3) 75 (1997), no. 2, 436-480.
6. P. J. Cameron and J. J. Seidel, Quadratic forms over GF (2), Indag. Math. 35 (1973), 1-8.
7. P. G. Casazza and M. Fickus, Fourier transforms of finite chirps, EURASIP J. Appl. Signal Process. (2006), Frames and overcomplete representations in signal processing, communications, and information theory, Art. ID 70204, 1-7.
8. P. Delsarte, J. M. Goethals, and J. J. Seidel, Spherical codes and designs, Geometriae Dedicata 6 (1977), no. 3, 363-388.
9. D. Han and D. Larson. Frames, bases and group representations, Memoirs American Math. Soc. 147 (2000), no. 697.
10. R. Holmes and V. I. Paulsen, Optimal frames for erasures, Lin. Alg. Appl. 377 (2004),
31-51.
11. S. D. Howard, A. R. Calderbank, and W. Moran, The finite Heisenberg-Weyl groups in radar and communications, EURASIP J. Appl. Signal Process. (2006), no. Frames and overcomplete representations in signal processing, communications, and information theory, Art. ID 85685, 1-12.
12. R. F. Streater and A. S. Wightman. PCT, Spin and Statistics and All That. Princeton University Press, Landmarks in Mathematics and Physics, 2000.
13. A. Bandeira, J. Cahill, D. Mixon, A. Nelson. Saving phase: Injectivity and stability for phase retrieval. Applied and Computational Harmonic Analysis (ACHA). 2014. V. 37. I. 1. PP. 106-125.
