О геометрии решения приближенных уравнений и их симметриях

ISSN 2074-1863 Уфимский математический журнал. Том 9. № 2 (2017). С. 40-55.

УДК 517.911

О ГЕОМЕТРИИ РЕШЕНИЯ ПРИБЛИЖЕННЫХ УРАВНЕНИЙ И ИХ СИММЕТРИЯХ

В.В. ГОРБАЦЕВИЧ

Аннотация. Статья посвящена разработке геометрического подхода к теории приближенных уравнений (в том числе ОДУ и УРЧП) и их симметрий. Вводятся дуальные алгебры Ли, многообразия над дуальными числами и дуальные группы Ли. Описываются некоторые конструкции, применяемые к этим объектам. На основе этих построений показывается, как можно сформулировать основные понятия и методы в теории приближенных уравнений и их симметрий. Доказательства многих общих утверждений при этом получают почти автоматически из классических рассуждений, в отличие от использовавшихся до сих пор методов изучения приближенных уравнений.

Ключевые слова: приближенное уравнение, приближенная алгебра Ли, дуальные числа, дуальная алгебра Ли, многообразие над дуальными числами

Mathematics Subject Classification: 53С15, 34С30, 37J15, 51Н30

Введение

В статье рассматриваются понятия, возникающие при изучении приближенных симметрий различных функциональных и дифференциальных объектов, зависящих от малого параметра. Это направление было инициировано статьей [1] и затем продолжено во многих других работах. Основная цель статьи — не разбор новых примеров решения уравнений или вычисления симметрий новых уравнений, а построение геометрического подхода к изучению этих уравнений и симметрий и развитие соответствующей методологии. Основой здесь для автора был принцип, который пропагандировал А. Гротендик: «Для того чтобы решить задачу, нужно вложить ее в такую математическую среду, в которой ее решение станет очевидным». Статья посвящена как раз описанию такой среды, в которой приближенные уравнения и их симметрии формулируются и решаются наиболее естественным образом. В частности, будет доказано весьма общее утверждение — принцип соответствия, который можно сформулировать так:

Пусть F(х,и,е) = 0 — некоторое уравнение (возможно, дифференциальное — обыкновенное или с частными производными), где и — гладкая функция (скалярная или векторная) независимой переменной х (скалярной или векторной), а е — некоторый параметр (обычно предполагаемый близким к 0). Рассмотрим линеаризацию F(х,и) ~ fo(x,u) + ef\(x,u) левой части этого уравнения и продолжим ее па область дуальных (т.е. принадлежащих области дуальных чисел а + eb, где a,b g R, е2 = 0) значений аргумента. Получим уравнение F(х,й) = 0. Пусть и(х) — его точное решение. Тогда при ограничении на вещественные значения аргумента х = х мы получим приближенное решение исходного уравнения с точностью до о(е).

Совсем общая формулировка основного методологического результата статьи совершенно элементарна: две функции Ф(*,е) и Ф(*, е) равны с точностью до о(е) тогда и только тогда, когда равны соответствующие им функции дуального аргумента Ф(*) и Ф(*),

V.V. Gorbatsevich, On geometry of solutions to approximate equations and their

symmetries.

© горбацевич в.в. 2017.

Поступила Ц июля 2016 г.

В применении к решениям дифференциального уравнения мы получает приведенную выше формулировку принципа соответствия.

Используемые в этих формулировках понятия будут подробно описаны ниже. Общий принцип соответствия применим, в частности, к решению дифференциальных уравнений, зависящих от параметра, а также к вычислению симметрий этих уравнений,

В §1 будут сформулированы исходные понятия, связанные с методами решения приближенных уравнений и вычислением их симметрий. Рассматриваются варианты определения приближенных алгебр Ли - одного из важных понятий дальнейшего изложения, В §2 рассматриваются понятия симметрии и приближенной симметрий уравнений, В §3 дается определение приближенной алгебры Ли, которое связывается с понятием алгебр Ли над алгеброй дуальных чисел И2.

В §4 вводится основное для данной статьи понятие ^-многообразия, дающее геометрическую основу для изучения приближенных уравнений и их симметрий. Вводятся также обычные сопутствующие понятию многообразия объекты - касательные векторы и поля и др.

В §5 доказаны некоторые результаты общей теории дифференциальных уравнений с дуальным аргументом.

Параграф §6 посвящен краткому обсуждению теории Ли в Д2-еитуации, И, наконец, в §7 рассмотрено применение всего вышеизложенного к вопросу о методах решения и анализа приближенных уравнений. Дана геометрическая интерпретация одного широко распространенного метода решения приближенных уравнений и указано обобщение этого метода, В качестве иллюстрации рассмотрен несложный пример решения задачи Коши для приближенного ОДУ второго порядка, иллюстрирующий полезность перехода в Д2-облаеть,

1. Мотивировки и исходные понятия

Как известно, Софус Ли начал свои исследования по теории непрерывных групп преобразований, пытаясь перенести на дифференциальные уравнения методы теории Галуа, При этом ему понадобилось ввести понятие группы симметрий дифференциального уравнения и некоторых других функциональных или дифференциальных объектов. Затем было введено понятие, говоря современным языком, группы Ли преобразований и соответствующей алгебры Ли векторных полей, И только потом в математику были введены общие понятия группы Ли (не обязательно группы преобразований) и абстрактной алгебры Ли,

Первоначальное понятие группы симметрий (фактически речь обычно шла об однопа-раметричеекой группе преобразований), сохраняющей функцию или дифференциальное уравнение, было довольно подробно изучено и многократно были рассмотрены различные его применения (в теории дифференциальных уравнений, в геометрии, теоретической физике и др.). Но при практических применениях стали возникать определенные трудности при использовании понятия симметрии. Большинство дифференциальных уравнений в физике, механике и других областях применений математики являются приближенными, в них входят некоторые параметры (обычно малые), Равенство нулю некоторого такого параметра соответствует определенной «классической» ситуации. Но реальные физические и другие объекты описываются более точно при некоторых ненулевых, хотя и малых, значениях этих параметров. Типичные примеры такого рода параметра — параметр к (постоянная Планка) или же величина 1/с, где с — скорость света. При с ^ ж мы имеем 1 /с ^ 0 и переход от релятивистской теории к классической физике. Квантовая механика при к ^ 0 переходит в классическую. Такого рода малые параметры вызвали в математике понятие квантования, под которым, в самой общей ситуации, понимается некоторая деформация,

