Вестник Самарского университета. Естественнонаучная серия. 2016. № 3-4
УДК 512.8, 517.98, 621.391
63
С.Я. Новиков, М.Е. Федина1
ВОССТАНОВЛЕНИЕ СИГНАЛА ПО МОДУЛЯМ
ИЗМЕРЕНИЙ2
Рассмотрены вопросы выбора векторов в конечномерных вещественных и комплексных пространствах так, чтобы по модулям измерений скалярных произведений этих векторов с неизвестным вектором оказалось возможным восстановить его.
Исследуется инъективность нелинейных отображений.
Ключевые слова: пространство, сигнал, системы векторов, инъектив-ность отображения, скалярные произведения.
1. Введение
Пусть V — пространство с комплексным скалярным произведением и орто-нормированным базисом (ОНБ) {иг}М=1. Каждый "сигнал" V £ V единственным образом представляется суммой
Комплексные коэффициенты представления {V, щ) дают возможность полного восстановления сигнала и часто понимаются как "измерения" сигнала. Реальные измерения получаются вещественными, и зазор между {V,щ) и \^,щ)\ оказывается непреодолимым при восстановлении сигнала.
Представляя сигнал в различных базисах, мы получаем разностороннюю информацию о нем. Так, переходя от представления по ортам в СМ к представлению в базисе Фурье, получаем частотные характеристики сигнала, дающие широкие возможности для его цифровой обработки.
Последние годы значительное количество работ [1, 2, 3, 4] посвящено решению такой задачи:
Существуют ли системы Ф = такие, для которых возможно восстано-
вить V по набору чисел
В классе ОНБ такая задача неразрешима.
!© Новиков С.Я., Федина М.Е., 2016
Новиков Сергей Яковлевич (nvks@ssau.ru), кафедра теории вероятностей и математической статистики, Самарский университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
Федина Мария Ефимовна (phedina@ssau.ru), кафедра безопасности информационных систем, Самарский университет, 443011, Российская Федерация, г. Самара, ул. Акад. Павлова, 1.
2Работа выполнена при финансовой поддержке Министерства образования и науки РФ в рамках базовой части государственного задания, проект № 204.
М
г=1
В англоязычной литературе сформулированная выше задача и связанные с ней вопросы составили раздел прикладных исследований под названием "PHASE RETRIEVAL (возвращение, воспроизведение фазы)".
Пионерской работой в этом направлении оказалась работа [1], в которой была доказана теоретическая возможность точного восстановления сигнала (с точностью до унимодулярного множителя), если в качестве системы представления использовать избыточные системы.
Дальнейшее развитие этого направления идет по двум ветвям:
1) поиск таких систем "измерительных" векторов Ф = {фг}гЕ1, которые позволят восстановить произвольный сигнал v £ V по набору вещественных чисел
2) обоснование возможности восстановления произвольного сигнала с большой вероятностью для относительно небольшого числа "случайно выбранных" измерительных векторов.
Второй подход весьма близок теории сжатого зондирования (в английских статьях "Compressed sensing", перевод предложен Б.С. Кашиным).
В данной работе будут изложены результаты, полученные в первом направлении, и сформулированы некоторые нерешенные задачи. То, что касается второго, будет предметом отдельной статьи.
Основная проблема, поставленная в [3], до сих пор далека от окончательного решения:
Найти необходимые и достаточные условия на систему векторов представления Ф = {фi}N=1 (т.н. "измерительных векторов"), которые обеспечивают инъ-ективность и устойчивость отображения "измерения амплитуды" сигнала x
(A(x))(n) := \{х,фп)\\
Естественной областью определения этого отображения является пространство V = RM или CM. Заметим, однако, что A(x) = A(y), если y = cx для некоторого скаляра c с единичным модулем. Поэтому, строго говоря, отображение A : V ^ RN не может быть инъективным. В дальнейшем изложении будем считать областью определения A либо RM/{±1}, либо CM/T, где T — окружность единичного радиуса на комплексной плоскости. При таком соглашении вопросы об инъективности и устойчивости отображения A являются содержательными.
2. Инъективность в RM
Восстановить сигнал x невозможно, если отображение A не является инъек-тивным.
Следующее определение конкретизирует определение инъективности применительно к поставленной задаче.
Определение 1. Набор векторов Ф = {фi}N=i в RM (или CM) обеспечивает воспроизведение фазы (ВФ), если для любых x,y £ RM (или CM), таких, что \{x, фг)\ = \(y, фг)\ для всех i = 1,..., N, имеем x = cy, где c = ±1 для RM (и c £ T для CM, где T — единичная окружность на комплексной плоскости).
Важным для анализа инъективности оказалось следующее свойство, которое мы назвали свойством альтернативной полноты, заменяя тем самым английское сочетание "complement property".
Определение 2. Набор векторов Ф = {фп}1^,= 1 в ММ (СМ) назовем альтернативно полным (АП) если для любого Б С {1,... ,К} либо {фп}пев, либо {фп}пева полно в ММ (СМ).
Для пространства ШМ свойство альтернативной полноты эквивалентно инъек-тивности отображения А.
Теорема 1. Пусть Ф = {фп}^=1 С ММ и отображение А : ЕМ/{±1} ^ определено равенствами (А(х))(п) := \{х,фп)\2, п = 1,...,М. Тогда А инъективно тогда только тогда, когда Ф альтернативно полно.
Доказательство.
Необходимость. Предположим, что Ф £ (АП). Следовательно, найдется Б С {1, ...,М} такое, что ни {фп}п^Б, ни {фп}п^Бс не полно в ММ. Выбираем ненулевые векторы и, и £ ММ так, что {и,фп) = 0 для всех п £ Б и {и,фп) = 0 для всех п £ Бс. Для каждого п имеем
\{и ± и, фп)\2 = \ {и,фп)\2 ± 2{и,фп){и,фп) + \ {и,фп)\2 = \{и, фп)\2 + \ {и,фп)\2.
Отсюда следует, что \{и + и,фп)\2 = \{и — и,фп)\2 для каждого п, и А(и + + и) = А(и — и). Более того, так как и и и ненулевые, по предположению, то и и + и = ±(и — и). Таким образом, инъективности отображения А нет.
Достаточность. Предположим, что А не инъективно. Это означает, что существуют векторы х,у £ ММ такие, что х = ±у и А(х) = А(у). Обозначим Б := = {п : {х,фп) = —{у,фп)}. Имеем: {х + у,фп) = 0 для каждого п £ Б, и {х —
— у, фп) = 0 для каждого п £ Бс. По предположению, х = ±у, поэтому х + у = 0 и х — у = 0. Таким образом, и {фп}п^Б, и {фп}пеБс не полны в ММ. •
Таким образом оказалось, что система, которая допускает восстановление фазы, с необходимостью оказывается фреймом пространства ШМ. Из дальнейшего станет ясным, что аналогичная картина наблюдается и в СМ.
Напомним, что система {фп}^=\ называется фреймом (каркасом) пространства ММ (СМ), если существуют константы 0 < А < В такие, что Ух £ ММ (Ух £ СМ)
N
А||х||2 ^\{х,фг)\2 < В\\х\\2.
г=1
В конечномерном пространстве понятие фрейма эквивалентно понятию полноты системы, т. е. равенству врап{фп}^^=1 = КМ (вра^фп}^^ = СМ).
3. Инъективность в См
Подробнее теория фреймов изложена в [5, 6]. В пространстве СМ не найдено практически применимой характеризации инъективности.
Прежде чем доказывать известный ныне критерий (теорема 4 ниже), рассмотрим следующее вспомогательное утверждение. Оно позволяет увидеть аналогию между теоремами 2 и 4. Если а = (а1,..., аМ)Т £ СМ и Ь = (Ь1,..., Ьм)т £ СМ, то скалярное (внутреннее) произведение {а,Ь) можно представить в виде Ь*а, где Ь*
— сопряженный вектор-строка Ь* = (Ь1,... ,Ьм). Внешнее произведение аЬ* дает М х М-матрицу с элементами
а1 Ь1 а1 Ь2 . . . а1 ЬМ а2Ь1 а2Ь2 .. . а2Ьм
\ ам Ь1 аМ Ь2 . . . аМ ЬМ )
Лемма 1. Пусть Ф = {фп}^\ С RM. Рассмотрим множество {фпф*пи^=1 векторов в RM. Альтернативная полнота Ф эквивалентна равенству 8рап{фпфПм^=1 = RM для любого и £ RM \ {0}.
Доказательство. Очевидно, что {фпф*пи^=1 = {{и,фп)фп}^1. Необходимость. Предположим противное. Пусть существует вектор и £ RM \ {0} такой, что
span{(u, фп= RM. (**)
Пусть S1 = {n : (и, фп) = 0}. Ясно, что
8рап{фп}пей1 = RM. (*)
Рассмотрим два случая:
1) число элементов в S1 больше, чем N — M. Тогда число элементов в SC меньше M. Вместе со (*) последнее противоречит АП Ф;
2) число элементов в S1 меньше или равно N — M. Из АП Ф и (*) следует,
что
span^^es^ = RM. (*)
Однако это противоречит (**).
Достаточность. Предположим противное. Пусть существует S1 С {1,..., N} такое, что
span^^esi = RM и span^^esj = RM. (* * *)
Выберем и ортогональным 8рап{^„}„ея1 • Тогда span{(u, фп)^}^,^ = RM• Однако это противоречит (***). •
Теорема 2. Пусть Ф = {фп}^\_ С CM и отображение A : CM/T ^ RN определено соотношениями (A(x))(n) := Kx^n)!2, n = 1,...,N. Векторы {фпф*пи}^1 будем рассматривать как векторы пространства R2M. Пусть S(и) := spanR{фnф*nu}N=1. Эквивалентны следующие утверждения:
(а) A инъективно.
(б) dim S(u) > 2M — 1 для каждого и £ CM \ {0}.
(в) S(u) = spanR{rn}^ для каждого и £ CM \ {0}.
Доказательство. (а)^(в): Предположим, что отображение A инъективно. Требуется доказать, что линейная оболочка векторов {ф^п^^г, является ортогональным дополнением вектора in. Ортогональность в данном случае понимается в смысле вещественного скалярного произведения (a,b)-R = Re (a,b). Заметим, что
\(и ± у,фп)\ = \{и,фп)\ ± 2Re(u,фn)(фn,v) + ,
поэтому
\(и + у,фи)\2 - \(и - у,фи)\2 = 4Re(u, фи){фи, V) =
= 4ае ф*пи(фп,ги) = 4(фпф*пи,ю)ж^ (1)
Если V £ Б(и)^, то правая часть (1) равна нулю. В силу инъективности А получаем существование комплексного числа ш, ^ = 1 такого, что и + V = ш(и — — v)• Так как и = 0, то ш = — 1 Имеем
_ 1 — ш _ (1 — ш)(1 + ш) _ 21тш .
V — ~ и — т~ г^ и — т~ ¡77 %и,
1+ ш \1 + ш\2 \1+ ш\2
т.е. Б(и)^ С врапк{ш}
Для доказательства включения spanR{rn} С S(u)^ возьмем v = aiu с некото-
1+ai 1 — ai
рым а £ М, положим ш := и заметим, что |ш| = 1. Имеем:
u + v = u + aiu = (1 + ai)u = —+—i (u — aiu) = w(u — v).
1 — —i
Следовательно, левая часть (1) равна нулю, и v £ S(u)^. Заметим, что в доказательстве последнего включения инъективность A не требовалась. Это замечание будет использоваться в дальнейшем.
(б)^ (в): из (в) сразу следует (б). Для доказательства обратной импликации заметим сначала, что iu ортогонален каждому из фп ф*п'а :
ltfn4>*nu, iu)R = Re{^n^*nu, iu) = Re(u, фп){фи, iu) = —Re il(u^n)l2 = 0.
Это означает, что spanR{iu} С S(u)^. В силу (б), dim S(u)^ ^ 1. Из последних двух соотношений получаем (в).
(в)^ (а): предположим, что A(x) = фПх12 = 1фПу12 = A(y). Если x = y, инъективность получена. Если x — y = 0, то применим (в) к u = x — y. Следующие ниже преобразования покажут нам, что x + y £ S(x — y)^. Действительно,
(VnV*n(x — y),x + y)R = Re{<fn<f*n(x — y),x + y) = Re(x + y)*^n^*ri(x — y) =
= Re {\фПx]l2 — x*<fin<fi*ny + y*<fn<f*nx — Wry?) = = Re( x* ФпфПу + x*ФпФПу) = 0.
Итак, x + y £ S(x — y)^ = spanR{i(x — y)}, то есть существует a £ М такое, что x + y = ai(x — y), отсюда y = x, что означает x = y modT. •
Следующая теорема показывает, что инъективность отображения "измерений амплитуды" (A(x))(n) = l(x, фп)12, n = 1,..., N в комплексном пространстве эквивалентна в широком смысле (с точностью до вещественного множителя) инъек-тивности другого отображения
(B(x))(n) := arg ((x, Фп)2), n = 1,..., N,
которое измеряет только фазы восстанавливаемого сигнала, то есть является отображением "измерений фазы". Будем считать, по определению, что arg ((x^n)2) = = 0, если l(x^n)l = 0.
Теорема 3. Пусть Ф = {фп}^=1 С CM и отображение A : CM/T ^ MN определено соотношениями (A(x))(n) := l(x^n)l2, n = 1,...,N. Отображение A инъек-тивно тогда и только тогда, когда отображение (B(x))(n) := arg ((x,фп)2') , n = = 1,...,N инъективно по модулю М \ {0}, то есть B(x) = B(y) ^ x = y modМ \ {0}.
Доказательство. По теореме 2 инъективность A эквивалентна следующему утверждению:
Vx £ CM \{0}, spanЛ{фnф*nX}N=l =spanR{ix}±. (2)
Как было отмечено в доказательстве теоремы 2, включение spanR{ix} С (spanR^^nx}N=1) справедливо для любого x £ CM \ {0}. Поэтому (2) эквивалентно следующей импликации:
Vx £ CM \ {0}, iy) =0 Vn =1,...,N =^ y = 0
или y = x modМ \ {0}. (3)
Докажем, что Ке(ф„ф„х, iy) = Im{x, ф„)(у, ф„) = 0 тогда и только тогда, когда arg ((х,ф„)2) = arg (у(у,ф„)2) . После этого оказывается возможным заменить левую часть импликации (3) на равенство arg ((х,ф„)2) = arg {(у,ф„)2) и теорема будет доказана. Положим для краткости a := (х,ф„) и b := (у,ф„). Теперь остается доказать, что Im ab = 0 тогда и только тогда, когда имеет место хотя бы одно из трех равенств: arg(a2) = arg(b2), a = 0 или b = 0.
Если a или b равно нулю, равенство очевидно. В оставшихся случаях имеем 2arg(a) = arg(a2) = arg(b2) = 2arg(b), отсюда 2(arg(a) — arg(b)) кратно 2п, и arg(ab) = arg(a) — arg(b) кратно п. Следовательно, ab € R и Im alb = 0.
Пусть Im ab = 0. Запишем полярное представление чисел a = твгв и b = = se1^. Из условия, rs вт(в—ф) = 0. Такое возможно когда либо т, либо s равны нулю. Если ни одна из этих возможностей не имеет места, то sin в cos ф = cos в sin ф. Отсюда получаем, что если в кратно п/2, то и ф кратно п/2. Поэтому arg(a2) = = 2arg(a) = п = 2arg(b) = arg(b2). В оставшемся случае делим обе части равенства на cos в cos ф и получаем равенство tan в = tan ф, из которого вытекает в = ф mod п и arg(a2) = 2arg(a) = 2arg(b) = arg(b2). •
Лемма 2. Пусть Ф = {^„}N=1 С CM и отображение A : CM/T ^ RN определено соотношениями (A(x))(n) := \(х,ф„)\2, п = 1,...,N. Если A инъективно, то отображение C : CM/{±1} ^ CN, определенное равенствами (C(х)) (n) := = (х, ф„)2, n = 1,... ,N также инъективно.
Доказательство. Пусть A инъективно. Справедливы следующие утверждения (первое по определению, второе — по теореме 3):
1) если Уп = 1,..., N, \(х, ф„)\2 = \(у, ф„)\2, то у = х mod T;
2) если Уп = 1,..., N, arg(^, ф„)2) = arg(^, ф„)2), то у = х mod R \ {0}.
Пусть (х, ф„)2 = (х, ф„)2 для всех п = 1,..., N. Выполнены условия утверждений (1) и (2). Из равенства у = х mod T вытекает, что х = 0 тогда и только тогда, когда у = 0. Если х = 0 и у = 0, то существует ш € T П R \ {0} = {±1} такое, что у = шх. В любом случае у = х mod{±1}, то есть C инъективно. •
Теорема 4. Пусть Ф = {фп}^,=1 С См и отображение А : См/Т ^ RN определено соотношениями (А(х))(п) := \(х,фп)\2, п = Если А инъективно, то Ф удовлетворяет свойству альтернативной полноты.
Интересно сравнить теорему 4 с теоремой 2 и понять, почему в комплексном случае не проходит доказательство первой части теоремы 2. Для комплексного поля равенство теоремы 2 принимает вид
\(и ± v,фn)\ = \(и,фп)\ ± 2Re(u,фn)(v,фn) + \(и,фп)\ = \(и,фп)\ + \(ю,фп)\ •
В проведенном доказательстве теоремы 2 было показано, что отрицание свойства (АП) приводит к равенству A(u+v) = А(и—V) и при этом u+v = ±(u—v), что противоречит инъективности А на RM/{±1}^ Однако, может иметь место, например, такое равенство и + V = г(и — V), и инъективность А на См/Т не доказана.
Доказательство теоремы 4. В силу леммы 2, инъективность А влечет инъективность С на См/{±1} • Поэтому достаточно доказать, что если С инъективно на См/{±1}, то Ф £ (АП) Будем рассуждать как в доказательстве теоремы 2. Предположим что Ф £ (АП) Тогда существует Б С {1,_,Ж} такое, что ни {фп}п^з, ни {фп}п^зс не полны в См• Следовательно, найдутся ненулевые векторы и, V £ См такие, что (и, фп) = 0 для всех п £ Б и (V, фп) = 0 для всех
n £ Sc. Для каждого n имеем
(и ± У,фп)2 = (и,фп)2 ± 2(и,фп)(у,фп) + (v, фп)2 = (и, фп) 2 + (v, фп)2.
Так как (и + ги,фп)2 = (и — ги,фп)2 для каждого n, то С(и + v) = С(и — v). Однако, так как и и v ненулевые, по предположению, то u+v = ±(u—v), что противоречит инъективности С на CM/{±1}. •
Свойство альтернативной полноты недостаточно для инъективности отображения в пространстве CM. Достаточно рассмотреть векторы (1, 0), (0, 1) и (1, 1). Они удовлетворяют (АП), но A((1,i)) = (1,1,2) = A((1; —i)), хотя (1,i) ф (1, — —г) mod T; более того, векторы с вещественными координатами вообще не могут давать инъективные измерения в комплексном пространстве, так как они не различают векторы с комплексно сопряженными координатами.
4. Системы измерительных векторов
Поиск хороших достаточных условий для инъективности в комплексном пространстве остается актуальным. В вещественном пространстве найдены примеры хороших "измерительных векторов", например, так называемые "полные спарки": множество {фп}^=1 С Rм называется "полным спарком", если каждое его подмножество из М векторов полно в Rм •
Лемма 3. Всякий полный спарк {фп}^=\ в RM с N ^ 2М — 1 удовлетворяет свойству альтернативной полноты.
Доказательство. Предположим противное: существует Б С {1, 2, • • •, N} такое, что ни {фп}п£в, ни {фп}п£вс не полны в RM• По определению полного спарка, отсюда следует, что \Б \ < М и \Бс\ < М, то есть N < 2М — 1, что противоречит условию. •
В силу теоремы 2 всякий полный спарк представляет собой хорошее множество для инъективности отображения в вещественном пространстве. С другой стороны, найдены эффективные методы построения полных спарков [7]. Оставаясь в рамках детерминированной модели, естественно возникает вопрос об определении для пространства размерности М наименьшего числа N*(М) "измерительных векторов", обеспечивающих инъективность отображения А•
Для Rм число N*(М) может быть найдено с использованием теоремы 2.
Предложение 1. Если Ф = {фп}^=\ в RM и N < 2М — 2, то отображение А : RM/{±1}^ RN не является инъективным. Если N = 2М — 1, то отображение А инъективно тогда и только тогда, когда Ф — полный спарк.
Доказательство. Если N ^ 2М — 2, то множество {1, 2, • • •, N} можно разбить на подмножества Б и Бс так, чтобы мощность каждого не превосходила М — 1 Ни одно из множеств {фп}пез, {фп}пеза не может быть полным.
Если N = 2М — 1, и Ф = {фп}^=\ — полный спарк, то инъективность А следует из леммы 3 и теоремы 2. Обратно, если А инъективно, то Ф = {фп}^=\ является альтернативно полным семейством. Возьмем произвольное подмножество Б С {1, 2, • • •, N} с \Б\ = М• Тогда \Бс\ = М — 1 и {фп}п^зс не может быть полным. Следовательно, {фп}пев полно, и Ф — полный спарк. •
Таким образом, N*(М) = 2М — 1 для RM•
Переходим к рассмотрению ситуации в См •
В этом случае поставленный вопрос имеет историю. Весьма долго представлялась правдоподобной гипотеза о равенстве N*(M) = 3M — 2. Эта гипотеза была опровергнута в [8], в этой работе было доказано, что N*(M) ^ 4M — 2а(М — 1) — 3, где a(M — 1) ^ log*(M) обозначает количество единиц в бинарном представлении числа M — 1 .
С другой стороны, в [1] доказано, что N*(M) ^ 4M — 2.
В настоящее время правдоподобной представляется
(4M — 4)-гипотеза. Пусть Ф = {фп^=1 Ç CM и отображение А : CM/T ^ RN
определено соотношениями (А(ж))(п) := \(х,фп)\2, n = 1,...,N.
Если M ^ 2, то:
1) если N < 4M — 4, то А не является инъективным;
2) если N ^ 4M—4, то А может быть инъективным для некоторых фреймов.
Ниже будет дано обоснование этой гипотезы.
Определим M -мерное пространство HM xM комплексных самосопряженных M х M матриц. Это пространство над полем вещественных чисел. Для заданного набора {фп|N=i Ç CM определим оператор матричного анализа A : HMxM ^ RN равенствами (AH)(n) = (Н,фпфП)нБ, здесь (-,cdot)HS обозначает скалярное произведение Гильберта-Шмидта, индуцирующее матричную норму Фробениуса [9]. Подробнее,
M
(H,G)hs = J2 hmngmn = Tr[G*H].
m,n=1
Заметим, что A — линейный оператор, и
(Axx*)(n) = (xx*,VnV*n)HS = ЩфпфПххх*] = = Tr[f*nxx* фп] = vnxx*Pn = \ {x,^n)\2 = (A(x))(n).
Полученное равенство показывает, что x mod T может быть расширен до xx*, при таком расширении процесс измерения амплитуды становится линейным за счет увеличения размерности пространства, в котором проводятся измерения.
Посмотрим, как отражается инъективность при таком расширении.
Лемма 4. Отображение А не является инъективным тогда и только тогда, когда в ядре оператора A существует матрица ранга 1 или 2.
Доказательство. Необходимость. Предположим, что А не инъективно, т. е. существуют x,y G CM/T такие, что x = y и А(x) = А(у). Отсюда следует, что Axx* = Ayy*, и xx* — уу* G ker A. Ранг матрицы xx* — уу* равен 1 или 2 в зависимости от зависимости (независимости) векторов x и у.
Достаточность. Сначала предположим, что ядро оператора A содержит матрицу H ранга 1. Тогда существует x G CM такой, что H = xx* и (А^))(п) = = (Axx*)(n) = 0 = (А(0))(п). Но x ф 0modT, и А оказывается не инъективным.
Теперь предположим существование (самосопряженной) матрицы H ранга 2 в ядре оператора A. По спектральной теореме существуют ортонормированные векторы ui, u* G CM и ненулевые вещественные числа Ai ^ A* такие, что H = = Aiuiu* + A*u*u*. Так как H G ker A, то для каждого n имеем:
0 = (H, фпфп)нs = {Aiuiu* + A*u*u*, фпфп)нs =
= Ai\(ui^n)\2 + A*\(u*, фп)\*. (4)
Положим x := |Ai|l/2ui и у := |A2|l/2u2, заметим, что уф x mod T, так как они ненулевые и ортогональные. Для завершения доказательства покажем, что A(x) = = A(y). Если A1 и A2 имеют одинаковые знаки, то по (4), \(х,фп)\'2 + \(у,фп)\2 =0 для каждого n, следовательно, \(х,фп)\'2 =0= \(у,фп)\2.
Если Ai > 0 > A2, то хх* — уу* = Aiuiu* + A2u2u* = H £ ker A, следовательно, A(x) = Axx* = Ауу* = A^). •
Теорема 5. (4М — 4)-гипотеза верна для М = 2^
Доказательство. (а) Пусть N < 4М — 4 = 4^ Так как А определяет линейное отображение из 4-мерного вещественного пространства в ^мерное вещественное пространство, ядро А с необходимостью нетривиально, так как, по теореме о ранге и дефекте,
ё1ш1ш(А) + ё1шкег(А) = 4^
Кроме того каждый ненулевой элемент ядра имеет ранг 1 или 2, и, по лемме 4, А не будет инъективным.
(б) Рассмотрим матрицу, образованную 16 вещественными переменными:
Ф^)
xi + ix2 x5 + ix<6 x<9 + ixio xi3 + ix 14 x3 + ix4 x7 + ixg xil + ixi2 xi5 + ixi6
(5)
Если обозначить п-й столбец Ф(х) через фп(х), то получаем, что А инъектив-но именно тогда, когда х £ R16 порождает базис {фп(х)фп(х)*}п=1 в пространстве 2 х 2 самосопряженных матриц. В самом деле, в этом случае гг* однозначно определяется через Агг* = {(гг*, фп(х)фп(х)*)не}п=1 = А(г), что определит и г с точностью до унимодулярного множителя. Пусть А(х) обозначает 4х 4 матричное представление оператора матричного анализа, в котором п-я строка состоит из координат фп(х)фп(х)* в некотором базисе пространства Н2х2, например,
10 00 1 01 1 0i
01 , 01 , 72 10 , 72 —i 0
(6)
Тогда V = {х : Re det А(х) = 1ш detA(x) = 0} оказывается вещественным алгебраическим многообразием в R16, и замечаем, что А инъективно, если х £ Vс• Так как Vс открыто в топологии Зарисского, оно либо пусто, либо плотно в евклидовой топологии. Но Vс = выбирая х так, что
Ф^)
получаем x £ Vе (см. теорему 4.1 в [5]).
Для доказательства аналогичного результата для М = 3 в [3] предложен проверочный тест.
Тест Кахилла (СаЫ11) проверки инъективности отображения А при М = 3 :
1) выбрать измерительные векторы {фп}^,=1 С С3;
2) определить оператор матричного анализа А : Н3х3 ^ RN так, чтобы АН = = {(Н,фпф*п)нв }?=1;
3) если dimker(A) =0, то А инъективно, а, следовательно, и А инъективно;
4) если dimker(A) > 0, то выбрать Н £ ker(A), Н = 0;
4 а) если dimker(A) = 1, и det(H) =0, то А инъективно; 4 б) во всех остальных случаях А неинъективно.
Теорема 6. Для M = 3 тест Кахилла корректно определяет инъективность
A.
Доказательство. Если A инъективно, то A(x) = Axx* = Ayy* = A(y), то есть, y = x mod T.
Далее, предположим, что A имеет одномерное ядро. По лемме 4, A инъективно тогда и только тогда, когда ядро A порождено матрицей полного ранга 3.
Наконец, если размерность ядра не меньше 2, то существуют линейно независимые (ненулевые) матрицы A и B в ядре. Если det(A) = 0, то её ранг должен быть 1 или 2, и, по лемме 4, A не инъективно.
Если det(A) = 0, рассмотрим функцию
f : t ^ det(Acost + B sint) t e [0; п].
Так как f (0) = det(A) и f (п) = det(—A) = — det(A), по теореме о промежуточных значениях существует to e [0, п] такое, что f (to) = 0, то есть матрица A cos to + +B sin to вырождена. Эта матрица ненулевая, так как A и B линейно независимы, поэтому она имеет ранг 1 или 2. По лемме 9, A неинъективно. •
Теорема 7. (4M — 4)-гипотеза верна для M = 3.
Доказательство. (a) Если N < 4M — 4 = 8, то, по теореме о ранге и дефекте (см. доказательство теоремы 5), размерность ядра оператора матричного анализа A : H3x3 ^ RN не меньше 2, поэтому (тест Кахилла), A не является инъектив-ным.
(б) Рассмотрим вещественную 3 х 8-матрицу, подобную (5). Отображение A инъективно, если x e R48 порождает семейство {фп(х)}п=1 С C3, которое удовлетворяет тесту Кахилла. Для этого, согласно теореме о ранге и дефекте, ядро оператора матричного анализа должно иметь размерность не больше 1 и порождаться несингулярной матрицей. Для представления оператора A(x) матричного анализа в виде 8 х 9-матрицы используется ортонормированный базис пространства H3x3, подобный (6). Элементы этой матрицы, которую будем обозначать A(x), являются полиномами от x. Рассмотрим матрицу
B(x, y) =
T
y1
A(x)
где у € К9, и пусть и(х) обозначает вектор из миноров первой строки матрицы В(х,у). Тогда (у,и(х)) = 6.е1(В(х,у)). Отсюда следует, что и(х) лежит в ядре А(х), так как каждая строка А(х) ортогональна и(х).
Вспомогательное утверждение. и(х) = 0 тогда и только тогда, когда размерность ядра оператора А(х) не меньше 2, то есть строки А(х) линейно зависимы.
(^) Координаты вектора и(х) являются знакопеременными определителями 8 х 8 подматриц матрицы А(х), которые с необходимостью будут нулевыми в силу линейной зависимости строк.
(^) Имеем 0 = {у, 0) = (у,и(х)) = det(B(х,у)) для всех у € К9. Это означает, что даже если у = 0 и ортогонален строкам А(х), строки В(х,у) линейно зависимы а, следовательно, и строки А(х) должны быть линейно зависимыми. о
Используем доказанное вспомогательное утверждение для завершения доказательства теоремы. Определим матрицу и(х) € Н3х3 с координатами вектора и(х) в выбранном ранее базисе. Элементы матрицы и(х) являются полиномами от х. Отображение А инъективно тогда и только тогда, когда det и(х) = 0. Для обоснования рассмотрим три случая:
1) и(х) = 0, т. е. и(х) = 0, что эквивалентно dimker(A(x)) ^ 2^ По тесту Кахилла, А не инъективно;
2) ker(U(х)) =0, но det и(х) = 0^ По тесту Кахилла, А не инъективно;
3) ker(U(х)) =0, и det и(х) = 0^ По тесту Кахилла, А инъективно. Определим вещественное алгебраическое многообразие V = {х : det и(х) =
= 0} С R48 • Заметим, что А инъективно, только если х £ Vс• Множество Vс открыто в топологии Зарисского, поэтому оно либо пусто, либо всюду плотно в евклидовой топологии. Первая возможность исключается, так как, например, следующая матрица проходит тест Кахилла:
Ф
2 1
10 0 0
1
— 10 0 1 1 —1 —2 2 0 1—11—1 2i i -1
При таком выборе координат ядро оператора матричного анализа А одномерно и порождается невырожденной матрицей. Следовательно, соответствующее отображение A инъективно. • Теорема 7 опровергает высказанную ранее т. н. гипотезу Райта (Wright) [10]
о том, что
N*(M) < 3M — 2.
Авторы благодарны И.Я. Новикову за полезные обсуждения.
Литература
[1] Balan R., Casazza P., Edidin D. On signal reconstruction without noizy phase Applied and Computational Harmonic Analysis (ACHA). 2006. V. 20. I. 3. PP. 345-356.
[2] Bodmann B. G., Hammen N. Stable phase retrieval with low-redundancy frames. Available online: arXiv:1302.5487.
[3] Bandeira A., Cahill J., Mixon D., Nelson A. Saving phase: Injectivity and stability for phase retrieval. Applied and Computational Harmonic Analysis (ACHA). 2014. V. 37. I. 1. PP. 106-125.
[4] Fickus M., Mixon D.G., Nelson A.A., Wang Ya. Phase retrieval from very few measurements. Linear Algebra and Its Applications. 2014. V. 449. PP. 475-499.
[5] Новиков И.Я., Протасов В.Ю., Скопина М.А. Теория всплесков. М.: Физматлит. 2005.
[6] Новиков С.Я., Лихобабенко М.А. Фреймы конечномерных пространств. Самара. Самарский госуниверситет. 2013.
[7] Alexeev B., Cahill J., Mixon D. Full spark frames. Journal of Fourier Analysis and Applications. 2012. V. 18. I. 6. PP. 1167-1194.
[8] Heinosaary T., Mazzarella I., Wolf M.M. Quantum tomography ... Available online: arXiv:1206.1405.
[9] Хорн Р., Джонсон Ч. Матричный анализ. Москва. Мир. 1989.
[10] Vogt A. Position and momentum distributions do not determine the quantum mechanical state. In: Mathematical Foundations of Quantum Theory (Marlow A. R.,ed.). Academic Press. New York. 1978. PP. 365-372.
References
[1] Balan R., Casazza P., Edidin D. On signal reconstruction without noizy phase Applied and Computational Harmonic Analysis (ACHA). 2006. V. 20. I. 3. PP. 345-356.
[2] Bodmann B. G., Hammen N. Stable phase retrieval with low-redundancy frames. Available online: arXiv:1302.5487.
[3] Bandeira A., Cahill J., Mixon D., Nelson A. Saving phase: Injectivity and stability for phase retrieval. Applied and Computational Harmonic Analysis (ACHA). 2014. V. 37. I. 1. PP. 106-125.
[4] Fickus M., Mixon D. G., Nelson A. A., Wang Ya. Phase retrieval from very few measurements. Linear Algebra and Its Applications. 2014. V. 449. PP. 475-499.
[5] Novikov I.Ya., Protasov V. Yu., Skopina M. A. Wavelet theory Moscow, Physmatlit. 2005.
[6] Novikov S. Ya., Likhobabenko M. A. Frames in finitedimensional spaces. Samara. Samara University. 2013.
[7] Alexeev B., Cahill J., Mixon D. Full spark frames. Journal of Fourier Analysis and Applications. 2012. V. 18. I. 6. PP. 1167-1194.
[8] Heinosaary T., Mazzarella I., Wolf M. M. Quantum tomography ... Available online: arXiv:1206.1405.
[9] Horn R., Johnson Ch. Matrix Analysis. Moscow, Mir. 1989.
[10] Vogt A. Position and momentum distributions do not determine the quantum mechanical state. In: Mathematical Foundations of Quantum Theory (Marlow A. R.,ed.). Academic Press. New York. 1978. PP. 365-372.
S.Y. Novikov, M.E. Fedina3
RESTORING THE SIGNAL BY MODULES OF MEASUREMENT
The questions of the selection of vectors in finite-dimensional real and complex spaces so that the modules of scalar products of these vectors with the unknown vector were possible to restore it, are observed.
Also, the injectivity of nonlinear mappings is studied.
Key words: space, signal, systems of vectorsB, injectivity of mappings, scalar products.
Статья поступила в редакцию 15/VT/2016. The article received 15/VT/2016.
3Novikov Sergey Yakovlevich (nvks@ssau.ru), the Dept. of Theory of Probability and Mathematical Statistics, Samara University, Samara, 443011, Russian Federation.
Fedina Maria Efimovna (phedina@ssau.ru), the Dept. of Security of Information Systems, Samara University, Samara, 443011, Russian Federation.