НАУЧНОЕ ИЗДАНИЕ МГТУ ИМ. Н. Э. БАУМАНА
НАУКА и ОБРАЗОВАНИЕ
Эл № ФС77 - 48211. Государственная регистрация №0421200025. ISSN 1994-0408
электронный научно-технический журнал
Классификация и конструирование
обратимых линейных дифференциальных операторов
на одномерном многообразии
# 07, июль 2014
Б01:10.7463/0714.0718107
Четвериков В. Н.1,а
УДК 517.977
1Россия, МГТУ им. Н.Э. Баумана
Мтеи^епкоу. vl@yandex.ru
Исследуется задача описания обратимых линейных дифференциальных операторов с одной независимой переменной, актуальная в связи с задачами преобразования и классификации систем управления, в частности, с задачей плоскостности. Применяемый подход основан на сопоставлении каждому обратимому линейному дифференциальному оператору таблицы чисел и описании этих таблиц на наглядном элементарно-геометрическом языке. Приводятся математические структуры, которые необходимо задать для однозначного определения оператора, и алгоритм его построения. Обсуждаются способы применения и обобщения результатов работы.
Ключевые слова: обратимые линейные дифференциальные операторы; преобразования систем управления; спектральные последовательности цепных комплексов
Введение
Обратимые дифференциальные операторы возникают при решении многих математических задач [1, §2.3]. Среди них — задачи преобразования и классификации систем дифференциальных уравнений (или, более общо, дифференциальных объектов) с использованием С-преобразований. В классической теории систем рассматриваются преобразования, при которых зависимые и независимые переменные одной системы функционально выражаются через зависимые и независимые переменные другой системы. Они называются преобразованиями Ли [2, гл. 4], их теория построена достаточно подробно. Однако существуют обратимые преобразования, при которых переменные одной системы зависят не только от переменных второй системы, но и от производных зависимых переменных по независимым. Такие преобразования называются С-преобразованиями [2, гл. 6]. В случае одной зависимой переменной любое С-преобразование есть преобразование Ли. Однако это не так в случае нескольких зависимых переменных. Более того, как показывают примеры, в этом случае С-преобразований значительно больше, чем преобразований Ли. Отсутствие
удобного описания С-преобразований сдерживает развитие методов их применения в задачах классификации дифференциальных объектов. С-преобразования линейных систем интерпретируются как обратимые линейные дифференциальные операторы. В случае нелинейных систем линеаризация С-преобразования интерпретируется как обратимый линейный дифференциальный оператор [3]. Поэтому исследования обратимых линейных дифференциальных операторов следует рассматривать как первый шаг в описании С-преобразований как линейных, так и нелинейных систем.
Проблема описания С-преобразований стала наиболее острой в последние 20 лет в связи с развитием теории плоских систем управления. Плоскими называют системы, эквивалентные линейным управляемым системам относительно группы С-преобразований (см. [4]). Методы управления, развитые для линейных систем, переносятся на плоские системы. Многочисленные работы (см. ссылки в [5]) демонстрируют, что такие системы описывают многие природные явления и процессы. Поэтому задача проверки плоскостности систем является актуальной.
В работах [6, 3] получен критерий плоскостности систем на языке обратимых линейных дифференциальных операторов. А именно, система плоская, если существует обратимый линейный дифференциальный оператор, отображающий столбец 1-форм, определяемый системой, в столбец точных 1-форм.
В известной монографии [1, §2.3] можно найти и другие (иногда неожиданные) применения обратимых дифференциальных операторов, а также результаты, характеризующие, насколько общим является свойство обратимости для дифференциальных операторов.
Данная работа продолжает исследования, начатые в [7, 3]. Здесь мы рассматриваем обратимые линейные дифференциальные операторы на одномерном многообразии, т.е. с одной независимой переменной. Применяемый подход аналогичен подходу, используемому ранее в [8], и основан на сопоставлении каждому обратимому линейному дифференциальному оператору наборов целых чисел и описании этих наборов на наглядном элементарно-геометрическом языке. Таким образом, обратимому оператору ставится в соответствие элементарно-геометрическая модель, которую мы называем ^-схемой квадратов. При этом одному классу принадлежат те, и только те обратимые операторы, которые имеют одинаковые наборы чисел. Обратимый оператор определяется этим набором неоднозначно. Математические структуры, которые необходимо задать для его однозначного определения, и алгоритм его конструирования, т.е. построения по ^-схеме квадратов и по этим математическим структурам, также приводятся.
Полученное описание обратимых линейных дифференциальных операторов может быть использовано следующим образом. Условия на обратимый оператор, которые следуют из постановки решаемой задачи, переформулируются на языке элементарно-геометрической модели. Ищется модель, удовлетворяющая этим условиям (или доказывается, что такой модели нет). Искомый оператор восстанавливается по этой модели.
Статья организована следующим образом. В разд. 1 определяются обратимые линейные дифференциальные операторы и приводятся примеры таких операторов. Элементарно-гео-
метрические модели обратимых операторов определяются в разд. 2. Основные результаты формулируются в разд. 3, а доказываются в разд. 4. В заключении, в частности, формулируются пути возможных обобщений достигнутых результатов.
1. Обратимые линейные дифференциальные операторы
С целью возможных обобщений в будущем приведем наиболее общее — алгебраическое — определение линейных дифференциальных операторов (подробное изложение алгебраической теории таких операторов см. в [2, гл. 1]). Пусть А — коммутативная алгебра над полем к характеристики нуль.
Рассмотрим два А-модуля Р и Q и группу Нотк (Р, ф) к-гомомофизмов из Р в ф. В этой группе введем две А-модульные структуры, определенные умножениями
(аД)(р) = аД(р), (а+Д)(р) = Д(ар), а е А, р е Р, Д е Нотк(Р,ф). (1)
ак?
Обозначим ¿а(Д) = а+Д — аД, а е А. Также введем обозначение ¿ао,,,,ак = ¿ао о ... о
ао, ..., ак е А.
Гомоморфизм Д е Нотк(Р, ф) называют линейным дифференциальным оператором над алгеброй А порядка не выше к, если ¿ао,...,ак (Д) = 0 для любого набора а0, ..., ак элементов из А.
Линейный дифференциальный оператор из А в А (т.е. Р = ф = А) называют скалярным.
Скажем, что линейный дифференциальный оператор Д имеет порядок к, если он порядка не выше к, но не является линейным дифференциальным оператором Д порядка не выше к — 1. Порядок линейного дифференциального оператора Д будем обозначать огё Д.
В случае, когда к = К, А является алгеброй С^(М) гладких (бесконечно дифференцируемых) функций на многообразии М, а Р и ф — модулями гладких сечений двух локально тривиальных векторных расслоений £, ( над М (о теории расслоений см. [9, гл. 2 и 3]), сформулированное определение эквивалентно обычному определению линейного дифференциального оператора. В частности, скалярный оператор на одномерном многообразии М представляет собой конечную сумму производных по координате Ь на М с коэффициентами из А:
Д = ак(Ь)^ + ак- 1(Ь)+ ... + ао(Ь), ак(Ь) ф 0, к = огёД.
В многомерном же случае линейный дифференциальный оператор можно представить как матрицу скалярных операторов.
Обозначим через В1Й+(Р, ф) множество линейных дифференциальных операторов из Р в ф степени не выше к, рассматриваемое как А-модуль относительно второго (правого) умножения, определенного в (1). Заметим, что если Р, ф, Я — А-модули, Д е В1й+(Р, ф), V е (ф, Я), то V о Д е В1й++г(Р, Я). Это утверждение вытекает из формулы
8а(Ч о Д) = О Д + V о §а (Д)
и того факта, что отображение 8а уменьшает порядок оператора на единицу.
Отметим также, что в общем случае композиция линейных дифференциальных операторов некоммутативна, даже если операторы скалярные.
Дифференциальный оператор А е Нотк(Р, Q) называют (двусторонне) обратимым, если существует такой дифференциальный оператор А-1 е Нотк( Р, Q), что оператор А-1 о А есть тождественное отображение модуля Р, а оператор А о А-1 — тождественное отображение модуля Q. В этом случае оператор А-1 называют обратным к А.
Отметим, что данное определение симметрично относительно замены А на А-1. Поэтому оператор А-1 также обратим, а А — обратный к нему оператор.
В этой работе рассматриваются обратимые линейные дифференциальные операторы на одномерном многообразии М (с одной независимой переменной). Рассуждения локальны, поэтому можно считать, что М = М — действительная прямая, а Р и Q — модули векторных функций на М. Нетрудно показать, что если существует обратимый линейный дифференциальный оператор из Р в Q, то Р и Q имеют одинаковую размерность, которую мы обозначаем через т, т.е. Р = Q = Ат, где А = С~(М).
Приведем два примера обратимых операторов рассматриваемого типа. В обоих примерах а, Ь, с — произвольные скалярные линейные дифференциальные операторы.
Пример 1. В трехмерном случае, т.е. т = 3, операторы А и А-1, заданные матрицами
/ 0 с А ( -а 0 1 \
А
b 10, А-1 = ba 1 -b \ 1 ac a ) \ 1 — cba —c cb )
являются взаимно обратными, в чем легко убедиться непосредственным перемножением матриц.
Пример 2. В двумерном случае операторы А и А-1, заданные матрицами
А ( 1 + bc —b\ А-1 =( 1 + ba b
\ —a — c — abc 1 + ab I la + c + cba 1 + cb
также являются взаимно обратными.
2. Схемы квадратов
Пусть M — одномерное многообразие, A = Crx(M), P, Q — A-модули гладких сечений двух локально тривиальных векторных расслоений, А: P ^ Q — дифференциальный оператор порядка l. Рассмотрим следующие A-модули:
Gp = Diff+(A,Q), Fk = {a G Diff++1(A, Q): a = А о в, в G Diff+(A, P)} , p, k > 0.
Так как Gp и Fk — подмодули модуля Diff+(A, Q), где s = max(p, k + l), то их пересечение Fk П Gp есть также подмодуль модуля Diff+(A, Q). Модули Gp, Fk являются модулями гладких сечений локально тривиальных векторных расслоений над M. Обозначим соответствующие расслоения через и (k. Слои этих расслоений над точкой t G M состоят из
линейных дифференциальных операторов в точке t и имеют постоянную размерность, не зависящую от t. Пересечение модулей Fk П Gp является модулем гладких сечений пересечения векторных расслоений £p и (k. Слои пересечения векторных расслоений могут иметь разную размерность, зависящую от точки t G M. Такой объект называют векторным расслоением с особенностями. Под размерностью модуля R гладких сечений векторного расслоения с особенностями над многообразием M мы понимаем целочисленную функцию, которая точке t G M ставит в соответствие размерность слоя этого расслоения над t. Обозначим ее через dim R.
Пусть оператор A: P ^ Q обратим, l = ord А, а L — порядок обратного оператора. Точку t G M будем называть d-регулярной точкой обратимого оператора А, если в некоторой окрестности этой точки пересечение слоев £p и (k имеет постоянную размерность для любых p = 0, 1, ..., l и k = 0, 1, ..., L. Таким образом, в окрестности d-регулярной точки функция dim(Fk П Gp) постоянна для указанных k и p. Далее будем рассматривать обратимые операторы только в окрестности d-регулярных точек.
Каждый линейный дифференциальный оператор порождает последовательность целых неотрицательных чисел dk,p = dim(Fk П Gp), p, k > 0. В дальнейшем анализе будем относить обратимые линейные дифференциальные операторы к одному классу, если у этих операторов соответствующие последовательности совпадают.
Последовательность {dk,p} однозначно определяет последовательность чисел
Pk,p = Kk,p - Kk-I,p-1, к, p > 0, (2)
где
Kk,p dk,p dk-i,p dk,p-1 + dk—i,p—i
(полагаем, что dk,p = Kk,p = 0, если к < 0 или p < 0).
Нетрудно доказать, что по последовательности {pk,p} последовательность {dk,p} восстанавливается однозначно: каждое число dk>p есть сумма чисел pkl,pi с k1 < к и p1 < p, умноженных на некоторые целочисленные коэффициенты.
Отметим, что с ростом индексов k и p числа dk>p возрастают, при этом pk,p = 0 для больших k, p (см. ниже).
Вернемся к примеру 1, полагая, что а — оператор первого порядка, а b и c — операторы второго порядка. Для оператора А найдем числа dk,p, а затем вычислим величины pk,p по формулам (2). В результате получим:
P0,1 = P0,2 = P0,3 = P2,0 = P4,0 = P5,0 = 1, P5,1 = P4,2 = P2,3 = —1,
а остальные числа Pk,p равны нулю.
Аналогично можно найти последовательность {Pk,p} для оператора А из примера 2. Ее ненулевые значения таковы:
P0,1 = P0,2 = P0,3 = P2,0 = P4,0 = P5,0 = 11, P5,1 = P4,2 = P2,3 = — 1.
рис.1
рис.2
Последовательность рк,р полезно изображать графически, приписав каждой точке с координатами (к, р) числорк,р. Нарис. 1 представлена последовательность {рк,р} для обратимого дифференциального оператора Д из примера 1: черные целочисленные точки отвечаю ненулевым значениям последовательности, белые — нулевым. Возможные значения ненулевых членов последовательности, как показано на рисунке, в данном случае равны 1 и — 1.
На рис. 2 аналогичным образом изображена последовательность {рк,р} для обратимого оператора из примера 2.
Рис. 1. Значения рк,р для оператора из примера 1
Рис. 2. Значения рк,р для оператора из примера 2
Отметим следующие особенности расположения ненулевых членов последовательности {рк,р}:
• все ненулевые значения последовательности лежат на осях координат, а также на вертикальных и горизонтальных прямых, проходящих через точки со значениями — 1;
• на этих прямых можно построить систему квадратов (на рисунках они выделены серым) так, что все точки с ненулевым значением будут расположены на сторонах этих квадратов, причем каждая точка со значением — 1 будет расположена в верхнем левом углу квадрата, так что количество квадратов совпадает с количеством значений — 1 в последовательности {рк,р};
• сумма значений точек, расположенных на каждой из осей координат, равна размерности т модулей Р и ф (в данном случай порядку матрицы, изображающей оператор);
• сумма значений точек, расположенных на вертикальной или горизонтальной прямой, проходящей через точку со значением —1, равна нулю.
Отметим также, что размеры квадратов на рис. 1 и 2 совпали с заданными порядками скалярных операторов а, Ь, с.
Оказывается, что эти особенности расположения чисел рк,р — общие для всех обратимых линейных дифференциальных операторов на одномерном многообразии. Для описания таких последовательностей введем некоторые понятия.
Рассмотрим конечный набор квадратов, расположенных в первой четверти координатной плоскости, у которых вершины имеют целочисленные координаты, а стороны параллельны осям координат. Верхнему правому (т.е. дальнему от начала координат) углу каждого квадрата припишем значение — 1, а нижнему левому (т.е. ближнему от начала координат) углу — значение 1. Если есть точки, являющиеся вершинами нескольких квадратов, значения, назначенные по каждому квадрату, в этих точках суммируем. Остальным точкам первой четверти с целочисленными координатами припишем нулевые значения. Получим последовательность чисел рк,р, к > 0, р > 0. Предположим, что существует такая последовательность {ак,р} целых неотрицательных чисел, что последовательности чисел рк,р = рк,р + ак,р удовлетворяет следующим условиям:
^ те
XРк,р = 0 (р > 0), ^Рк,р = 0 (к > 0); (3)
к=0 р=0 те те
X рМ = X] ро,р = т; (4)
к=0 р=0
р-1
X Р,р > *к,р, X РкЛ > **,р (к > 0, р > 0). (5)
¿=0 ¿=0
Здесь гк р — количество квадратов с верхним правым углом в точке (к, р).
Указанный конечный набор квадратов будем называть ^-схемой квадратов, а соответствующую этому набору последовательность {рк,р} — т-таблицей ^-схемы квадратов.
Отметим, что ^-схема квадратов определяет последовательность {рк,р} неоднозначно. Во-первых, можно добавить произвольное неотрицательное число п к значению в начале координат, что не нарушит условий (3)-(5). Но при этом сумма значений по осям координат увеличится на п. В результате из т-таблицы получаем (т+п)-таблицу. Во-вторых, есть примеры ^-схем квадратов, которым соответствует несколько т-таблиц. Это возможно, когда есть несколько вариантов выбора последовательности {ак,р}, используемой для получения последовательности {рк,р} (в рассмотренных выше примерах 1 и 2 ^-схема квадратов имеет единственную т-таблицу, где т = 3 в примере 1 и т = 2 в примере 2).
3. Основные результаты
Теорема 1. Пусть М — одномерное многообразие, А = Сте(М), Р = Q = Ат. Тогда:
а) для всякого обратимого линейного дифференциального оператора А € Нотк(Р, Q) последовательность {рк,р}, построенная согласно формулам (2), в окрестности ^-регулярной точки совпадает с т-таблицей некоторой ^-схемы квадратов;
б) для всякой ^-схемы квадратов существует такой обратимый линейный дифференциальный оператор А € Нотк(Р, Q), для которого последовательность {рк,р} совпадает с т-таблицей данной ^-схемы квадратов.
Отметим, что обратимый оператор определяется ^-схемой квадратов неоднозначно.
Покажем, как построить обратимый оператор по заданной схеме и какие структуры для этого еще нужно задать. Пусть задана ^-схема квадратов и для нее построена последовательность чисел рк,р, удовлетворяющих условиям (3)-(5). Обозначим через Z множество верхних правых углов квадратов ^-схемы, а через В — множество нижних левых углов этих квадратов и точек плоскости (к, р), в которых ак,р = 0, причем точка (к, р) включается в В ак,р раз. Из равенств (3)-(4) следует, что сумма положительных чисел ак,р равна т. Поэтому множество В имеет на т элементов больше, чем множество Z. Обозначим через Z0,p подмножество элементов Z с координатами (к, р), а через В0,р — аналогичное подмножество в В.
Построим три инъективных отображения ф, ' из множества Z в множество В. Отображение ф ставит в соответствие дальнему углу квадрата ^-схемы его ближний угол. Отображение ^ элементу г е Z с координатами (к, р) ставит в соответствие некоторый элемент из В с координатами (к — в, р), где 5 > 0. Аналогично отображение ' элементу из Z с координатами (к, р) ставит в соответствие элемент из В с координатами (к, р — в), в > 0. Существование таких отображений доказано ниже (см. разд. 4).
Поставим в соответствие каждому квадрату ^-схемы скалярный дифференциальный оператор, порядок которого равен длине стороны квадрата. Этот оператор можно выбрать произвольно, но его порядок определен квадратом однозначно. Между квадратами ^-схемы и элементами множества Z имеется взаимно-однозначное соответствие, так что каждый квадрат можно считать элементом множества Z. Через П обозначим дифференциальный оператор, соответствующий квадрату г е Z.
Каждому элементу Ь е В поставим в соответствие дифференциальный оператор Уь, который построим индуктивно по второй координате р элемента Ь. В случае р = 0 оператор Уь построим следующим образом. Модуль С0 = Б1й+(А, Q) изоморфен Q. Пусть е1, .. ., ет — какой-либо базис модуля Q (он существует, так как по условию Q = Ат). Обе координаты любого элемента из Z положительны. Поэтому число рк,0 совпадает с количеством элементов в подмножестве В0О. Из первого равенства (4) следует, что множество В имеет т элементов с нулевой второй координатой. Обозначим их через Ь1, . .., Ьт. Аналогично доказывается, что множество В имеет т элементов с нулевой первой координатой. Обозначим эти элементы через Ь1, . .., Ь^. Отметим, что элементы множества В0О обозначены двумя разными способами. Однако это не является препятствием для дальнейших рассуждений. Для каждого г = 1, ..., т элементу Ьг поставим в соответствие дифференциальный оператор , который каждому а е А ставит соответствие аег.
Рассмотрим теперь элемент Ь е В со второй координатой р > 0. Ниже (см. разд. 4) доказано, что определен элемент ^-1(Ь) из Z, а элементы (ф о ^-1)(Ь) и ('^ъ о ^-1)(Ь) принадлежат В и их вторые координаты меньше р. По предположению индукции им уже поставлены в соответствие дифференциальные операторы. Определим оператор Уь равенством
Таким образом, каждому элементу Ь Е В поставлен в соответствие дифференциальный оператор Уь: А ^
Пусть /1, . .., /т — какой-либо базис модуля Р. Заданной ^-схеме соответствует обратимый линейный дифференциальный оператор А: Р ^ определенный формулой
А(а/ + ... + От/т) = УЬ1 (01) + ... + УЬ1 (ат),
а1,
ат Е А.
(7)
4. Доказательство теоремы 1
Введем алгебраические структуры, необходимые для описания модулей П Ср. Рассмотрим фактормодули С^/С^ь г > 0 (здесь С— = 0). Обозначим через [а]* элемент фактормодуля С^/С^ь соответствующий дифференциальному оператору а € С*. Пусть
£ — координата на М в окрестности рассматриваемой ^-регулярной точки. Для оператора
й
а € С* обозначим через Да композицию а с производной по ¿: Да = а о —. Тогда Да Е С*+1 и для всех г > 0 определены отображения Д: С* ^ С*+1. Эти отображения индуцируют отображения Д: Ср/Ср-1 ^ Ср+1 /Ср согласно формуле
Д[а]г
d'
а о —
Р > 0.
(8)
р+1
Покажем, что отображения (8) являются гомоморфизмами модулей. Оператор умножения на функцию а Е А будем также обозначать через а. Если а Е Ср, то [а]р+1 = 0. Поэтому для любой функции а Е А справедливо равенство
Д(а+[а]р) = а+Д[а]р,
поскольку
Д(а+[а]р) = Д[а о а]р =
d'
а о а— dt
р+1
d ^а
ао-оа — а о — dí dí
d
а о — dt
р+1 dаx +
р+1
— ( ^ ) [а]р+1 = а+Д[а]р.
Гоморфизм (8) для каждого р > 0 является изоморфизмом, так как фактормодуль Ср/Ср-1
состоит из элементов вида
йр+1
а
(Р_
а Е С0, а гомоморфизм (8) отображает такой элемент
в элемент
а
^р+1
фактормодуля Ср+1/Ср. Таким образом, для любого р фактормодуль
Ср+1/Ср изморфен С0/С-1, а следовательно,
Для описания модулей П Ср используем теорию цепных комплексов и их спектральных последовательностей (см. [10, гл. 4, §1 и гл. 9, §1]). При этом для удобства все используемые понятия далее формулируются, а необходимые факты доказываются. Изоморфизмы (8) определяют точные комплексы
0
Б
С.
р+1
Ср-1
Ср
0, р > 0.
(9)
р
Предположим, существует такое натуральное число Ь, что Ор С Ер+ь для р > 0. Модули Ор/Ор-1 имеют фильтрацию
Ор Ор П + Ор -1 Ор П Ек + Ор-1 Ор П Е0 + Ор-1
Ор-1 Ор-1 Ор-1 Ор-1
Так как композиция с производной по Ь оператора из Ек преставляет собой оператор из Ек+ь то изоморфизмы (8) сохраняют построенные фильтрации:
Ор П Рк + Ор-1 ^ Ор+1 П Ек+1 + Ор р > о (10)
О р - 1 О р
А это значит, мы имеем фильтрацию комплексов (9) и их спектральные последовательности (см. [10, гл. 9, §1, теорема 2]). Элементы Е0 этих последовательностей состоят из модулей
Е>0 _ Ор П рк + Ор-1
' Ор П Ек-1 + Ор-1'
Точнее, комплексу (9) с номером р соответствуют две строки Ек0р и Е^ комплексу (9) с номером р +1 — две строки, но Е<0+1р и Е°+2,р и т.д. При этом вторая строка предыдущего комплекса совпадает с первой строкой последующего комплекса, поэтому эти строки обозначены одинаково.
Отметим, что мы используем отличную от обычной нумерацию элементов Е0, потому что для дальнейших рассуждений удобно иметь нумерацию, симметричную относительно р, к. Для доказательства симметричности Е^0р используем следующий известный факт.
Теорема 2 (Нетер об изоморфизме [10, Введение, §4]). Пусть X и У — некоторые подмодули А-модуля Z, и пусть X + У — подмодуль, порожденный множеством X и У.
Тогда вложение X С X + У переводит X П У в У и индуцирует изоморфизм А-модулей
X X + У
X П У и У .
Обозначим X _ Ор П Ек и У _ Ор П Ек-1 + Ор-1. Так как Ек-1 С Ек, то X + У _ Ор П Ек + Ор-1, X П У _ Ор П Ек-1 + Ор-1 П Ек. Применяя теорему Нетер, получаем симметричную относительно индексов к и р формулу
Е0 р -Т7Ор П Ек-ет, к, р > 0. (11)
Ор П Гк-1 + Ор-1 П рк
Используя известные факты о размерностях пространств (модулей), получаем
dim E0kp = dim(Gp П Fk) - dim(Gp П Ffc_i + Gp_i П Fk) =
= dk,p - dim(Gp П Fk_i) - dim(Gp_i П Fk) + dim(Gp П Fk_i П Gp_i П Fk) =
= dk,p — dk_i ,p — dk ,p_i + dk_i ,p_i = Kk,p. (12)
Отметим, что до сих пор мы никак не использовали обратимость оператора А. Предположим теперь, что А — обратимый оператор и A_i: Q ^ P — обратный к нему оператор
некоторого порядка L. Тогда Gp С Fp+L для p > 0, так как если а Е Gp, то а = А о в, где в = А-1 о а Е Diff++L(A, P). Таким образом, применимы рассуждения, приведенные выше. Более того, оператор А можно заменить оператором А-1 и повторить эти рассуждения. Поскольку модуль Diff+(A, P) изоморфен модулю Fk, а модуль
А-1 о Diff+(A, Q) = {А-1 о в: в Е Diff+(A,Q)} —
модулю Gp, получим те же элементы E0,p только с заменой p на k и k на p. Согласно общей теории [10, гл. 9, §1], определены комплексы
0 — E0p Efc0+1;p+1 0, k, p > 0, (13)
которые представляют собой факторизацию комплексов (9). А именно, если элемент £ Е Ekp соответствует оператору а Е Gp П Fk в представлении (11), то d0(£) есть элемент, соответствующий оператору Да Е Gp+1 П Fk+1. Заметим, что замена А на А-1 приводит к тем же комплексам (13).
Докажем, что последовательность pkp, вычисленная для обратимого оператора А по формулам (2), удовлетворяет равенствам (3), (4). При k > p + L имеем Gp С Fp+L С Fk-1, поэтому Gp П Fk-1 = Gp. Из формул (11) получаем E0,p = 0 при k > p + L. Заменив А на А-1, получим E°;p = 0 при p > k + /. Учитывая (12) и (2), заключаем, что
Pfc,p = 0, (14)
если k > p + L или p > k + /. Кроме того, при p > 0 имеем
p+L p+L
X dim E0,p = X (dim(Gp П Ffc + Gp-1) - dim(Gp П Ffc-1 + Gp-1) fc=0 fc=0
= dim(Gp П Fp+L + Gp-1) - dim(Gp П F-1 + Gp-1) =
= dim Gp - dim Gp-1 = dim(Gp/Gp-1) = m,
так как модуль Gp/Gp-1 изморфен Q, а dim Q = m. Отсюда с учетом (12) заключаем, что
p+L
X Kfc,p = m. fc=0
Комбинируя это равенство с (2) и (14), находим
оо
I 0, р> 0;
[ш, р = °.
Таким образом, доказаны первые равенства в (3) и (4). Вторые равенства доказываются аналогично заменой А на А-1.
Неравенства (5) докажем после того, как определим d-схему квадратов для обратимого оператора А. Сейчас докажем неравенства
kp
^ Php > 0, ^ Pk,3 > 0, k,p > 0. (16)
j=o j=o
Гомоморфизм (10) есть ограничение изоморфизма (8) на подмодуль, а значит, есть мономорфизм. Обозначим
_ Gp П Fk + Gp_i Gk,p = ^ . Gp_i
Так как (10) — мономорфизм, то
dim Gk,p < dim Gk+i,p+i. (17)
С другой стороны, Gk_i>p есть подмодуль модуля Gk,p, а фактормодуль Gk>p/Gk_i>p по определению совпадает с Е^О^. Поэтому из формулы (12) следует равенство
dim Gk,p - dim Gk_i,p = dim E°kp = Kkp.
Так как F_i = 0, то E0p = G0,p. Из полученных равенств следует соотношение
s s
dim Gs,p = dim Go,p + J^(dim Gkp - dim Gk_i,p) = Xk,p. (18)
k=i k=0
Объединяя это равенство с (2) и (17), получаем
s s s
^Т Pk,p = ^2 Kk,p Kk_i,p_i = dim Gsp - dim Gs_i,p_i > 0. k=o k=o k=o
Итак, доказано первое равенство (16). Второе равенство доказывается аналогично.
Докажем, что pk,p = 0 вне прямоугольника 0 < k < L, 0 < p < l. Пусть pk,p = 0 для некоторого p > l. Выберем минимальное k0, такое, что pk0,p0 = 0 при p0 > l. Тогда
ko
Y1 Pk,po = Pko,po и из первого неравенства в (16) получаем pko,po > 0. Из соотношений (14)
k=0
следует, что k0 > p0 - l > 0. Поэтому (см. (3))
те
^ Pko,p = 0. (19)
p=0
p0_i p0 С другой стороны, из (16) следует, что ^ Pko,p > 0. Значит, ^ Pko,p > Pko,po > 0. Учи-
p=0 p=0
те
тывая (19), получаем Pko,p < 0. Таким образом, существует такое натуральное pi, что
p=p0+i
k0 k0_i pi > p0 и Pkopi < 0. Еще раз используя (16), получаем Pkpi > 0, а значит, Pkpi > 0.
k=0 ' k=0 Следовательно, существуют такие k < k0 и pi > p0 > l, что Pk,pi = 0. Это противоречит
минимальности k0 • Поэтому pk,p = 0 для всех p > l. Для k > L это равенство доказывается аналогично.
Построим d-схему квадратов для обратимого оператора А. Модули гомологий комплексов (13) с фиксированным p и всеми k > 0 в совокупности составляют член E1 [10, гл. 9, §1] комплекса (9). Однако в отличии от E0 элементы E1, соответствующие комплексам вида (9) с разными номерами p, не совпадают. Модуль гомологий комплекса (13) в члене E°,p обозначим через Zk,p, а в члене E°+1 +1 — через Bk+1,p+1, т.е.
E0
Zfc,p = ker do, Bfc+1,p+1 = coim do = Л+У"!. (20)
do (Efc,pj
Из сравнения размерностей образа и прообраза гомоморфизма (13) следует, что
dim Ek,p — dim ker d0 = dim E°+1p+1 — dim coim d0. Учитывая (12) и (20), получаем равенство Kk p — dim Zk p = Kk+1p+1 — dim Bk+1,p+1, а значит,
pfc,p = dim Bfc,p — dim Zfc-1,p-b (21)
В формуле (21) размерности модулей Bkp и Zk-1p-1 понимаются как целочисленные функции точки t G M. По определению в некоторой окрестности U d-регулярной точки t0 G M левая часть равенства (21) постоянна, так как выражается через постоянные функции (см. (2)). Фиксируем k < L и p < l, для которых dim Zk p = 0 в U. Число s = min dim Zk p определено, так как функция dim Zk p может принимать только конечное количество значений. Пусть dim Zk p = s в точке t1 G U. Тогда образ линейного отображения (13) в точке t1 имеет размерность j = Kk p — s. Поэтому существуют j сечений из E0p, которые линейно независимы в точке t1 . В виду непрерывности этих сечений они линейно независимы и в любой точке некоторой окрестности U1 точки t1, а значит, dim Zk p < Kk p — j = s в U1. Так как s есть минимальное значение dim Zk p в U1, то dim Zk p = s, а значит, модуль Zk p имеет постоянную размерность в U1.
Точка t1 также является d-регулярной точкой, так как лежит в U. Заменив t0 и U на t1 и U1, повторим наши рассуждения для других k < L и p < /. Так как таких чисел k и p конечное число, в результате в окрестности произвольной d-регулярной точки t0 G M мы найдем точку t* и ее окрестность U*, в которой dim Zkp будет постоянно для всех k < L и p < /. Из формулы (21) следует, что в этой окрестности размерности соответствующих модулей Bk,p тоже постоянны. Далее будем рассматривать d-регулярную точку t* G M и докажем первое утверждение теоремы в окрестности U* этой точки. Так как числа (2) не зависят от выбора точки из U, а U* С U, то из первого утверждения теоремы для U* вытекает, что это утверждение верно для всей окрестности U.
Формула (21) объясняет как строить искомую d-схему квадратов: каждый базисный элемент модуля Zk p при 0 < k < L, 0 < p < l определяет квадрат с верхним правым углом
в точке (к + 1, р +1). Более точно, будем рассматривать диагонали вида
р0 р0 р0 Ek_s,p_s, ..., Ek_i,p_i, Ek,p,
где ZkíP _ 0, 0 < к < Ь, 0 < р < I, и соответствующие гомоморфизмы вида (13). Пусть Е0-зр-8 — первый ненулевой элемент на этой диагонали. Рассмотрим цепочку гомоморфизмов
п ¿0 т-,0 ¿0 ¿0 „0 ¿0 т-10 /ООЧ
0 —> Ек-з,р-з —> ... —> Ек-1 ,р-1 —> Ек,р. (22)
Так как Ек0-3-1 р-8-1 _ 0 по выбору в, то Е^ _ Бк-зр>-3. Рассмотрим композиции Л Е0_ л р_л ^ Е0 р гомоморфизмов (22), ] _ 1, ..., в, и их образы в Zk р:
Zk,p 3 Zk,p П d0(Ek_i,p_i) 3 Zk,p П d0(Ek_2,p_2) ... 3 Zk,p П d0(Ek_s,p_s). (23)
Выберем какой-либо базис модуля Zk>p П d0(E0_sp_s) из цепочки (23), дополним его до базиса модуля Zk>p П d0_i(Ek_s+ip_s+i) и так до базиса модуля Zkp. Каждый базисный элемент £ из Zkp П d%0(Ek0_i,p_i) есть образ d0(n), где п G E(0_ip_i. По построению £ G Zkp П П d0+i(E0_i_i p_i_i), поэтому п G d0(Ek_i_i p_i_i), а значит, п определяет некоторый ненулевой элемент Bk_i,p_i. Базисному элементу £ поставим в соответствие квадрат, противоположные углы которого (k - i,p - i) и (k + 1,p +1). Рассмотрев все базисные элементы модулей Zk,p при 0 < k < L, 0 < p < l, получим набор квадратов. Этот набор конечен, так как и количество рассматриваемых значений k, p, и размерности dim Zk,p конечны. В модуле Bk p каждому квадрату с нижним левым углом в точке (k,p) соответствует ненулевой элемент п. По построению набор таких элементов (соответствующих разным квадратам) линейно независим. Поэтому dim Bk p не меньше, чем количество квадратов с нижним левым углом в точке (k,p). Разница между dim Bkp и количеством указанных квадратов есть то натуральное число ak,p, которое добавляется при построении m-таблицы d-схемы квадратов.
Докажем первое неравенство (5). Пусть k > 0, p > 0, s = dim Zkp, а H — подмодуль модуля Gk,p, состоящий из таких элементов а, что Da G Gk,p+i. Ограничение мономорфизма (10) на H есть мономорфизм из H в Gk,p+i, а значит, dim Gkp+i > dim H. Кроме того, H 3 Gk_i,p и Zk,p = H/Gk_ i^ так как Zk,p ядро гомоморфизма (13), а Ek,p = Gk,p /Gk E0+i,p+i = Gk+i,p+i/Gk,p+i. Используя (2) и (18), получаем
k k k ^T Pj,p+i = ^2 Kj,p+i Kj_i,p = dim Gk,p+i - dim Gk_i,p > dim H - dim Gk_i,p = s. j=0 j=0 j=0
Так как по построению d-схемы квадратов s есть количество квадратов с верхним правым углом в точке (k + 1,p + 1), то доказано первое неравенство (5). Второе неравенство доказывается аналогично.
Таким образом, с одной стороны, набор (21) есть набор чисел (2) оператора А, а с другой стороны, он совпадает с m-таблицей построенной d-схемы квадратов (соотношения (3)-(4) для нее мы доказали ранее). Первое утверждение теоремы 1 доказано.
Перейдем к доказательству второго утверждения теоремы. Пусть задана ^-схема квадратов, для которой построена т-таблица чисел рк,р, к, р > 0, удовлетворяющая условиям (3)-(5). Обозначим через гк,р и Ьк,р количество элементов в и Вк,р соответственно. Из определения т-таблицы ^-схемы квадратов следует, что
рк,р = Ьк,р — ¿к,р. (24)
Обозначим через Ь — максимальное значение первой координаты элементов Z и В, а через I — максимальное значение второй координаты этих элементов. Из формулы (24) следует, что рк,р = 0 вне прямоугольника 0 < к < Ь, 0 < р < /.
Докажем существование отображения Пусть р > 0 и Z0,p = 0 для некоторого к > 0.
Выберем минимальное к, такое, что Z0,p = 0. Из формулы (24) следует, что р^,р = Ь.,-,р при
к-1
< к, а значит, ^ р^-,р — количество элементов в множестве В = В0р и В0р и ... и В0-1 р.
¿=о ' , , ,
Поэтому из первого неравенства (5) получаем, что в множестве В больше, чем гк,р элементов. Пусть ^ отображает элементы Zk р в какие-то (любые, но различные) 5 элементов множества В. Рассмотрим множество Z0+1p. Сумма
к к Ху р^',р ^ у Ь.?,р — ^к,р
¿=0 .7=0
равна количеству элементов множества В и Вк,р, которые не лежат в образе ). Из
первого неравенства (5) следует, что таких элементов больше, чем элементов в Zk+1p.
Аналогично предыдущему определим инъективное отображение из Z0+1p в множество
£
(В и В° ) \ ). Увеличивая далее к, получим инъективное отображение ^ из и Z0,p
т , , к=1 , £
в и В°,р. Далее, из первого равенства (3) и формулы (24) следует, что в указанных мно-
к=0 ,
жествах одинаковое количество элементов, а значит, инъекция ^ есть биекция. Перебирая р = 1, ..., /, получаем определение причем ^ взаимно однозначно отображает элементы Z на элементы В с ненулевой второй координатой.
Аналогично определяется взаимно однозначное отображение ^ из Z в подмножество элементов В с ненулевой первой координатой, причем ^ не меняет первую координату и уменьшает вторую.
Покажем, что индуктивное определение операторов Уь (см. разд. 3) корректно. Заметим, что при р > 0 и Ь € В°,р определен элемент (р-1 (Ь) из Z, так как биективно отображает Z на элементы В с ненулевой второй координатой. Отображения ф и р^ уменьшают вторую координату. Поэтому вторые координаты элементов (ф о р-1)(Ь) и (р^ о р-1)(Ь) из В меньше р, а значит, применимо предположение индукции и формула (6), которая определяет Уь.
Используя индукцию по р, покажем, что при Ь € В°,р порядок оператора Уь равен р. Для р = 0 это доказано при определении Уь. Если р > 0 и 5 — длина стороны квадрата, соответствующего элементу р-1(Ь) € Z, то скалярный оператор по определению имеет
порядок а вторая координата элемента (ф о р- )(Ь) есть р — Поэтому первое слагаемое
правой части формулы (6) имеет порядокр, поскольку при композиции скалярного оператора с матричным оператором порядки операторов складываются. Порядок второго слагаемого правой части формулы (6) меньше р, так как вторая координата элемента (р^ о р-1)(Ь) меньше р. Следовательно, порядок оператора (6) равен р, т.е. V € Ср \ Ср-1.
Докажем обратимость оператора А, определенного формулой (7). Обратный оператор А-1 строится аналогично построению А с заменой Q на Р, Р на р^, на , на р^, и с использованием выбранных базисов: е1, .. ., ет в модуле Q и /1, .. ., /т в модуле Р. А именно, каждому элементу Ь € В поставим в соответствие дифференциальный оператор V из А в Р. При этом используем индукцию по к (первая координата элемента Ь). Для к = 0 элементу Ь1, г = 1, ..., т, соответствует оператор УЬ1, отображающий а в а/. Для к > 0 верна формула
^ = ^о^Рхь) о — (25)
Оператор А-1 определяется формулой
А-1(«1в1 + ... + атет) = УЬ1 (а1) + ... + УЬт (ат), аь...,ат € А. (26)
Покажем, что операторы А и А-1 взаимно обратные. Для этого докажем сначала формулу связи двух типов построенных операторов:
V = А о Ь € В. (27)
Используем индукцию по первой координате к элемента Ь. Из определения УЬ1 и формулы (7)
г
имеем
(А о УЬх )(а) = А(а/г) = УЬх (а), а € А, г = 1, ...,т.
Таким образом, доказана формула (27) для к = 0. В случае к > 0 из формулы (25) следует, что
А о ^ = А о ^о^-1)(Ь) о ^-1(Ь) — А о ^о^-1)(Ь).
Первая координата элементов (ф о р-1)(Ь) и (р^ о р-1 )(Ь) меньше к. По предположению индукции формула (27) для этих элементов справедлива. Следовательно,
А о = о П^-1(Ь) — ^(^^о^—1)(Ь).
Обозначим Ь = (р^ о р-1)(Ь). Тогда р-1(Ь) = р-1(Ь) и (ф о р-1)(Ь) = (ф о р-1)(Ь), а значит, используя (6) для элемента Ь, получаем
А о ^ = ^о^Х*) о п^-1 (Ь) — V;, = ^ о^—1(Ь) = V*.
Формула (27) доказана. Комбинируя ее с определением А-1 (см. (26)) и определением операторов , ..., , получаем
(А о А-1)(а1в1 + ... + атет) = (А о ^)Ы + ... + (А о )(ат) =
= ^ (а1) + ... + (ат) = а^1 + ... + атет.
Таким образом, композиция А о А-1 есть тождественный оператор. Тождественность композиции А-1 о А доказывается аналогично.
Докажем, что набор чисел (2) оператора А есть т-таблица заданной ^-схемы квадратов. Доказательство основано на последовательном построении базисов модулей Ср/Ср-1, Ск,р, Ек,р, Бк,р, Zk,p, к,р > 0, а затем сравнении формулы (21) для чисел рк,р оператора А с формулой (24) для т-таблицы заданной ^-схемы квадратов. Отметим сначала, что если (к,р) — координаты элемента Ь е Б, то по построению УЬ € В1Г:£+(А, Q) _ Ор, а УЬ е Б1й+(А, Р). Используя равенство (27), получаем УЬ е Рк П Ср, а значит, [УЬ]р е Скр.
Используя индукцию по р, докажем, что базисом Ор/Ор-1 при р > 0 является набор таких элементов
Ор-р У]и, ьр е Б<0.,р., рг < р, г _ 1, ..., т, (28)
что для любого г _ 1, ..., т не существует элемента г из Z0o , для которого ф(г) _ Ьр и р0 < р. Отметим, что если элемент из Бк0 ,р попадает в наборы (28) для нескольких р, то он имеет несколько обозначений. Верхний индекс в обозначении Ьр указывает на номер фактормодуля.
При р _ 0 имеем С0/С-1 ~ С0, так как С-1 _ 0. Кроме того, из определения С0 _ _ Б1Й:+(А, Q) следует, что 5а(а) _ 0 для любых а е С0 и а е А, а значит,
£«(а)(1) _ (а+а)(1) — (аа)(1) _ а(а) — аа(1) _ 0.
Из равенства а (а) _ аа(1) следует, что элемент а однозначно определяется значением а(1) е Q. Раскладывая элемент а(1) по выбранному базису е1, . .., ет модуля Q, получаем:
а(а) _ аа(1) _ а(а^1 + ... + атвт) _ (а+У^ + ... + ат^Ьт)(а),
где Ь1, . .., Ьт — элементы Б с нулевой второй координатой; УЬ1, .. ., УЬт — операторы из С0, определенные выше. Таким образом, а _ а+УЬ1 + ... + атУЬт, а значит, [УЬ1 ]0, .. ., [Уьт]0 — базис модуля 6^/^-1.
Пусть (28) — базис рассматриваемого типа модуля Ср/Ср-1. Так как отображение (8) является изоморфизмом модулей, то набор элементов
Вр-р1+1[^К ]р1, ..., £р-рт+1[Уьт ]рт (29)
есть базис модуля Ср+1 /Ср. Предположим, существует такой элемент г из Z0o , р0 < р +1, что ф(г) _ Ьро для некоторого г0 от 1 до т. Тогда по предположению индукции р0 > р, т.е. р0 _ р +1. Рассмотрим квадрат ^-схемы, соответствующий элементу г. Вторая координата верхнего правого угла есть вторая координата г, т.е. р + 1, а вторая координата его нижнего левого угла есть вторая координата Ь1ро, т.е. рг0. Поэтому сторона квадрата равна р + 1 — р^0, а значит, порядок оператора естьр+1 — 'р%0). Формулу (6) для элемента Ь _ перепишем в виде
_ о П2 — У^(г). (30)
По определению р^, элемент Ь = р^(г) имеет ту же вторую координату, что и г, т.е. р + 1. Поэтому У^(,г) € Ср+1. По определению р^ вторая координата элемента р^ (г) меньшер +1, а значит, У^(¿) € Ср и [У^(г)]р+1 = 0. Пусть а0 — коэффициент при старшей производной оператора . Тогда а0 = 0 и из (30) следует, что
УЫ*)]Р+1 = [Уьро ◦ пг]р+1 - [У^(г)]Р+1 = аоДр+1-Рг0 Уо]Рго. (31)
го го п
В наборе (29) заменим элемент Др-ко+1[Уьр ]Ко на элемент [У^^(г)]Р+1, т.е. положим ЬРп+1 = = Р^(г).
Если существует еще такой номер г1, что ЬР1 = ^(г) для некоторого г из , р0 < р +1, то с элементом набора (29) с номером г1 поступим аналогично. Для остальных номеров г = 1, ..., т положим ЬР+1 = ЬР. Получим базис модуля СР+1/СР вида (28). Докажем, что при к > 0, р > 0 модуль Ск,р имеет базис
ДР-Р*[У^, Ь € В^, р* < р, к, < к -р + р*, (32)
где рассматриваются лишь те Ь, для который не существует такого элемента г из ^0о , что = Ь и р0 < р. Из определения Ск,р и теоремы Нетер для X = Ср П Рк и У = СР-1 следует, что
_ Ср П Рк + СР-1 _ Ср П Рк СА:,р = _
Gp_i GP-i n Fk
Так как Fk С Gk+Z, то Gk+Z П Fk = Fk. Используя эти факты и равенство G-1 = 0, получаем
k+Z k+Z
X dim Gk,p = X (dim(Gp n Fk) - dim(Gp_i n Fk p=0 p=0
= dim(Gk+Z П Fk) - dim(G_1 n Fk) = dim Fk = dimDiff+(A, P) = m(k + 1).
Обозначим через gk p количество элементов в наборе (32). Этот набор есть часть набора (28), а следовательно, он линейно независим в каждой точке t е M .С другой стороны, элементы (32) лежат в Gk,p. Поэтому gk p < dim Gk,p и
k+Z k+Z
Xgk,p dimGk,p = m(k + 1). (33)
p=0 p=0
Покажем, что первая сумма в соотношении (33) также равна m(k +1). Для этого заметим, что наборы (32) для разныхp не пересекаются, а указанная сумма равна количеству элементов в наборах (32) для p = 0, 1, ..., k + /. Рассмотрим элементы
Ds[Vb]pt, b е , k* < k, s = 0, 1, ...,k - k*, p = 0, 1, ...,k + l. (34)
Элемент Ds[Vb]pt не входит в набор (32), только если существует элемент z из , для которого -0(z) = b и p0 < p = p* + s. Из равенства -0(z) = b и определения ф следует, что p* — k* = p0 — k0, и поэтому s > k0 — k*. Кроме того, из неравенства s < k — k*
получаем условие k0 = p0 — p* + k* < s + k* < k. Таким образом, каждому элементу z G Z°o , k0 < k, p0 < p, соответствуют k — k0 + 1 элементов набора (34), которые при s = k0 — k*, k0 — k* + 1, ..., k — k* не входят в наборы (32). Всего количество элементов в наборах (34) при p > 0 равно
k+l k
— k* + i)bkt,pt,
Pt=0 k.=0
количество элементов в (34), которые не входят в наборы (32), равно
k+l k
— k0 + 1)zk0,P0.
Po=0 ko=0
Учитывая это и меняя во второй сумме индексы k0,p0 на k*,p*, получаем
k+l k+l k k+l k
J29k,P — k* + 1)bk*,P* — — k* + 1) Zk»,P» =
p=0 Pt =0 k* =0 Pt=0 k*=0
k+l k
= — k* + 1)(bk»,Pt — Zkt,Pt).
Pt=0 kt=0
Используя формулу (24) для элементов m-таблицы заданной d-схемы квадратов, а также равенства (3)-(4), получаем
k+l k+l k k k+l
Y^ 9k,P = YY(k — k* + 1)Pkt,Pt = Y(k — k* +1) pk*P* = m(k+1).
P=0 Pt=0 kt=0 kt=0 Pt=0
Следовательно, соотношение (33) возможно только в случае равенства. Поэтому для всех k > 0, p > 0, имеем gk,P = dim Gk,P, а значит, модуль Gk,P имеет базис (32).
Сравнивая базис (32) модуля Gk,P и аналогичный базис модуля Gk-i,P, получаем, что элементы набора (32), для которых
k < k 1 - p + p*, образуют базис его подмодуля Gk-1,P С Gk,P. Остальные элементы в (32) определяют базис фактормодуля E°P = = Gk,P/Gk-l,P, т.е. базис модуля E°,P определяют элементы
DP-Pt [V6]Pt, b G Б".*., p* < p, k* = k — p + p*, (35)
где для каждого b не существует элемента z из Z°o , для которого ^>(z) = b и p0 < p.
Базисный элемент модуля E0P, соответствующий элементу DP-Pt [Vb]Pt, лежит в образе гомоморфизма d0, если p > p*, и определяет элемент Bk,P, если p = p*. Из определения отображения ф следует, что если b G Bk P, z G Z°o,Po и ф^) = b, то p0 > p. Поэтому для любого элемента b G B°,P не существует такого элемента z из Z°o , что ф(z) = b и p0 < p. Таким образом, для всех k > 0, p > 0 набор элементов [Vb]P, b G B0P, определяет базис модуля Bk,P, а значит, dim Bk,P = bk,P.
Наконец, гомоморфизм (13) отображает базисный элемент модуля Ekp, соответствующий элементу Dp_p* [Vb]pt, в базисный элемент модуля E0+1 p+1, соответствующий элементу Dp-p*+1[Vb]pt, если не существует элемента z е Z°+1 +1, для которого ф^) = b. Если такой элемент z существует, то Dp+1-p* [Vb]p<, = 1 [V^h(z)]p+1 (см. (31)). Из определения следует, что первая координата ph(z) е B меньше k + 1, а значит, элемент [V^h(z)]p+1 соответствует нулевому элементу в Ej°+1 p+1. Таким образом, гомоморфизм (13) отображает часть базиса (35) модуля p в часть базиса модуля Ej°+1 p+1, а оставшуюся часть базиса (35) в нуль. Следовательно, размерность ядра этого гомоморфизма равна количеству элементов во второй части этого базиса. Между элементами этой части базиса (35) и элементами множества Zk0+1 +1 установлено взаимно однозначное соответствие. Поэтому dim Zk p = zk+1 ,p+1.
Из доказанных равенств следует, что правые части формул (21) и (24) совпадают, а значит, набор чисел (2) оператора А есть m-таблица данной d-схемы квадратов. Теорема 1 полностью доказана.
Заключение
Как следует из доказанной теоремы, каждый обратимый линейный дифференциальный оператор однозначно определяет ^-схему квадратов. Из доказательства этой же теоремы вытекает, что такой оператор однозначно определяется ^-схемой квадратов, выбором базисов модулей ^ и Р, выбором скалярных дифференциальны« операторов П, г € Z, и выбором отображений р^, р^ из Z в В, причем р^, не меняет вторую координату элемента г € Z и уменьшает первую, а р^, наоборот, не меняет первую координату и уменьшает вторую (процедура конструирования обратимого оператора приведена в разд. 3).
Отображения р^ и р^ имеют простую интерпретацию на языке спектральный последовательностей (см. [10, гл. 9, §1]). Например, если р^(г) = Ь, ф(г) = Ь0, (к,р) — координаты г, (к — з,р) — координаты Ь, а (к — /,р — /) — координаты Ь0, то оператор Уь определяет элемент Е|-5р, оператор Д1-1УЬо —элемент Е|-1р-1, а оператор Е|-1р-1 ^ отображает
второй из указанных элементов в первый.
Модуль всех линейных дифференциальных операторов из А в Q имеет систему вложенных подмодулей
... с Ср = ЫЙ+(А, Q) С Ср+1 = ЫЙ+^А, Q) С ... (36)
Такой модуль называют фильтрованным, а семейство его подмодулей (36) — фильтрацией. Определяющим свойством фильтрации (36) является точность последовательностей (9). Теорема 1 и ее доказательство основаны на исследовании, как преобразуются фильтрации (36) под действием обратимого оператора. При этом свойство точности последовательностей (9) играет ключевую роль. Модули линейных дифференциальных операторов в частных производных также являются фильтрованными, а их фильтрации определяются свойством точности последовательностей, которые называются в [2, п. 1.2.16] последовательностями Спен-
сера для символов. Но в отличие от последовательностей (9) они имеют больше двух ненулевых членов, и поэтому модули Ek0,P для них нельзя расположить на плоскости, а следует рассматривать более сложную трехмерную конструкцию. Тем не менее, если удастся найти аналог d-схемы квадратов для этой трехмерной конструкции, можно ожидать обобщения теоремы 1 на случай операторов в частных производных.
По-видимому, дифференциальные операторы с запаздыванием и разностные операторы также допускают аналогичные конструкции, и результаты этой работы могут быть обобщены и на такие обратимые операторы.
Работа выполнена при финансовой поддержке гранта Президента РФ поддержки ведущих научных школ НШ-53.2014.1 и грантов РФФИ 13-07-00736 и 14-01-00424.
Список литературы
1. Громов М. Дифференциальные соотношения с частными производными. М.: Мир, 1990. 536 с.
2. Виноградов А.М., Красильщик И.С., Лычагин В.В. Введение в геометрию нелинейных дифференциальных уравнений. М.: Наука, 1986. 336 с.
3. Четвериков В.Н. Метод линеаризации для решения задач плоскостности и поиска оператора совместности // Дифференциальные уравнения. 2006. Т. 42, № 10. С. 1518-1527.
4. Fliess M., Levine J., Martin Ph., Rouchon P. A Lie — Backlund approach to equivalence and flatness of nonlinear systems // IEEE Trans. Automatic Control. 1999. Vol. 44, no. 5. P. 922-937. DOI: 10.1109/9.763209
5. Martin Ph., Murray R., Rouchon P. Flat systems // Proc. of the 4th European Control Conf. Plenary lectures and Mini-courses. Brussels, 1997. P. 211-264.
6. Pomet J.-B. A differential geometric setting for dynamic equivalence and dynamic linearization // Geometry in Nonlinear Control and Differential Inclusions / B. Jakubczyk, W. Respondek, T. Rzezuchowski, Eds. Warsaw: Banach Center Publications, 1995. P. 319-339.
7. Chetverikov V.N. Invertible linear differential operators on two-dimensional manifolds Vienna, 1993. 16 p. (Preprint of the Erwin Schrodinger International Institute for Mathematical Physics; no. 55).
8. Четвериков В.Н. Управляемость плоских систем // Дифференциальные уравнения. 2007. Т. 43, № 11. С. 1518-1527.
9. Хьюзмоллер Д. Расслоенные пространства. М.: Мир, 1970. 442 с. 10. Спеньер Э. Алгебраическая топология. М.: Мир, 1971. 680 с.
SCIENTIFIC PERIODICAL OF THE BAUMAN MSTU
SCIENCE and EDUCATION
EL № FS77 - 48211. №0421200025. ISSN 1994-0408
electronic scientific and technical journal
Classification and Construction
of Invertible Linear Differential Operators
on a One-Dimensional Manifold
# 07, July 2014
DOI: 10.7463/0714.0718107
Chetverikov V. N.
1Bauman Moscow State Technical University, Moscow, Russian Federation
achetverikov. vl@yandex.ru
Keywords: invertible linear differential operators, transformations of control systems, spectral sequences of chain complexes
Invertible linear differential operators with one independent variable are investigated. The problem of description of such operators is important, because it is connected with transformations and the classification of control systems, in particular, with the flatness problem.
Each invertible linear differential operator represents a square matrix of scalar differential operators. Its product with an operator-column is an operator-column whose order does not exceed the sum of orders of initial operators. The operators-columns, the product with which leads to order fall, i.e. the order of the product is less than sum of orders of factors, are interesting for the description of invertible operators. In this paper the classification of invertible operators is based on dimensions dk,p of intersections of modules Gp and Fk for various k and p, where Gp is the module of all operators-columns of order not above p, and Fk is the module of compositions of the invertible operator with all operators-columns of order not above k. The invertible operators that have identical sets of numbers dk,p form one class.
In the paper the general properties of tables of numbers dk,p for invertible operators are investigated. A correspondence between invertible operators and elementary-geometrical models which in the paper are named by d-schemes of squares is constructed. The invertible operator is ambiguously defined by its d-scheme of squares. The mathematical structure that must be set for its unique definition and an algorithm for the construction of the invertible operator are offered.
In the proof of the main result, methods of the theory of chain complexes and their spectral sequences are used. In the paper all necessary concepts of this theory are formulated and the corresponding facts are proved.
Results of the paper can be used for solving problems in which invertible linear differential operators are arisen. Namely, it is necessary to formulate the conditions on the invertible operator
which follow from the problem statement, in language of d-schemes of squares. Further the d-scheme, satisfying to these conditions, is searched and the required operator is restored on this scheme.
A discussion of possible generalizations of offered methods on the cases of partial-differential operators, delay-differential operators and difference operators concludes the paper.
References
1. Gromov M. Partial differential relations. Springer Berlin Heidelberg, 1986. 363 p. (Ser. Ergebnisse derMathematik undihrer Grenzgebiete; vol. 9). DOI: 10.1007/978-3-662-02267-2 (Russ. ed.: Gromov M. Differentsial'nye sootnosheniia s chastnymi proizvodnymi. Moscow, MirPubl., 1990. 536 p.).
2. Krasil'shchik S., Lychagin V.V., Vinogradov A.M. Geometry ofJet Spaces and Nonlinear Partial Differential Equations. Gordon and Breach, New York, 1986.441 p. (Russ. ed.: Vinogradov A.M., Krasil'shchik I.S., Lychagin V.V. Vvedenie v geometriiu nelineinykh differentsial'nykh uravnenii. Moscow, Nauka Publ., 1986. 336 p.).
3. Chetverikov V.N. Linearization method for the flatness problem and for finding the compatibility operator. Differentsial'nye uravneniia, 2006, vol.42, no. 10, pp. 1518-1527. (English translation: Differential Equations, 2006, vol.42, no. 10, pp. 1479-1489. DOI: 10.1134/S0012266106100120).
4. Fliess M., Levine J., Martin Ph., Rouchon P. A Lie — Backlund approach to equivalence and flatness of nonlinear systems. IEEE Trans. on Automatic Control, 1999, vol. 44, no. 5, pp. 922-937. DOI: 10.1109/9.763209
5. Martin Ph., Murray R., Rouchon P. Flat systems. Proc. of the 4th European Control Conf. Plenary lectures and Mini-courses. Brussels, 1997, pp. 211-264.
6. Pomet J.-B. A differential geometric setting for dynamic equivalence and dynamic linearization. In: JakubczykB., Respondek W., Rzezuchowski T., eds. Geometry in Nonlinear Control and Differential Inclusions. Warsaw, Banach Center Publications, 1995, pp. 319-339.
7. Chetverikov V.N. Invertible linear differential operators on two-dimensional manifolds. Preprint no. 55. Vienna, Erwin Schrodinger International Institute for Mathematical Physics, 1993. 16 p.
8. Chetverikov V.N. Controllability of flat systems. Differentsial'nye uravneniya, 2007, vol. 43, no. 11, pp. 1518-1527. (English translation: Differential Equations, 2007, vol.43, no. 11, pp. 1558-1568. DOI: 10.1134/S0012266107110110).
9. Husemoller D. Fibre Bundles. McGraw Hill, 1966. 376 p. (Russ. ed.: Husemoller D. Rass-loennyeprostranstva. Moscow, MirPubl., 1970. 442 p.).
10. SpanierE.H. Algebraic Topology. McGraw-Hill. 1966. 528 p. (Russ. ed.: SpanierE.H. Alge-braicheskaia topologiia. Moscow, MirPubl., 1971. 680 p.).