Вычисление собственных чисел эллиптических дифференциальных операторов с помощью теории регуляризованных рядов Текст научной статьи по специальности «Математика»

УДК 517.5

DOI: 10.14529/mmph160205

ВЫЧИСЛЕНИЕ СОБСТВЕННЫХ ЧИСЕЛ ЭЛЛИПТИЧЕСКИХ ДИФФЕРЕНЦИАЛЬНЫХ ОПЕРАТОРОВ С ПОМОЩЬЮ ТЕОРИИ РЕГУЛЯРИЗОВАННЫХ РЯДОВ

С.И. Кадченко, О.А. Торшина

Магнитогорский государственный технический университет, г. Магнитогорск, Российская Федерация E-mail: kadchenko@masu.ru

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

Ключевые слова: дифференциальные операторы; спектральная теория; ре-гуляризованные следы; теория возмущений; собственные числа.

Рассмотрим на проективной плоскости F дифференциальный оператор

Т = -Д = --

1

sinddd

. в d 1 sin в---

12 Л

дв sin2 в dg2

действующий в гильбертовом пространстве Н. Собственные числа оператора Лапласа-Бельтрами образуют систему ортонормированных сферических функций. Обозначим

Хп = n(n +1) (n = 0,- собственные числа, vn = 2n +1 - кратность собственного числа Хп оператора Т; vn i (i = 0,2n) - собственные функции, ln ={Л\Л = Лп + n +1 + ip, - ^ < p < - прямые на комплексной плоскости, ¡nn t - собственные числа оператора T + P

i - n(n +1)| < const.

Определим собственные числа дифференциального оператора с помощью теории регуляризованных рядов. Рассмотрим четырех поправки теории возмущений. Первая поправка теории возмущений [1]

2n о ■ 1

X(Pvm,vm )=-пflp(Msin0d^ = 0.

i=1

2n

Вторая поправка теории возмущений ап (р) к сумме представляется формулой [3]:

i=0

« (p )=-2~Sp

j- i

Я

(T-ЯЕ) 1P2 ] (T-ЯЕ) 1 йЯ

= I

«

k ,n

1, k Фп\\ Яп

2k 2n

«

,k = II ((, U)((, %) •

i=0 j=0

F

Вычислим вторую поправку теории возмущений. Для этого введем так называемое Л преобразование f (a) функции p(0,p) [2]. Обозначим T(a) - пересечение конуса со сферой,

ф (в, в', a) - якобиан перехода в, p, в', p к координатам в, p, в', a.

a' \ л. sina

ф(в,в ,а) = ±. -—

J(-cosa + cos (0-0'))(cosa- cos (в + 0))

При этом знак «плюс» берется, если рр < p, знак «минус», если рр > p. Тогда [4]

а,«=

f (а) = JJ p (в,р) sin OdOdp JJ p (в, р) sin в

F T (а)

(2k + l)(2n + 1)П 2 2

sin а

4п2

■J f (а)ЁLi (P2k (cosa))XL (P2n (cosa))ssinada.

0 1=1 1=1 Очевидно, что функция /(а) нечетная, то есть /(а) = —/(п — а), 0 <а<п. Докажем, что / (0) = / (п) = 0. Для этого введем функцию [6]

Z(a) =

/ ч-2 „ . П

(sin a) , 0 <a< —

t ■ \-2 П

-(sin a) , — <a<n.

Если f (0) = f (п) ф 0, то J f (a)g(a)sinada = С другой стороны [5]

о

П

J f (a)g(a)sinada <

о

< JJ JJ P (в(Р)Р (в', P')sin в sin Od0d0'dpdp'

F F

sin a

< const

JJJJ

d0d0'dpdp'

< const

sin a

JJJJ

ф(в, 0',a)d0d0'dpda

F F

d0d0'dpda

sin a

F F

= const JJJJ- ----—

F F J(-cosa + cos (в - - cos (в + в))

<

В результате, /(0) = /(п) = 0. Докажем, что функция /(а) почти всюду дважды дифференцируема и /"(а)е Ц [0,п]. Запишем равенство [8]

J f (a)g(a)da = ±JJJJ p (0,p)p (в',p)g(a)sinвsinвx

dвdpd0da.

F F 1

cosa + cos (в-в'))(cosa- cos (в + в))

Далее имеем [7]

x

J f (a)g' (a)da =±JJJJ p (0,p)p (в',p')g'(a)sinвsinв'x

dвdpdв'da =

F F 1

<j(-cosa + cos (в-в)))

= ± JJ JJ dp{p (0,p)p (в^^И^тв'

F F

cosa-cos (в + в'))

,' p)} dp da

х . dddtydO'da =

J(- cos a + cos (в - в')) (cos a - cos (в + в))

= ЦЦp (в, ф) рф (в',р)с(а)твsin вх

F F

х . 1 dвdфdв'da±

J(-cosa + cos (в-в/))(cosa- cos (в + в))

±2ÍÍÍÍР(в,ф)р(в',фУ{-Sinag((a) -1 sinfínввва±

(-cosa + cos (в-в)) (cosa-cos (в + в))

F F

±1ÍÍÍÍР(в,Ф)Р(в,ф)-(-*тадв в)) (-1 sin^in^^Va

2 FF (-cosa + cos (в-в)) (cosa-cos (в + в))

Произведем замены: в первом интеграле в-в' = aсhф, в + в' = aсhф', а во втором в-в = a^ф, в + в' = ashфchф', в третьем интеграле в-в = aсhф, в + в' = ashфchф'. В результате по формуле Фубини получим f'(a)e Ц [0,п], где производная взята в смысле С. Л. Соболева [9]. Далее

п

Í f (a)g"(a)da = ÍÍÍÍ р (вф)рф(в',ф)х 0 FF

sin a

х-.----Y7----тc/(a)c(a)slnвsln в'dвdфdв'da +

(-cosa + cos (в-в ))(cosa-cos (в + в')) '

±2ÍÍÍÍР(в,Ф)р(в,Ф)( (-\-—)slnвsшв>dвdфdв>da±

2 ff (-cosa + cos (в-в)) (cosa-cos (в + в))

± 1 ÍÍÍÍ р (в,ф)р (в'фГ (-1-—щslnвslnв>dвdфdв>da.

2 ff (-cosa + cos (в-в)) (cosa-cos (в + в))

По аналогии в формуле для f''(a) получим

2

cosa sin a

1/О /О ' 1/О ^/О *

(-cosa + cos (в-в')) (cosa-cos (в + в)) (-cosa + cos (в-в')) (cosa-cos (в + в'))

Сделаем замены в + в' = аскф, в - в' = a4 sh фскф', -в + в = 2аскф, в-в = оО sh30ch0'. Воспользовавшись теоремой Фубини, получим, что f"(a)e Ц [0,п] и существует почти всюду [11].

Вычислим [10]

2n t

Ц (P2n (cosa)) + Ц (P2n (cosa)) = X

m=0

'2 - k2 . , _dL _ k2

. nm nm

sin- + cos- х

2 2

2n 2 2n 2n 7 2

х П p2n (cosa) = X п Pk (cosa) =

k=0,kФm m - k m=1 k=0,k*m m - k d 2 - k 2

V ГГ dф2 (n -s)! P(s), в P(s), в ,

= X П -+ \P2n (cosв ^P2n (cosacoss(ф-ф) =

m=1 k=0,k*m m - k s=0 ds (n + S)!

s 2 - к2

= 2ly((=+S)^2(:)(cOS^')^2(;) (cosacos s {<P-V')t П VtT =P2n

s=0ds (= + S )! m=1 к=0,к * mm - к

Тогда

ак ,= =

(2к +1)(2= +1)

4п

2 2

f f (a)ZLt (P2k (cosa))ZL (P2= (cosa))sinada =

i=1

2

Z

i=1

(2k +1)(2= +1)

4n2

f f (a)P2k (cosa)P2 = (cosa)sinada.

В результате вторая поправка теории возмущений

a= (p ) = -(2= Z l^^lf f (a)P2 к (cosa)P2 = (cosa)sin ad0C =

4П к=1, к * = \Лк -Л=\ 0

2к +1

(2= +1) £ 2к +1

f + f

0 Я-£

f (a)P2к (cosa)P2= (cosa)sin ada +

4п2 к=1,к*= К -A

я-е

+ f f (о:)Р2к (cosa)P2= (cosa)sin ada к

Воспользовавшись [5]

Р2к (cosa) =

cos

{(2к +12 )a-n¡ 4}

(sin a)

12

(2/nf + 0(1)

2к12 202 O (1)

(sina) 2к3/2 ' (sin a) 2к5/2 "

+

+ sin{(2к + 32)a-n/4} 0(1) +

Получим:

, ч 2т +1 ^ 2к +1

a= (pZ

4п2 к=1,к*= \Лк

f + f

0 я-е

f (a^inaP^ (cosa)P2 = (cosar )da +

+

f f (a)sin(a)\

cos

{(2к + 1)a-n 4}

sin{(2к + 32)a-n¡4} 0(1)

(sin a))2 0(1)

(2/nf O(1)

(2/ к )1/2 (2/ к )3/2 cos {(2к +1/2 )a-n 4}

(sin a) (2к )3/2 (sina)5/2 (2к )5/2

sin{(2к + 3/2)a-n/4} 0(1)

(sin a)2 0(1)

+

(2/nf0(1)

(2/ к )1/2 (2/ к )3/2

+

ff p (0»n0

,F

ein = 1+0 (_£i+01 e+0

(sin a) (2к )3/2 (sina)5/2 (2к )5/2

где е - положительное число

f (a) = ff p (0,p)sin0d0dp

F

С учетом асимптотических формул [4] получим a= (p ) = 0 (e2ln=) + 0

)da>,

=

+01 £|+0(£=, 4 ,

= ) v = ) 4n =

n

я-е

3/2

sin a

eln = | eln =

+ 0 —— I +

1

— f f (a)ctg(n-a)da + 0

= ) v =

eln =

n

+

+0^)+0 (^ ||+0 (=£)++0 (-)+0 (£W |+0 (

Выбирая в качестве £ = —ij- , найдем

n3'4

1 П-£

ап(p) =—— [ f (a)ctg(n-a)da+O 4пп

ln n

П—£

Учитывая, что j f (a)ctg (n - a)da = 0, получим

' an(p) = O(^j.

Вычислим третью поправку теории возмущений

ßn (p ) =-П- sp

Ьт

j - j

ln ln-1

3

P (T -AE) 1

3 1

i A =-Sp

-ni

j - j

£ (2n + 1)(2t +1)(2' + 1)Ax

(An -A)(A-Л)(Л, -A xf¿1 jjjjjjpC-^pCJ)pO")(cosa)Pn(cosß)Pn(cosy)x

F F F

X sin в sin в sin 6"d (pd pd p"d 6'd в", cos a = cos в cos в + sin в sin в' cos (p — p), cos в = cos в cos в" + sin в sin в" cos (p — p), cos у = cos в cos в" + sin в sin в" cos (p — p"). Если в шестикратном интеграле сделать замену p — n + p, в —П — в, то Pn (cos(n — a)) = Pn (—cosa) = (—1)" Pn (cosa), P¡ (cos(л — у)) = P¡ (—cos/) = (—1)lp (cos/). В результате

вп (P) = 0.

Для оценки четвертой поправки докажем теорему.

Теорема 1. Для оператора Лапласа-Бельтрами с дважды непрерывно дифференцируемым потенциалом на ln = {Л \ Л = Лп + п +1 + ip, — <^< p < справедливо равенство

(T — ЛЕ)—1

2 < Oil

2 |e| + n

где Ле ln и e = Im(Л)

Доказательство. Имеем

Ч—1 р ~

n—1

(T -AE)-1 = £

■= £

+

+

2 k=0 (Ak-An -n -1) + e2 k=o (Ak-An -n -1) + e2 (n +1) + <

+ £

k=n+1 (Лк —Л" — n — 1) + e Оценим каждое из слагаемых в отдельности

V,

2 2 1 2 3

n—1

A = £

k=0 (Ak -An -n -1)2 + e2 (2k + 1)dk 0 (k (k +1)-n (n +1)-n -1)2 + e2 (2k + 1)dk =

n—1

= const j

n—1

= const j

0

((n +1)2 -k(k +1)) + e2

п( п—1) = СО^ |

ёт

0(1)

« \( , 1\2 Г , 2 |е + п о 1(п +1) — т I + е

Оценим

Л =-

2п +1

0(1)п 0(1)

(п +1)+ е2 (п +1)+ е2 (п + |е|) \е\ + п

Здесь была использована эквивалентность норм |а| + |Ь| и 2 + |Ь|2 . Далее

2к +1

Аз = I

= I

к=п+1 (Лк — Яп — п — 1) + е к=п+1

(п +1)2 — к (к + 1)

+ е

Оценивая А3 аналогично А1, получим

Аз <

0(1) \е\ + п

Теорема доказана.

Учитывая результаты тео

эемы 1, рассмотрим четвертую поправку теории возмущений.

1

2П

Р (Т — ЯЕ) 1 )4 ёЯ

<

< 0(1) | ||(Т — ЯЕ)1 (Т — ЯЕ) 1 ёЯ< 0(1) |

, =0 [ЛI.

+ я)3 V п2

Теорема 2. Пусть р - потенциал, удовлетворяющий условию Липшица, тогда для собственных чисел оператора Т + Р верна оценка

2п 1 п—£ [ п Л

— п(п + 1)(2п +1)+—I /(а)сх%а(п — а)ёа = 01 1. ;=о 4п п ; V п' )

Интерес к подобного рода задачам все время возрастает в связи с широкой областью их применения [12-14].

Литература

1. Дубровский, В.В. Проблема решения задач на собственные значения для дифференциальных операторов со сложным вхождением спектрального параметра / В.В. Дубровский, О.А. Торшина // Новые математические методы. Электромагнитные волны и электронные системы. - 2002. - Т. 7, по. 9. - С. 4-10.

2. Дубровский, В.В. Формула первого регуляризованного следа для дифференциального оператора Лапласа-Бельтрами / В.В. Дубровский, О.А. Торшина // Дифференциальные уравнения и их приложения. - 2002. - № 1. - С. 9-19.

3. Торшина, О.А. Алгоритм вычисления регуляризованного следа оператора Лапласа -Бельтрами с потенциалом на проективной плоскости / О.А. Торшина // Вестник МаГУ. Математика. - 2003. - Вып. 4. - С. 183-215.

4. Торшина, О.А. Следы дискретных операторов с частными производными / О.А. Торшина // Альманах современной науки и образования. Научно-теоретический тематический журнал. -2012. - № 4(59). - С. 220-222.

5. Торшина, О.А. Формула первого регуляризованного следа Лапласа-Бохнера с потенциалом на проективной плоскости / О.А. Торшина // Воронежская зимняя математическая школа. -2004. - С. 104-105.

6. Торшина, О.А. Формула первого регуляризованного следа оператора Лапласа-Бельтрами с негладким потенциалом на проективной плоскости / О.А. Торшина // Вестник Самарского государственного технического университета. Математика. - 2006. - С. 32-40.

7. Торшина, О.А. Формула регуляризованного следа дифференциального оператора со сложным вхождением спектрального параметра / О.А. Торшина // Общие проблемы управления и их приложения. - 2003. - С. 467-468.

8. Торшина, О.А. Численный метод вычисления поправок теории возмущений / О.А. Торшина // Альманах современной науки и образования. - 2013. - № 12. - С. 168-170.

9. A new method for approximate evaluation of the first eigenvalues in the spectral problem of hy-drodynamic stability of poiseuille flow in a circular pipe / V.V. Dubrovskii, S.I. Kadchenko, V.F. Kravchenko, V.A. Sadovnichii // Doklady Mathematics. - 2001. - Т. 64, № 2. - С. 165-168.

10. A new method for the evaluation of the first eigenvalues in the spectral problem of hydrody-namic stability of viscous fluid flow between two rotating cylinders / V.V. Dubrovskii, S.I. Kadchenko, V.F. Kravchenko, V.A. Sadovnichii // Doklady Mathematics. - 2001. - Т. 64, № 3. - С. 425-429.

11. A new metod for approximate evaluation of the first eigenvalues in the Orr-Sommerfelg eigenvalue problem / V.V. Dubrovskii, S.I. Kadchenko, V.F. Kravchenko, V.A. Sadovnichii // Doklady Mathematics. - 2001. - Т. 63, № 3. - С. 355-358.

12. Computation of the first eigenvalues of a discrete operator / V.V. Dubrovskii, S.I. Kadchenko, V.F. Kravchenko, V.A. Sadovnichii // Электромагнитные волны и электронные системы. - 1998. -Т. 3, № 2. - С. 4.

13. Computation of the first eigenvalues of the hydrodynamic stability problem for a viscous fluid flow between two rotating cylinders / V.A. Sadovnichii, V.V. Dubrovskii, S.I. Kadchenko, V.F. Kravchenko // Differential Equations. - 2000. - Т. 36. - Вып. 6. - С. 819-824.

14. Evaluation of eigenvalues of the problem of hydrodynamic stability of viscous liquid flow between two rotating cylinders at small Reynolds numbers / V.A. Sadovnichii, V.V. Dubrovskii, S.I. Kadchenko, V.F. Kravchenko // Doklady Mathematics. - 1998. - Т. 58, № 3. - С. 483-486.

Поступила в редакцию 10 декабря 2015 г.

DOI: 10.14529/mmph160205

CALCULATION OF EIGENVALUES OF ELLIPTIC DIFFERENTIAL OPERATORS USING THE THEORY OF REGULARIZED SERIES

S.I. Kadchenko, O.A. Torshina

Magnitogorsk State Technical University of G.I. Nosova, Magnitogorsk, Russian Federation E-mail: kadchenko@masu.ru

The study of the spectral properties of perturbed differential operators is one of the significant problems of the spectral theory. In order to solve this problem it is necessary to determine the asymptotic behavior of the spectrum. But when investigating the asymptotic behavior, the improvement of remainder term is often impossible. Moreover, even the separation of the second term of the asymptotics from the remainder term is impossible. As a consequence it is necessary to come over to the study of deeper spectrum structure. A standard research tool is the derivation of formulas for regularized traces. The author makes a calculation of four amendments of the perturbation theory with the help of the theory of regularized series, followed by the access to the eigenvalues of elliptic differential operators with potential on a projective plane. In this case the projective plane is identified with the sphere by comparing opposite points and poles puncturing.

Keywords: differential operators; spectral theory; regularized traces; perturbation theory; eigenvalues.

References

1. Dubrovskiy V.V., Torshina O.A. Problema resheniya zadach na sobstvennye znacheniya dlya dif-ferentsial'nykh operatorov so slozhnym vkhozhdeniem spektral'nogo parametra [Problem solving eigenvalue problems for differential operators with a complex spectral parameter]. Novye matematicheskie

metody. Elektromagnitnye volny i elektronnye sistemy [New mathematical methods. An electromagnetic waves and electronic systems], 2002, no. 9, Vol. 7, pp. 4-10. (in Russ.).

2. Dubrovskiy V.V., Torshina O.A. Formula pervogo regulyarizovannogo sleda dlya different-sial'nogo operatora Laplasa-Bel'trami (The formula for the first regularized trace of the differential Lap-lace-Beltrami). Differentsial'nye uravneniya i ikh prilozheniya [Differential equations and their applications], 2002, no. 1, pp. 9-19. (in Russ.).

3. Torshina O.A. Algoritm vychisleniya regulyarizovannogo sleda operatora Laplasa-Bel'trami s potentsialom na proektivnoy ploskosti [The algorithm for calculation of the regularized trace for the Laplace-Beltrami operator with potential on the projective plane]. Vestnik MaGU. Matematika, 2003, Vol. 4, pp. 183-215. (in Russ.)

4. Torshina O.A. Sledy diskretnykh operatorov s chastnymi proizvodnymi [The traces of discrete operators with partial derivatives]. Al'manakh sovremennoy nauki i obrazovaniya [Almanac of modern science and education], 2012, no. 4(59), pp. 220-222. (in Russ.).

5. Torshina O.A. Formula pervogo regulyarizovannogo sleda Laplasa-Bokhnera s potentsialom na proektivnoy ploskosti (The formula for the first regularized trace of the Laplace-Bochner with potential on the projective plane). Voronezhskaya zimnyaya matematicheskaya shkola [Voronezh Winter Mathematical School], 2004, pp. 104-105. (in Russ.).

6. Torshina O.A. Formula pervogo regulyarizovannogo sleda operatora Laplasa-Bel'trami s neglad-kim potentsialom na proektivnoy ploskosti [The formula for the first regularized trace for the Laplace-Beltrami operator with potential on the non-smooth projective plane]. Vestnik Samarskogo gosu-darstvennogo tekhnicheskogo universiteta. Matematika [Vestnik of Samara State Technical University. Mathematics], 2006, pp. 32-40. (in Russ.).

7. Torshina O.A. Formula regulyarizovannogo sleda differentsial'nogo operatora so slozhnym vkhozhdeniem spektral'nogo parametra [The formula of the regularized trace of a differential operator with complex spectral parameter]. Obshchie problemy upravleniya i ikh prilozheniya [General management problems and their applications], 2003, pp. 467-468. (in Russ.).

8. Torshina O.A. Chislennyy metod vychisleniya popravok teorii vozmushcheniy [A numerical method for calculating corrections of perturbation theory]. Al'manakh sovremennoy nauki i obrazovaniya [Almanac of modern science and education], 2013, no. 12, pp. 168-170. (in Russ.).

9. Dubrovskii V.V., Kadchenko S.I., Kravchenko V.F., Sadovnichii V.A. A new method for approximate evaluation of the first eigenvalues in the spectral problem of hydrodynamic stability of poiseuille flow in a circular pipe. DokladyMathematics, 2001, Vol. 64, no. 2, pp. 165-168.

10. Dubrovskii V.V., Kadchenko S.I., Kravchenko V.F., Sadovnichii V.A. A new method for the evaluation of the first eigenvalues in the spectral problem of hydrodynamic stability of viscous fluid flow between two rotating cylinders. Doklady Mathematics, 2001, Vol. 64, no. 3, pp. 425-429.

11. Dubrovskii V.V., Kadchenko S.I., Kravchenko V.F., Sadovnichii V.A. A new metod for approximate evaluation of the first eigenvalues in the Orr-Sommerfelg eigenvalue problem. Doklady Mathematics, 2001, Vol. 63, no. 3, pp. 355-358.

12. Dubrovskii V.V., Kadchenko S.I., Kravchenko V.F., Sadovnichii V.A. Computation of the first eigenvalues of a discrete operator. Elektromagnitnye volny i elektronnye sistemy [Electromagnetic Waves and Electronic Systems], 1998, Vol. 3, no. 2, p. 4. (in Russ.).

13. Sadovnichii V.A., Dubrovskii V.V., Kadchenko S.I., Kravchenko V.F. Computation of the first eigenvalues of the hydrodynamic stability problem for a viscous fluid flow between two rotating cylinders. Differential Equations, 2000, Vol. 36, Issue 6, pp 819-824. DOI: 10.1007/bf02754405

14. Sadovnichii V.A., Dubrovskii V.V., Kadchenko S.I., Kravchenko V.F. Evaluation of eigenvalues of the problem of hydrodynamic stability of viscous liquid flow between two rotating cylinders at small Reynolds numbers. Doklady Mathematics, 1998, Vol. 58, no. 3, pp. 483-486.

Received December 10, 2015