Аннотация:
Алгебра эйконалов E метрического графа Ω – это операторная C∗-алгебра, определяемая динамической системой, которая описывает распространение волн, порожденных источниками в граничных вершинах Ω. В работе описывается каноническая блочная форма алгебры E произвольного компактного связного метрического графа. Переход к этой форме равносилен построению функциональной модели, реализующей E в виде алгебры непрерывных матричнозначных функций на ее спектре ˆE. Результаты предполагается использовать в обратной задаче, состоящей в реконструкции графа по спектральным и динамическим граничным данным.
Библиография: 28 наименований.
Ключевые слова:
динамическая система на метрическом графе, достижимые множества, C∗-алгебра эйконалов, каноническая форма.
Работа выполнена при поддержке РФФИ (грант № 20-01-00627A), Фонда Volkswagen Stiftung, Математического института им. Леонарда Эйлера (соглашение № 075-15-2019-1620) и при частичной поддержке Конкурса “Молодая математика России”.
Поступило в редакцию: 22.04.2021 Исправленный вариант: 09.10.2021
Существует подход к обратным задачам математической физики – метод граничного управления (BC-метод) [1]. Подход имеет выражено междисциплинарный характер: он основан на связях обратных задач с теорией систем и теорией управления, использует асимптотические методы, функциональный анализ, теорию операторов и др. Разработана алгебраическая версия BC-метода, основанная на связях с банаховыми алгебрами, давшая новое решение задачи реконструкции риманова многообразия по граничным данным [1]–[3]. Наша перспективная цель – применение этой версии к обратным задачам на графах. Данная работа – шаг в этом направлении.
Алгебраическая версия основана на фундаментальном факте: топологическое пространство характеризуется адекватной алгеброй. Как пример, компактное хаусдорфово пространство Ω с точностью до гомеоморфизма определяется алгеброй непрерывных функций A=C(Ω) (И. М. Гельфанд, 1943). При этом спектр алгебры – множество ˆA ее неприводимых представлений, снабженное адекватной топологией, – гомеоморфен пространству: ˆA≅Ω. Как следствие, располагая любым представлением A′ алгебры A и находя его спектр ^A′≅ˆA≅Ω, мы получаем гомеоморфную копию пространства Ω. По этой схеме решается задача реконструкции: из данных обратной задачи извлекается представление A′ и находится его спектр ^A′, который и доставляет решение задачи – гомеоморфную копию подлежащего восстановлению многообразия Ω. Содержательная часть подхода заключается в нахождении алгебры A′ по известным данным. В качестве последней используется алгебра эйконаловE, определяемая динамической системой, которая описывает распространение волн в Ω.
Варианты BC-метода для обратных задач на графах предложены в [4]–[6]. Версия, использующая алгебру эйконалов, инициирована в [7] и дополнена в [8]. Наша работа развивает этот подход. Его общее направление – изучение связей между свойствами алгебры E (блочной структурой, алгебраическими инвариантами, представлениями) и геометрией графа. Перспективная цель – реконструкция графа по его граничным данным.
1.2. Алгебра эйконалов
Для наглядности и без потери общности граф Ω можно представлять как связный компактный граф в R3, состоящий из гладких кривых (ребер) {e1,…,el}=E, скрепленных во внутренних вершинах1[x]1Любой метрический граф допускает такую реализацию.{v1,…,vm}=V. Имеются граничные вершины {γ1,…,γn}=Γ, из которых выходит по одному ребру. Метрика (внутреннее расстояние) в Ω индуцирована евклидовой метрикой из R3.
Ребра графа “материальны”: вдоль них распространяются колебания (волны), инициированные точечными источниками (управлениями), которые помещены в граничных вершинах. Волны движутся от границы с единичной скоростью, постепенно заполняя граф. Процесс описывается динамической системой
utt−Δu=0в H,0<t<T,u|t=0=ut|t=0=0в Ω,u=fв Γ×[0,T],
где H=L2(Ω), Δ – лапласиан, определенный на гладких функциях, удовлетворяющих условиям сшивания (Кирхгофа) во внутренних вершинах; f=f(γ,t) – граничное управление класса L2(Γ×[0,T])=:FT; u=uf(x,t) – решение (волна), uf(⋅,t)∈H при 0⩽t⩽T.
Имеется возможность управлять волнами не со всей границы, а с ее части Σ⊂Γ: в этом случае используются управления класса
FTΣ:={f∈FT∣suppf⊂Σ×[0,T]}=⨁∑γ∈ΣFTγ.
Каждой граничной вершине сопоставлено семейство достижимых множествUtγ:={uf(⋅,t)∣f∈FTγ}, 0⩽t⩽T, и соответствующих проекторов Ptγ в H на Ut. Оператор ETγ:=∫T0tdPtγ называется эйконалом, отвечающим вершине γ. Эйконалы суть самосопряженные операторы – элементы алгебры ограниченных операторов B(H).
Для C∗-алгебры A и множества S⊂A через ∨S обозначим C∗-алгебру, порожденную этим множеством, т. е. минимальную C∗-подалгебру в A, содержащую S. Алгебра эйконалов, отвечающая выделенному семейству граничных вершин Σ⊂Γ, есть операторная C∗-алгебра
ETΣ:=∨{ETγ∣γ∈Σ}⊂B(H).
1.3. Результаты и комментарии
Как установлено в [7], алгебра ETΣ имеет блочную структуру: она изометрически изоморфна некоторой подалгебре алгебры ⨁Jj=1C([0,ϵj];Mmj) и отличается от последней наличием связей между блоками. Для простейших графов (трехлучевых звезд) характер этих связей и их эволюция с изменением T рассмотрены в [8].
Главный результат данной работы – каноническое блочное представление алгебры эйконалов. Это представление, во-первых, выделено отсутствием связей между блоками и, во-вторых, инвариантно: с точностью до тривиальных преобразований (перестановки блоков, замены параметризации и т. п.) оно определяется по любому представлению алгебры (1.1). Последнее и подкрепляет надежду на полезность ETΣ в обратных задачах, данные которых определяют ее с точностью до изометрии.
Разумеется, об эффективности подхода в полной мере можно будет судить по конкретным приложениям в обратных задачах. В данной работе их нет, но приведенные выше результаты представляются важным шагом в этом направлении.
Алгебра эйконалов относится к классу C∗-алгебр с конечномерными представлениями разных размерностей [9], [10]. Приведение к канонической форме равносильно построению функциональной модели, которая реализует ETΣ в виде алгебры непрерывных матрично-значных функций на ее спектре. Она относится к типу моделей, описанных в работе Н. Б. Васильева [9].
Определение (1.1) фактически воспроизводит определение соответствующих алгебр, используемых в [1]–[3] для реконструкции многообразий. Собственно, успех такого применения и мотивировал попытку перенести подход на задачи на графах. Препятствием оказалась некоммутативность алгебры ETΣ. Эта проблема встречалась и раньше – в обратной задаче электродинамики [3], но там она снималась факторизацией по идеалу компактных операторов, с помощью которой дело сводилось к коммутативной C(Ω;R). Некоммутативность ETΣ неустранима, что существенно осложняет ее изучение.
Обратные задачи на графах – вполне актуальная тема. Различные постановки и подходы содержатся в работах С. А. Авдонина, П. Б. Курасова, М. Новачик, А. С. и В. С. Михайловых, П. А. Кучмента, В. А. Юрко. В [11]–[17] BC-метод, а также другие подходы используются для решения динамических и спектральных обратных задач для различных классов графов. Работы В. А. Юрко и его учеников используют спектральный подход к обратным задачам для дифференциальных операторов на графах [18]–[20]. Упомянем содержательный обзор П. А. Кучмента и Г. М. Берколайко по всей тематике квантовых графов, в том числе, и обратным задачам на них [21].
Весьма объемная вводная часть работы фактически повторяет соответствующие разделы из [7] и [8]. Это неизбежно, поскольку изложение относящихся к ETΣ фактов и результатов требует основательной подготовки. Техническая часть довольно сложна, поскольку работа имеет дело с максимально общим объектом – произвольным компактным связным метрическим графом, имеющим границу. Для понимания может быть полезен разбор простых примеров, проведенный в [7], [8] с иллюстрациями. Также рекомендуем работу [22], в которой на примере простейшего графа детально описана процедура приведения ETΣ к канонической форме.
Данная работа рассчитана на специалистов по C∗-алгебрам со вкусом к приложениям и/или специалистов в области математической физики, разделяющих мнение о полезности абстракций. Алгебра эйконалов – сложный и богатый свойствами объект, достойный всестороннего исследования. Для обратных задач отдельный интерес представляют возможные связи между ее алгебраическими инвариантами и геометрией графа. Как пример, упомянем гипотезу о соответствии кластеров в спектре ETΣ внутренним вершинам Ω. Очень интересен вопрос о том, как влияет на структуру спектра наличие циклов в графе.
§ 2. Волны на графе
2.1. Граф
Пусть Ω=E∪W есть связный компактный граф в R3 с ребрами {e1,…,eL}=E и вершинами {w1,…,wM}=W. Ребра суть гладкие2[x]2Всюду в работе гладкий означает C∞-гладкий. кривые, концами которых служат вершины. Удобно считать ребра открытыми, не относя к ним их концы. Вершина w и ребро eинцидентны (мы пишем w≺e), если w является концом e. Вершины {γ1,…,γN}=Γ, каждая из которых инцидентна единственному ребру, называются граничными; вершины {v1,…,vM−N}=V=W∖Γ суть внутренние.
Число μ(w) ребер, инцидентных вершине w называется ее валентностью; для γ∈Γ имеем μ(γ)=1. Дополнительно предполагаем, что вершин с μ(w)=2 нет, так что во внутренних вершинах выполнено μ(v)⩾3.
Граф снабжен метрикой (внутренним расстоянием) τ, индуцированной евклидовой метрикой из R3. При этом τ(a,b) есть минимум длин кусочно гладких кривых, лежащих в Ω и соединяющих точки a и b. Для множества A⊂Ω через
ΩrA:={x∈Ω∣τ(x,A)<r},r>0,
обозначается его метрическая окрестность радиуса r.
Каждое ребро e параметризовано длиной τ, отсчитываемой от одного из его концов. Для функции y на графе знак производной по длине dy/dτ зависит от выбора конца, но вторая производная d2y/dτ2 от этого выбора не зависит. Для вершины w и инцидентного ей ребра e определяется производная
[dydτ]+e(w):=lime∋x→wdydτ(x)
по направлению, исходящему из w. Величина
Fw[y]:=∑e≻w[dydτ]+e(w)
называется потоком функции y через вершину w.
Метрика τ на графе определяет (вещественное) гильбертово пространство H=L2(Ω) со скалярным произведением:
(y,u)H:=∫Ωyudτ=∑e∈E∫eyudτ.
Пусть C(Ω) есть пространство непрерывных функций с нормой ‖y‖=supΩ|y(⋅)|. Функцию y отнесем к классу Соболева H2(Ω), если y∈C(Ω) и на каждом ребре выполнено dy/dτ,d2y/dτ2∈L2(e).
Определим класс Кирхгофа
K:={y∈H2(Ω)∣Fv[y]=0,v∈V}.
Оператор Лапласа на графе вводится определением
Δ:H→H,DomΔ=K,(Δy)|e=d2ydτ2,e∈E.
Он плотно задан и замкнут.
2.2. Волны
Начально-краевая задача, описывающая распространение волн в графе, имеет вид
utt−Δu=0в H,0<t<T,
u(⋅,t)∈Kпри 0⩽t⩽T,
u|t=0=ut|t=0=0в Ω,
u=fна Γ×[0,T].
Здесь T>0 – финальный момент времени, f=f(γ,t) – граничное управление, u=uf(x,t) – решение (волна). При гладком (по t) управлении f, исчезающем вблизи t=0, задача имеет единственное классическое решение uf.
Из определения (2.1) на каждом ребре e решение uf удовлетворяет уравнению однородной струны utt−uττ=0. Отсюда видно, что волны распространяются от границы внутрь Ω с единичной скоростью. Как следствие, если управление действует с части границы Σ⊆Γ, т. е. выполнено suppf⊂Σ×[0,T], имеем соотношение
suppuf(⋅,t)⊂¯ΩtΣ,t>0.
Ниже вводятся (обобщенные) решения задачи (2.2)–(2.5) для управлений класса L2(Γ×[0,T]). Их определение требует некоторой подготовки.
Через δx обозначим меру Дирака, т. е. функционал на C(Ω), принимающий значения по правилу ⟨δx,y⟩=y(x). Пусть δ=δ(t) есть дельта-функция Дирака. Ближайшая цель – определить и описать фундаментальное решение задачи (2.2)–(2.5) с T=∞, отвечающее управлению f=δγδ(t), которое действует мгновенно из граничной вершины γ.
Для описания решения uδγδ удобно использовать следующий формализм “динамики импульсов”.
0. Под импульсом понимается мера aδx; постоянная a≠0 называется его амплитудой.
1. Каждый импульс aδx(t) движется по ребру со скоростью 1 в одном из двух возможных направлений, так что |˙x(t)|=1 при x(t)∈e.
2 (принцип суперпозиции). Импульсы движутся независимо друг от друга. Если в момент t имеется несколько импульсов a1δx(t),…,apδx(t), расположенных в точке x(t)∈Ω∖Γ, то они складываются, образуя импульс [a1+⋯+ap]δx(t).
3 (прохождение через внутреннюю вершину). Перемещаясь вдоль ребра e и проходя через внутреннюю вершину v, импульс aδx(t) делится на μ(v) импульсов: один отраженный и μ(v)−1 прошедших. Отраженный импульс движется вдоль e в обратном направлении и имеет амплитуду (2−μ(v)/μ(v))a. Каждый из прошедших импульсов движется вдоль своего (инцидентного v) ребра в сторону от v и имеет амплитуду (2/μ(v))a.
Таким образом, общая амплитуда составляет
2−μ(v)μ(v)a+[μ(v)−1]2μ(v)a=a,
что соответствует закону сохранения токов Кирхгофа Fv[y]=0.
4 (отражение от границы). Как только импульс aδx(t) достигает вершины γ∈Γ, он мгновенно меняет свое направление на противоположное и инвертируется: изменяет амплитуду с a на −a.
Приняв эти правила, можно описать решение uδγδ следующим образом (напомним, что τ – это расстояние в Ω):
(∗) при 0⩽t⩽τ(γ,V) имеем uδγδ=δx(t), где x(t) – точка ребра e≻γ такая, что τ(x(t),γ)=t. Таким образом, при малых временах uδγδ – это уединенный импульс единичной амплитуды, вошедший в граф из вершины γ и движущийся по e с единичной скоростью;
(∗∗) дальнейшая эволюция при временах t>τ(γ,V) определяется правилами 1–4.
Нетрудно убедиться, что такое описание является вполне детерминированным. В каждый момент времени t⩾0 решение uδγδ представляет собой набор конечного числа импульсов, движущихся в Ω. С физической точки зрения эта картина описывает, например, распространение острых сигналов (скачков напряжения) в электросети – графе, составленном из проводов.
Таким образом, фундаментальное решение есть некоторое пространственно-временное распределение в Ω×{t⩾0}. Структура его такова, что для управлений вида f=δγφ c φ∈L2[0,T] свертка по времени
uf(x,t):=[uδγδ(x,⋅)∗φ](t),x∈Ω,0⩽t⩽T,
корректно определена. Кроме того, можно показать, что uf∈C([0,T];H) и, если φ является гладким и исчезает около t=0, то uf совпадает с классическим решением задачи (2.2)–(2.5).
С этого момента функция uf, определяемая соотношением (2.7), рассматривается как (обобщенное) решение для управлений указанного вида. В более общем случае для f∈L2(Γ×[0,T]) вида f=∑γ∈Σfγ с fγ=δγφγ полагаем
uf(x,t):=∑γ∈Σufγ(x,t),x∈Ω,0⩽t⩽T.
Нетрудно показать, что и для обобщенного решения соотношение (2.6) остается верным. Оно показывает, что метрическая окрестность ΩTΣ есть часть графа, захваченная волнами, идущими от Σ, к моменту t=T.
2.3. Гидра
Здесь мы вводим пространственно-временной граф, который используется для эффективного описания волн.
Фиксируем граничную вершину γ. Рассматривая фундаментальное решение как пространственно-временное распределение, определим множество
Hγ:=suppuδγδ⊂ΩׯR+,
которое будем называть гидрой [6]. По существу это пространственно-временной граф, образованный траекториями импульсов в ходе эволюции, описываемой правилами 1–4 и (∗), (∗∗) (рис. 1)3[x]3Иллюстрации заимствованы из работы [7].. Его ребра суть характеристики волнового уравнения (2.2).
На гидре зададим функцию (амплитуду) a(⋅) по следующему правилу:
– для точки h∈Hγ такой, что π(h)=x∈Ω∖Γ и ρ(h)=t>0, имеем uδγδ(⋅,t)=aδx(⋅) и определяем a(h):=a;
– для h∈Hγ такой, что π(h)∈Γ и ρ(h)>0, полагаем a(h):=0;
– для h∈Hγ такой, что π(h)=γ и ρ(h)=0, полагаем a(h):=1.
Как видно, амплитуда является кусочно постоянной функцией, определенной на всей гидре Hγ (рис. 2). Уточним, что в точках самопересечения p согласно правилу 2 эволюции импульсов имеем a(p)=−4/9+1/3=−1/9. Точки вида p суть вершины Hγ, которые проектируются в Ω∖[V∪Γ].
Приведем представление, ради которого введена гидра (см. [7]). Используя запись h=(x,t)∈Hγ и a(h)=a(x,t), для управления f=δγφ(t) с φ∈L2[0,T] согласно (2.7) получим
uf(x,T)=∑t∈ρ(π−1(x))a(x,t)φ(T−t),x∈Ω.
В общем случае, когда управление f=∑γ∈Σδγφγ(t) действует из нескольких граничных вершин, согласно (2.8) имеем
uf(x,T)=∑γ∈Σ∑t∈ρ(π−1(x))aγ(x,t)φγ(T−t),x∈Ω,
где aγ суть амплитуды на гидрах HTγ.
Приведенные представления вполне эффективны: по ним можно вычислять значения волн. Однако для предстоящего анализа алгебры эйконалов потребуется их модификация, к описанию которой мы переходим. Модификация использует разбиение графа Ω на части (семейства), согласованное со структурой гидр. Во всех подробностях и с графическими иллюстрациями оно описано в [7].
2.4. Разбиение Π
В дальнейшем мы имеем дело с усеченными гидрами
HTγ:=Hγ∩{Ω×[0,T]}.
Введем одно понятие общего характера. Пусть дано множество X и рефлексивное симметричное (но, вообще говоря, не транзитивное!) бинарное отношение ∼0 на нем. Элементы, связанные этим отношением, будем называть соседними. Соседство определяет эквивалентность по следующему правилу. Мы говорим, что x и y эквивалентны (и пишем x∼y), если в X имеется конечный набор элементов x1,…,xn такой, что выполнено x∼0x1∼0⋯∼0xn∼0y.
Класс эквивалентности [x] элемента x∈X по введенному выше отношению ∼ можно описать конструктивно. Введем операцию ext, расширяющую подмножество B⊂X по правилу
B↦extB:=⋃b∈B{x∈X∣x∼0b};
обозначим ext1B:=extB и extjB:=extextj−1B, j⩾2. Нетрудно убедиться в справедливости представления
[x]=⋃j⩾1extj{x}.
Если множество X конечно, то последовательность extj на некотором шаге стабилизируется: ext1{x}⊂⋯⊂extN{x}=extN+1{x}=⋯=[x]. Этот случай встретится в работе.
Рассмотрим эквивалентность такого вида на гидре. Скажем, что точки h,h′∈HTγ суть соседние (h∼0h′), если выполнено хотя бы одно из условий: π(h)=π(h′) или ρ(h)=ρ(h′). Через γ∼ обозначим отношение эквивалентности, порожденное таким соседством. Класс эквивалентности
L[h]:={h′∈HTγ∣h′γ∼h}
назовем решеткой. Можно показать, что этот класс состоит из конечного числа точек. Для подмножества B⊂HTγ определим решетку
L[B]:=⋃h∈BL[h].
Заметим, что операция B↦L[B] обладает следующими свойствами:
Первые три свойства показывают, что она является топологическим замыканием (по Куратовскому).
С каждой точкой x∈¯ΩTγ свяжем множество
Λ[x]:=π(L[π−1(x)])⊂¯ΩTγ
(замыкание в метрике Ω), которое назовем множеством определенности точки x. Это множество конечно. Нетрудно проверить, что отношение
x∼x′⟺Λ[x]=Λ[x′]
есть эквивалентность, а операция A↦Λ[A]:=⋃x∈AΛ[x]=π(L[π−1(A)]) – топологическое замыкание. Множества A=Λ[A] мы называем Λ-замкнутыми.
На полной гидре точку h∈Hγ назовем угловой, если π(h)∈V∪Γ или h является точкой самопересечения (как p на рис. 2). Последние суть вершины гидры валентности 4, которые проектируются в Ω∖[V∪Γ].
На усеченной гидре HTγ, помимо содержащихся в ней угловых точек полной гидры, отнесем к угловым и точки множества ρ−1(T). Через CornHTγ обозначим множество всех угловых точек усеченной гидры.
Решетка L[CornHTγ] разделяет гидру на конечное число открытых пространственно-временных интервалов. На каждом интервале амплитуда a принимает постоянное значение.
Точки, составляющие конечное множество
Θ:=π(L[CornHTγ])⊂¯ΩTγ,
называются критическими. Оставшиеся точки x∈¯ΩTγ∖Θ назовем регулярными. Критические точки делят ¯ΩTγ на части. Множество регулярных точек
Π:=¯ΩTγ∖Θ
представляет собой совокупность конечного числа открытых интервалов, каждый из которых принадлежит некоторому ребру e. Таким образом, ¯ΩTγ=Π∪Θ есть разбиение части графа Ω, захваченной волнами, определяемое структурой гидры HTγ.
Пусть ω=(c,c′)⊂Π есть максимальный интервал, состоящий из регулярных точек. Максимальность означает, что концы интервала c и c′ суть критические точки, так что расширение ω сохранением регулярности внутренних точек невозможно. Несложно убедиться, что множество
Φ:=Λ(ω)=π(L[π−1(ω)])
состоит из максимальных интервалов ω1,…,ωm одинаковой длины:
Φ=m⋃k=1ωk,diamω1=⋯=diamωm=τ(c,c′)=:ϵΦ,
где τ – расстояние на графе. Будем говорить, что интервалы ωk являются клетками семейства Φ.
Сравнивая определения (2.13) и (2.14), приходим к представлению
Φ=⋃x∈ωΛ[x],
где ω – любая из клеток Φ.
Пусть ω′⊂Π есть максимальный интервал, не лежащий в семействе Φ. Он определяет еще одно семейство Φ′=Λ[ω′], состоящее из клеток и так далее. В результате множество Π оказывается конечным объединением непересекающихся семейств Φ1,…,ΦJ, каждое из которых состоит из непересекающихся клеток:
Π=J⋃j=1Φj=J⋃j=1mj⋃k=1ωjk,
где mj – количество клеток в Φj.
В пару к множеству определенности (2.13) каждому x∈¯ΩTγ∖Γ сопоставим множество
Ξ[x]:=ρ(L[π−1(x)])⊂[0,T].
При этом для x≠x′ выполнено либо Ξ[x]=Ξ[x′], либо Ξ[x]∩Ξ[x′]=∅. Также определим Ξ[B]:=⋃x∈BΞ[x].
Пусть Φ=⋃mΦk=1ωk⊂Π есть семейство. Легко проверить, что множество
Ψ:=Ξ[Φ]=nΦ⋃i=1ψi⊂[0,T]
состоит из временных интервалов ψi:=(ti,˜ti) таких, что 0⩽t1<˜t1⩽t2<˜t2⩽⋯⩽tnΦ<˜tnΦ⩽T и имеющих одинаковую длину ˜ti−ti=ϵΦ. Множество Ψ также будем называть семейством, состоящим из временных клеток ψi.
В дальнейшем используются функции τi:Φ→[0,T], связанные с разбиением графа на семейства. Они вводятся следующим образом4[x]4Приводимое определение отличается от принятого в [7]. Причину изменений мы укажем позже.. Пусть x∈Φ, положим
τi(x):=ψi∩ρ(L[π−1(x)]),i=1,…,nΦ.
Так как для любого xk∈Λ[x] выполнено L[π−1(x)]=L[π−1(xk)], то эти функции постоянны на множествах определенности: τi(x)=τi(xk). Из определения следует, что
τi(x)≠τi′(x)приi≠i′,x∈Φ.
С изменением x внутри клетки ω=(c,c′)⊂Φ изменяется множество Λ[x]. При этом значение τi(x) заметает клетку ψi=(ti,˜ti)⊂Ψ и, как нетрудно видеть из определения (2.18), справедливо одно из двух представлений
τi(x)=ti+τ(x,c)илиτi(x)=˜ti−τ(x,c).
Имея ввиду это представление, можно сказать, что функции τi зависят от x∈ω линейно.
Свойство (2.20) позволяет доопределить функции τi в критических точках: если x∈ω=(c,c′), x→c, то τi(c)=ti или τi(c)=˜ti в зависимости от того, какое из представлений (2.20) имеет место.
Для каждого семейства Φ⊂Π набор функций τi свой и при необходимости мы отмечаем это обозначением τiΦ. Добавим, что в силу (2.19) и (2.20) равенство τiΦ(x)=τi′Φ′(x) для разных Φ и Φ′ возможно лишь в критических точках x.
Разбиение графа на семейства мотивировано, в частности, тем, что волны uf зависят от управлений f локально в следующем смысле. Как видно из (2.9), значения uf(⋅,T)|Φ определяются значениями f|Ξ[Φ]. Более того, условия
suppf⊂Ξ[Φ]иsuppuf(⋅,T)⊂Φ
являются равносильными.
Модификация представлений (2.9) и (2.10), о которой говорилось в конце предыдущего пункта, состоит в эффективном представлении волн на множествах определенности. Оно использует функции (векторы) на Λ[x], к описанию которых мы переходим.
2.5. Амплитудные векторы
Текущие рассмотрения по-прежнему относятся к одной граничной вершине γ∈Γ, фиксированному T>0 и соответствующему разбиению Π.
Пусть A={x1,…,xm}⊂¯ΩTγ есть конечное Λ-замкнутое множество и
Ξ[A]=ρ(L[π−1(A)])={t1,…,tn},
где t1<⋯<tn. Множеству A сопоставим функции (амплитудные векторы) αi:A→R, i=1,…,n, по правилу
Набор {α1,…,αn} назовем α-набором над множеством A. Для каждого x∈¯ΩTγ множество определенности Λ[x] является Λ-замкнутым и над ним имеется свой α-набор.
есть время заполнения всего графа Ω волнами, идущими (с единичной скоростью) от вершины γ. Справедлива лемма 1.
Лемма 1. Пусть x есть точка в ¯ΩTγ, Λ[x] – ее множество определенности, и пусть α1,…,αn есть α-набор над Λ[x]. Если T<Tγ, то векторы αi линейно независимы.
Доказательство. Фиксируем какой-либо момент времени t∗=ti∈Ξ[x] с i≠1. В решетке L[π−1(x)] обязательно найдется точка (x∗,t∗) такая, что x∗ принадлежит границе ¯Ωt∗γ∖Ωt∗γ области Ωt∗[γ], захваченной волнами к моменту t=t∗. Эта граница непуста в силу t∗⩽T<Tγ; при этом x∗∈Λ[x], t∗=τ(x∗,γ) и, очевидно, a(x∗,t∗)=αi(x∗)≠0.
В то же время, во всех точках (x∗,t) с t∈Ξ[x∗]=Ξ[x], t<t∗, выполнено a(x∗,t)=0, поскольку такие (x∗,t) не лежат на гидре HTγ. Последнее соответствует простому факту: при указанных временах t волны от γ не успевают дойти до точки x∗ (см. (2.6)).
Таким образом, в точках (x∗,t) с t=tk∈Ξ[x], ti<t∗ (k=1,…,i−1) имеем a(x∗,tk)=αk(x∗)=0 и, в то же время, a(x∗,t∗)=αi(x∗)≠0. Это исключает линейную зависимость αi от α1,…,αi−1. В силу произвольности i получаем линейную независимость всего набора α1,…,αn. Лемма доказана.
Следствие 1. При T<Tγ справедливо равенство Ξ[¯ΩTγ]=[0,T].
В самом деле, из доказательства леммы 1 видно, что условие T<Tγ обеспечивает непустоту множества ρ−1(t)⊂HTγ при 0⩽t⩽T, что равносильно равенству Ξ[¯ΩTγ]=[0,T].
Вернемся к представлению (2.9). В терминах амплитудных векторов оно может быть записано в форме
представляющей волну не только в точке x, но и на всем множестве определенности Λ[x]. Более того, согласно (2.15), варьируя точку x внутри клетки ω⊂Φ, мы представляем волну uf(⋅,T) на всем семействе Φ.
Заключительный шаг к модификации исходного представления (2.9) состоит в переходе в (2.22) к более удобной системе амплитудных векторов.
2.6. β-представление волн
Пусть опять A=Λ[A]={x1,…,xm}⊂¯ΩTγ. Введем пространство l2(A) функций (векторов) на A со скалярным произведением
⟨f,g⟩=∑x∈Af(x)g(x)=m∑k=1f(xk)g(xk).
Оно содержит подпространство
A[A]:=span{α1,…,αn},dimA[A]⩽n,
определяемое α-набором над A. Используя процедуру Грама–Шмидта, перейдем в A[A] к набору
βi:={α1‖α1‖,если i=1,αi−∑i−1j=1⟨αi,βj⟩βj‖αi−∑i−1j=1⟨αi,βj⟩βj‖,если i⩾2 и αi∉span{α1,…,αi−1},0,если αi∈span{α1,…,αi−1},
где βi – вектор с компонентами βi(x1),…,βi(xm). Для ненулевых элементов набора выполнено ⟨βi,βj⟩=δij, а их линейная оболочка, очевидно, совпадает с A[A]. Совокупность векторов {β1,…,βn} назовем β-набором над множеством A.
Для каждого x∈¯ΩTγ над множеством определенности Λ[x] имеется свой β-набор. Согласно лемме 1 при T<Tγ векторы αi линейно независимы. Как следствие, при таких T все векторы βi суть ненулевые и выполнено dimA[Λ[x]]=n⩽m.
С использованием β-набора над множеством A=Λ[x] представление (2.22) приобретает окончательный вид
где ρ – матрица перехода, связывающая наборы α и β.
2.7. Гидра HTΣ
Введенные выше понятия и объекты относились к одной граничной вершине γ. В дальнейшем, указывая на это при необходимости, мы используем обозначения γ∼, Lγ, Λγ и т. п.
Для набора граничных вершин Σ⊆Γ определим пространственно-временной граф
HTΣ:=⋃γ∈ΣHTγ⊂Ω×[0,T].
На нем задаются аналоги объектов, введенных ранее для отдельных гидр HTγ. Опишем их.
Проекции из HTΣ в ¯ΩTΣ и в [0,T] суть π((x,t)):=x и ρ((x,t)):=t; под π−1 и ρ−1 понимаются полные прообразы в HTΣ.
Соседство h∼0h′ на HTΣ по определению означает, что π(x)=π(x′) и/или ρ(x)=ρ(x′). Соседство определяет эквивалентность hΣ∼h′. Через LΣ[h] обозначается класс эквивалентности (решетка) точки h∈HTΣ. Операция HTΣ⊃B↦LΣ[B] есть (топологическое) замыкание.
Множество ΛΣ[x]:=π(LΣ[π−1(x)]) назовем множеством определенности точки x∈¯ΩTΣ. Отметим очевидное вложение Λγ[x]⊂ΛΣ[x] для γ∈Σ. Операция ¯ΩTΣ⊃A↦ΛΣ[A] есть (топологическое) замыкание.
Множество угловых точек CornHTΣ состоит их всех угловых точек гидр HTγ⊂HTΣ плюс точки (трансверсального) пересечения ребер разных гидр HTγ. Критические точки в ¯ΩTΣ суть ΘΣ:=π(LΣ[CornHTΣ]), регулярные суть ΠΣ:=¯ΩTΣ∖ΘΣ. Имеет место разбиение
ΠΣ=J⋃j=1Φj=J⋃j=1mj⋃k=1ωjk,diamωjk=ϵj:=ϵΦ,
на семейства и клетки, вполне аналогичное разбиению (2.16). Для каждого семейства Φ⊂ΠΣ множество
ΞΣ[Φ]:=ρ(LΣ[π−1(Φ)])=nΦ⋃i=1ψi⊂[0,T]
состоит из временных интервалов ψi=(ti,˜ti) таких, что 0⩽t1<˜t1⩽t2<˜t2⩽⋯⩽tnΦ<˜tnΦ⩽T, все интервалы имеют одинаковую длину ˜ti−ti=ϵΦ. На каждом семействе определен набор функций τi:
τiΦ(x):=ψi∩ρ(LΣ[π−1(x)]),x∈Φ,i=1,…,nΦ.
Эти функции постоянны на множествах определенности: τiΦ(x)=τiΦ(x′) для регулярных x,x′∈ΛΣ[x]. Для них выполнено (2.19) и справедливо представление (2.20). Последнее позволяет доопределить функции τiΦ в критических точках (концах клеток ω) по непрерывности.
Пусть γ∈Σ, Φ⊂ΠΣ и x∈Φ. Нетрудно показать, что множество ΛΣ[x]∩ΩTγ является Λγ-замкнутым (в ΩTγ). Следовательно, на нем имеется β-набор векторов β1γΦ,…,βnΦγΦ. Условимся считать их заданными на всем ΛΣ[x], продолжив с ΛΣ[x]∩ΩTγ на ΛΣ[x] нулем.
Повторив построение для всех γ∈Σ, получим совокупность β-наборов
{β1γΦ,…,βnΦγΦ∣γ∈Σ},
где βiγΦ – вектор с компонентами (βiγΦ)1,…,(βiγΦ)mΦ и mΦ=♯ΛΣ[x]. Каждый из наборов ортонормирован в l2(ΛΣ[x]). Соответственно на каждом семействе Φ⊂ΠΣ заданы функции βiγΦ(⋅), принимающие постоянные значения (βiγΦ)k на клетках ωk⊂Φ.
Для данного семейства Φ и разных вершин γ∈Σ набор {β1γΦ,…,βnΦγΦ} содержит одно и то же число векторов, равное nΦ, а функции τiΦ суть одни и те же. Тем не менее, примем по определению
nγΦ:=nΦ,τiγΦ(x):=τiΦ(x),x∈Φ,γ∈Σ.
Это на первый взгляд избыточное обозначение (индекс γ) окажется удобным в дальнейших рассмотрениях.
§ 3. Эйконалы
3.1. Достижимые множества и проекторы
Здесь мы рассматриваем задачу (2.2)–(2.5) как динамическую систему и снабжаем ее атрибутами теории управления – пространствами и операторами.
Пространство управлений FT:=L2(Γ×[0,T]) со скалярным произведением
(f,g)FT=∑γ∈Γ∫T0f(γ,t)g(γ,t)dt
называется внешним пространством системы (2.2)–(2.5). Оно содержит подпространства управлений, действующих с отдельных граничных вершин γ∈Γ:
FTγ:={f∈FT|suppf⊂{γ}×[0,T]}.
Каждое управление f∈FTγ имеет вид f=δγφ с некоторым φ∈L2[0,T]. Подмножеству граничных вершин Σ⊆Γ соответствует подпространство
FTΣ:=∑γ∈Σ⊕FTγ
действующих с него управлений (слагаемые ортогональны в FT).
Пространство H=L2(Ω) называется внутренним; волны uf(⋅,t) суть его элементы, зависящие от времени. Для множества B⊂Ω определим подпространство H⟨B⟩:={y∈H∣suppy⊂¯B} функций, локализованных в B.
Множество волн
Usγ:={uf(⋅,s)∣f∈FTγ}⊂H,0⩽s⩽T,
называется достижимым (с вершины γ к моменту t=s). Из представлений (2.9) и (2.10) видно, что Usγ суть (замкнутые) подпространства в H. С ростом s они расширяются: Usγ⊂Us′γ при s<s′.
Отмеченная в (2.21) локальность соответствия “управление–волна” ведет к разложению по семействам
UTγ=∑Φ⊂ΠΣ⊕UTγ⟨Φ⟩,
где подпространство UTγ⟨Φ⟩⊂H⟨Φ⟩ состоит из волн uf(⋅,T)∈UTγ, локализованных в Φ∩ΩTγ. Ортогональность слагаемых – следствие дизъюнктности семейств: Φj∩Φk=∅ для j≠k.
Фиксируем граничную вершину γ∈Σ; пусть PTγ есть проектор в H на подпространство UTγ. Обсудим его свойства и опишем его действие.
где PTγ⟨Φ⟩ суть проекторы в H на UTγ⟨Φ⟩. Таким образом, проектор PTγ приводится подпространствами UTγ⟨Φ⟩ и его характеризация сводится к описанию действия проекторов PTγ⟨Φ⟩.
Как показано в [7], проекторы PTγ⟨Φ⟩ выражаются через векторы (2.25) и отвечающие им функции βiγΦ(⋅) следующим образом:
где y∈H – произвольная функция на графе5[x]5В точном смысле (3.3) есть представление “почти всюду” в Ω. Для y∈C(Ω) оно верно “всюду”.. Это представление выводится из представления (2.23), для получения которого и вводились векторы βi.
Как видно из (3.2) и (3.3), если проектируемая на UTγ функция y∈H непрерывна и такова, что y|Λγ[x]≡0, то и (PTγy)|Λγ[x]≡0. Другими словами, значения проекции PTγy на множестве Λγ[x] вполне определяются значениями y на Λγ[x]. Именно в этом состоит повод называть Λγ[x], а с ними и ΛΣ[x], множествами определенности.
Напомним, что векторы βiγΦ суть элементы подпространства A[ΛΣ[x]]⊂l2(ΛΣ[x]). Из сказанного выше следует, что проектор PTγ определяет в l2(ΛΣ[x]) оператор pγΦ[x], проектирующий на A[Λγ[x]]⊂A[ΛΣ[x]]. Пусть χ1,…,χmΦ есть (стандартный) базис в l2(ΛΣ[x]), состоящий из индикаторов точек множества ΛΣ[x]; в нем (βiγΦ)k=⟨βiγΦ,χk⟩. Согласно (3.3) в этом базисе матрица проектора pγΦ[x] принимает вид
Семейство {Psγ∣0⩽s⩽T} проекторов в H на достижимые множества Usγ определяет оператор эйконала (коротко – эйконал)
ETγ:H→H,ETγ:=∫T0sdPsγ.
Из определения следует, что ETγ есть ограниченный самосопряженный положительный оператор. Как и проектор PTγ, эйконал ETγ приводится подпространствами H⟨Φ⟩: для него справедливы соотношения ETγH⟨Φ⟩⊂H⟨Φ⟩ и разложение
ETγ=∑Φ⊂ΠETγ⟨Φ⟩,
где ETγ⟨Φ⟩:=ETγ|Φ есть сужение ETγ, действующее в H⟨Φ⟩.
Как показано в [7], справедливо согласованное с (3.3) представление
в котором y∈H произвольно, а функции τiγΦ даются определениями6[x]6В соответствующем представлении в [7] использовались некорректно определенные функции τiΦ (см. сноску 4). Представление (3.5) с функциями, определенными в (2.18), (2.26), исправляет эту ошибку.(2.24) и (2.26).
Оператор ETγ индуцирует в l2(ΛΣ[x]) оператор eγΦ[x]. Согласно (3.4) и (3.5) в базисе χ1,…,χmΦ его матрица имеет вид
Подпространство AγΦ[x]=span{β1γΦ,…,βnγΦγΦ}⊂l2(ΛΣ[x]) приводит эту матрицу, а ее ненулевой блок в базисе β1γΦ,…,βnγΦγΦ есть diag{τiγΦ(x)}nγΦi=1. Как видно из (2.19), все ее собственные значения τiγΦ(x) внутри клеток семейства различны; при изменении x в ω⊂Φ они заметают интервалы (ti,˜ti)⊂Ξ[Φ]. Согласно следующему далее предложению 1 при T<Tγ объединение (по всем семействам Φ⊂Π) всех сегментов [ti,˜ti] совпадает с [0,T].
Из приведенных рассмотрений легко усматриваются общие свойства эйконала как оператора в H (см., например, [23]).
Предложение 1. Пусть для эйконала ETγ выполнено ¯RanETγ=UTγ и KerETγ=H⊖UTγ. Эйконал приводится частями достижимого множества: ETγUTγ⟨Φ⟩⊂UTγ⟨Φ⟩, Φ⊂ΠΣ. При T<Tγ оператор ETγ|UTγ имеет простой абсолютно непрерывный спектр, заполняющий сегмент [0,T].
Замечание 1. Из представления (3.5) следует, что при временах T>Tγ спектр ETγ|UTγ есть объединение сегментов [0,T0]∪[T1,T2]∪⋯∪[TN−1,TN], где Tγ⩽T0<T1<⋯<TN⩽T, а каждый сегмент состоит из областей значений функций7[x]7Однако точное описание спектра – открытый вопрос. Есть гипотеза, что он всегда исчерпывается интервалом [0,T0] с достаточно большим T0. Вопрос связан с тонкими деталями строения достижимых множеств UTγ метрического графа.τiγΦ (замыканий временных клеток ψiγΦ, см. (2.17)).
3.3. Параметризация
Выберем семейство Φ=⋃mΦk=1ωk⊂ΠΣ; пусть ω=(c,c′)⊂Φ – одна из его клеток. Напомним, что все клетки имеют одинаковую длину ϵΦ=τ(c,c′). Для x∈ω введем обозначение x=x(r), если τ(c,x)=r. Наряду с x, множество определенности также оказывается параметризованным: ΛΣ[x(r)]={xk(r)}mΦk=1. При изменении r в интервале (0,ϵΦ) точки xk(r) непрерывно меняют положение и заметают клетки ωk. Таким образом, семейство Φ параметризовано.
Параметризованы по r и все элементы представлений (3.3) и (3.5): векторы
и функции τiγΦ(r):=τiγΦ(x(r)). Последние согласно (2.20) принимают значения
τiγΦ(r)=tiΦ+rилиτiγΦ(r)=˜tiΦ−r=(tiΦ+ϵΦ)−r.
Отметим, что параметризаций семейства Φ суть две: принятая выше с r=τ(x,c) и отвечающая параметру r=τ(x,c′). Они вполне равноправны. В дальнейшем будем считать, что каждое семейство Φ⊂ΠΣ параметризовано каким-либо одним из двух способов.
Параметризация определяет матричные представления функций и операторов на графе.
Пусть Φ⊂ΠΣ есть параметризованное семейство, y∈H – функция на графе, x=x(r)∈ΛΣ[x(r)]={xk(r)}mΦk=1⊂Φ, 0<r<ϵΦ. Легко проверить, что отображение
где τiγΦ(r) имеют вид (3.7). При изменении параметра r матрица BγΦ не меняется в силу постоянства ее элементов в клетках ω1,…,ωmΦ семейства Φ. Вместе с ней постоянна и матрица B∗γΦBγΦ. Последняя в силу того, что столбцы βiγΦ составляют ортонормированный набор, является проектором в RmΦ на подпространство
Aγ[Φ]:=span{β1γΦ,…,βnγΦγΦ}=[B∗γΦBγΦ]RmΦ.
Матрица-проектор B∗γΦBγΦ допускает разложение в сумму попарно ортогональных одномерных проекторов
B∗γΦBγΦ=nγΦ∑i=1PiγΦ,PiγΦ:=⟨⋅,βiγΦ⟩βiγΦ,
где ⟨⋅,⋅⟩ – стандартное скалярное произведение в RmΦ. В соответствии с (3.5) имеем представление
Опишем параметризацию пространств и операторов, отвечающую разбиению ΠΣ в целом. Аналог разложения (3.1) приобретает вид
UTΣ=⊕∑Φ⊂ΠΣUTΣ⟨Φ⟩,UTΣ⟨Φ⟩:=span{UTγ⟨Φ⟩∣γ∈Σ},
причем каждое из подпространств-слагаемых приводит все эйконалы одновременно:
ETγUTΣ⟨Φ⟩⊂UTΣ⟨Φ⟩,γ∈Σ.
Используя параметризации в семействах, имеем
UΦH⟨Φ⟩=L2([0,ϵΦ];RmΦ),UΦUTΣ[Φ]=L2([0,ϵΦ]);AΣ[Φ]),
где AΣ[Φ]:=span{Aγ[Φ]∣γ∈Σ}, а части ETγ⟨Φ⟩ эйконалов умножают элементы L2([0,ϵΦ];RmΦ) на матрицы-функции (3.9).
Соглашение 1. Условимся об обозначениях. Для пространств S1,…,Sn сумма S=⨁∑jSj есть пространство наборов s={s1,…,sn}, sj∈Sj (с покомпонентными операциями). Для операторов A1,…,An, Aj∈EndSj, сумма A=⨁∑jAj∈EndS есть оператор, действующий по правилу As:={A1s1,…,Ansn}. Для матриц M1,…,Mn, Mj∈Mϰj, сумма M=⨁∑jMj∈Mϰ1+⋯+ϰn есть блочно-диагональная матрица с блоками M1,…,Mn (мы также пишем [M]j=Mj). Для алгебр {A1,…,An}A=⨁jAj есть прямая сумма алгебр-слагаемых (мы также пишем [A]j=Aj и называем Aj блоками).
Параметризация всего ΠΣ реализуется оператором U:=⨁∑Φ⊂ΠΣUΦ:
где \dot D_{\gamma\Phi}(\,{\cdot}\,):=D_{\gamma\Phi}(\,{\cdot}\,)+ I, I – единичная матрица соответствующей размерности и \dot\tau_{\gamma\Phi}^i(r) :=\tau_{\gamma\Phi}^{i}(r)+1. Аналог предложения 1 имеет следующий вид.
Предложение 2. Для оператора \dot E^T_{\gamma} выполнено {\operatorname{Ran}\dot E^T_{\gamma}}=\mathscr{U}^T_\gamma, \operatorname{Ker}\dot E^T_{\gamma}=\mathscr{H}\ominus\mathscr{U}^T_\gamma и \dot E^T_{\gamma}\mathscr{U}^T_\gamma\langle\Phi\rangle\subset\mathscr{U}^T_\gamma \langle\Phi\rangle, \Phi\subset\Pi. При T<T_\gamma он имеет собственное значение 0 бесконечной кратности и простой абсолютно непрерывный спектр, заполняющий сегмент [1,T+1]. При T>T_\gamma согласно замечанию 1 выполнено
где T_\gamma\leqslant T_0<T_1<\dots <T_N\leqslant T, а сегменты суть объединения областей значений функций \dot\tau^i_{\gamma\Phi} (смещенных на 1 временных клеток \dot\psi^i_{\gamma\Phi}).
§ 4. Алгебра эйконалов
4.1. Определения и общие факты
Напомним, что C^*-алгеброй называется банахова алгебра с инволюцией [24], [25]. Таковыми, в частности, являются алгебры ограниченных операторов \mathfrak{B}(\mathscr{H}) в гильбертовом пространстве \mathscr{H}, в которых роль инволюции играет операторное сопряжение. Запись {\mathfrak A}\cong{\mathfrak B} будет означать, что C^*-алгебры {\mathfrak A} и {\mathfrak B} связаны изометрическим *-изоморфизмом (далее коротко – изометричны). Для множества S\subset \mathfrak{A} через \vee S обозначается минимальная C^*-(под)алгебра в \mathfrak{A}, содержащая S.
Под \mathbb{M}^n понимается алгебра вещественных (n\times n)-матриц, рассматриваемых как операторы в \mathbb R^n и снабженных соответствующей (операторной) нормой. Она неприводима.
Через C([a,b],\mathbb{M}^n) обозначается алгебра непрерывных \mathbb{M}^n-значных функций с нормой \|c\|=\sup_{a\leqslant t\leqslant b} \|c(t)\|_{\mathbb{M}^n}. Тем же символом мы обозначаем операторную (под)алгебру в \mathfrak{B}(L_2([a,b];\mathbb R^n)), элементы которой умножают квадратично-суммируемые \mathbb R^n-значные функции на функции из C([a,b],\mathbb{M}^n). Соответствие c\mapsto c\cdot устанавливает изометрию этих алгебр.
C^*-подалгебру \mathfrak A\subset \mathbb{M}^n условимся считать неприводимой, если выполнено \mathfrak A\cong \mathbb{M}^k, где k\leqslant n. Такая алгебра в подходящем базисе в \mathbb R^n принимает блочно-диагональную форму и состоит из двух блоков, один из которых есть \mathbb{M}^k, а второй (если имеется) – нулевой.
Приведем сводку известных результатов8[x]8В работе [8] в аналогичной сводке “On matrix algebras” допущена ошибка: утверждение 3 неверно. Однако после соответствующих исправлений все результаты работы сохранили силу..
Предложение 3. Любая C^*-подалгебра алгебры \mathbb{M}^n изометрична прямой сумме \bigoplus_{k}\mathbb{M}^{n_k}, где \sum_{k}{n_k}\leqslant n.
Предложение 4 (см. [24]). Пусть \mathfrak{P}\subset\mathbb{M}^n и C^*-подалгебра \mathfrak{A}\subset C([a,b];\mathfrak{P}) такова, что для любых t,t'\in[a,b] и p,p'\in\mathfrak{P} найдется элемент u\in\mathfrak{A}, для которого выполнено u(t)=p, u(t')=p'. Тогда \mathfrak{A}= C([a,b];\mathfrak{P}).
Представление C^*-алгебры \mathfrak{A} – это гомоморфизм \pi\colon \mathfrak{A}\to\mathfrak{B}(H), где H – гильбертово пространство. Эквивалентность представлений \pi\sim\pi' означает, что \iota\,\pi(a)=\pi'(a)\iota, a\in\mathfrak{A}, где \iota\colon H\to H' – изометрия пространств представлений. Представление неприводимо, если операторы \pi(\mathfrak{A}) не имеют общего ненулевого инвариантного подпространства в H.
Спектром C^*-алгебры \mathfrak{A} называется множество \widehat{\mathfrak{A}} классов эквивалентности ее неприводимых представлений. Класс эквивалентности (точку спектра), отвечающий представлению \pi, будем обозначать через \widehat\pi. Спектр снабжен канонической топологией Джекобсона [24], [25].
Изометрия алгебр \mathrm{u}\colon \mathfrak{A}\to\mathfrak{B} определяет соответствие представлений
неприводимы; их классы эквивалентности исчерпывают спектр алгебры C([a,b],\mathbb{M}^n). Для любого неприводимого представления \pi алгебры C([a,b],\mathbb{M}^n) существует единственная точка t\in[a,b] такая, что \pi\sim\pi_t.
В случае \mathbb M_a=\mathbb M_b=\mathbb M^n имеем \dot C([a,b]; \mathbb M^n)= C([a,b]; \mathbb M^n). Алгебры \dot C([a,b]; \mathbb M^n) мы будем называть стандартными. Спектр стандартной алгебры состоит из классов \widehat\pi_t, t\in(a,b), неприводимых представлений вида (4.3) и представлений \widehat\pi_a, \widehat\pi_b, которые могут оказаться приводимыми. Если, например, n_a\geqslant 2, то \pi_a распадается на неприводимые представления
где [\,{\cdots}\,]^k есть k-й блок блочно-диагональной матрицы в представлениях (4.4). В этом случае мы говорим, что \widehat\pi_a^1,\dots,\widehat\pi_a^{n_a} образуют кластер в спектре стандартной алгебры. Этот термин мотивирован тем, что они неотделимы друг от друга в топологии Джекобсона. Аналогичный кластер может иметься и на правом конце t=b. В то же время, все \widehat\pi_t с разными t\in(a,b) отделимы друг от друга и от кластеров (см. [7], [8]). Спектр алгебры C([a,b]; \mathbb M^n) кластеров не содержит.
Центральным объектом работы является алгебра эйконалов графа \Omega
Переход от \mathfrak{E}^T_{\gamma} к \dot{\mathfrak{E}}^T_{\gamma} состоит в добавлении проектора P^T_\gamma, играющего роль единицы в \dot{\mathfrak{E}}^T_{\gamma} (см. (3.11)), а алгебра \mathfrak{E}^T_{\gamma} оказывается подалгеброй в \dot{\mathfrak{E}}^T_{\gamma}.
Согласно функциональному исчислению самосопряженных операторов и в силу ортогональности проекторов P^i_{\gamma \Phi} в (3.12) имеем
для \varphi\in C(\sigma_{\mathrm{ac}}(\dot E^T_\gamma)). Задаваемое первым равенством соответствие \varphi(\dot E^T_\gamma)\leftrightarrow\varphi есть изометрия алгебр \dot{\mathfrak{E}}^T_\gamma и C(\sigma_{\mathrm{ac}}(\dot E^T_\gamma)).
Соглашение 2. Далее всюду, если не оговорено противное, мы имеем дело только со смещенными эйконалами и опускаем (\,\dot{}\,) в относящихся к ним обозначениях: \dot E^T_\gamma\equiv E^T_\gamma, \dot\tau^k_{\gamma l}\equiv\tau^k_{\gamma l}, \dot\psi^k_{\gamma l}\equiv\psi^k_{\gamma l}, \dot{\mathfrak{E}}^T_{\Sigma}\equiv\mathfrak{E}^T_{\Sigma} и так далее.
здесь нулевые внедиагональные блоки опущены, \Pi_\Sigma=\Phi^1\cup\dots \cup\Phi^J. Обозначение (\,\cdot_j) подчеркивает, что аргументы r_j\in[0,\epsilon_j] у функций \tau^i_{\gamma \Phi^j} суть разные в соответствии с представлением (4.7).
Определим наборы проекторов \mathbb P_{\Phi^j}:=\{P^i_{\gamma \Phi^j}\mid i=1,\dots,n_{\gamma \Phi^j};\, \gamma\in\Sigma\}. Тогда с привлечением алгебр
Как можно предвидеть из (4.8) и (4.10), описание структуры алгебры эйконалов сведется к установлению связи между ее блоками [U\mathfrak{E}^T_\Sigma U^{-1}]_j(\,\cdot_j), отвечающими разным семействам \Phi^j\subset\Pi_\Sigma. Именно эти связи отличают U\mathfrak{E}^T_\Sigma U^{-1} от алгебры в правой части (4.10), у которой блоки вполне независимы. Приводимая ниже лемма – шаг в изучении связей между блоками алгебры U\mathfrak{E}^T_\Sigma{ U^{-1}}.
где O_k – нулевые матрицы соответствующей размерности. Как и P^i_{\gamma \Phi^j}, эти проекторы попарно ортогональны. Определим также “точки” \mathbf{r}:=\{r_1,\dots ,r_J\} с координатами r_j\in[0,\epsilon_j], матрицы
Лемма 2. Фиксируем \gamma, i, j. Пусть точки \mathbf{r} и \mathbf{r}' таковы, что их координаты удовлетворяют условиям r_j\in(0,\epsilon_j) и {r_j}\ne {r_j'}. Тогда найдется элемент e\in U\mathfrak{E}^T_\Sigma U^{-1}, для которого выполнено e(\mathbf{r})=\mathcal{P}^i_{\gamma j} и e(\mathbf{r}')=O.
Доказательство. В силу попарной ортогональности проекторов в (4.11) при \mathscr{U}\in\mathbb N имеем
число \tau^i_{\gamma j}(r_j) встречается один раз. Выберем полином так, чтобы он обращался в 1 при t=\tau^i_{\gamma j}(r_j) и в 0 во всех остальных точках набора. Для e:=q(UE^T_\gamma U^{-1}), очевидно, имеем e(\mathbf{r})= \mathcal P^i_{\gamma j}, e(\mathbf{r}')=O. Лемма доказана.
Следствие 2. Если координаты точки \mathbf{r} не принимают крайних значений, т. е. r_j\notin\{0,\epsilon_j\} для всех j, то выполнено
Если координаты таковы, что для фиксированного j выполнено 0<a_j\leqslant r_j\leqslant b_j<\epsilon_j, а остальные координаты произвольны, то справедливо соотношение
Первое соотношение – прямое следствие утверждения леммы; второе легко выводится из первого с использованием предложения 4. Равенство (4.13) указывает на отсутствие связей между блоками [U\mathfrak{E}^T_\Sigma U^{-1}]_j при принятых ограничениях на координаты.
Если координаты точки \mathbf{r} принимают крайние значения, то равенство (4.12), вообще говоря, нарушается: в матричной алгебре [U\mathfrak{E}^T_\Sigma U^{-1}](\mathbf{r}) появляются связи между ее блоками. Именно они отличают левую и правую части во вложении (4.10): у состоящей из стандартных алгебр правой части такие связи отсутствуют. Поясним это на примерах.
Из определения и свойств функций \tau^i_{\gamma \Phi^j} (см. (2.19)) следует, что равенства \tau_{\gamma \Phi^j}^i(r_j)=\tau_{\gamma\Phi^{j'}}^{i'}(r_{j'}) возможны только при крайних значениях координат, т. е. при r_j\in\{0,\epsilon_j\} и r_{j'}\in\{0,\epsilon_{j'}\}. Пусть функции \tau^i_{\gamma\Phi^j} таковы, что выполнено \tau^i_{\gamma\Phi^j}(\epsilon_j)=\tau^{i+1}_{\gamma\Phi^j}(\epsilon_j)=\tau. В этом случае j-й блок эйконала в (4.8) примет вид
Поэтому в семейство образующих алгебры [U\mathfrak{E}^T_\Sigma U^{-1}](\mathbf{r}) проекторы P^i_{\gamma \Phi^j} и P^{i+1}_{\gamma \Phi^j} войдут не порознь, а в составе суммы P^i_{\gamma \Phi^j}+ P^{i+1}_{\gamma \Phi^j} и, таким образом, число образующих уменьшится на 1. Как следствие, возможно (и в содержательных случаях происходит) нарушение (4.12): вместо равенства гарантировано лишь вложение [U\mathfrak{E}^T_\Sigma U^{-1}](\mathbf{r})\subset\bigoplus_{j=1}^J \mathfrak{P}_{\Phi^j}.
Аналогичным образом могут связываться проекторы, входящие в разные блоки. Пусть \tau^i_{\gamma \Phi^j} таковы, что для каких-то разных j и j' выполнено \tau^i_{\gamma\Phi^j}(\epsilon_j)=\tau^{i'}_{\gamma\Phi^{j'}}(\epsilon_{j'})=\tau. В этом случае в (4.8) блоки с номерами j и j' примут вид
и окажутся связанными (через проекторы P^i_{\gamma \Phi^j} и P^{i'}_{\gamma \Phi^{j'}}), что также может привести к уменьшению числа образующих алгебры [U\mathfrak{E}^T_\Sigma U^{-1}](\mathbf{r}).
4.3. Приводимость
Алгебры \mathfrak{P}_{\Phi^j}, определяемые наборами проекторов \mathbb P_{\Phi^j} в (4.9), вообще говоря, приводимы: по предложению 3 имеем
Тем самым, приведение алгебры U\mathfrak{E}^T_\Sigma{ U^{-1}} к сумме неприводимых блоков в (4.14) сведется к адекватной группировке проекторов внутри каждого из наборов \mathbb P_{\Phi^j}. При группировке удобно перейти к новой нумерации.
Первым шагом является переход к сквозной нумерации алгебр, входящих в (4.14):
где L:={q_1+\dots+q_{J}}. Аналогично (т. е. формальной заменой \mathfrak{P} на \mathbb P в (4.17)) перейдем от наборов \mathbb P_{\Phi^j}^k к наборам \mathbb P_{l}, где l=1,\dots, L. Каждый набор \mathbb P_{l}, в свою очередь, можно разбить на поднаборы, отвечающие отдельным вершинам \gamma\in\Sigma:
где \varepsilon_1=\dots =\varepsilon_{q_1}=\epsilon_1; \varepsilon_{q_1+1}=\dots =\varepsilon_{q_1+q_2}=\epsilon_2; \dots, а представление эйконалов (3.10) и соотношение (4.12) примут согласованный с (4.18) вид
где координаты точки \mathbf{r} = \{r_1,\dots,r_{L}\} не принимают крайних значений.
Отметим, что именно при смене нумерации становится необходимым использовать индекс \gamma для функций \tau^i_\Phi (см. (2.26)). Это связано с тем, что исходные равные функции \tau_{\gamma \Phi^j}^i=\tau_{\gamma'\Phi^j}^i для фиксированного семейства \Phi^j могут соответствовать проекторам P_{\gamma \Phi^j}^i и P_{\gamma'\Phi^j}^i, которые после смены нумерации оказываются в разных блоках \mathfrak{P}_l и \mathfrak{P}_{l'}.
Как показывают следствие 2 и последующие комментарии, связи между блоками алгебры U\mathfrak{E}^T_\Sigma{ U^{-1}} возможны только на границах интервалов [0,\varepsilon_j]. Для детального описания этих связей удобно использовать следующий формализм.
Фиксируем вершину \gamma\in\Sigma и рассмотрим набор (k,l,r_l), отвечающий значению \tau_{\gamma l}^k(r_l) функции \tau_{\gamma l}^k. Скажем, что наборы (k,l,r_l) и (k',l',r_{l'}) связаны и будем писать (k,l,r_l)\leftrightarrow (k',l',r_{l'}), если выполнено \tau_{\gamma l}^k(r_l)=\tau_{\gamma l'}^{k'}(r_{l'}). По свойствам (2.19) и (2.20) такие равенства возможны только при крайних значениях параметров r_l\in\{0,\varepsilon_l\} и r_{l'}\in\{0,\varepsilon_{l'}\}; в приводимом ниже предложении фигурируют только такие значения. Эти же свойства ведут к следующим свойствам связи \leftrightarrow.
Предложение 6. Для заданного набора (k,l,r_l) выполняется одно и только одно из следующих условий:
1. Если (k,l,r_l)\in \mathbf 1, то для любого \widetilde{\mathbf{r}} с координатой \widetilde r_{l}\neq r_l существует элемент e\in U\mathfrak{E}^T_\Sigma U^{-1} такой, что выполнено
2. Если (k,l,r_l)\in \mathbf{2} и (k,l,r_l)\leftrightarrow(k',l,r_l), то для любого \widetilde{\mathbf{r}} с \widetilde r_{l}\neq r_l существует элемент e\in U\mathfrak{E}^T_\Sigma U^{-1} такой, что выполнено
3. Если (k,l,r_l)\,{\in}\, \mathbf 3 и (k,l,r_l)\,{\leftrightarrow}\,(k',l',r_{l'}), то для \mathbf{r}\,{=}\,\{r_1,\dots, r_l,\dots,r_{l'},\dots,r_\mathcal L\} и любого \widetilde{\mathbf{r}} с координатами \widetilde r_{j}\neq r_j, \widetilde r_{j'}\neq r_{j'} существует элемент e\in U\mathfrak{E}^T_\Sigma U^{-1} такой, что выполнено
Доказательство вполне аналогично доказательству леммы 2: как и последнее, оно сводится к выбору подходящего полинома q. Причина “склейки” проекторов в суммы та же, что и в рассмотренных ниже примерах следствия 2.
Как отмечалось, отличие алгебры U\mathfrak{E}^T_\Sigma U^{-1} от \bigoplus_{l=1}^L C([0,\varepsilon_l];\mathfrak{P}_l) состоит в возможных связях между блоками матричной алгебры (U\mathfrak{E}^T_\Sigma U^{-1})(\mathbf{r}), которые могут появиться, когда координаты r_l принимают крайние значения. Для изучения связей будет использована граничная алгебра
где \mathbf{0}=\{0,\dots,0\} и \boldsymbol{\varepsilon}=\{\varepsilon_1,\dots,\varepsilon_L\}. Для их описания сначала рассмотрим результаты общего характера, касающиеся структуры матричных алгебр типа \bigoplus_{l=1}^L \mathfrak{P}_l.
Пусть в гильбертовом пространстве \mathscr{G} со скалярным произведением \langle\,{\cdot}\,,{\cdot}\,\rangle задан набор одномерных проекторов P^1,\dots,P^n: P^i=\langle\,{\cdot}\,,\beta^i\rangle\beta^i, где \|\beta^i\|=1; обозначим
На наборе B:=\{\beta^1,\dots,\beta^n\} введем отношение \beta^i\sim_0\beta^{i'}, если \langle\beta^i,\beta^{i'}\rangle\neq 0. Оно определяет эквивалентность \beta^i\sim\beta^{i'}, если найдутся векторы \beta^{i_1},\dots,\beta^{i_k} такие, что \beta^i \sim_0\beta^{i_1}\sim_0 \cdots\sim_0\beta^{i_k}\sim_0 \beta^{i'}. Перенесем эту эквивалентность на проекторы, приняв P^i\sim P^{i'}, если \beta^i\sim\beta^{i'}.
Пусть B=B_1\,{\cup}\,{\cdots}\,{\cup}\, B_q есть разбиение на классы эквивалентности; обозначим \mathscr A_k:=\operatorname{span}B_k. Из определения отношения \sim следует \mathscr A_k\perp \mathscr A_l при k\neq l. Таким образом, имеем разложение \mathscr A=\mathscr A_1\oplus\dots\oplus \mathscr A_q, которое, очевидно, приводит все проекторы P^i.
Предложение 7. Алгебра \mathfrak P приводится подпространствами \mathscr A_k, причем справедливо разложение
где \mathfrak P|_{\mathscr{A}_k}\cong\mathbb{M}^{\varkappa_k}, \varkappa_k=\operatorname{dim}\mathscr A_k.
Приводимость очевидна; равенство для размерности следует из того, что любой \beta^i\in\mathscr A_k является циклическим в \mathscr A_k для части \mathfrak{P}|_{\mathscr A_k}. Для приведения \mathfrak P к виду (4.21) можно воспользоваться процедурой (2.11), (2.12).
Следующие рассмотрения моделируют ситуацию, которая встретится при изучении алгебры эйконалов. А именно, в абстрактной форме обсуждаются возможные связи между ее блоками в представлении (4.8).
Пусть имеются три гильбертовых пространства \mathscr{G}_k, k=1,2,3, в каждом из которых задан набор одномерных проекторов P^1_k,\dots,P^{n_k}_k\colon P^i_k=\langle\,{\cdot}\,,\beta^i_k\rangle\beta^i_k, \|\beta^i_k\|=1, где \beta^i_k суть векторы из наборов
где O_k – нулевой оператор, действующий в k-й компоненте пространства \mathscr{G}_1\oplus\mathscr{G}_2\oplus\mathscr{G}_3. Об алгебрах \mathfrak{P}_k в (4.22) будем говорить как о блоках алгебры \mathfrak{P}. Заметим, что в текущих рассмотрениях, грубо говоря, алгебры \mathfrak{P}_1 и \mathfrak{P}_2 отвечают какой-то паре выделенных блоков в (4.8), а \mathfrak{P}_3 – “все остальное”.
Скажем, что алгебра \mathfrak{Q}\subset\mathfrak{P}разделяет (не связывает) блоки \mathfrak{P}_1 и \mathfrak{P}_2 в (4.22), если вместе с любым элементом q_1\oplus q_2\oplus q_3\in\mathfrak{Q} в ней содержатся элементы q_1\oplus O_2\oplus q_3' и O_1\oplus q_2\oplus q_3'', где q_3', q_3'' – некоторые элементы \mathfrak{P}_3. В противном случае будем говорить, что она эти блоки связывает. Аналогично определяется связь (или ее отсутствие) для любой пары блоков в (4.22).
Отметим очевидный факт: если алгебра \mathfrak{Q} допускает систему образующих, каждая из которых имеет вид либо q_1\oplus O_2\oplus q_3', либо O_1\oplus q_2\oplus q_3'', то она разделяет блоки \mathfrak{P}_1 и \mathfrak{P}_2.
есть полный набор образующих алгебры \mathfrak{P}. Зададим на нем отображение (инволюцию) \mathcal{T}\colon \mathbb P\to\mathbb P такое, что если \mathcal{T}(\mathcal{P}_k^i)=\mathcal{P}_{k'}^{i'}, то \mathcal{T}(\mathcal{P}_{k'}^{i'})=\mathcal{P}_{k}^i, а также выполнено одно (и только одно) из условий:
Это отображение определяет разбиение множества \mathbb{P} на пары \{\mathcal P,\mathcal T(\mathcal P)\}, причем компоненты в каждой из пар либо совпадают, либо ортогональны друг другу. Видно, что такие \mathcal{T} найдутся, причем во множественном числе. Однако в алгебре эйконалов отображение \mathcal{T} окажется вполне конкретным и будет определяться значениями функций \tau^k_{\gamma l}(r) при r=0 и r=\varepsilon.
По данному отображению определим проекторы \mathcal{Q}_k^i:
среди которых могут быть совпадающие. Если \mathcal{T} не тождественно, то за счет совпадений их общее число, очевидно, будет меньше n_1+n_2+n_3. Образуем алгебру
Она определяется отображением \mathcal T. В то же время, разным допустимым \mathcal T и \mathcal T' может отвечать одна и та же алгебра \mathfrak{Q}. Этот факт используется ниже в доказательстве теоремы 1.
Далее обсуждается вопрос о том, при каких условиях введенная таким образом алгебра \mathfrak{Q} разделяет или связывает блоки \mathfrak{P}_1 и \mathfrak{P}_2.
Разобьем определенный в (4.24) набор проекторов \mathbb{P} на части \mathbb{P}_k:=\{\mathcal{P}_k^i\mid i=1,\dots,{n}_k\}, k=1,2,3. Каждой части сопоставим матрицы
Следующий результат о связи блоков используется ниже при изучении структуры алгебры эйконалов.
Теорема 1. Пусть блоки \mathfrak{P}_1 и \mathfrak{P}_2 алгебры (4.22) таковы, что каждый из отвечающих им наборов B_1 и B_2 является классом эквивалентности по отношению \sim. Тогда алгебра \mathfrak{Q} может связывать \mathfrak{P}_1 и \mathfrak{P}_2 только в том случае, если
Из вида проекторов (4.23) легко усмотреть, что связывать блоки \mathfrak{P}_1 и \mathfrak{P}_2 алгебра \mathfrak{Q} может только при условии \mathbb{Q}_{12}\neq \varnothing; в противном случае она их разделяет.
2. Предположим, что алгебра \mathfrak{Q} допускает представление
Тогда она разделяет блоки \mathfrak{P}_1 и \mathfrak{P}_2. В самом деле, в этом случае в алгебру \mathfrak{Q} вместе с элементами \mathcal Q^i_1=\mathcal{P}_1^i+\mathcal T(\mathcal{P}_1^i) и \mathcal Q^i_2=\mathcal{P}_2^i+\mathcal T(\mathcal{P}_2^i) порознь войдут все проекторы \mathcal{P}_1^i={P}^i_1\oplus O_2\oplus O_3 и \mathcal{P}_2^i=O_1\oplus{P}^i_2\oplus O_3. Вводя их в число образующих алгебры {\mathfrak{Q}} вместо элементов из \mathbb Q_{12}, легко убедиться в том, что разделение имеет место.
Здесь поясним следующее. Определение отображения \mathcal T исключает одновременное присутствие элементов \mathcal Q^i_k=\mathcal P^i_k+\mathcal T(\mathcal P^i_k) и \mathcal P^i_k в списке образующих \{\mathcal{Q}_k^i\mid i=1,\dots,{n}_k;\,k=1,2,3\} алгебры \mathfrak Q (см. (4.25)). Тем не менее, содержаться в самой алгебре они вполне могут. Это соответствует возможности заменить \mathcal T на другое отображение \mathcal T'\colon \mathbb{P}\to\mathbb P так, чтобы отвечающая ему алгебра \mathfrak Q' совпала с исходной \mathfrak Q.
то для алгебры \mathfrak{Q} справедливо представление (4.27) и, следовательно, она разделяет блоки \mathfrak{P}_1 и \mathfrak{P}_2. Последующие рассмотрения в силу равнозначности \mathfrak{P}_1 и \mathfrak{P}_2 вполне аналогичны для частей \mathbb{Q}_{11}\cup\mathbb{Q}_{13} и \mathbb{Q}_{22}\cup\mathbb{Q}_{23}. Проведем их для случая \mathbb{Q}_{11}\cup\mathbb{Q}_{13}\neq\varnothing. При этом имеем \mathbb{Q}_{12}\ne \varnothing, ибо в противном случае блоки \mathfrak{P}_1 и \mathfrak{P}_2 заведомо разделены.
Каждый элемент \mathcal{Q}_1^i\in\mathbb{Q}_{11}\cup\mathbb{Q}_{12}\cup\mathbb{Q}_{13} имеет вид
Каждому вектору \beta_1^i из набора B_1 отвечает проектор \mathcal{P}_1^i=P_1^i \oplus O_2\oplus{O}_3=\langle\,{\cdot}\,,\beta^i_1\rangle\beta^i_1\oplus O_2 \oplus O_3, который, в свою очередь, определяет соответствующий проектор \mathcal{Q}_1^i вида (4.28). Это позволяет определить отображение \mathbf{b}\colon B_1\to \mathbb{Q}_{11}\cup\mathbb{Q}_{12}\cup\mathbb{Q}_{13} по правилу
Заметим, что если \mathcal{Q}_1^i=\mathcal P^i_1+\mathcal T(\mathcal P^i_1)=\mathcal P^i_1+\mathcal P^{i'}_1\in\mathbb{Q}_{11}, то имеются два вектора \beta_1^i, \beta_1^{i'} таких, что \mathbf{b}(\beta_1^i)=\mathbf{b}(\beta_1^{i'})=\mathcal{Q}_1^i=\mathcal{Q}_1^{i'}. Далее под \mathbf{b}^{-1}(\,\cdot\,) понимается полный прообраз.
По условиям леммы B_1 есть класс эквивалентности по отношению \sim. Из условия \mathbb{Q}_{11}\cup\mathbb{Q}_{13}\neq \varnothing следует, что найдется пара векторов \beta_1^i\in\mathbf{b}^{-1}(\mathbb{Q}_{12}) и \beta_1^{i'}\in\mathbf{b}^{-1}(\mathbb{Q}_{11}\cup\mathbb{Q}_{13}) таких, что \langle\beta_1^i,\beta_1^{i'}\rangle\neq 0. В самом деле, отсутствие такой пары означало бы, что
что невозможно по определению эквивалентности \stackrel{*}{\sim}.
Выбранная пара векторов определяет проекторы \mathcal{Q}_1^i= \mathcal P^i_1+\mathcal T(\mathcal P^i_1)=\mathbf{b}(\beta^i_1)\in\mathbb Q_{12} и \mathcal{Q}_1^{i'}= \mathbf{b}(\beta^{i'}_1)\in\mathbb{Q}_{11}\cup\mathbb{Q}_{13}, а с ними и элемент
По выбору вектора \beta^i_1 имеем \mathcal{Q}_1^i=\mathbf{b}(\beta^i_1)\in\mathbb{Q}_{12}. Поэтому (\mathcal{Q}_1^i)_3=O_3, а значит, и (\widetilde{\mathcal{Q}}_1^i)_3=O_3. Аналогично, если \mathcal{Q}_1^{i'}\in\mathbb{Q}_{11}\cup\mathbb{Q}_{13}, то (\mathcal{Q}_1^{i'})_2=O_2, а значит, и (\widetilde{\mathcal{Q}}_1^i)_2=O_2. Таким образом, ненулевой может быть только компонента (\widetilde{\mathcal{Q}}_1^i)_1. В этой ситуации имеются две возможности.
1) Пусть \mathcal{Q}_1^{i'}\in\mathbb{Q}_{13}. Тогда (\mathcal{Q}_1^{i'})_1= P_1^{i'} и (\mathcal{Q}_1^i)_1= P_1^i, а (\widetilde{\mathcal{Q}}_1^i)_1 имеет вид
2) Пусть теперь \mathcal{Q}_1^{i'}\in\mathbb{Q}_{11}. Если \mathcal{Q}_1^{i'}=\mathcal{P}_1^{i'} , то рассмотрения, вполне аналогичные приведенным выше, ведут к тому же равенству (4.29). Рассмотрим случай, когда \mathcal{Q}_1^{i'}=\mathcal{P}_1^{i'} + \mathcal{T}(\mathcal{P}_1^{i'}), причем \mathcal{T}(\mathcal{P}_1^{i'})= \mathcal{P}_1^j\in\mathbb{P}_1 и \mathcal{P}_1^j\mathcal{P}_1^{i'}=O_1. Проектору \mathcal{P}_1^j соответствует вектор \beta_1^j\in B_1. Тогда (\widetilde{\mathcal{Q}}_1^i)_1 имеет вид
Сопоставляя (4.29) с (4.30), приходим к соотношению
\begin{equation*}
\widetilde{\mathcal{Q}}_1^i = c \mathcal{P}_1^i,\qquad c\geqslant\langle\beta_1^i,\beta_1^{i'}\rangle^2>0.
\end{equation*}
\notag
Это значит, что в алгебру \mathfrak{Q} по-отдельности входят проекторы \mathcal{P}_1^i\in\mathbb{P}_1 и \mathcal{T}(\mathcal{P}_1^i)\in\mathbb{P}_2. Об этом результате скажем, что проектор \mathcal Q_1^i=\mathcal{P}_1^i+\mathcal{T}(\mathcal{P}_1^i) распался на независимые (в алгебре {\mathfrak{Q}}) одномерные части \mathcal{P}_1^i и \mathcal{T}(\mathcal{P}_1^i).
Далее рассмотрим отображение \mathcal T'\colon {\mathbb P}\to\mathbb P, отличающееся от \mathcal{T} значениями только на двух проекторах \mathcal{P}_1^i и \mathcal{T}(\mathcal{P}_1^i), и положим
При этом алгебры \mathfrak{Q} и \mathfrak{Q}', определенные отображениями \mathcal{T} и \mathcal{T}', очевидно совпадают, а для \mathfrak{Q}' справедливо свое представление вида (4.26):
Таким образом, распад \mathcal Q^i_1 привел к тому, что часть \mathbb Q_{12}, отвечающая за связь блоков \mathfrak{P}_1 и \mathfrak{P}_2, уменьшилась на один проектор.
Повторяя рассмотрения для части \mathbb{Q}_{12}'\subset\mathbb{Q}_{12}, убедимся, что и из нее можно изъять очередной проектор, не меняя при этом алгебры \mathfrak{Q}. Продолжение процедуры за конечное число шагов приведет к распаду всех проекторов, содержавшихся в \mathbb{Q}_{12} и, как следствие, к представлению (4.27).
Итак, показано, что условие \mathbb{Q}_{11}\cup\mathbb{Q}_{13}\cup\mathbb{Q}_{22}\cup\mathbb{Q}_{23}= \varnothing необходимо для того, чтобы алгебра связывала блоки \mathfrak{P}_1 и \mathfrak{P}_2. Заметим, что условие \mathbb{Q}_{11}\cup\mathbb{Q}_{13}=\varnothing эквивалентно тому, что \mathcal{T}(\mathbb{P}_1)\subset\mathbb{P}_2, а условие \mathbb{Q}_{22}\cup\mathbb{Q}_{23}=\varnothing – тому, что \mathcal{T}(\mathbb{P}_2)\subset\mathbb{P}_1. Поскольку отображение \mathcal T инволютивно, отсюда получаем, что \mathbb{Q}_{11}\cup\mathbb{Q}_{13}\cup\mathbb{Q}_{22}\cup\mathbb{Q}_{23}= \varnothing выполняется тогда и только тогда, когда \mathcal{T}(\mathbb{P}_1)=\mathbb{P}_2.
4. С этого момента будем считать, что условие \mathcal{T}(\mathbb{P}_1)=\mathbb{P}_2 выполнено, а значит, для алгебры \mathfrak{Q}, связывающей блоки \mathfrak{P}_1 и \mathfrak{P}_2, справедливо представление
Отсюда видно, что \mathfrak{Q}=\mathfrak{Q}_{12}\oplus\mathfrak{Q}_{3} , где \mathfrak{Q}_{12}:=\vee \mathbb{Q}_{12}, \mathfrak{Q}_{3}:=\vee \mathbb{Q}_{33}, причем алгебра \mathfrak{Q} разделяет блоки \mathfrak{Q}_{12} и \mathfrak{Q}_3.
Далее мы уточним структуру алгебры \mathfrak{Q}_{12}. С этой целью удобно использовать матричную запись:
через векторы \beta_k^i, отвечающие проекторам P_1^i и P_2^{i'}. Теперь предположим, что {G}(\mathbb{P}_1)\neq{G}(\mathcal{T}(\mathbb{P}_1)). Отметим, что это возможно, только если в каждом из наборов B_1 и B_2 больше одного элемента. В этом предположении в алгебре \mathfrak{Q}_{12} найдутся такие \mathcal{Q}^i и \mathcal{Q}^j, что |\langle\beta_1^i,\beta_1^j\rangle|\neq|\langle\beta_2^{i'},\beta_2^{j'}\rangle|. Для произведения \mathcal{Q}^i\mathcal{Q}^j\mathcal{Q}^i\in\mathfrak{Q}_{12} получаем представление
Из них с учетом неравенства |\langle\beta_1^i,\beta_1^j\rangle|\neq|\langle\beta_2^{i'},\beta_2^{j'}\rangle| заключаем, что алгебра \mathfrak{Q}_{12} содержит каждый из проекторов
5. Таким образом, условие {G}(\mathbb{P}_1)={G}(\mathcal{T}(\mathbb{P}_1)) также необходимо для того, чтобы алгебра \mathfrak{Q} связывала блоки \mathfrak{P}_1 и \mathfrak{P}_2. При его выполнении справедливо представление
Укажем на важное обстоятельство. Если условия теоремы 1 выполнены, то наличие связи между блоками \mathfrak{P}_1 и \mathfrak{P}_2 исключает связь какого-либо из этих блоков с блоком \mathfrak{P}_3. Это следует из факта о разделении блоков \mathfrak{Q}_{12} и \mathfrak{Q}_3, отмеченного курсивом в начале п. 4 доказательства.
Следующий результат Д. В. Корикова [26], позволяет уточнить структуру алгебры \mathfrak{Q}_{12} из теоремы 1. Пусть выполнено равенство \mathcal{T}(\mathbb{P}_1)=\mathbb{P}_2 и соответственно n_1=n_2=:n. Выберем согласованную нумерацию в этих наборах:
Прямое применение критерия унитарной эквивалентности семейств подпространств из работы [26] приводит к следующему результату.
Лемма 4. Пусть наборы \mathbb{P}_1 и \mathbb{P}_2 удовлетворяют условиям теоремы 1. Тогда отображение \mathcal{T} продолжается с образующих \mathcal{P}_k^i до изометрии алгебр \mathcal{I}\colon \mathfrak{P}_1\to\mathfrak{P}_2, а алгебра \mathfrak{Q}_{12} имеет вид
Также несложно показать, что нарушение хотя бы одного из равенств леммы 4 приводит к тому, что алгебра \mathfrak{Q}_{12}не связывает блоки \mathfrak{P}_1 и \mathfrak{P}_2. Тем самым, одновременное выполнение условий теоремы 1 и леммы 4 гарантирует, что блоки связаны и существует изометрия \mathcal{I}, а невыполнение хотя бы одного из условий приводит к разделению соответствующих блоков алгеброй \mathfrak{Q}.
Рассмотрим более общую ситуацию. Пусть имеется N гильбертовых пространств \mathscr{G}_k, k=1,\dots,N. В каждом \mathscr{G}_k задан свой набор одномерных проекторов
определяемый набором векторов B_k:=\{\beta^i_k\mid i=1,\dots,{n}_k\}\subset \mathscr{G}_k, причем каждый B_k является классом эквивалентности по отношению \sim. Проекторы порождают алгебры
где {O}_k – нулевой оператор, действующий в k-й компоненте пространства \bigoplus\sum_{k=1}^N\mathscr{G}_k. О слагаемых \mathfrak{P}_k в (4.31) будем говорить как о блоках алгебры \mathfrak{P}_*.
алгебры \mathfrak P зададим отображение \mathcal{T}\colon\mathbb{P}\to\mathbb{P} такое, что если \mathcal{T}(\mathcal{P}_k^i)=\mathcal{P}_{k'}^{i'}, то \mathcal{T}(\mathcal{P}_{k'}^{i'})=\mathcal{P}_{k}^i, а также выполнено одно (и только одно) из условий:
и опишем ее структуру, используя результаты теоремы 1 и леммы 4.
Из блоков, составляющих алгебру \mathfrak{P} в (4.31), составим всевозможные пары \{\mathfrak{P}_k, \mathfrak{P}_{k'}\} с k\ne k' и отберем те из них, в которых компоненты связаны через алгебру \mathfrak{Q} (в том же смысле, что и \mathfrak{P}_1, \mathfrak{P}_2 в теореме 1 и лемме 4). Такой отбор однозначен, поскольку, как отмечалось после доказательства теоремы, каждый из \mathfrak{P}_k может быть связан лишь с одним \mathfrak{P}_{k'}. Перенумеруем блоки в (4.31), выделяя пары связанных блоков и независимые блоки:
таков, что \mathcal{T}(\mathcal{P}_{2N_1+j}^i)\notin\mathbb{P}_{2N_1+j}, то в алгебру \mathfrak{Q} по-отдельности входят проекторы \mathcal{P}_{2N_1+j}^i и \mathcal{T}(\mathcal{P}_{2N_1+j}^i). Это позволяет заменить \mathcal T|_{\mathbb P'} новым отображением \mathcal T'\colon {\mathbb P'}\to{\mathbb P'}, которое определяется на \mathcal{P}_k^i\in\mathbb{P}' по правилам
определяет ту же алгебру \mathfrak{Q}, что и \mathcal T, и в дополнение к (4.34) приводится независимыми блоками \widetilde{\mathcal T}(\mathbb P_{2N+j})=\mathbb P_{2N+j}. Внутри этих блоков оно либо действует тождественно, либо сопоставляет проектору ортогональный к нему.
Отображение \mathcal T определяет алгебру {\mathfrak{Q}} через проекторы-образующие Q^i_k согласно (4.32). Вполне аналогично, отображение \widetilde{\mathcal T} задает соответствующие проекторы \widetilde Q^i_k, генерирующие ту же алгебру. Особенность последних заключается в их форме: по построению имеем
т. е. все двумерные \widetilde{\mathcal{Q}}^i_k суть суммы проекторов, входящих в один и тот же блок \mathfrak{P}_k.
Обозначим \mathbb{Q}_k:=\{\widetilde{\mathcal{Q}}_k^i\mid i=1,\dots,n_k\}. Следующий результат резюмирует проведенные выше рассмотрения. Напомним, что блоки пронумерованы согласно (4.33).
Предложение 8. Для алгебры \mathfrak{Q} справедливо представление
где \mathfrak{Q}_k^{\mathrm{I}}:=\{A\oplus \mathcal TA\mid A\in\mathfrak{P}_{2k-1}\}\subset\mathfrak{P}_{2k-1}\oplus\mathfrak{P}_{2k} и \mathfrak Q^{\mathrm{II}}_k:=\vee \mathbb Q_k\subset \mathfrak P_k.
4.5. Каноническое представление
Используем установленные выше результаты для описания структуры алгебры эйконалов. Оно сведется к некоторому каноническому представлению.
Группировка проекторов в (4.16), проведенная в соответствии с разложением (4.15), есть разбиение наборов \mathbb P_{\Phi^j} на классы эквивалентности \mathbb P_{\Phi^j}^k (они же суть классы \mathbb P_{l}) по отношению \sim. Она мотивирована предложением 7 и подготавливает применение теоремы 1.
Роль алгебры \mathfrak P из теоремы 1 (см. (4.31)) играет алгебра
Зададим отображение \mathcal{T}\colon {\mathbb P}^\partial\to{\mathbb P}^\partial с помощью введенного ранее формализма, определяющего связи между наборами (k,l,r_l):
\begin{equation}
\mathcal{T}(\mathcal P_{\gamma l}^{k r_{l}}):= \begin{cases} \mathcal{P}_{\gamma l'}^{k' r_{l'}}, &\text{если } (k',l',r_{l'})\leftrightarrow(k,l,r_l), \\ \mathcal{P}_{\gamma l}^{k r_l}, &\text{если } (k,l,r_l)\text{ не связан ни с каким } (k',l',r_{l'}). \end{cases}
\end{equation}
\tag{4.36}
Несложно проверить, что набор \mathbb P^\partial и отображение \mathcal T удовлетворяют всем условиям теоремы 1, а определяемая ими согласно (4.25) алгебра \mathfrak Q совпадает с граничной алгеброй (4.20):
При приведении алгебры эйконалов к канонической форме будет использовано соединение ее блоков, связанных через граничную алгебру. Опишем эту конструкцию.
Пусть имеются две стандартные алгебры \mathfrak{A}=\dot C([0,\varepsilon];\mathfrak{P}) и \mathfrak{B}=\dot C([0,\varepsilon'];\mathfrak{P}') такие, что
Алгебры такого вида появятся при переходе к каноническому представлению в ситуации, когда в граничной алгебре \partial (U\mathfrak{E}^T_\Sigma U^{-1})\,{=}\,\mathfrak Q в представлении (4.35) имеются блоки типа \mathfrak{Q}_{k}^I, связывающие граничные значения \mathfrak{A}(\varepsilon) и \mathfrak{B}(0) какой-то пары блоков \mathfrak{A} и \mathfrak{B} алгебры эйконалов.
Для a\in\mathfrak{A} и b\in\mathfrak{B} таких, что a(\varepsilon)=\mathcal I b(0), определим элемент a\sqcup b\in C([0,\varepsilon+\varepsilon'],\mathfrak{P}) по правилу
Видно, что \mathfrak{A}\sqcup\mathfrak{B} есть подалгебра в C([0,\varepsilon+\varepsilon'];\mathfrak{P}), являющаяся стандартной алгеброй, для которой справедливо представление
\begin{equation*}
\mathfrak{A}\sqcup\mathfrak{B}=\{c\in C([0,\varepsilon+\varepsilon'];\mathfrak{P})\mid c(0)\in\mathfrak{A}(0),\, c(\varepsilon+\varepsilon')\in\mathcal I [\mathfrak{B}(\varepsilon')]\},
\end{equation*}
\notag
Заметим, что алгебры \mathfrak{A}^{\oplus}\mathfrak{B} и \mathfrak{A}\sqcup\mathfrak{B} изометрически изоморфны. Резюмируя рассмотрения, условимся говорить, что алгебры \mathfrak{A} и \mathfrak{B} допускают соединение через концыr=\varepsilon и r'=0. Отметим при этом, что вид алгебр \mathfrak{A}(0) и \mathfrak{B}(\varepsilon') на возможность соединения и его результат влияния не оказывают.
Очевидным образом изменяя определение (4.37), можно ввести соединение \mathfrak{A}\sqcup\mathfrak{B}\subset C([0,\varepsilon+\varepsilon'];\mathfrak{P})
\begin{equation}
\begin{aligned} \, &\text{через концы }r=\varepsilon \text{ и } r'=\varepsilon', \text{ если } \mathfrak{A}(\varepsilon)=\mathcal I[\mathfrak{B}(\varepsilon')]; \\ &\text{через концы }r=0 \text{ и } r'=\varepsilon', \text{ если } \mathfrak{A}(0)=\mathcal I[\mathfrak{B}(\varepsilon')]; \\ &\text{через концы } r=0 \text{ и } r'=0, \text{ если } \mathfrak{A}(0)=\mathcal I[\mathfrak{B}(0)]. \end{aligned}
\end{equation}
\tag{4.38}
На элементах алгебры \mathfrak{A}=\dot C([0,\varepsilon],\mathfrak{P}) определим транспозициюt\colon a\mapsto a^t, a^t(r):=a(\varepsilon-r), r\in[0,\varepsilon], и примем \mathfrak{A}^t:=\{a^t\mid a\in\mathfrak{A}\}. Изометрия {\mathcal M}\colon \mathfrak{P}\to\mathfrak{P} определяет преобразование \check{\mathcal M}\colon C([0,\varepsilon],\mathfrak{P})\to C([0,\varepsilon],\mathfrak{P}), (\check{\mathcal M}a)(r):=\mathcal M [a(r)], r\in[0,\varepsilon]. Для стандартной алгебры \mathfrak{A} примем \check{\mathcal M}\mathfrak{A}:=\{\check{\mathcal M} a\mid a\in\mathfrak{A}\}\subset C([0,\varepsilon],\mathfrak{P}).
Алгебры \mathfrak{A}^t и \mathcal M\mathfrak{A} также являются стандартными; легко видеть, что \check{\mathcal M}\mathfrak{A}\cong{\mathfrak{A}}^t\cong{\mathfrak{A}} и выполнено
Пусть блоки [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_l и [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_{l'} таковы, что их граничные значения [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_l(\varepsilon_l) и [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_{l'}(0) образуют один блок типа {\mathfrak{Q}}_k^I в граничной алгебре \partial(\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}):=\partial (U\mathfrak{E}^T_\Sigma U^{-1}) = \mathfrak Q (см. предложение 8). В этом случае отображение \mathcal T (см. (4.36)) определяет изометрию \mathcal I\colon \mathfrak{P}_{l'}\to\mathfrak{P}_l:
\begin{equation*}
\mathcal I (P^{k'}_{\gamma l'}) = P^{k}_{\gamma l },\quad \text{если}\quad \mathcal T(\mathcal P^{k' 0}_{\gamma l'})=\mathcal P^{k \varepsilon_l}_{\gamma l }.
\end{equation*}
\notag
Как следствие, корректно определена алгебра [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_l\sqcup[\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_{l'} – соединение этих блоков через концы r_l=\varepsilon_l и r_{l'}=0. Несложно проверить, что элементы
По свойствам отображения \mathcal T из существования связи между блоками следует, что n_{\gamma l}=n_{\gamma l'} для любого \gamma\in\Sigma, а соединения [\mathbf{U}_0 E^T_{\gamma}]_l\sqcup [\mathbf{U}_0 E^T_{\gamma}]_{l'} могут быть представлены в следующем виде:
Функции [\tau^k_{\gamma l}\sqcup\tau^{k'}_{\gamma l'}] являются непрерывными в силу равенства \tau^k_{\gamma l}(\varepsilon_{l})=\tau^{k'}_{\gamma l'}(0), определяющего связь между наборами (k,l,\varepsilon_{l}) и (k',l',0), которая, в свою очередь, определяет отображение \mathcal T и изоморфизм \mathcal I. Напомним, что каждая \tau^k_{\gamma l} есть линейная функция одного из двух видов:
\begin{equation*}
\text{либо }\tau^k_{\gamma l}(r)= t^k_{\gamma l}+r, \text{ либо } \tau^k_{\gamma l}(r) =t^k_{\gamma l }-r,
\end{equation*}
\notag
где t^k_{\gamma l }=\mathrm{const}\geqslant 0.
Пусть для определенности \tau^k_{\gamma l}(r)= t^k_{\gamma l}+r. Тогда из условия \tau^k_{\gamma l}(\varepsilon_{l})=\tau^{k'}_{\gamma l'}(0) и того, что равенство \tau^k_{\gamma l}(r_{l})=\tau^{k'}_{\gamma l'}(r_{l'}) возможно только в случае граничных значений параметров r_l и r_{l'} (в нашем случае r_l=\varepsilon_l, r_{l'}=0), следует, что \tau^{k'}_{\gamma l'}(r)= t^{k'}_{\gamma l'}+r, причем t^{k'}_{\gamma l'}= t^{k}_{\gamma l}+\varepsilon_l. Это приводит к следующим равенствам:
Аналогичные рассуждения справедливы и в случае \tau^k_{\gamma l}(r)=t^k_{\gamma l}-r. Таким образом, соединение функций \tau^k_{\gamma l} и \tau^{k'}_{\gamma l'} является линейной функцией того же вида, что и сами \tau^k_{\gamma l} и \tau^{k'}_{\gamma l'}. Из этого следует, что соединение блоков [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_l\sqcup[\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_{l'} есть стандартная алгебра с образующими того же вида, что и у исходных блоков.
В рассмотренном случае граничные значения [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_l(\varepsilon_l) и [\mathbf{U}_0 \mathfrak{E}^T_{\Sigma}]_{l'}(0) связаны. Вполне аналогично рассматриваются случаи других возможных связей между граничными значениями, допускающих соединение блоков (см. (4.38)).
Нетрудно видеть, что отображение \mathbf{U}_1 есть изометрия: изометричность алгебр \mathbf{U}_0\mathfrak{E}^T_{\Sigma} и \mathbf{U}_1\mathfrak{E}^T_{\Sigma} вполне аналогична изометричности \mathfrak{A}^{\oplus}\mathfrak{B} и \mathfrak{A}\sqcup\mathfrak{B}. Изометрия \mathbf{U}_1 переводит алгебру эйконалов в алгебру того же вида, что и \mathbf{U}_0\mathfrak{E}^T_{\Sigma}, но с меньшим на 1 количеством блоков; при этом в граничной алгебре \partial(\mathbf{U}_1\mathfrak{E}^T_{\Sigma}) блоков типа {\mathfrak{Q}}_{k}^I становится также на 1 меньше.
Последовательно заменяя аналогичным образом все пары связанных блоков на их соединения, придем к некоторому изометрическому изоморфизму \mathbf{U}_N\colon \mathfrak{E}^T_{\Sigma}\to\mathbf{U}_N\mathfrak{E}^T_{\Sigma} такому, что в граничной алгебре \partial(\mathbf{U}_N\mathfrak{E}^T_{\Sigma}) уже не содержится ни одного блока типа \mathfrak{Q}_{k}^I. Таким образом, отображение \mathbf{U}_N переводит алгебру эйконалов в алгебру той же структуры, что и \mathbf{U}_0\mathfrak{E}^T_{\Sigma}, но уже с независимыми блоками. Отметим также, что в процессе перехода к представлению \mathbf{U}_N\mathfrak{E}^T_{\Sigma} несколько (больше двух) блоков исходного представления \mathbf{U}_0\mathfrak{E}^T_{\Sigma} могут соединиться в один новый блок.
Шаг 3. На заключительном шаге еще раз используется неприводимость алгебр \mathfrak{P}_l. Для каждой из них подберем преобразование, реализующее изометрию \mathfrak{P}_l\,{\cong}\,\mathbb M^{\varkappa_l}. Как следствие, определится изометрия \mathbf{U}\colon \mathfrak{E}^T_{\Sigma}\,{\to}\, \bigoplus_{l=1}^\mathcal L \dot C([0,\zeta_l];\mathbb{M}^{\varkappa_l}), заданная на образующих равенствами
где \widetilde P_{\gamma l}^k\in\mathbb{M}^{\varkappa_l} суть одномерные (матричные) проекторы, попарно ортогональные для фиксированной вершины \gamma\in\Sigma, а \widetilde\tau_{\gamma l}^k – линейные функции одного из двух видов: либо \widetilde\tau^k_{\gamma l}(r)=\widetilde t^{\,k}_{\gamma l}+r, либо \widetilde\tau^k_{\gamma l}(r)=\widetilde t^{\,k}_{\gamma l}-r, r\in[0,\zeta_l], где \widetilde t^{\,k}_{\gamma l }\geqslant0 суть постоянные, а каждое \zeta_l есть сумма какого-то набора длин \varepsilon_k.
Таким образом, последовательными преобразованиями исходного параметрического представления (4.6) мы получаем представление такой же структуры, но уже с независимыми блоками, являющимися стандартными алгебрами. Сформулируем итоговый результат.
Теорема 2. Существует изометрия \mathbf U, доставляющая алгебре \mathfrak{E}_{\Sigma}^T и ее образующим-эйконалам представление
В нем \widetilde\tau_{\gamma l}^k – линейные функции от r_l\in[0,\zeta_l] такие, что |d\widetilde \tau_{\gamma l}^k/dr_l|=1, а их области значений суть сегменты длины \zeta_l, которые могут иметь разве что общие концы. \widetilde P_{\gamma l}^k\in \mathbb{M}^{\varkappa_l} суть проекторы, попарно ортогональные для каждого \gamma и такие, что выполнено \vee\{\widetilde P_{\gamma l}^k\mid k=1,\dots ,s_{\gamma l};\, \gamma\in\Sigma\}=\mathbb{M}^{\varkappa_l}.
Представление такого вида будем называть каноническим. Оно не единственно, но можно показать, что любые два таких представления отличаются друг от друга нумерацией блоков [\mathbf{U}\mathfrak{E}^T_{\Sigma}]_l, их транспозицией [\mathbf{U}\mathfrak{E}^T_{\Sigma}]_l\to[\mathbf{U}\mathfrak{E}^T_{\Sigma}]_l^t и изометриями [\mathbf{U}\mathfrak{E}^T_{\Sigma}]_l\to\check{\mathcal M}[\mathbf{U}\mathfrak{E}^T_{\Sigma}]_l. Неоднозначность, связанная с транспозицией, очевидно, отвечает двум направлениям изменения аргумента каждой из функций \widetilde\tau_{\gamma l}^k на интервале [0,\delta_l] (двум возможным параметризациям l-го блока).
В качестве комментария к теореме 2 отметим следующее. Исходное параметрическое представление алгебры эйконалов (4.8) отвечало разбиению части графа \Omega^T_\Sigma на семейства \Phi^1,\dots,\Phi^J. Мы полагаем, что переход к каноническому представлению соответствует некоторому новому разбиению. Это не доказано, но подкрепляется известными примерами [7], [8].
§ 5. Приведение к канонической форме
5.1. Каноническое представление и спектры
Главное достоинство канонического представления состоит в том, что из него легко извлекается полная информация о спектре алгебры \mathfrak{E}^T_\Sigma и определяется ряд ее инвариантов. Перейдем к их описанию.
С учетом отмеченной после теоремы 2 неоднозначности канонического представления будем считать, что нумерация и параметризация его блоков фиксированы. Само представление (4.39) перепишем в обновленных удобных для дальнейшего обозначениях (заменив \zeta_l на \varepsilon_l и s_{{\gamma l}} на n_{\gamma l} , а также убрав (\,\widetilde{\ }\,)):
суть временные клетки, соответствующие каноническому представлению. Правая часть представления (5.1) для оператора \mathbf{U} E_{\gamma}^T определяет его спектр и в силу изометричности \mathbf{U} спектр эйконала:
где каждый из сегментов [T^\gamma_{i-1},T^\gamma_i], в свою очередь, покрывается клетками \overline{\psi_{\gamma l}^k}, которые либо не пересекаются, либо имеют общие концы. По предложению 2 то же верно и для клеток \overline{\psi^i_{\gamma\Phi}}, связанных с параметрическим представлением (3.12). Сравнивая разложения (5.2) и (3.13), заключаем, что каноническому представлению соответствует новая (каноническая) нарезка спектра эйконала на временные клетки. Можно показать, что каждая новая клетка состоит из старых, т. е. переходу к каноническому представлению отвечает укрупнение временных клеток.
Как отмечалось в п. 4.1, спектр стандартной алгебры \dot C([a,b];\mathbb M^n) состоит из (классов эквивалентности) неприводимых представлений, отвечающих внутренним точкам сегмента [a,b] и кластеров (если таковые имеются) \{\widehat\pi^1_a,\dots,\widehat\pi^{n_a}_a\} и \{\widehat\pi^1_b,\dots,\widehat\pi^{n_b}_b\}, примыкающих к его концам. Такой спектр, оснащенный топологией Джекобсона, гомеоморфен пространству, которое естественно назвать сегментом с расщепленными концами. Оно описывается следующей конструкцией (см., например, [27]).
Рассмотрим n_a полусегментов [a,b) и n_b полусегментов (a,b] с топологией из \mathbb R. Отождествим внутренние точки всех полусегментов, имеющие одинаковые координаты. Образовавшееся факторпространство \mathscr{S}_{[a,b]} состоит из части S_{(a,b)}, гомеоморфной (a,b) и двух наборов попарно неотделимых точек (двух кластеров) K_a и K_b, которые состоят из n_a и n_b точек и отвечают концам a и b соответственно. Каждый из кластеров неотделим от S_{(a,b)}. Часть S_{(a,b)}:=\operatorname{int}\mathscr{S}_{[a,b]} есть множество внутренних точек, каждая из которых имеет окрестность, гомеоморфную (открытому) интервалу вещественной оси.
По первому из представлений (5.1), спектр алгебры \mathbf U \mathfrak{E}_{\Sigma}^T гомеоморфен дизъюнктному объединению сегментов
Каждый сегмент \mathscr{S}_{[0,\varepsilon_l]} характеризуется как максимальная компактная линейно-связная компонента пространства \mathscr{S}^T_\Sigma, а часть S_{(0,\varepsilon_l)}=\operatorname{int}\mathscr{S}_{[0,\varepsilon_l]} – как множество его внутренних точек.
Спектры алгебр \mathfrak{E}_{\Sigma}^T и \mathbf U \mathfrak{E}_{\Sigma}^T связаны гомеоморфизмом \mathbf U_*\colon \widehat{\mathfrak{E}_{\Sigma}^T}\to\widehat{\mathbf U \mathfrak{E}_{\Sigma}^T} (см. (4.2)). Следовательно, \widehat{\mathfrak{E}_{\Sigma}^T} гомеоморфен пространству \mathscr{S}^T_\Sigma и допускает представление
где \eta\colon\mathscr{S}^T_\Sigma\to\widehat{\mathfrak{E}_{\Sigma}^T} – гомеоморфизм, S_l=\eta(S_{(0,\varepsilon_l)}), \mathscr{K}^l_0=\eta(K^l_0), \mathscr{K}^l_{\varepsilon_l}=\eta(K^l_{\varepsilon_l}). Оно определяет разбиение спектра на максимальные компактные линейно-связные компоненты и, как таковое, имеет инвариантный топологический смысл. Последнее относится и к разбиению сегментов \mathscr{S}_l на внутренние точки S_l и кластеры \mathscr{K}^l_0, \mathscr{K}^l_{\varepsilon_l}. Множество внутренних точек спектра есть \operatorname{int}\widehat{\mathfrak{E}_{\Sigma}^T}=S_1\cup\dots \cup S_\mathcal L.
5.2. Координаты
Напомним, что алгебры \mathfrak{E}_{\gamma}^T=\vee\{{E}^T_{\gamma}\}\subset\mathfrak{E}_{\Sigma}^T, отвечающие отдельным эйконалам, мы называем парциальными. Через \pi и \widehat\pi обозначаются неприводимое представление и его класс эквивалентности.
Соответствие \varphi({E}^T_{\gamma})\leftrightarrow\varphi есть изометрия алгебр \mathfrak{E}_{\gamma}^T и C(\sigma_{\mathrm{ac}}( E^T_\gamma)), определяемая первым равенством в (4.5) (см. [23]). Каждая \mathfrak{E}_{\gamma}^T есть коммутативная подалгебра в \mathfrak{E}_{\Sigma}^T. Ее спектр (множество характеров) \widehat{\mathfrak{E}_{\gamma}^T} исчерпывается мерами Дирака:
при этом \widehat\delta_t=\{\delta_t\} (см. [25], [28]). Таким образом, каждому характеру отвечает точка (число) t из объединения сегментов (5.2), которую мы будем рассматривать как его \gamma-координату. Отметим также, что любое матричное приводимое представление алгебры \mathfrak{E}_{\gamma}^T имеет вид
Числа t_1,\dots,t_p, однозначно определяющие класс \widehat\rho, назовем его \gamma-координатами. Следующий шаг – координатизация спектра алгебры \mathfrak{E}_{\Sigma}^T.
Если коммутативная алгебра имеет конечное число образующих, то они составляют координаты на ее спектре (см. [28; гл. III, теорема 6]). Здесь используется адекватный аналог этого факта для некоммутативной C^*-алгебры, образующие которой суть самосопряженные операторы с простым спектром. А именно, каждому \widehat\pi\in\widehat{\mathfrak{E}_{\Sigma}^T} мы сопоставим набор \{\widehat\pi|_{\mathfrak{E}_{\gamma}^T}\mid \gamma\in\Sigma\} сужений на парциальные алгебры. Каждый элемент набора есть представление вида (5.4), уже снабженное \gamma-координатами t^1_{\gamma},t^2_{\gamma},\dots . Соответствие
доставляет требуемую координатизацию спектра. Уточним некоторые детали.
Соответствие (5.5), вообще говоря, не является инъективным, но оказывается таковым при сужении на \operatorname{int}\widehat{\mathfrak{E}_{\Sigma}^T}. Поскольку множество \widehat{\mathfrak{E}_{\Sigma}^T}\setminus\operatorname{int}\widehat{\mathfrak{E}_{\Sigma}^T}, состоящее из точек, образующих кластеры, конечно (так что, почти все точки спектра суть внутренние), то употребление термина координаты представляется оправданным.
Если \widehat\pi\in S_l\subset\mathscr{S}_l (см. (5.3)), то \mathbf U_*\widehat\pi есть внутренняя точка спектра \widehat{\mathbf U_*\mathfrak{E}_{\Sigma}^T} и ей, по первому из представлений (5.1), соответствует определенное значение параметра r\in (0,\varepsilon_l). В этом случае точку \widehat\pi обозначим через \widehat\pi_r и выясним, каковы ее \gamma-координаты.
Пусть \varphi\in C(\sigma_{\mathrm{ac}}( E^T_\gamma)). Из эквивалентных представлений, составляющих \widehat\pi_r, выберем \pi_r из условия
Сопоставляя начало и конец этой выкладки, заключаем, что представление \pi_r|_{\mathfrak{E}^T_\gamma} эквивалентно приводимому представлению \bigoplus\sum_{k=1}^{n_{\gamma l}}\delta_{t^k_{\gamma l}} алгебры \mathfrak{E}^T_\gamma, а числа \tau_{\gamma l}^k(r) суть \gamma-координаты точки \widehat\pi_r\in S_l\subset\operatorname{int}\widehat{\mathfrak{E}^T_\Sigma}. Они определены корректно, поскольку выбор другого \pi\ne \pi_r, \pi\in\widehat\pi_r приведет лишь к замене в (5.6) проекторов P^k_{\gamma l} на унитарно-эквивалентные им проекторы.
Как видно из (5.7), при вариациях точки \widehat\pi в S_l ее координаты t_{\gamma l}^k(\widehat\pi) заметают временные клетки \psi_{\gamma l}^k\subset\sigma_{\mathrm{ac}}( E^T_\gamma), и выполнено
что указывает на инвариантный характер этих клеток, а значит, и разбиения (5.2). Мы имеем ввиду, что они однозначно определяются самой алгеброй эйконалов; точнее – структурой (5.3) ее спектра. Длины клеток \varepsilon_l, l=1,\dots,\mathcal L, суть числовые инварианты алгебры {\mathfrak{E}^T_\Sigma}.
При фиксированном k соответствие \operatorname{int}\psi^k_{\gamma l}\ni t^k_{\gamma l}(\widehat\pi)\leftarrow\widehat\pi\in S_l биективно и определяет естественную параметризацию сегмента9[x]9Точнее, одну из двух возможных параметризаций.\mathscr{S}_l. Для его внутренних точек примем
где \tau^k_{\gamma l}(0) – левый конец клетки \psi^k_{\gamma l}. Затем распространим параметризацию на “концы” \mathscr{S}_l\setminus S_l по непрерывности, принимая для них r=0 и r=\varepsilon_l соответственно. Из (5.7) следует, что для внутренних точек выполнено r(\widehat{\pi}_r)=r.
Выбрав параметризацию описанным выше способом, мы, очевидно, параметризуем остальные клетки \psi^{k'}_{\gamma l}, k'\ne k, и определяем функции \tau^k_{\gamma}(r), r\in(0,\varepsilon_l), k=1,\dots,n_{\gamma l}, согласно (5.7). Они также являются инвариантами алгебры эйконалов.
5.3. Приведение
Каноническое представление (5.1) было получено “переформатированием” параметрического представления (4.8). Покажем, как прийти к нему, отправляясь от самой алгебры \mathfrak{E}^T_\Sigma и используя ее инварианты.
Итак, в нашем распоряжении эйконалы E^T_\gamma, \gamma\in\Sigma, и образованная ими алгебра \mathfrak{E}^T_\Sigma.
Шаг 1. Найдем спектр \widehat{ \mathfrak{E}^T_\Sigma}, оснастим топологией (Джекобсона) и выделим в нем сегменты и их компоненты согласно (5.3).
Шаг 2. Найдем спектры \sigma_{\mathrm{ac}}( E^T_\gamma) и введем на \operatorname{int}\widehat{ \mathfrak{E}^T_\Sigma}\gamma-координаты. Определим клетки \psi^k_{\gamma l} по (5.9). Параметризуем сегменты по (5.10).
Шаг 3. Для каждого l выберем представления \pi_r\colon \widehat{\mathfrak{E}^T_\Sigma}\to\mathbb M^{\varkappa_l}, \pi_r\,{\in}\,\widehat\pi_r\,{\in}\, S_l, r\in(0,\varepsilon_l), так, чтобы для всех e\in \mathfrak{E}^T_\Sigma\mathbb M^{\varkappa_l}-значные функции \pi_r(e) были непрерывны по r\in(0,\varepsilon_l) и продолжались по непрерывности10[x]10Такой выбор возможен уже в силу существования канонического представления. на [0,\varepsilon_l].
Шаг 4. Определим собственные значения \tau^k_{\gamma l}(r) и собственные проекторы P^k_{\gamma l} матриц \pi_r(E^T_\gamma) (см. (5.8)). По ним зададим отображение
и продолжим его с образующих-эйконалов до изометрии алгебр {\mathbf U}\colon \mathfrak{E}^T_\Sigma\to\bigoplus_{l=1}^\mathcal L \dot C([0,\varepsilon_l];\mathbb M^{\varkappa_l}).
Отображение \mathbf U приводит \mathfrak{E}^T_\Sigma к канонической форме (точнее к одной из ее версий: см. комментарий после теоремы 2).
Все полученные результаты относятся к смещенным эйконалам \dot E^T_\gamma (см. соглашение 2). Их переформулировка для исходных E^T_\gamma=\dot E^T_\gamma-P^T_\gamma очевидна и фактически сводится к замене в (5.1) функций \tau^k_{\gamma l} на \tau^k_{\gamma l}-1.
Комментарии
Отображение {\mathbf V}\colon \mathfrak{E}^T_\Sigma\to\bigoplus_{l=1}^\mathcal L \dot C(\mathscr{S}_l;\mathbb M^{\varkappa_l}), задаваемое на образующих соотношением
реализует элементы алгебры \mathfrak{E}^T_\Sigma как матрично-значные функции на ее спектре и, тем самым, доставляет инвариантную функциональную модель алгебры эйконалов. К этой модели, со всеми ее атрибутами – клетками \psi^k_{\gamma l}, функциями \tau^k_{\gamma l}, длинами \varepsilon_l, размерностями \varkappa_l – можно перейти, отправляясь от любой изометричеcкой копии алгебры \mathfrak{E}^T_\Sigma и используя процедуру Шаг 1–4. Это важно для обратной задачи, поскольку ее данные определяют одну из таких копий. Вопрос в том, в какой мере эти атрибуты определяют строение графа \Omega.
Следующее наблюдение, возможно, окажется полезным в обратной задаче. На точках спектра \widehat{\mathfrak{E}^T_\Sigma} введем отношение \widehat\pi\stackrel{\gamma}{\sim_0}\widehat\pi', если среди \gamma-координат этих точек имеются совпадающие, т. е. выполнено \tau^k_{\gamma l}(\widehat\pi)=\tau^{k'}_{\gamma l'}(\widehat\pi'). Затем определим эквивалентность \widehat\pi\sim\widehat\pi', если найдутся точки \widehat\pi_1,\dots,\widehat\pi_p\in\widehat{\mathfrak{E}^T_\Sigma} и вершины \gamma_1,\dots,\gamma_{p+1}\in\Sigma такие, что \widehat\pi\stackrel{\gamma_1}\sim_0\widehat\pi_1\stackrel{\gamma_2}\sim_0\cdots \stackrel{\gamma_p}\sim_0\widehat\pi_p\stackrel{\gamma_{p+1}}{\sim_0}\widehat\pi'. Можно показать, что при факторизации спектра по отношению \widehat\pi\sim\widehat\pi' отождествляются только точки, входящие в кластеры. При этом факторизация превращает \widehat{\mathfrak{E}^T_\Sigma} в пространство, гомеоморфное некоторому графу. В примерах прослеживается, что получившийся граф гомеоморфен заполненной волнами области \Omega^T_\Sigma, факторизованной по некоторому отношению, имеющему простой геометрический смысл.
В известных примерах [6], [8] кластеры появляются, когда внутренняя вершина графа \Omega перекрывается волнами, идущими, по крайней мере, от двух граничных вершин. Мы полагаем, что это общий факт. Также интересен вопрос: можно ли характеризовать наличие циклов в \Omega^T_\Sigma в терминах алгебры \mathfrak{E}^T_\Sigma (см. [6])? Вопрос открыт.
Список литературы
1.
М. И. Белишев, “Граничное управление и томография римановых многообразий (BC-метод)”, УМН, 72:4(436) (2017), 3–66; англ. пер.: M. I. Belishev, “Boundary control and tomography of Riemannian manifolds (the BC-method)”, Russian Math. Surveys, 72:4 (2017), 581–644
2.
M. Belishev, “Geometrization of rings as a method for solving inverse problems”, Sobolev spaces in mathematics, Int. Math. Ser. (N. Y.), III, Applications in mathematical physics, Springer, New York, 2009, 5–24
3.
M. I. Belishev, M. N. Demchenko, “Elements of noncommutative geometry in inverse problems on manifolds”, J. Geom. Phys., 78 (2014), 29–47
4.
M. I. Belishev, “Boundary spectral inverse problem on a class of graphs (trees) by the BC-method”, Inverse Problems, 20:3 (2004), 647–672
5.
M. I. Belishev, A. F. Vakulenko, “Inverse problems on graphs: recovering the tree of strings by the BC-method”, J. Inverse Ill-Posed Probl., 14:1 (2006), 29–46
6.
M. I. Belishev, N. Wada, “On revealing graph cycles via boundary measurements”, Inverse Problems, 25:10 (2009), 105011, 21 pp.
7.
M. I. Belishev, N. Wada, “A C^*-algebra associated with dynamics on a graph of strings”, J. Math. Soc. Japan, 67:3 (2015), 1239–1274
8.
M. I. Belishev, A. V. Kaplun, “Eikonal algebra on a graph of simple structure”, Eurasian J. Math. Comput. Appl., 6:3 (2018), 4–33
9.
Н. Б. Васильев, “C^*-алгебры с конечномерными неприводимыми представлениями”, УМН, 21:1(127) (1966), 135–154; англ. пер.: N. B. Vasil'ev, “C^*-algebras with finite-dimensional irreducible representations”, Russian Math. Surveys, 21:1 (1987), 137–155
10.
P. Niemiec, “Models for subhomogeneous C^*-algebras”, Colloq. Math., 166:1 (2021), 75–106; arXiv: 1310.5595
11.
S. Avdonin, P. Kurasov, “Inverse problems for quantum trees”, Inverse Probl. Imaging, 2:1 (2008), 1–21
12.
S. Avdonin, P. Kurasov, M. Nowaczyk, “Inverse problems for quantum trees II: recovering matching conditions for star graphs”, Inverse Probl. Imaging, 4:4 (2010), 579–598
13.
S. Avdonin, G. Leugering, V. Mikhaylov, “On an inverse problem for tree-like networks of elastic strings”, ZAMM Z. Angew. Math. Mech., 90:2 (2010), 136–150
14.
P. Kurasov, M. Nowaczyk, “Inverse spectral problem for quantum graphs”, J. Phys. A, 38:22 (2005), 4901–4915
15.
P. Kurasov, M. Nowaczyk, “Geometric properties of quantum graphs and vertex scattering matrices”, Opuscula Math., 30:3 (2010), 295–309
16.
P. Kurasov, “Graph Laplacians and topology”, Ark. Mat., 46:1 (2008), 95–111
17.
M. Nowaczyk, “Inverse spectral problem for quantum graphs with rationally dependent edges”, Operator theory, analysis and mathematical physics, Oper. Theory Adv. Appl., 174, Birkhäuser, Basel, 2007, 105–116
18.
В. А. Юрко, “О восстановлении операторов Штурма–Лиувилля на графах”, Матем. заметки, 79:4 (2006), 619–630; англ. пер.: V. A. Yurko, “On recovering Sturm–Liouville operators on graphs”, Math. Notes, 79:4 (2006), 572–582
19.
V. A. Yurko, “Inverse spectral problems for differential operators on arbitrary compact graphs”, J. Inverse Ill-Posed Probl., 18:3 (2010), 245–261
20.
V. A. Yurko, “An inverse problem for higher order differential operators on star-type graphs”, Inverse Problems, 23:3 (2007), 893–903
21.
G. Berkolaiko, P. Kuchment, Introduction to quantum graphs, Math. Surveys Monogr., 186, Amer. Math. Soc., Providence, RI, 2013, xiv+270 pp.
22.
А. В. Каплун, “Каноническое представление алгебры эйконалов трехлучевого графа”, Математические вопросы теории распространения волн. 51, Зап. науч. сем. ПОМИ, 506, ПОМИ, СПб., 2021, 57–78
23.
М. Ш. Бирман, М. З. Соломяк, Спектральная теория самосопряженных операторов в гильбертовом пространстве, 2-е изд., испр. и доп., Лань, СПб.–М.–Краснодар, 2010, 464 с.; англ. пер. 1-го изд.: M. S. Birman, M. Z. Solomjak, Spectral theory of self-adjoint operators in Hilbert space, Math. Appl. (Soviet Ser.), 5, D. Reidel Publ. Co., Dordrecht, 1987, xv+301 с.
24.
Ж. Диксмье, C^*-алгебры и их представления, Наука, М., 1974, 399 с. ; пер. с фр.: J. Dixmier, Les C^*-algèbres et leurs représentations, Cahiers Scientifiques, XXIX, 2ème éd., Gauthier-Villars, Paris, 1969, xv+390 pp.
25.
Дж. Мерфи, C^*-алгебры и теория операторов, Факториал, М., 1997, 336 с.; пер. с англ.: G. J. Murphy, C^*-algebras and operator theory, Academic Press, Inc., Boston, MA, 1990, x+286 с.
W. Arveson, An invitation to C^*-algebras, Grad. Texts in Math., 39, Springer-Verlag, New York–Heidelberg, 1976, x+106 pp.
28.
М. А. Наймарк, Нормированные кольца, 2-е изд., Наука, М., 1968, 664 с. ; англ. пер. 1-го изд.: M. A. Naimark, Normed rings, Wolters-Noordhoff Publ., Groningen, 1970, xvi+572 с.
Образец цитирования:
М. И. Белишев, А. В. Каплун, “Каноническое представление C^*-алгебры эйконалов метрического графа”, Изв. РАН. Сер. матем., 86:4 (2022), 3–50; Izv. Math., 86:4 (2022), 621–666
Д. В. Кориков, “О представлениях алгебры гармонических эйконалов”, Математические вопросы теории распространения волн. 54, Зап. научн. сем. ПОМИ, 533, ПОМИ, СПб., 2024, 124–139
М. И. Белишев, А. В. Каплун, “Канонические формы алгебры эйконалов метрического графа и его геометрия”, Краевые задачи математической физики и смежные вопросы теории функций. 50, Зап. научн. сем. ПОМИ, 519, ПОМИ, СПб., 2022, 35–66