зависящая от параметра. При таких «квантованиях» группы симметрии меняются, причем некоторые из симметрии «разрушаются». Это указывало на необходимость изучать такие симметрии, которые не исчезают при ненулевых значениях параметров, а просто как-то «деформируются». Так появились методы изучения «приближенной» симметрии. Имеется несколько подходов к понятию приближенной симметрии дифференциального уравнения, В данной работе мы основываемся на том из них, который был введен в [1] и потом был довольно подробно исследован рядом автором (см, например обзоры [2], [3]), Однако наш подход будет несколько отличен от принятого этими авторами — он будет более геометричен и основан на понятии многообразия над алгеброй дуальных чисел,

В дальнейшем будет использоваться стандартное обозначение $ ~ У для эквивс1-лентноети двух функций по переменной е: запись f (х,е) ~ д(х,е) обозначает условие $(х, е) — д(х, е) = о(е) при е ^ 0, Можно рассматривать и о(ер) при произвольном натуральном р, но здесь мы этого обычно использовать не будем (указывая иногда на возможные в этом случае обобщения).

Начнем с предварительных соображений и определений (введенных и первоначально изученных в [1]), Пусть имеется уравнение Р(х,е) = 0, где Р — некоторая функция (гладкая или даже аналитическая), зависящая от переменной х (возможно, векторной — ж = (х\,х2 ...)) и малого параметра е. Значению е = 0 будет соответствовать уравнение ¡0(х) = 0 (где /о(х) = Р(х, 0)), которое мы будем называть точным, а исходное — приближенным. Аналогично можно рассматривать и зависящие от скалярного параметра системы уравнений (когда Р — это вектор-значное отображение), а также дифференциальные уравнения, например Р (х, и,и',... ,е) = 0 ил и Р (х, у, и, их ,иу,... ,е) = 0,

Теперь рассмотрим понятие однопараметричеекой приближенной группы преобразований и ее инфинитезимальный аналог — приближенное векторное поле (или, что эквивалентно, приближенный дифференциальный оператор, соответствующий этой группе). Рассмотрения будут локальными, поэтому можно предполагать, что все объекты определены в некоторой окрестности начала координат евклидова пространства И", При этом исходить мы будем из классических определений теории групп преобразований,

Однопараметрическая (локальная) группа преобразований — это такое преобразование х' = ф(х, а), зависящее от а € И (возможно, что а пробегает только некоторый симметричный интервал в локальной ситуации), для которого при а,Ь € И выполняются два тождества:

ф

ф(х, 0) = х.

Отметим, что обратное преобразование здесь всегда существует — оно соответствует значению параметра — а. Векторное поле X, соответствующее этой однопараметричеекой группе преобразований (его еще иногда называют инфинитезимальным преобразованием или генератором), определяется так: оно имеет вид X = X(х)где

Л

X (х) = — (ф(х,а))а=о. аа

Инвариантность функции Р(х) относительно однопараметричеекой группы преобразований определяется естественным образом:

ф

для любого допустимого а € И,

Это условие эквивалентно, как известно еще со времен С, Ли, обращению в нуль результата применения соответствующего дифференциального оператора к этой функции — т.е. условию X • Р = 0,

Рассмотрим теперь ситуацию, когда имеется малый параметр е. Здесь повторяются те же конструкции, которые были приведены выше, только появляется дополнительная переменная — параметр е — и равенства заменяются на «приближенные равенства», В частности, приближенная однопараметрическая (локальная) группа преобразований — это такое преобразование х' & ф(х, а), где а € И (или же а пробегает некоторый интервал в локальной ситуации), для которого при а,Ь € И выполняются тождества

ф(ф(х, а),Ь) & ф(х, а + Ь) ф(х, 0) & х.

В этих соотношениях подразумевается наличие малого параметра е.

Условие приближенной инвариантности функции Р (когда ф(Р(х),а) & Р(х) для любого а € И принимает вид X • Р & 0, Здееь X = X(х,е) — векторное поле (гладкое или аналитическое, в зависимости от рассматриваемой нами ситуации), зависящее от параметра е. При этом X(х,е) можно рассматривать и как однопараметричеекое семейство (зависящее от параметра е) «точных» векторных полей.

Для дифференциальных уравнений понятия симметрии и приближенной симметрии определяются естественным образом, при этом используется естественное продолжение преобразования координат на производные.

Когда мы переходим от однопараметричееких групп преобразований дифференциальных уравнений к полным группам спмметрнй (являющихся группами Ли, возможно, бесконечномерными), соответствующие инфинитезимальные преобразования образуют алгебры Ли векторных полей. При наличии параметра е получаемые семейства векторных полей (соответствующих приближенным снмметрпям) не образуют, вообще говоря, алгебру Ли, так как аксиомы алгебры Ли (в частности, условие Якоби) выполняются только приближенно. Соответствующие инфинитезимальные объекты естественно называть приближенными алгебрами Ли, о них будет рассказано подробнее ниже.

Рассмотрение функций и векторных полей с точностью до о(е) эквивалентно их линеаризации по параметру е. В частности, функция Р(х, е) заменяется на функцию /0(х) + е/1(х), где ¿о, — некоторые гладкие или аналитические функции. Для приближенного векторного поля X(х,е) имеем

X(х, е) & Х0(х) + еХ1(х),

где Х0, Х1 — обычные векторные поля. Коммутирование подобных приближенных векторных полей происходит так:

[X, У] = [*о + , уо + 6^1] = [л0, Уо] + е([л0, + [Х1, ЭД).

При этом билинейность коммутирования и кососимметричность сохраняются, тождество Якоби тоже выполняется. Тем самым мы получаем некоторую алгебру Ли, Но если рассматривать обычный коммутатор линеаризованных векторных полей, то его результат уже не будет, вообще говоря, линейным по е. Поэтому указанное выше коммутирование уже не является операцией коммутирования векторных полей, оно производится немного по-другому. Мы приходим к новому понятию, в котором в единый объект объединяются все алгебры Ли, зависящие от параметра е. Этот объект X(е) называется (следуя [4]) приближенной алгеброй Ли, и сейчас мы рассмотрим несколько различных подходов к рассмотрению таких объектов.

Вот первый подход. Пусть Ф — некоторая алгебра Ли векторных полей (или их ростков) в некоторой окрестности точки х € И", Рассмотрим множество векторных полей вида X + бУ, где Х,У € Ф, а е — некоторый параметр. Это множество обозначим Ф(б) или представим его в виде Ф + бФ, Вообще говоря, это множество (ростков) векторных полей не образует алгебру Ли, так как операция коммутирования может приводить к векторным полям, не лежащим в Ф(б) (ибо может появится елагаемое с б2), Но если мы определим на Ф(е) операцию коммутирования по-другому (как это сделано выше для приближенных векторных полей), то получим, как нетрудно убедиться, алгебру Ли (но это уже не будет алгебра Ли векторных полей). Такие алгебры Ли называются в указанных выше статьях приближенными.

Кроме уже построенных таким способом приближенных алгебр Ли, полезно рассматривать и некоторые их подалгебры — но не все, а только те, которые выдерживают умножение на е. Именно такие подалгебры Ли соответствуют группам Ли приближенных еиммет-рий дифференциальных уравнений, изучение которых и инициировало введение понятия приближенной алгебры Ли,

Отметим также, что рассмотрение выражений сходного вида а + ша°, где а, а0 — трехмерные векторы, а ш — символ, для которого ш2 = 0, используется в механике как объект винтового исчисления (см, например [5]),

Рассмотрим теперь и другую конструкцию. Пусть Ф — произвольная вещественная алгебра Ли, Рассмотрим двумерную алгебру И2 (алгебру дуальных чисел вида а + еЪ, иногда

И

е, причем б2 = 0. Отметим, что алгебра И2 (коммутативная, ассоциативная и с единицей) имеет точное матричное представление, матрицы которого имеют вид (0 ьа), где а,Ь € И,

И

тривиальна, и поля комплексных чисел С, группа автоморфизмов которого изоморфна (с единственным нетривиальным автоморфизмом — комплексным сопряжением), алгебра И2 имеет одномерную группу автоморфизмов — она состоит из преобразований вида а + еЬ ^ а + еаЬ при произвольных ненулевых а € И,

Положим Ф = Ф  И2 — это тензорное произведение алгебр Ф и И2. Можно записать элементы ф в виде а + еЪ, где а,Ь € Ф . Операция коммутирования в алгебре ф задается естественным образом (коммутирование по первому и умножение по второму тензорному сомножителям), при этом получаем структуру алгебры Ли, Фактически ф можно рассматривать как полупрямую сумму Ф+а^ Ф подалгебры, изоморфной Ф, и абелева идеала, соответствующее присоединенному действию а^ адгебры Ли Ф на себе как на векторном пространстве. Ясно, что если Ф — алгебра Ли векторных полей, то мы такой конструкцией получаем приближенную алгебру Ли из предыдущей конструкции.

Эту конструкцию можно обобщать. Например, рассмотрим алгебру вида Ир = И [б]/ < бр >, где И [б] — алгебра полиномов от б, а < ер > идеал, порожденный элементом ер. Получаем ассоциативную и коммутативную алгебру с единицей размерности р над И, она иногда называется алгеброй плюральных чисел. При р = 2 эта алгебра уже была введена выше. Рассмотрим теперь тензорное произведение Ф  с естественной операцией коммутирования — полученную алгебру Ли можно рассматривать как приближенную алгебру Ли при использовании приближений порядка р с помощью о(бр). Она представляется в виде полупрямой суммы подалгебры, изоморфной Ф и нильпотентного идеала вида бФ + б2  Ф    + ер-1  Ф,

2. О симметриях уравнений

Приступим к более подробному изучению приближенных еимметрий уравнений. Основная наша цель — построение геометрического подхода и методики изучения приближенных

групп симметрии дифференциальных уравнений (ОДУ и УРЧП), Начнем же с приближенного решения самих уравнений. Будем вести рассмотрение на простейшем примере — рассмотрим обыкновенное дифференциальное уравнение первого порядка и задачу Коши для него на некотором отрезке (а,Ь):

и' = Р

'и\х=хо и .

Здесь и(х) — неизвестная функция. Уравнение при б = 0 будем называть точным. Рассмотрим группу (бесконечномерную группу Ли) локальных диффеоморфизмов интервала изменения переменной х (окрестности точки жо). Она действует на искомые функции и и на дифференциальное уравнение в целом. Стационарные подгруппы этого действия при каждом б — это и есть группы (локальные) Ли симметрий рассматриваемых дифференциальных уравнений при данных е. Предположим для простоты, что эти группы симметрий конечномерны (так бывает довольно часто для важных и интересных дифференциальных уравнений). Тогда, как известно, при малых значениях б размерность (конечная!) стационарной подгруппы не больше, чем у стационарной подгруппы при б = 0 (т.е. при деформации «точки» размерность ее стационарной подгруппы не возрастает). Это явление при изучении симметрий дифференциальных уравнений изучается очень подробно. Если размерность группы симметрий при малых б постоянна, то говорят об устойчивых снмметрпях и о возможности продолжения симметрии уравнения, соответствующего е = 0 (т.е. точного уравнения). Если же размерность группы симметрий при малых е = 0 падает, то говорят, что некоторые симметрии точного уравнения не допускают продолжения, и потому даже для точного уравнения их можно и не рассматривать (т.е. они неустойчиво зависят от параметра). Полезно отметить, что свойство устойчивости группы симметрий для разных, хотя и близких между собой, уравнений может проявляться по-разному. Подчеркнем, что выше речь шла о группах симметрий, понимаемых в обычном (точном) смысле этого слова. На самом деле рассматриваются и приближенные симметрии, к которым указанные выше соображения о размерности стационарных подгрупп напрямую применять пока не получается. Так что приведенные выше соображения о характере изменения размерности групп симметрий следует понимать в основном как наводящие соображения.

Для изучения дифференциальных уравнений, зависящих от параметра, и нахождения приближенных групп Ли их симметрий было предложено два метода (имеющих довольно точные географические локализации).

Один из них (развиваемый изначально в уфимской математической школе) основан на линеаризации уравнений и дифференциальных операторов симметрии (через рассмотрение их с точностью до о(е)-, иногда используются и приближения более высоких порядков). Здесь возникают принципиально новые математические объекты, которые были названы приближенными группами преобразований (с соответствующими им приближенными алгебрами Ли), Были доказаны (копируя классические доказательства в «точном» случае) аналог точной теоремы Коши о существовании приближенного решения для приближенной задачи Коши и другие аналоги классических результатов. Активно изучаются приближенные симметрии получающихся приближенных уравнений.

Второй метод был создан в киевской математической школе, и он более классичен — рассматривается разложение неизвестной функции по параметру, У получившегося приближенного уравнения рассматриваются точные (классические) симметрии.

При применении двух этих методов к конкретным уравнениям они довольно часто дают различные результаты. Имеется несколько работ (например [6]), в которых на целом ряде дифференциальных уравнений из различных разделов физики и механики, для которых

можно вычислить точные группы симметрии при произвольных значениях параметра, выясняется, какой из этих двух методов дает более точные результаты, В [6] отдается предпочтение второму методу, который, однако, в целом получил среди специалистов меньшее распространение. При этом следует отметить, что группы спмметрнй дифференциальных уравнений — весьма тонкий объект, и многое зависит от природы изучаемого уравнения и от выбора малого параметра. Поэтому «статистический анализ» при сравнении эффективности двух указанных методов представляется не очень уместным.

Как уже было отмечено выше, в данной статье мы будет изучать именно первый метод и те новые математические объекты, которые с ним связаны (второй метод в этом смысле намного беднее),

3. Приближенные и дуальные алгебры Ли

Выше было отмечено, что можно ввести понятие приближенной алгебры Ли векторных полей — здесь обычные условия, задающие алгебры Ли, заменяются на их приближенные аналоги. Но более естественно перейти к более общим конструкциям, что мы и сделаем.

Приближенными алгебрами Ли мы назовем поначалу алгебры Ли вида Ь  И2 _ тензорные произведения обычных алгебр Ли Ь и двумерной алгебры И2 =< 1,е > (алгебры дуальных чисел). Отметим, что такие алгебры Ли можно рассматривать и как Д2-модули. Более того, приближенные алгебры Ли симметрий приближенных дифференциальных уравнений тоже имеют структуру Д2-модулей, Это следует из такого простого факта: если векторное поле X является инфинитезимальной симметрией дифференциального уравнения, то и векторное поле еХ — тоже. Тем самым алгебра Ли инфинитезимальных симметрий инвариантна относительно умножений на элементы алгебры И2.

Кроме алгебр Ли вида можно (и нужно!) рассматривать и их подалгебры. Но если

мы интересуется алгебрами Ли приближенных симметрий, то естественно рассматривать только такие подалгебры, которые являются Д2-модулями, Именно их и имеет смысл называть приближенными алгебрами Ли (в том понимании, какое возникает при изучении приближенных дифференциальных уравнений и их симметрий), Тем самым мы приходим к необходимости рассматривать класс всех алгебр Ли, определенных над 02. Их мы и будем называть приближенными алгебрами Ли, Итак, вот окончательное определение:

Определение 1. Приближенная алгебра Ли — это алгебра, Ли над алгеброй Б2 дуальных чисел. Другими словами, это алгебра Ли, которая дополнительно обладает структурой Б2-модуля, согласованной с операцией коммутирования.

При этом структура Д2-модуля задается линейным оператором, соответствующим умножению на е. Мы будем обозначать этот оператор через £. Естественно изучать и соответствующие (инвариантные относительно £) подалгебры Ли, идеалы, гомоморфизмы и т.п. Отметим, что как Д2-модуль алгебра Ли может быть и тривиальной (когда оператор £ — нулевой). Тем самым обычные алгебры Ли, для которых как раз £ = 0, можно рассматривать и как Д2-алгебры Ли, Однако в дальнейшем мы будем в основном полагать, что £ = 0, Более того, иногда полезно рассматривать только «невырожденные» Д2-структуры,

£

полагается, естественно, четномерной). Именно таковы алгебры Ли вида Ь  И2.

Иод базисом Д2-алгебры Ли можно понимать минимальный набор векторов, порождающих эту алгебру Ли как Д2-модуль, Говоря другими словами, это такой минимальный набор векторов, которые вместе с их произведениями на е порождают алгебру Ли как векторное пространство над И, Такого рода базисы были введены в [4] и там они названы существенными. Отметим, что только для алгебр Ли над И2 имеет смысл понятие существенного базиса. Число элементов в нем — это некоторый не совсем стандартный аналог размерности рассматриваемой алгебры Ли над И2. Например, подалгебра Ли в

Ь  И2, порожденная элементем вида еХ — одномерна над И, а порожденная элементом X € Ь — двумерна,

В качестве примеров Д2-алгебр Ли можно брать алгебры Ли вида Ь  И2 для произвольных вещественных алгебр Ли, В частности, если алгебра Ли Ь линейна, то при тензорном умножении ее па И2 просто матричные элементы этой алгебры Ли берутся лежащими в И2 (е соответствующими алгебре Ли Ь соотношениями). Например, если Ь = в1(п, И), то получаем алгебру Ли в1(п, И2), состоящую из матриц порядка п с элементами из И2, причем след этих матриц равен 0, Отметим, что если алгебра Ли Ь полупроета, то Ь  И2 уже полупростой никогда не будет (все алгебры Ли вида Ь  И2 имеют ненулевой радикал),

4. Гладкие многообразия с ^-структурой

Введя алгебры Ли над И2, естественно на этом не остановиться и рассматривать гладкие многообразия, имеющие дополнительную Д2-етруктуру, Это позволит рассматривать приближенные дифференциальные уравнения и многие связанные с ними понятия (в частности, группы приближенных еимметрий) более естественно, используя геометрические методы.

Отметим, что при изучении приближенных уравнений нам можно было бы и не выходить за пределы обычных многообразий, что и делалось до сих пор. Используя отношение можно ввести понятие приближенной алгебры Ли векторных полей — это будет такой набор векторных полей, который по модулю о(е) замкнут относительно коммутирования и удовлетворяет обычным по форме аксиомам алгебры Ли, но выполненным по модулю о(е). Но такой подход представляется искусственным, и мы его развивать далее не будем. Переходим теперь к описанию понятия ^-многообразий и связанных с ним понятий (касательные пространства, отображения, группы Ли и др.). Отметим, что соответствующие геометрические построения были подробно развиты в рамках дифференциальной геометрии в казанской математической школе [7],

Рассмотрим понятие многообразия над алгеброй дуальных чисел (аналогичны построения и над алгебрами плюральных чисел, но мы здесь на этом останавливаться не будем). Многообразия определяются по стандартной схеме — карты, диффеоморфные областям в Д2, склеиваются с помощью функций перехода, которые сохраняют Д2-етруктуру, Дадим более подробное описание.

Рассмотрим пространство Д2 — это свободный Д2-модуль над алгеброй И2, имеющий п образующих. Как векторное проетранетво над И от изоморфен И2", Но, в отличие от пространства И2", его элементы допускают умножение на элементы а + еЬ дуальной алгебры И2, что задает — например, с помощью оператора £ — на нем Д2-етруктуру,

Гладкое отображение / : и ^ V области и С в область V С Б™ называется гладким ^-отображением (иногда употребляют термины «голоморфное» или «аналитическое» отображение по аналогии с комплексным случаем), если его дифференциал сохраняет естественную Д2-структуру на касательных пространствах (отождествляемых с Д2 и И™ соответственно). Другими словами, оператор умножения на е должен быть перестановочен с дифференциалом.

Вкратце рассмотрим условие Д2-гладкоети или «голоморфности» для функций одной переменной (в общем случае это — условия Шефферса, впервые изучившего гладкие отображения над алгебрами и которые являются аналогом условий Коши-Римана), Условие голоморфности для функции Р(г) = f (х, у) + ед(х, у) дуальной перемеиной г = х + еу можно записать в виде условия на дифференциал, который имеет вид ¿Р = (А + еВ)  ¿г (где А = А(г),В = В (г) € И при фиксированном значении аргумента г). Раскрывая скобки и приравнивая вещественные и «мнимые» части, получаем:

гх = А,д'х = В,ГУ = 0,д'у = А.

Отсюда следует, что f (х,у) = ф(х) — произвольная функция от ж, а д(х,у) = у  фх + ф(х), где ф(х) — еще одна произвольная функция от х. Тем самым получаем общий вид Д2-гладких функций одной переменной

Р(х + еу) = ф(х) + е(у  фх + ф(х)).

Отметим, что если у = 0 т.е. если аргумент вещественен, то значения Д2-гладкой функции имеют вид Р(х) = ф(х) + еф(х), т.е, Р(х) является произвольной гладкой вектор-функцией от х. Отсюда следует, что функция дуального аргумента однозначно определяется своими значениями при вещественных значениях аргумента. Если рассматривать Р(х + еу) как вектор-функцию двух переменных, то ее матрица Якоби имеет вид

( Ф'(х) 0 \

ууф'(х)+ф' (х) ф'(х) ) .

Ясно, что эта матрица перестановочна с матрицей (10) оператора £ на И2, что соответствует определению дифференцируемой функции дуального аргумента.

Так как в качестве компонент Д2-гладкой функции могут фигурировать произвольные гладкие вещественные функции, то ясно, что Д2-гладкие функции далеко не всегда будут иметь разложения в степенные ряды (которые можно ввести, используя понятие мультипликативной нормы на И2, указанной ниже). Поэтому название «голоморфные» для Д2-функций представляется не очень удачным, так как оно по ассоциации с ТФКП наводит на мысль об аналитичности, которой в Д2-еитуации может и не быть.

Аналогично можно вывести и общий вид Д2-гладких функций нескольких переменных; об этом можно прочесть, например, в [7],

Отметим, что с помощью разложения по формуле Тейлора можно строить — причем совершенно естественно — продолжение обычных гладких функций в дуальную область значений аргумента. Для функции и(х) такое продолжение будем об означать и(х). Напомним, что рассмотрение выражений вида а + ша0, где а, а0 — трехмерные векторы, а ш — символ, для которого ш2 = 0 используется в механике как объект винтового исчисления (см,[5]), Именно из винтового исчисления были позаимствованы и многие формулы при изучении функций дуального аргумента.

Переходим теперь непосредственно к определению ^-многообразий. Топологическое пространство М называется гладким ^-многообразием, если оно хауедорфово, удовлетворяет 2-й аксиоме счетности и имеет атлас с картами в Д2, функции перехода которого являются Д2-гладкими функциями. Но можно подойти и по-другому: пусть М — обычное четномерное гладкое многообразие. Оно называется ^-многообразием (или многообразием, имеющем Д2-структуру), если для некоторого его атласа (набора карт) соответствующие функции перехода являются Д2-гладкими, Так или иначе, для ^-многообразий их касательные проетранетва имеют Д2-етруктуру, Отображения в этой категории — это гладкие отображения, сохраняющие Д2-структуру (в частности, дифференциалы которых являются И2-отображениями), Используя понятие оператора £, соответствующего Д2-етруктуре, условие гладкости отображения f : М ^ N записываем в виде $(£м) = £м ((#).

Отметим, что для произвольного гладкого Д2-многообразия М пространство Т(М) его касательного расслоения (о касательных векторах см, ниже) имеет естественную Д2-етруктуру (подробнее о Др-структурах па касательных пространств высших порядков см, например [7]), В частности, ^-многообразие можно рассматривать как гладкое многообразие, моделью для которого служит касательное расслоение над произвольным евклидовым пространством.

Отметим также, что для задании Д2-топологии мы можем пользоваться обычной метрикой в евклидовом пространстве. Однако далее нам потребуется существование на И2

мультипликативной нормы. Сам факт существования мультипликативной нормы доказан, как известно, для произвольных банаховых алгебр и потому верен и для двумерной алгебры И2. Однако нам будет удобно иметь в своем распоряжении мультипликативную норму на И2, заданную явно.

Для произвольного дуального числа а + еЬ положим ||а + еЪ\\ = 1тах(| а |, | Ь |),

Лемма 1. Норма ||а + еЪ\\ = 1 таж(| а |, | Ь |) на Б2 является мультипликативной, т.е. \1ХУ\1 < |М||Ы| для ЛК)бых х,у € Б2.

С помощью указанной нормы можно построить «стандартное» дифференциальное исчисление для функций над 02. Вводятся понятие предела (по норме), производной, неопределенного и определенного интегралов (последний - как предел интегральных сумм). Все привычные свойства сохраняются. Далее, вводятся функциональные ряды, степенные ряды и т.д. Далее здесь мы не будем далее развивать эти перспективы, однако в другом месте они будут нами использованы.

Мы привели стандартное определение многообразия над алгеброй И2. Соответствующую Д2-структуру будем называть невырожденной. Но оказывается, что естественно рассматривать и более общий случай. Это необходимо, например, при рассмотрении фактор-пространств групп Ли и однородных пространств. При этом можно использовать понятие е-етруктуры па многообразии — так будем называть поле линейных операторов в касательных пространствах (некоторые геометры предпочитают говорить о поле аффиноров) на многообразии, квадрат которых равен нулю. Тем самым задается и поле Р образов этих операторов, а также поле их ядер — мы получаем два распределения на многообразии, При этом возникает вопрос об интегрируемости этих распределений, на чем мы здесь подробно останавливаться не будем.

Наличие на многообразии М Д2-етруктуры эквивалентно тому, что на М существует

£

£

распределение (половинной размерности в невырожденном случае) на ^-многообразии, которое (в силу самого определения ^-многообразия) интегрируемо и тем самым определяет слоение (называемое обычно каноническим). При этом все И2-диффеоморфизмы (т.е. взаимно однозначные и взаимно Д2-гладкие отображения) сохраняют это слоение. Оно является важным инвариантом в категории ^-многообразий.

Под Д2-гладкой кривой на многообразии М понимается Д2-гладкое отображение 7 : и ^ М, где и С И2 — некоторое открытое подмножество в И2 (можно в качестве и брать некоторый круг относительно введенной выше нормы).

Естественным образом вводится понятие касательного вектора - например, определяя их как дифференцирования ростков Д2-гладких функций в данной точке или через отношение эквивалентности Д2-кривых в данной точке. Получаем касательное пространство ТХ(М), которое естественным образом имеет структуру свободного Д2-модуля,

Естественным образом на ^-многообразии определяются далее Д2-дифференциальные формы и другие объекты из стандартного набора дифференциальной геометрии. Напомним, что при этом умножение на е € И2 в касательном пространстве мы рассматриваем

£

Далее, естественно вводится касательное расслоение произвольного ^-многообразия М (которое можно рассматривать и как касательное расслоение многообразия М, если игнорировать его Д2-етруктуру), Тогда Д2-векторным полем на М называется Д2-гладкое сечение этого расслоения. Множество всех Д2-векторных полей на многообразии образует

алгебру Ли над И2. Произвольные Д2-векторпые поля на многообразии М можно рассматривать как инфинитезимальные преобразования этого многообразия, сохраняющие на нем Д2-етруктуру,

5. О дифференциальных уравнениях с дуальным аргументом

Здесь мы рассмотрим вопрос о задаче Коши в Д2-случае,

Теорема 1. Пусть функция /(х, и) является 02-гладкой в некоторой окрестности, точки (х0,й0). Тогда, задача, Коши

(1и

Тх = 1 (х,и) = й0

имеет решение в некоторой, окрестности, точки х0, причем, единственное.

Доказательство

Задача Коши может быть переписана в следующем виде:

г-х

и(х) = й0 + / {(х,и)йх.

</хо

Это соотношение рассматривается в Д2-области, о понятии интеграла в которой уже было сказано выше. Далее можно следовать, например, схеме доказательства из книги В,И, Арнольда [8], которое основано на отображении Пикара и применении процесса итераций (или, что здесь одно и то же, теоремы о неподвижной точке). Все рассуждения из этой книги автоматически переносятся в нашу ситуацию, если только вместо абсолютной величины числа использовать введенную выше мультипликативную норму для дуальных чисел. Нет необходимости полностью воспроизводить здесь все рассуждения из [8],

Можно продолжить изучение дифференциальных уравнений дуального аргумента, следуя, например, тому, как это сделано в книге [8] для вещественного аргумента. Вот несколько результатов такого рода. Доказательства их здесь не приводятся, так как они буквально копируют доказательства из [8] (с заменой абсолютной величины вещественного числа на норму дуального числа).

Теорема 2.

— = / (х,и,а) ах

й\х=хо = ^0

в некоторой, окрестности, точки х0 имеет, причем, единственное, решение, гладко зависящее от па,ра,м,етра, а.

Если исходные данные задачи, Коши гладко зависят от параметра, а, то и решение этой задачи, Коши тоже гладко зависит от указанного параметра.

Справедлива, конечно, и теорема о зависимости (непрерывной, гладкой) решения задачи Коши от начальных условий (которые можно рассматривать и как параметры).

Справедливы также и теорема о максимальном продолжении решений (вплоть до границы области, где задано уравнение). Доказательство основано на Теореме 1, в силу которого решение существует в окрестности каждой допустимой областью определения точке. Затем берется объединение всех таких окрестностей и выделяется связная его компонента,

содержащая точку, заданную начальным условием задачи Коши, Решение в этой области является, очевидно, максимально продолженным.

Верен, конечно, и тот факт, что для линейного дифференциального уравнения решения определены там, где определены правые части уравнений. Другими словами, практически все утверждения классической теории обыкновенных дифференциальных уравнений почти автоматически в рамках изложенного выше геометрического подхода переносятся на случай приближенных уравнений и их приближенных решений. Так же автоматически получаем и критерий приближенной инвариантности уравнений, сформулированный выше,

6. О Д2-теории Ли

Следующий наш шаг — введение понятия D2-rpynn Ли. Группы Ли G над D2 — это такие группы Ли, которые имеют Д2-структуру на группе Ли G как па многообразии и на отображениях, входящих в стандартное определение группы Ли (т.е. отображение группового умножения и отображение обратного элемента должны быть Д2-гладкими), Касательное пространство Д2-группы Ли в единичном ее элементе естественным образом имеет структуру Д2-адгебры Ли (т.е. алгебры Ли, являющейся Д2-модулем, о чем уже было сказано выше). Наконец, естественным стандартным образом определяется понятие Д2-дейетвия Д2-группы Ли па ^-многообразии. Это действие задается Д2-гладким отображением G х M ^ M при условии выполнения обычных соотношений, задающих действие группы на множестве, В частности, мы имеем стандартное определение однопа-раметричеекой подгруппы в категории Д2-объектов (причем параметр подгруппы принимает значения из D2, хотя иногда можно ограничиться и одними только вещественными значениями параметра),

Если па ^-многообразии задана Д2-гладкая однопараметричеекая группа преобразований, то касательные векторы ее орбит дают нам Д2-векторное поле (поле скоростей). Докажем, что верно и обратное — векторному полю на ^-многообразии M соответствует локальная однопараметричеекая (над D2 и двумерная над R) группа преобразований этого многообразия М. Это даст нам аналог классической теоремы С, Ли:

Предложение 1. Пусть X(z) - гладкое D2-eeKmopnoe поле, заданное в некотором открытом, подмножестве U С D2. Тогда, существует локальная однопараметричеекая группа преобразований некоторого открытого подмножества W С DV;, для, которой соответствующее касательное векторное поле совпадает с X.

Если, же X - векторное поле на, компактном 02-м,ногообра,зии М, то ем,у соответствует однопараметричеекая группа преобразований многообразия М.

Доказательство получается стандартным процессом интегрирования соответствующего дифференциального уравнения (отметим, что параметр в этом уравнении пробегает алгебру D2 или некоторое открытое подмножество в ней). Основой в локальном случае является теорема о существовании и единственности решения задачи Коши (см,§5 выше), Что касается глобального случая, то здесь с учетом уже рассмотренного локального случая применимы стандартные рассуждения, используемые в классической теории Ли,

Далее, стандартными рассуждением можно доказать, что, как и в классическом случае, в категории Д2-групп и алгебр Ли имеет место со ответствие Д2-подадгебр Л и и D2-подгрупп Ли (локальных). Мы не будем стараться описать в данной статье все множество возможностей, которые возникают при рассмотрении категории Д2-групп Ли и связанных с ними объектов. Поэтому мы ограничиваемся только перечислением некоторых важных шагов на пути изучения этих довольно новых объектов.

Естественно рассматриваются замкнутые Д2-подгруппы Ли H С Си соответствующие им однородные пространства. Но фактор-пространство G/H далеко не всегда обладает

естественной Дг-етруктурой, Точнее, ото обладает некоторой e-структурой (которая переносится с помощью экспоненциального отображения), которая индуцируется естественной Д2-структуры на фактор-пространстве соответствующих алгебр Ли, индуцированной оператором При этом получается не обязательно певырожденная Д2-етруктура, именно поэтому нам приходится рассматривать и вырожденные Д2-структуры на многообразиях. Например, пусть Н — связная подгруппа Ли некоторой Д2-группы Ли, соответствующая абелевой подалгебре Ли 1т£ (подробнее об этой подалгебре Ли говорится в другой статье автора), и пусть эта подгруппа замкнута (так будет, например, если G — одноевязная ниль-потентная группа Ли), Тогда на фактор-пространстве G/H индуцируется тривиальная D2-структура (для которой оператор £ нулевой). При этом на самой под группе Ли Н тоже Д2-етруктура тривиальна. Но если па некоторой замкнутой подгруппе Ли F индуцируется невырожденная Д2-структура (что бывает hg очвнь часто), то на фактор-пространстве G/F мы все же получаем полноценную невырожденную Д2-етруктуру,

7. О приближенных дифференциальных уравнениях

Данные выше определения позволят нам переформулировать, следуя указанному выше принципу А. Гротендика, методику анализа дифференциальных уравнений с параметром в терминах приближенности и о(е) в новом, более естественном, геометрическом оформлении, Более того, мы получаем обоснование принципа соответствия, сформулированного в самом начале статьи. Конкретизируем теперь эту методику.

Пусть задано некоторое дифференциальное уравнение, зависящее от параметра:

du

dx = f

и некоторое начальное условие: ulx=xo = u°

В принципе, даже независимую переменную можно было бы предполагать зависящей от того же параметра е (см, ниже). Мы не будем вдаваться здесь в другие возможные обобщения.

Будем функцию f и решение и рассматривать с точностью до о(е). В [1] и последующих работах был изложен метод вычисление приближенных решений и групп симметрии соответствующих дифференциальных уравнений.

Рассмотрим этот метод для указанного выше уравнения. Положим и(х, е) ~ и0(х) + eui(x)-это будет приближенная неизвестная функция. Далее, имеем:

f (х, и, е) ~ fo(%, ио(х) + 6Ul(x)) + tfl(x, щ(х) + tul(x)). Преобразуя это выражение по модулю о(е), получим приближенное дифференциальное уравнение:

и'0(х) + tul (х) = fo(x,uo(x)) + e(fo)'u(x,uo(x))ui(x) + efi(x,uo(x)) + 0. Приравнивая подобные члены (относительно е), получаем систему уравнений:

и'о = fo

u'i(x) = Ш'и(х,ио(х)))щ(х) + fi(x,Uo(x)). Первое уравнение — это исходное точное уравнение (соответствующее параметру е = 0), Второе уравнение решается уже в предположении, что из первого уравнения функция и0(х) найдена, В результате получаем приближенное решение

и(х) ~ щ(х) + tui(x).

Учет начального условия в задаче Коши делается стандартным способом, и потому мы на этом останавливаться не будем.

Далее находятся симметрии рассматриваемого нами дифференциального уравнения с параметром е. Для этого используются «приближенные» векторные поля

X(х, б) « Х0(х) + еХ1(х).

Условие приближенной инвариантности дифференциального уравнения имеет вид

X(х, б)  /(х, и(х)) « 0.

Применяя тот же метод, что при решении выше дифференциального уравнения, можно получить метод последовательного нахождения векторных полей (точнее, их компонент) ^ и Этим мы вкратце изложили метод, впервые предложенный в [1].

А вот теперь мы рассмотрим модификацию этого метода, используя изложенные выше геометрические конструкции. Мы получим более короткий, концептуальный подход.

От уравнения и' = $(х,и,е) перейдем к его дуальному продолжению:

и' = £(х, и).

Начальное условие запишем в виде

= & = + .

Объекты с тильдой можно рассматривать как точные на многообразии над алгеброй И2, Здесь «приближенное» дифференциальное уравнение уже записывается как точное. То же касается и задачи Коши, Поэтому, применяя к многообразию над И2 классические результаты о дифференциальных уравнениях, приходим к концептуальному доказательству основных общих результатов теории приближенных решений дифференциальных уравнений, В частности, доказана:

Теорема 3.

Для, любой гладкой функции /(х, и, б) и гладкой функции щ(б) задача, Коши

и' = /

и1х=хо = Щ

всегда, имеет в некоторой окрестности, точки х0 приближенное решение, причем, единственное.

Напомним, что приближенное решение - это такая функция, которая с точностью до о(б) дает решение исходного дифференциального уравнения с малым параметром.

Точно так же получаем и другие утверждения о приближенных решениях дифференциальных уравнений, В [1] и подобных статьях соответствующие утверждения сопровождались доказательства (копирующими классические, но все равно необходимыми). При нашем же подходе нужда проводить такие доказательства часто отпадает. Результаты сразу следуют их стандартных (не приближенных) утверждений, хотя и примененных к отличным от исходных многообразиям и другим Д2-объектам,

До сих пор мы находились в пределах уже подробно разработанной ранее теории работы с приближенными решениями и симметриями. Нами была только дана их геометрическая интерпретация. Выйдем теперь за эти пределы и переместимся в намного менее разработанную область. При изучении дифференциальных уравнений предлагается перейти к дуальным (в нашей терминологии) значениям не только для уравнения и неизвестной функции, но и для независимой переменной. Другими словами, полностью перейти в категории объектов над дуальными числами.

Отметим, что идея варьировать и зависимую переменную встречалась и раньше, но в другой ситуации. Например, в [9] дается обобщение метода из [1], в котором по параметру раскладывают и зависимую и независимую переменные. Мы рассматриваем еще одно обобщение метода из [1].

Мы опишем построенный на основе введенных выше геометрических конструкции некоторый новый метод.

Здесь линеаризации подвергаются не только уравнение (при нахождении приближенных его решений) и дифференциальный оператор — инфинитезимальный оператор симметрии уравнения, но и все остальные величины, входящие в уравнение — независимые и зависимые переменные. При этом возникают дополнительные функциональные параметры, которые можно использовать для упрощения процесса вычислений (обращая с их помощью в нуль некоторые слагаемые).

Рассмотрим этот обобщенный метод на простом примере. Пусть дана задача Коши для обыкновенного дифференциального уравнения первого порядка, зависящая от малого параметра е:

— = / (x,u,t), ах

и\х=хо = И°(б).

Рассмотрения будем вести с точностью до о(б).

Используя е-продолжение гладких функций, уже упоминавшееся выше, исходное дифференциальное уравнение перепишем в виде уравнения для Д2-функции относительно Д2-переменной с ^-аналитической правой частью:

— = / (x,u), ах

Применяя к этому Д2-уравнению известные методы теории дифференциальных уравнений, разработанные для вещественной и комплексной областей, мы будем получать новые результаты, касающиеся приближенных решений. Однако нужно отметить, что при этом мы фактически вводим две независимые переменные х,у (х = х + еу) и нужно (подробнее см, ниже) указать связь между ними (заметим, например, что в выражение для || входит Выбор указанной связи - это отдельный шаг описанного метода.

Покажем, как можно вычислять производную ^ функции дуального аргумента г = х + еу. Здесь нам понадобится формула для деления дуальных чисел. Отметим, что деление это удается провести не всегда.

Имеем = ^(1 — е^) (проверяется, например, прямым вычислением). При а = 0 такое

деление возможно и однозначно. Рассмотрим теперь производную ^ =  Используя

приведенную выше формулу для обратного числа, получим ^ = ^ (1 — б ^), Если мы зададим у = р(х) как функцию ж, то получим формулу и' = и' — ер'и', где штрих означает производную по х. Развитие и использование указанного нового метода мы предоставляем специалистам в этой области,

В заключение отметим, что переход к дуальным числам при решении и исследовании уравнений с малым параметром заметно упрощает даже самые рутинные вычисления. Например, пусть нам нужно решить такую задачу Коши для линейного ОДУ второго порядка:

и" = (1 + е)и — 2б

ulx=o = 1 - 2e,u'|ж=о = б Будем решать его в Д2-области, На самом деле нечто подобное можно делать и с использованием символа о(е), но при этом приходится постоянно следить за корректностью проводимых упрощений формул, В дуальной области корректность достигается автоматически, Нужно только использовать несколько простых формул: еа+еЬ = еа(1 + eb), sin(a + eb) = sin(a) + eb • cos(a), (1 + e)d = 1 + ed и т.п. В частности, решая приведенную выше задачу, для характеристических корней соответствующего однородного уравнения получим выражения \12 = ±(1 + е/2). Далее, для общего решения однородного уравнения получим (используя приведенную выше формулу для экспоненты):

Uo = Ciex + с2е-х + 2(Ciex - С2е-Х).

Методом подбора частного решения для нашей простой правой части находим частное решение и* = 2е, а затем и общее решение нашего дифференциального уравнения. Все вычисления — чисто арифметические. Наконец, два начальных условия дают «точные» в дуальной области, но приближенные в содержательной интерпретации выражения для решения поставленной задачи Коши (С1 = ,С2 = ). По аналогичной схеме можно рассмотреть и множество других примеров, в том числе и тех, которые были рассмотрены в многочисленных статьях, посвященных решению приближенных дифференциальных уравнений. Также подобными вычислениями решаются многие задачи по вычислению симметрий дифференциальных уравнений. При этом, рассматривая случаи е = 0 и е ^ 0, выделяются симметрии, которые разрушаются при ненулевых значениях параметра е. Но мы не будем входить здесь в соответствующие подробности, так как целью данной статьи было не проведение новых вычислений, а выработка нового, геометрического подхода к таким вычислениям.

СПИСОК ЛИТЕРАТУРЫ

Байков В.А., Газизов Р.К., Ибрагимов Н.Х. Приближенные симметрии // Мат. сборник. Т. 136. 1988. С. 436-450.

Байков В.А., Газизов Р.К., Ибрагимов Н.Х. Методы возмущений в групповом анализе // Итоги науки и техн. Сер. Соврем, пробл. мат. Нов. достиж. Т. 34. 1989. С. 85-147. N.H. Ibragimov CRC handbook of Lie group analysis of differential equations. Vol.3:New trends in theoretical development and computational methods, CRC Press, FL. 1996.

Байков В.А., Газизов P.K., Ибрагимов Н.Х. Приближенные группы преобразований // Дифференциальные уравнения. Т. 29.1993. С. 1712-1732. Дименберг В.Ф. Винтовое исчисление. М.: Наука. 1965.

М. Pakdemirli, М. Yurusoi, I.T. Dolarci Comparison of Approximate Symmetry Methods for Differential Equations // Acta Appl. Math. V. 80. 2004. C. 243-271.

Вишневский В.В., Широков А.П., Шурыгин В.В. Пространства над алгебрам,и. Казань : Изд-во КГУ. 1985.

Арнольд В.И. Обыкновенные дифференциальные уравнения. Ред. журнала «Регулярная и хаотическая механика». 2000.

Z.-Y. Zhang, X.-L. Yong, Y.-F. Chen A new method to obtain approximate symmetry of nonlinear evolution equation from perturbations // Chinese Phisics, В. V. V. 18. 2009. P. 2629-2633.

Владимир Витальевич Горбацевич, П-626 Проспект Мира, д. 116А, кв. 68, 129626, Москва, Россия E-mail: vgorvich@yandex.ru